i tried many times but i still can not figure out why matlab show such message,could somebody help me? i am puzzled.

-------------------------------------------------matlab message---------------------------------------------------------

Configuring Dynare …

[mex] Generalized QZ.

[mex] Sylvester equation solution.

[mex] Kronecker products.

[mex] Sparse kronecker products.

[mex] Bytecode evaluation.

[mex] k-order perturbation solver.

[mex] k-order solution simulation.

Starting Dynare (version 4.1.2).

Starting preprocessing of the model file …

Found 14 equation(s).

Evaluating expressions…done

Computing static model derivatives:

- order 1

Computing dynamic model derivatives: - order 1
- order 2

Processing outputs …done

Preprocessing completed.

Starting MATLAB/Octave computing.

??? Error: File: C:\Users\Thinkpad\Desktop\dynare\quant and price\monetary.m Line: 131 Column: 1

Missing variable or function.

Error in ==> dynare at 132

evalin(‘base’,fname) ;

-------------------------the mod file--------------------------------------------------------------------------------------------------------------

var c y k R i pai phai omigabar m z v n x r;

varexo eps1 eps2;

parameters h sigma alpha theta delta omiga beta bphai ingta psai l gamma la lb lc kk rhox;

h = 0.65;

sigma = 3.57;

alpha = 0.4;

theta = 11;

delta = 0.024;

omiga = 0.5884;

beta = 0.98;

bphai = 0.7533;

ingta = 0.9728;

psai = 0.60;

l = 2.0;

gamma = 6.16;

la = 0.62;

lb = 1.27;

lc = 0.9238;

kk = 0.66;

rhox = 0.85;

model(linear);

pai=omiga*pai(-1)/(1+beta*omiga^2)+(2-omiga-beta*omiga+beta*omiga^2)*pai(+1)/(1+beta*omiga^2)+(1-omiga)*(1-beta*omiga)*phai;
k(+1)=(1-delta) k+deltai;
c=h*c(-1)/(1+h)+c(+1)/(1+h)-(1-h)

*(R-pai(+1))/((1+h)*(1+ingta

*sigma);*

c=hc(-1)+c(+1)-(1-h)c=h

*r(+1)/sigma;*

y=(1+alpha(theta-1)/theta)y=(1+alpha

*c+(1-alpha)*(theta-1)*i/theta;*

y=z+alphan+(1-alpha)y=z+alpha

*k;*

omigabar/bphai=(1+ingtal)omigabar/bphai=(1+ingta

*psai*(omigabar(-1)+pai(-1))-(1+beta)*l)*l)

*psai*pai+(1+ingta*psai*beta*(omigabar(+1)+pai(+1))+(1-beta

*psai)*(1-psai)

*(ingta*n+sigma*(c-c(-1))/(1-h));

m=sigma

*c/((1-h)*c(-1)/((1-h)

*gamma)- sigma*h*gamma)-R/gamma;*

m=m(-1)-pai+v;

v=lav(-1)-lb

m=m(-1)-pai+v;

v=la

*pai(+1)-lc*y+x;

phai=alpha

*omigabar+(1-alpha)*

x=rhoxx(-1)+eps2;

*r-z;*

n=r+k-omigabar;

z=kkz(-1)+eps1;n=r+k-omigabar;

z=kk

x=rhox

end;

steady;

check;

shocks;

var eps1; stderr 0.01;

var eps2; stderr 0.01;

end;

stoch_simul;

monetary.mod (1.21 KB)