function [param, ss] = calib

ALPHA=0.7556;
BETA=0.99;
ETA=2;
XSS=1.1;
PHIP=0.5;
PHIPm=0.06;
RHO=0.8;
RHOa=0.9;
RHOr=0.98;
RHOx=0.95;
THETA=1.5;
KTILDE=0.0858;
cratio=0.97;
imratio=0.2;
xratio=0.23;
STDERR_AE=0.01;

a=1;
v=0.9;
rf=1.0175;
dp=1;
dpf=1.005;
q=1;
ge=1;
mu=dp;
l=0.4;
r=dp/BETA;
w=a/v;
yj=a*l;
y=yj*((1/6)/q)^((1-ALPHA)/ALPHA);
gdp=y/(1+imratio);
x=xratio*gdp;
yf=imratio*gdp/q;
c=cratio*gdp;
lamda=1/c;
xf=(q^THETA)/x;
m=PHIPm/(lamda*(1-BETA/dp));
bf=0.03*gdp/(q*(1-rf/dpf));
ca=q*bf*(1-1/dpf);
b=(q*bf*(1-rf/dpf)-m*(1-1/mu))/(1-r/dp);

param = [ALPHA, BETA, ETA, XSS, KTILDE, THETA, RHO, RHOa, RHOr, RHOx, PHIP, PHIPm, cratio, imratio, xratio, STDERR_AE, dpf]';
ss = [c, l, w, r, dp, lamda, m, v, yf, yj, q, ca, x, b, bf, y, mu, rf, ge, a, xf, gdp]';



