function [ys, params, check] = my_baseline_steadystate(ys, exo, M_, options_)

xid = 0.75 ;
xix = 0.75 ;
ximc	 = 0.75 ;
ximi = 0.75 ;
ximx = 0.75 ;
kappad = 0.5 ;
kappax = 0.5 ;
kappamc = 0.5 ;
kappami = 0.5 ;
kappamx = 0.5 ;
kappaw = 0 ;
work_cap_para = 0.5 ;
sigmaLScaled = 0.2 ;
b = 0.65 ;
SppScaled = 0.5 ;
sigmaa = 0.2 ;
rhoR = 0.85	;
rpi = 1.7 ;
rdeltapi = 0 ;
rq = 0 ;
ry = 0.125 ;
rdeltay = 0 ;
etai = 1.5 ;
etaf = 1.5 ;
etac = 1.5 ;
etax = 1.5 ;
phitildes = 1.25 ;
rhomuz = 0.5 ;
rhomupsi = 0.5 ;
rhoepsilon = 0.85 ;
rhoUpsilon = 0.85 ;
rhozetac = 0.85 ;
rhophitilde = 0.85 ;
rhog = 0.75 ;
rhosigma = 0.85 ;
rhogamma = 0.85 ;
rhoeta = 0.85 ;
rhosigma_a = 0.85 ;
a11 = 0.5 ;
a22 = 0.0 ;
a33 = 0.5 ;
a12 = 0 ;
a13 = 0 ;
a21 = 0 ;
a23 = 0 ;
a24 = 0 ;
a31 = 0 ;
a32 = 0 ;
a34 = 0 ;
c21 = 0 ;
c31 = 0 ;
c32 = 0 ;
c24 = 0 ;
c34 = 0 ;
lambdaw = 1.5 ;
xiw = 0.75 ;
alpha = 0.1996 ;
beta = 0.99943 ;
muzplus = 1.0050 ;
muzPercent = 0.50 ;
pibreve = 1.005 ;
pitarget = 1.005 ;
pistar = 1.005 ;
etaa = 0.0 ;
etag = 0.291 ;
phitildea = 0.01 ;
lambdad = 1.2 ;
lambdamx = 1.2 ;
lambdamc = 1.2 ;
lambdami = 1.2 ;
lambdax = 1.2 ;
omegac = 0.25 ;
omegai = 0.43 ;
omegax = 0.35 ;
tauc = 0.25 ;  
tauw = 0.35 ; 
tauy = 0.30 ;  
tauk = 0.25 ;  
taub = 0.0 ;
nuf = 1 ;
nustar = 1 ;
nux = 1 ;
thetaw = 1;
epsilon = 1 ;
Upsilon = 1 ;
zetac = 1 ;
zetah = 1 ;
taux = 1 ;
taumc = 1 ;
taumi = 1 ; 
taumx = 1 ;
taud = 1 ;
u = 1 ;
Stilde = 0 ;
Sprimetilde = 0 ;
rhotaux = 0 ;
rhotaumc = 0 ;
rhotaumi = 0 ;
rhotaumx = 0 ;
rhotaud = 0 ;
sig_muz = 0.15 ;
sig_mupsi = 0 ;
sig_epsilon = 0.5 ;
sig_upsilon = 0.5;
sig_zetac = 0.15 ;
sig_phitilde = 0.15 ;
sig_epsR = 0.15 ;
sig_g = 0.5 ;
sig_taud = 0.15 ;
sig_taux = 0.15 ;
sig_taumc = 0.15 ;
sig_taumi = 0.15 ;
sig_taumx = 0.15 ;
sig_gamma = 0.5 ;
sig_sigma = 0 ;
sig_ystar = 0.5 ;
sig_pistar = 0.5 ;
sig_Rstar = 0.5 ;
sig_zetah = 0.15 ;
sig_eta = 0 ;
sig_sigma_a = 0 ;
sig_sigmam = 0 ;
sig_tauy = 0 ;
sig_pitarget = 0 ;
phitilde = 0 ;
rhopi = 0.975 ;
rhotauy = 0.85 ;   
rhozetah = 0.85 ;
rhosigmam = 0.85 ; 
hours_target = 0.25 ; 
iy_target = 0.18 ;
xy_target = 0.44 ;
delta = 0.015 ;
varphi = 0.2775 ;

