def y_eq_func(ylag,pilag,v,s,slag,alpha,h,b,phi,gamma):
    """ equilibrium value for output

    Args:

        ylag (float): lagged output
        pilag (float): lagged inflation
        v (float): demand disturbance
        s (float): supply disturbance
        slag (float): lagged supply disturbance
        alpha (float): sensitivity of demand to real interest rate
        h (float): coefficient on inflation in Taylor rule
        b (float): coefficient on output in Taylor rule
        phi (float): degree of stickiness in inflation expectations
        gamma (float): effect of output on inflation in SRAS

    Returns:

        (float):  equilibrium value for output

    """

    return 1/(alpha*b+alpha*gamma*h+1)*(-pilag*alpha*h+alpha*gamma*h*phi*ylag+alpha*h*phi*slag-alpha*h*s+v)