hi
I really have problems with errors my model could please help me to complete it?
//=========================================================
// sim model(2018) NX model, 5 structural schoks
//=========================================================
// ENDOGENOUS VARIABLES
var
xt rr rf df gdp rst gdps qf cons ins pi pih pin delhn inss
pis ad adh adn as nxp nxpf ph pn ps gdpf pnp php qt epsic epsim epsin ;
//=========================================================
// EXOGENOUS INNOVATIONS
varexo
nuc num nul hepsilon nepsilon sepsilon;
//=========================================================
// DEEP AND STUCTURAL PARAMETERS
parameters
cbeta csigma ctau cpsi clambda cnu ceta calpha
cepsilon ckapa comegah comegan comegaf comegastarh comegastarf
crohc crohm crohl crohdh crohdn crohstara crohi cpsi1
cpsi2 ckapa1 ckapa2 ckapa3 ;
//=========================================================
// composit parameters
comega cphi ctheta psih psin psif psistarf f h g cgamma omegag m CH CZ;
//=========================================================
// fixed parameters
Cbeta = 0.98;
csigma = 1.57;
ctau= 2.1;
cpsi= 2.39;
Clambda=0.98;
Cnu0=0.43;
Ceta=1.56;
Calpha=0.54;
Cepsilon=4.33;
Ckapa=0.0007;
Comegah=0.68;
Comegan=0.2;
Comegaf=0.5;
Comegastarh=0.5;
Comegastarf=0.2;
Crohc=0.5;
Crohm=0.2;
Crohl=0.5;
Crohdh=0.2;
Crohdn=0.5;
Crohstara=0.2;
Crohi=0.5;
cpsi1= 0.5;
cpsi2 =0.5;
ckapa1 =0.2;
ckapa2 =0.2;
ckapa3= 0.5;
//=========================================================
// NX MODEL LINEARIZED EQUATIONS
model(Linear);
//=========================================================
//Composite paramters
# comega=cnu*(calpha-1)+ceta*(2-calpha) ;
# cphi= ceta*(2-calpha)-calpha*clambda*cnu;
# ctheta = calpha*(1-clambda);
# psih = (1-cbeta*comegah)*(1-comegah)/((1+cbeta)*comegah);
# psin = (1-cbeta*comegan)*(1-comegan)/((1+cbeta)*comegan);
# psif = (1-cbeta*comegaf)*(1-comegaf)/((1+cbeta)*comegaf);
# psistarf = (1-cbeta*comegastarf)*(1-comegastarf)/((1+cbeta)*comegastarf);
# f = (calpha*(1-clambda)*( ceta*(2-calpha)-calpha*clambda*cnu)*csigma+(1-ctheta)^2)/csigma;
# h = 1-clambda/(1-calpha*(1-clambda))-calpha*(1-clambda)*clambda*(cnu*(calpha-1)+ceta*(2-calpha));
# g = ckapa*calpha*(1-clambda)/(1-calpha*(1-clambda))*((calpha*(1-clambda)*( ceta*(2-calpha)-calpha*clambda*cnu)*csigma+(1-ctheta)^2)/csigma-ceta*(2-calpha)+calpha*clambda*cnu);
# cgamma = csigma*(ceta*(2-calpha)-calpha*clambda*cnu)*calpha*(1-clambda)/(1-calpha*(1-clambda))-calpha*(1-clambda);
# omegag = csigma*clambda*(cnu*(calpha-1)+ceta*(2-calpha))*(calpha*(1-clambda)/(1-calpha*(1-clambda)))-clambda;
# m = 1/(cbeta*(1+ckapa));
# CH = calpha*(ceta*(2-calpha)-calpha*clambda*cnu)*(calpha*(1-clambda)/(1-calpha*(1-clambda)))-calpha*(1-clambda);
# CZ =calpha*(cnu*(calpha-1)+ceta*(2-calpha))*clambda*(calpha*(1-clambda)/(1-calpha*(1-clambda)))-clambda;
//=========================================================
//LinearsEquations
xt=xt(+1)-f*(rr-rf)-g*df;
rf=1/f*(gdpf(+1)-gdpf)+(csigma*f-1)/(csigma*f)*rst-(1-ctheta)/(csigma*f)*(epsic(+1)-epsic);
gdpf=1/(csigma+cpsi*(1-ctheta))*(csigma*ctheta*gdps+cgamma*qf+(1+cpsi)*(1-ctheta)*ad+(1-ctheta)*(epsic-epsin));
gdp=(1-ctheta)*cons+ctheta*gdps-ctheta*clambda*comega*(1-ctheta)*delhn+ctheta*cphi*(1-ctheta)*qt;
rr=ins-(h*pih(+1)+(1-h)*pin(+1));
gdps=gdps(+1)-1/csigma*(inss-pis(+1));
rst=inss-pis(+1);
xt=gdp-gdpf;
ad=(1-clambda)*adh+clambda*adn;
pih=cbeta/(1+cbeta)*pih(+1)+1/(1+cbeta)*pih(-1)+psih*(cpsi+csigma/(1-ctheta))*xt-cgamma*psih/(1-ctheta)*(qt-qf)+omegag*psih/(1-ctheta)*delhn;
pin=cbeta/(1+cbeta)*pin(+1)+psin*(cpsi+csigma/(1-ctheta))*xt-psin*cgamma/(1-ctheta)*(qt-qf)+psin*(1+omegag/(1-ctheta))*delhn;
pis=cbeta/(1+cbeta)*pis(+1)+1/(1+cbeta)*pis(-1)+psistarf*((cpsi-csigma)*gdps-(1+cpsi)*as);
pi=(1-clambda/(1-ctheta))*pin+(clambda/(1-ctheta))*pin+(ctheta/(1-ctheta))*(qt-qt(-1));
qt=qt(+1)-csigma/(1+cgamma)*(xt(+1)-xt)-csigma*ctheta*clambda*comega/((1-ctheta)*(1+cgamma))*(pih(+1)-pin(+1))-ckapa*(1-ctheta)/((1+cgamma))*df-
(qf(+1)-qf);
qf=csigma/(1+cgamma)*(gdpf-gdps)-(1-ctheta)/(1+cgamma)*epsic;
df=m*df(-1)-m*(1-ctheta-clambda)/(1-ctheta)*pih-m*clambda/(1-ctheta)*pin-m*ctheta/(1-ctheta)*(qt-qt(-1)-qf+qf(-1))+m*(cbeta-1)*(nxp-nxpf);
nxp=1/((1-ctheta)*(1-calpha))*((1-ctheta-calpha)*xt+CZ*delhn+CH*(qt-qf))+nxpf;
nxpf= 1/((1-ctheta)*(1-calpha))*((1-ctheta-calpha)*gdpf+calpha*ctheta*gdps+CH*qf);
ins=crohi*ins(-1)+(1-crohi)*(cpsi1*pi+cpsi2*xt);
qt=qt(-1)-ckapa1*df(-1)+ckapa2*(df(-1)-df)-ckapa3*pi(-1);
pih=ph-ph(-1);
pin=pn-pn(-1);
pis=ps-ps(-1);
pnp=(ctheta-1+clambda)/(1-ctheta)*delhn-ctheta/(1-ctheta)*qt;
php=clambda/(1-ctheta)*delhn-ctheta/(1-ctheta)*qt;
delhn=ph-pn;
epsic=crohc*epsic(-1)+nuc;
epsim=crohm*epsim(-1)+num;
epsin=crohl*epsin(-1)+nul;
adh=crohdh*adh(-1)+hepsilon;
adn=crohdn*adn(-1)+nepsilon;
as=crohstara*as(-1)+sepsilon;
end;
//=========================================================
initval;
xt=0.33;
rr=0.02;
rf=0.005;
df=0.23;
gdp=0.50;
rst=0.005;
gdps=0.99;
qf=0.20;
cons=0.30;
ins=0.02;
pi=0.25;
pih=0.23;
pin=0.15;
delhn=0.01;
inss=0.01;
pis=0.11;
ad=0.10;
adh=0.10;
adn=0.10;
nxp=0.3;
nxpf=0.2;
ph=0.25;
pn=0.2;
ps=0.10;
epsic=0.010;
epsim=0.05;
epsin=0.02;
as=0.05;
gdpf=0.01;
pnp=0.01;
php=0.01;
qt=0.1;
end;
//=========================================================
check;
steady;
shocks;
var nuc; stderr 0.01;
var num; stderr 0.01;
var nul; stderr 0.01;
var hepsilon; stderr 0.01;
var nepsilon; stderr 0.01;
var sepsilon; stderr 0.01;
end;
stoch_simul(periods=30000, irf=30, graph_format =pdf);