nuf=work_cap_para;
nustar=work_cap_para;
nux=work_cap_para;

varkappad=1-kappad;
varkappax=1-kappax;
varkappamx=1-kappamx;
varkappami=1-kappami;
varkappamc=1-kappamc;
varkappaw=1-kappaw;

%some simple transformations
muz=1+muzPercent/100;
%mupsi=1+mupsiPercent/100;
mupsi=(muzplus/muz)^((1-alpha)/alpha);

sigmaL=sigmaLScaled*10;
Spp=SppScaled*10;

pic=pitarget;
pibar=pic;
pimx=pibar;
pimc=pibar;
pimi=pibar;
pii=pibar/mupsi;
piw=muzplus*pibar;
pix=pistar;

s=pic/pistar;
Rk=pibar*muzplus/beta;
R=(muzplus/beta-taub)/(1-taub)*pibar;
Rstar=(muzplus/beta-taub)/(1-taub)*pistar;
Rnustar=nustar*Rstar+1-nustar;
Rf=nuf*R+1-nuf;
Rx=nux*R+1-nux;

pitilded=pibar^kappad*pitarget^(1-kappad-varkappad)*pibreve^varkappad;

mc=1/lambdad*(1-beta*xid*(pitilded/pibar)^(lambdad/(1-lambdad)))/(1-beta*xid*(pitilded/pibar)^(1/(1-lambdad)))*((1-xid*(pitilded/pibar)^(1/(1-lambdad)))/(1-xid))^(1-lambdad);

phalo=((1-xid*(pitilded/pibar)^(1/(1-lambdad)))/(1-xid))^(1-lambdad)/((1-xid*(pitilded/pibar)^(lambdad/(1-lambdad)))/(1-xid))^((1-lambdad)/lambdad);

pitildex=pix^kappax*pix^(1-kappax-varkappax)*pibreve^varkappax;

mcx=1/lambdax*((1-xix*(pitildex/pix)^(1/(1-lambdax)))/(1-xix))^((1-lambdax))*(1-beta*xix*(pitildex/pix)^(lambdax/(1-lambdax)))/(1-beta*xix*(pitildex/pix)^(1/(1-lambdax)));

phalox=((1-xix*(pitildex/pix)^(1/(1-lambdax)))/(1-xix))^(1-lambdax)/((1-xix*(pitildex/pix)^(lambdax/(1-lambdax)))/(1-xix))^((1-lambdax)/lambdax);

pitildemc=pimc^kappamc*pitarget^(1-kappamc-varkappamc)*pibreve^varkappamc;

mcmc=1/lambdamc*((1-ximc*(pitildemc/pimc)^(1/(1-lambdamc)))/(1-ximc))^((1-lambdamc))*(1-beta*ximc*(pitildemc/pimc)^(lambdamc/(1-lambdamc)))/(1-beta*ximc*(pitildemc/pimc)^(1/(1-lambdamc)));

phalomc=((1-ximc*(pitildemc/pimc)^(1/(1-lambdamc)))/(1-ximc))^(1-lambdamc)/((1-ximc*(pitildemc/pimc)^(lambdamc/(1-lambdamc)))/(1-ximc))^((1-lambdamc)/lambdamc);

pitildemi=pimi^kappami*pitarget^(1-kappami-varkappami)*pibreve^varkappami;

mcmi=1/lambdami*((1-ximi*(pitildemi/pimi)^(1/(1-lambdami)))/(1-ximi))^((1-lambdami))*(1-beta*ximi*(pitildemi/pimi)^(lambdami/(1-lambdami)))/(1-beta*ximi*(pitildemi/pimi)^(1/(1-lambdami)));

phalomi=((1-ximi*(pitildemi/pimi)^(1/(1-lambdami)))/(1-ximi))^(1-lambdami)/((1-ximi*(pitildemi/pimi)^(lambdami/(1-lambdami)))/(1-ximi))^((1-lambdami)/lambdami);

pitildemx=pimx^kappamx*pitarget^(1-kappamx-varkappamx)*pibreve^varkappamx;

mcmx=1/lambdamx*((1-ximx*(pitildemx/pimx)^(1/(1-lambdamx)))/(1-ximx))^((1-lambdamx))*(1-beta*ximx*(pitildemx/pimx)^(lambdamx/(1-lambdamx)))/(1-beta*ximx*(pitildemx/pimx)^(1/(1-lambdamx)));

