function [equs]=my_baseline_steadystate_find_AL_delta_varphi(xxx,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)

AL=xxx(1)*10;
delta=xxx(2);
varphi=xxx(3);


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;
i=kbar*(1-(1-delta)/(muzplus*mupsi))/Upsilon;
ystar=1/alpha3*(alpha1*H^(-sigmaL)+alpha2*H);
x=px^(-etaf)*ystar;
whalo=((1-xiw*(pitildew/piw)^(1/(1-lambdaw)))/(1-xiw))^(1-lambdaw)/((1-xiw*(pitildew/piw)^(lambdaw/(1-lambdaw)))/(1-xiw))^((1-lambdaw)/lambdaw);
smallh=whalo^(lambdaw/(1-lambdaw))*H;

equ1=smallh-intensity_target;
equ2=i*pinvest/ybar-iy_target;
equ3=px^(-etaf)*ystar*px*varphi/ybar-xy_target;
equs=[equ1
      equ2
      equ3];