var k,c,w,r; 
parameters delta,TW,TR,TC,qw,alpha,cho,sigma,l; 
sigma= 0.5;   
alpha= 0.3;�@
cho= 0.01; 
delta=0.05;�@
TW=0.20;   
TR=0.05;   
TC=0.10;  
qw=0.1;     
l=100;       
cstar=(1-TR*alpha+(alpha-1)*(TW+qw)*k^alpha*l^(1-alpha))/(1-TR);
rstar=((1+cho)^(sigma)-1)/(1-TR);
kstar=l*((1+cho)^(sigma)-1/(alpha*(1-TR)))^(1/(alpha-1));
wstar=(1-alpha)*((1+cho)^(sigma)-1/(alpha*(1-TR)))^(1/(alpha-1))^alpha;
model; 
c(+1)=c*((1+r(+1)-TR*r(+1))/(1+cho))^(1/sigma);
k=k(-1)*(1-delta)+(1-TR)*r*k+(1-TW-qw)*w*l-(1-TR)*c;
r =alpha* (l/k)^(1-alpha);
w = (1-alpha)*(k/l)^alpha;
end;
initval;
k=kstar;
w=wstar;
c=cstar;
r=rstar;
end;
steady;