phalomx=((1-ximx*(pitildemx/pimx)^(1/(1-lambdamx)))/(1-ximx))^(1-lambdamx)/((1-ximx*(pitildemx/pimx)^(lambdamx/(1-lambdamx)))/(1-ximx))^((1-lambdamx)/lambdamx);

pitildew=pic^kappaw*pitarget^(1-kappaw-varkappaw)*pibreve^varkappaw*muzplus^thetaw;

whalo=((1-xiw*(pitildew/piw)^(1/(1-lambdaw)))/(1-xiw))^(1-lambdaw)/((1-xiw*(pitildew/piw)^(lambdaw/(1-lambdaw)))/(1-xiw))^((1-lambdaw)/lambdaw);

opts=optimset('Display','off','MaxIter',2000,'TolFun',1e-6,'TolX',1e-6);

xxx=[4      %AL
    0.025   %delta
    0.1];     %varphi

max_fsolve_tries=100;
fsolve_tries=1;
check=1; %Dynare steady state indicator (0 good, 1 bad)

while fsolve_tries<=max_fsolve_tries && check==1

    if fsolve_tries == 2
        if ~isempty(LastSuccessFsolve)
            xxx=LastSuccessFsolve;
        end

    elseif fsolve_tries>2
        ten=10;
        lowerbound=0.01;
        xxx(1)=4*(ten*rand*rand/2+lowerbound); %AL
        xxx(2)=0.025*(ten*rand*rand/2+lowerbound); %delta
        xxx(3)=0.2*(ten*rand*rand/2+lowerbound); %varphi
    end

    try
         [yyy,fval,exitflag]=fsolve(@my_baseline_steadystate_find_AL_delta_varphi,xxx,opts,pibar,mupsi, tauk,muzplus,R,alpha,Rf,epsilon,lambdad, ...
            etag,Rx,etax,omegax,etaf,omegai,etai,Upsilon,zetac,beta,b,tauc,zetah,sigmaL,omegac, ...
            etac,Rnustar,tauw,tauy,etaa,taud,s,phalo,pitildew,lambdamc,mc,piw,phalomi,phalomc, ...
            phalomx,Rstar,lambdami,lambdamx,phalox,lambdax,taumx,taumc,taumi,mcmi,mcmx,mcmc,taux,mcx, ...
            intensity_target,iy_target,xy_target,Rk,xiw,lambdaw);

        fsolve_tries=fsolve_tries+1;
 
        if  exitflag~=1
            ys=ones(size(ys))*1e+2;
            check=1;

        elseif exitflag==1
            AL=yyy(1)*10;
            delta=yyy(2);
            varphi=yyy(3);
            LastSuccessFsolve=yyy;

            pmx=taumx*Rnustar*varphi/mcmx;
            pmc=taumc*Rnustar*varphi/mcmc;
            pmi=taumi*Rnustar*varphi/mcmi;

            pc=(1-omegac+omegac*pmc^(1-etac))^(1/(1-etac));

            pinvest=(1-omegai+omegai*pmi^(1-etai))^(1/(1-etai));

            px=(taux*Rx/(mcx*varphi))*(omegax*pmx^(1-etax)+1-omegax)^(1/(1-etax));

            pkprime=pinvest/Upsilon;

            rkbar=(Rk*mupsi*pkprime/pibar-(1-delta)*pkprime-tauk*delta*mupsi*pkprime/pibar)/(1-tauk);

            wbar=(mc/(taud*(1/(1-alpha))^(1-alpha)*(1/alpha)^alpha*rkbar^alpha*Rf^(1-alpha)/epsilon))^(1/(1-alpha));

            kh=muzplus*((taud*mupsi^alpha*wbar*Rf)/((1-alpha)*mc))^(1/alpha);

            ThetacH=1/(zetah*AL)*((1-xiw*(pitildew/piw)^(1/(1-lambdaw)))/(1-xiw))^(1-lambdaw*(1+sigmaL))*(wbar*(1-tauy)/(lambdaw*(1+tauw)))*(1-beta*xiw*(pitildew/piw)^(lambdaw*(1+sigmaL)/(1-lambdaw)))*(zetac/(muzplus-b)*(muzplus-beta*b)/(pc*(1+tauc)))/(1-beta*xiw*(pitildew/piw)^(1/(1-lambdaw)));

            Thetay=phalo^(lambdad/(lambdad-1))*epsilon*(kh/(mupsi*muzplus))^alpha/(1+phalo^(lambdad/(lambdad-1))*(1/mc-phalo^(lambdad/(1-lambdad))));

            alpha1=omegac*(pc/pmc)^etac*ThetacH*phalomc^(lambdamc/(1-lambdamc));

            alpha2=(1-Rstar*s/(pibar*muzplus))*etaa*Thetay/varphi+omegai*(pinvest/pmi)^etai*kh*(1-(1-delta)/(muzplus*mupsi))/Upsilon*phalomi^(lambdami/(1-lambdami));

            alpha3=px^(1-etaf)/Rnustar-omegax*((omegax*pmx^(1-etax)+1-omegax)^(1/(1-etax))/pmx)^etax*phalox^(lambdax/(1-lambdax))*px^(-etaf)*phalomx^(lambdamx/(1-lambdamx));

            alpha4=(1-etag)*Thetay-pinvest^etai*kh*(1-(1-delta)/(muzplus*mupsi))/Upsilon*(1-omegai);

            alpha5=(1-omegac)*pc^etac*ThetacH;

            alpha6=(omegax*pmx^(1-etax)+1-omegax)^(etax/(1-etax))*(1-omegax)*phalox^(lambdax/(1-lambdax))*px^(-etaf);

            H=((alpha5+alpha6*alpha1/alpha3)/(alpha4-alpha6*alpha2/alpha3))^(1/(1+sigmaL));

            kbar=kh*H;
            ybar=Thetay*H;
            g=etag*ybar;
            q=varphi/pc;
            i=kbar*(1-(1-delta)/(muzplus*mupsi))/Upsilon;
            ystar=1/alpha3*(alpha1*H^(-sigmaL)+alpha2*H);
            x=px^(-etaf)*ystar;
            c=ThetacH*H^(-sigmaL);
            phi=ybar*(1/mc-phalo^(lambdad/(1-lambdad)));
            a=etaa*ybar;
            k=kbar*u;
            psizplus=(zetac*(muzplus-beta*b)/(muzplus-b))/(c*pc*(1+tauc));

            im=omegai*(pinvest/pmi)^(etai)*i;
            xm=omegax*(1/pmx*(omegax*pmx^(1-etax)+1-omegax)^(1/(1-etax)))^etax*phalox^(-lambdax/(lambdax-1))*px^(-etaf)*ystar;
            cm=omegac*(pc/pmc)^(etac)*c;

            Kd=lambdad*psizplus*ybar*mc/(1-beta*xid*(pitilded/pibar)^(lambdad/(1-lambdad)));
            Fd=psizplus*ybar/(1-beta*xid*(pitilded/pibar)^(1/(1-lambdad)));

            Kx=lambdax*psizplus*varphi*px*x*mcx/(1-beta*xix*(pitildex/pix)^(lambdax/(1-lambdax)));
            Fx=psizplus*varphi*px*x/(1-beta*xix*(pitildex/pix)^(1/(1-lambdax)));

            Kmc=lambdamc*psizplus*pmc*mcmc*cm/(1-beta*ximc*(pitildemc/pimc)^(lambdamc/(1-lambdamc)));
            Fmc=psizplus*pmc*cm/(1-beta*ximc*(pitildemc/pimc)^(1/(1-lambdamc)));

            Kmi=lambdami*psizplus*pmi*mcmi*im/(1-beta*ximi*(pitildemi/pimi)^(lambdami/(1-lambdami)));
            Fmi=psizplus*pmi*im/(1-beta*ximi*(pitildemi/pimi)^(1/(1-lambdami)));

            Kmx=lambdamx*psizplus*pmx*mcmx*xm/(1-beta*ximx*(pitildemx/pimx)^(lambdamx/(1-lambdamx)));
            Fmx=psizplus*pmx*xm/(1-beta*ximx*(pitildemx/pimx)^(1/(1-lambdamx)));

            Kw=(zetah*H^(1+sigmaL))/(1-beta*xiw*(pitildew/piw)^(lambdaw*(1+sigmaL)/(1-lambdaw)));
            Fw=(psizplus/lambdaw*H*(1-tauy)/(1+tauw))/(1-beta*xiw*(pitildew/piw)^(1+lambdaw/(1-lambdaw)));

            smallh=whalo^(lambdaw/(1-lambdaw))*H;
            sigmab=rkbar/pinvest;
            aofu=0.5*sigmab*sigmaa+sigmab*(1-sigmaa)+sigmab*(sigmaa/2-1);

            pitildedpi=pitilded/pibar;
            pitildexpi=pitildex/pix;
            pitildemcpi=pitildemc/pimc;
            pitildemipi=pitildemi/pimi;
            pitildemxpi=pitildemx/pimx;
            pitildewpi=pitildew/piw;

            data_pid=400*log(pibar);
            data_R=400*(R-1);
            data_wdiff=0 ;
            data_cdiff=0;
            data_idiff=0;
            data_qdiff=0;
            data_Hdiff=0;
            data_ydiff=0;
            data_xdiff=0;
            data_impdiff=0;
            data_pic=400*log(pic);
            data_pii=400*log(pii);
            data_ystardiff=0;
            data_pistar=400*log(pistar);
            data_Rstar=400*(Rstar-1);
            data_gdiff=0;
            data_Hhat=0;
            ylessd=ybar;
            Slev=1;
            nxdivGDP = 0;
            GDP12q_diff=12*100*log(muzplus);
            GDP4q_diff=4*100*log(muzplus);
            pic4q=data_pic;

            ys =[ log(mc); log(rkbar); log(wbar);  log(Rf); log(px); log(pix); log(k);
                log(Rnustar); log(H); log(pibar); log(mcx); log(pimc); log(mcmc);
                log(pimx); log(mcmx); log(pimi); log(mcmi); log(pic); log(pmc); log(pii); log(pmi);
                log(kbar);  log(i);  log(psizplus);  log(c);
                log(pkprime);  log(u);  log(pinvest);  log(s); a;  log(q); log(x); log(pc); log(R);
                log(ybar); log(Rx); log(pmx); log(Rk); Stilde; Sprimetilde;
                log(pitildedpi); log(Kd); log(Fd); log(phalo);
                log(pitildexpi); log(Kx); log(Fx); log(phalox);
                log(pitildemcpi); log(Kmc); log(Fmc); log(phalomc); log(cm);
                log(pitildemipi); log(Kmi); log(Fmi); log(phalomi); log(im);
                log(pitildemxpi); log(Kmx); log(Fmx); log(phalomx); log(xm);
                log(pitildewpi); log(Kw); log(Fw); log(whalo);
                log(smallh); aofu; log(piw);
                log(muzplus); log(epsilon);  log(mupsi);  log(pitarget);
                log(Upsilon);  log(zetac);  log(Rstar); phitilde;  log(tauy); log(zetah);
                log(ystar); log(pistar);  log(g);  log(muz);
                log(taud); log(taux); log(taumc); log(taumi); log(taumx);
                data_pid; data_R; data_wdiff; data_cdiff; data_idiff;
                data_qdiff; data_Hdiff; data_ydiff; data_xdiff; data_impdiff;
                data_pic; data_pii; data_ystardiff; data_pistar; data_Rstar;
                data_gdiff; log(ylessd); nxdivGDP; log(Slev); data_Hhat; GDP12q_diff;
                GDP4q_diff; pic4q];

           check=0;

            if check==0,
                vv=[rkbar,wbar,Rnustar,Rf,Rx,R,mc,mcx,mcmc,mcmi,mcmx,pibar,pix,pic,pii, ...
                    pimc,pimi,pc,px,pinvest,pmx,pmc,pmi,pkprime,kbar,sigmab,H,smallh,q,i,c,x, ...
                    psizplus,ybar,ystar,cm,xm,im,delta,varphi,AL];

                if min(vv) < 0
                    ys=ys*1e+2;
                    check=1;
                end
            end
        end
