function [y,yp,ys,ch,h,d,ks,wage,DIV_S,hs,kp,DIV_P,ps,Rs,Rp,PI_CB,k,i,utility,m,bg]=solve_star_nosb(lp,ls,hp,pp,PI,DELTA,THETA_S,THETA_P,ALPHA,Ap,ETA,GAMMA,TAU_S,TAU_P,KAPPA,SIGMA,RL,BETA_B,PHI,ALPHA_S,ALPHA_P,BETA,KAPPA_M,EPSILON,ZETAL)
R=1/BETA*PI;
pw=(EPSILON-1)/EPSILON;lambda_c=1-BETA_B*R/PI;

Rs=(ALPHA_S*ls+1-TAU_S-lambda_c*(1-GAMMA))/BETA_B*PI+RL*TAU_S;
Rp=(ALPHA_P*lp+1-TAU_P-lambda_c*(1-GAMMA*ZETAL))/BETA_B*PI+RL*TAU_P;
kp=lp*Rp/THETA_P/PI;
ks=ls*Rs/THETA_S/PI;
d=lp+ls-GAMMA*(ls+ZETAL*lp);

yp=Ap*kp^ALPHA*hp^(1-ALPHA);
wage=(1-ALPHA)*pp*yp/hp;
y=yp/(1-PHI)^SIGMA/(pp/pw)^(-SIGMA);
DIV_P=(-DELTA)*kp+pp*yp-wage*hp+(1-Rp/PI)*lp;
k=ks+kp;
ps=((pw^(1-SIGMA)-(1-PHI)^(SIGMA)*pp^(1-SIGMA))/(PHI^SIGMA))^(1/(1-SIGMA));
ys=(ps/pw)^(-SIGMA)*PHI^(SIGMA)*y;
hs=(1-ALPHA)*ps*ys/(wage);
i=DELTA*k;
DIV_S=(-DELTA)*ks+ps*ys-wage*hs+(1-Rs/PI)*ls;
PI_CB=(Rp*lp+Rs*ls-R*d-RL*(TAU_P*lp+TAU_S*ls))/PI-ALPHA_P/2*lp^2-ALPHA_S/2*ls^2-(1-TAU_S)*ls-(1-TAU_P)*lp+d;
ch=y-i;
h=(wage/ch/KAPPA)^(1/ETA);
lambda_h=1/ch;
m=KAPPA_M/(lambda_h*(1-BETA/PI));
bg=((TAU_P*lp+TAU_S*ls)*(RL/PI-1)+(1-1/PI)*m)/(1-R/PI);
utility=(log(ch)-KAPPA*h^(1+ETA)/(1+ETA)+KAPPA_M*log(m))/(1-BETA);