#
#  PROTO, Sensitivity Analysis Example
#
#  Equations of the model:
#
    eq1 :=           y = w*h + s  :
    eq2 :=   p*(1+t)*c = a*y      :
    eq3 :=         w*j = (1-a)*y  :
    eq4 :=         w*l = a*y - s  :
    eq6 :=           l = q/b      :
    eq7 :=           p = w/(b-1)  :
    eq8 :=       t*p*c = s        :
#
#  Some parameters and exogenous variables:
#
    b := 2.040:
    h := 100.0:
    p := 1.000:
#
#  For several values of parameter "a" solve 
#  the base run (tax=0.2), change the tax and then
#  solve the policy run (tax=0.3):
#
#  CASE 1: a = 0.285
#
    a := 0.285:
    t := 0.200:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
    t := 0.300:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
#
#  CASE 2: a = 0.255
#
    a := 0.255:
    t := 0.200:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
    t := 0.300:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
#
#  CASE 3: a = 0.315
#
    a := 0.315:
    t := 0.200:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
    t := 0.300:
    fsolve({eq1,eq2,eq3,eq4,eq6,eq7,eq8},{y,w,s,c,l,j,q});
#
#  Done
#		 
    quit;