end


M_.params(	1	)=	alpha       	;
M_.params(	2	)=	kappad      	;
M_.params(	3	)=	beta        	;
M_.params(	4	)=	rhopi       	;
M_.params(	5	)=	xid         	;
M_.params(	6	)=	xix         	;
M_.params(	7	)=	kappax      	;
M_.params(	8	)=	ximc        	;
M_.params(	9	)=	etaf        	;
M_.params(	10	)=	ximx        	;
M_.params(	11	)=	kappamc     	;
M_.params(	12	)=	kappamx     	;
M_.params(	13	)=	kappami     	;
M_.params(	14	)=	ximi        	;
M_.params(	15	)=	omegac      	;
M_.params(	16	)=	etac        	;
M_.params(	17	)=	omegai      	;
M_.params(	18	)=	etai        	;
M_.params(	19	)=	Spp         	;
M_.params(	20	)=	muzplus     	;
M_.params(	21	)=	mupsi       	;
M_.params(	22	)=	delta       	;
M_.params(	23	)=	b           	;
M_.params(  24  ) = sigmab      ;
M_.params(	25	)=	sigmaa      	;
M_.params(	26	)=	phitildea   	;
M_.params(	27	)=	a           	;
M_.params(	28	)=	phitildes   	;
M_.params(	29	)=	rhoR        	;
M_.params(	30	)=	rpi         	;
M_.params(	31	)=	ry          	;
M_.params(	32	)=	rq          	;
M_.params(	33	)=	rdeltapi    	;
M_.params(	34	)=	rdeltay     	;
M_.params(	35	)=	etax        	;
M_.params(	36	)=	omegax      	;
M_.params(	37	)=	phi         	;
M_.params(	38	)=	rhomupsi    	;
M_.params(	39	)=	epsilon     	;
M_.params(	40	)=	pibar       	;
M_.params(	41	)=	Upsilon     	;
M_.params(	42	)=	zetac       	;
M_.params(	43	)=	tauc        	;
M_.params(	44	)=	tauk        	;
M_.params(	45	)=	R           	;
M_.params(	46	)=	tauy        	;
M_.params(	47	)=	tauw        	;
M_.params(	48	)=	zetah       	;
M_.params(	49	)=	nustar      	;
M_.params(	50	)=	nuf         	;
M_.params(	51	)=	nux         	;
M_.params(	52	)=	etag        	;
M_.params(	53	)=	lambdad     	;
M_.params(	54	)=	lambdax     	;
M_.params(	55	)=	lambdamc    	;
M_.params(	56	)=	lambdami    	;
M_.params(	57	)=	lambdamx    	;
M_.params(	58	)=	lambdaw     	;
M_.params(	59	)=	wbar        	;
M_.params(	60	)=	pic         	;
M_.params(	61	)=	psizplus    	;
M_.params(	62	)=	H           	;
M_.params(	63	)=	q           	;
M_.params(	64	)=	pimi        	;
M_.params(	65	)=	mcmi        	;
M_.params(	66	)=	mc          	;
M_.params(	67	)=	pix         	;
M_.params(	68	)=	mcmx        	;
M_.params(	69	)=	mcx         	;
M_.params(	70	)=	pimx        	;
M_.params(	71	)=	muz         	;
M_.params(	72	)=	ybar        	;
M_.params(	73	)=	u           	;
M_.params(	74	)=	kappaw      	;
M_.params(	75	)=	AL          	;
M_.params(	76	)=	sigmaL      	;
M_.params(	77	)=	Rk          	;
M_.params(	78	)=	varphi      	;
M_.params(	79	)=	ystar       	;
M_.params(	80	)=	mcmc        	;
M_.params(	81	)=	px          	;
M_.params(	82	)=	pimc        	;
M_.params(	83	)=	pmc         	;
M_.params(	84	)=	pmi         	;
M_.params(	85	)=	kbar        	;
M_.params(	86	)=	i           	;
M_.params(	87	)=	c           	;
M_.params(	88	)=	pc          	;
M_.params(	89	)=	pkprime     	;
M_.params(	90	)=	pinvest     	;
M_.params(	91	)=	pmx         	;
M_.params(	92	)=	pii         	;
M_.params(	93	)=	pistar      	;
M_.params(	94	)=	Rstar       	;
M_.params(	95	)=	rhomuz      	;
M_.params(	96	)=	rhoepsilon  	;
M_.params(	97	)=	rhoUpsilon  	;
M_.params(	98	)=	rhozetac    	;
M_.params(	99	)=	rhozetah    	;
M_.params(	100	)=	rhog        	;
M_.params(	101	)=	rhotauy     	;
M_.params(	102	)=	rhophitilde 	;
M_.params(	103	)=	phitilde    	;
M_.params(	104	)=	sig_muz     	;
M_.params(	105	)=	sig_epsilon 	;
M_.params(	106	)=	sig_upsilon 	;
M_.params(	107	)=	sig_zetac   	;
M_.params(	108	)=	sig_zetah   	;
M_.params(	109	)=	sig_phitilde	;
M_.params(	110	)=	sig_epsR    	;
M_.params(	111	)=	sig_pitarget	;
M_.params(	112	)=	sig_g       	;
M_.params(	113	)=	sig_tauy    	;
M_.params(	114	)=	s           	;
M_.params(	115	)=	Stilde      	;
M_.params(	116	)=	Sprimetilde 	;
M_.params(	117	)=	xiw         	;
M_.params(	118	)=	rkbar       	;
M_.params(	119	)=	Rf          	;
M_.params(	120	)=	Rnustar     	;
M_.params(	121	)=	x           	;
M_.params(	122	)=	Rx          	;
M_.params(	123	)=	varkappad   	;
M_.params(	124	)=	varkappax   	;
M_.params(	125	)=	varkappamx  	;
M_.params(	126	)=	varkappami  	;
M_.params(	127	)=	varkappamc  	;
M_.params(	128	)=	varkappaw   	;
M_.params(	129	)=	thetaw      	;
M_.params(	130	)=	pibreve     	;
M_.params(	131	)=	pitildedpi  	;
M_.params(	132	)=	Kd          	;
M_.params(	133	)=	Fd          	;
M_.params(	134	)=	phalo       	;
M_.params(	135	)=	pitildexpi  	;
M_.params(	136	)=	Kx          	;
M_.params(	137	)=	Fx          	;
M_.params(	138	)=	phalox      	;
M_.params(	139	)=	pitildemcpi 	;
M_.params(	140	)=	Kmc         	;
M_.params(	141	)=	Fmc         	;
M_.params(	142	)=	phalomc     	;
M_.params(	143	)=	cm          	;
M_.params(	144	)=	pitildemipi 	;
M_.params(	145	)=	Kmi         	;
M_.params(	146	)=	Fmi         	;
M_.params(	147	)=	phalomi     	;
M_.params(	148	)=	im          	;
M_.params(	149	)=	pitildemxpi 	;
M_.params(	150	)=	Kmx         	;
M_.params(	151	)=	Fmx         	;
M_.params(	152	)=	phalomx     	;
M_.params(	153	)=	xm          	;
M_.params(	154	)=	pitildewpi  	;
M_.params(	155	)=	Kw          	;
M_.params(	156	)=	Fw          	;
M_.params(	157	)=	whalo       	;
M_.params(	158	)=	smallh      	;
M_.params(	159	)=	aofu        	;
M_.params(	160	)=	taud        	;
M_.params(	161	)=	taux        	;
M_.params(	162	)=	taumc       	;
M_.params(	163	)=	taumi       	;
M_.params(	164	)=	taumx       	;
M_.params(	165	)=	rhotaud     	;
M_.params(	166	)=	rhotaux     	;
M_.params(	167	)=	rhotaumc    	;
M_.params(	168	)=	rhotaumi    	;
M_.params(	169	)=	rhotaumx    	;
M_.params(	170	)=	piw         	;
M_.params(	171	)=	taub        	;
M_.params(	172	)=	pitarget    	;
M_.params(	173	)=	etaa        	;
M_.params(	174	)=	a11         	;
M_.params(	175	)=	a12         	;
M_.params(	176	)=	a13         	;
M_.params(	177	)=	a21         	;
M_.params(	178	)=	a22         	;
M_.params(	179	)=	a23         	;
M_.params(	180	)=	a24         	;
M_.params(	181	)=	a31         	;
M_.params(	182	)=	a32         	;
M_.params(	183	)=	a33         	;
M_.params(	184	)=	a34         	;
M_.params(	185	)=	sig_ystar   	;
M_.params(	186	)=	c21         	;
M_.params(	187	)=	sig_pistar  	;
M_.params(	188	)=	c24         	;
M_.params(	189	)=	c31         	;
M_.params(	190	)=	c32         	;
M_.params(	191	)=	sig_Rstar   	;
M_.params(	192	)=	c34         	;
M_.params(	193	)=	SppScaled   	;
M_.params(	194	)=	muzPercent  	;
M_.params(	195	)=	mupsiPercent	;
M_.params(	196	)=  intensity_target ;
M_.params(	197	)=  iy_target ;
M_.params(	198	)=  xy_target ;
M_.params(199)   =  ylessd ;

end