problem for solving DSGE model

[luckygift2002]

Dear All,

I have come aross a problem when solving my GSDE model. The error message shows in DYNARE version 3 is
??? Error using ==> dynare at 19
syntax error at line 70

- which is the command for end;

But if I am using version 4 , the error message would be
Warning: Divide by zero.
> In solve1 at 107
In dynare_solve at 94
In steady_ at 55
In steady at 27
In DSGE at 242
In dynare at 79
??? Error using ==> lnsrch1 at 41
Some element of Newton direction isn't finite. Jacobian maybe singular or there is a problem with
initial values

Error in ==> solve1 at 114
[x,f,fvec,check]=lnsrch1(xold,fold,g,p,stpmax,func,j1,j2,varargin{:});

Error in ==> dynare_solve at 94
[x,info]=solve1(func,x,j1(r(i):r(i+1)-1),j2(r(i):r(i+1)-1),jacobian_flag,varargin{:});

Error in ==> steady_ at 55
[oo_.steady_state,check] = dynare_solve([M_.fname '_static'],...

Error in ==> steady at 27
steady_;

Error in ==> DSGE at 242
steady;

Error in ==> dynare at 79
evalin('base',fname) ;


I have been trying to find the problems these days but I still cannot find any. Could anyone help me with this? I would be very appreciated for this. Many thank.

Jing