The below program generates an error message '??? Error: Unexpected MATLAB operator'
Can anyone help?
var k pi c n b g a;
varexo epsilon_g epsilon_a;
parameters alpha beta delta gamma_c gamma_c_tilde eta lambda lambda_p mu_p phi_b phi_g phi_pi psi rho rho_a rho_g theta theta_n theta_t omega sigma_tilde;
alpha = .33;
beta = .99;
delta = .025;
gamma_c = .7;
gamma_c_tilde = .1;
eta = 1;
lambda_p = .5;
mu_p = .2;
phi_b = .33;
phi_g = .1;
phi_pi = 1.5;
psi = .2;
rho = .01;
rho_a = .97;
rho_g = .9;
theta = .75;
theta_n = .27;
theta_t = .07;
sigma_tilde = 0.5;
model (linear);
k = (1-delta+(delta*alpha/(1-gamma_c_tilde)))*k(-1) + (delta*(1-alpha)/(1-gamma_c_tilde))*n(-1) - (delta*gamma_c/(1-gamma_c_tilde))*c(-1) - (delta/(1-gamma_c_tilde))*g(-1) + (delta/(1-gamma_c_tilde))*a(-1);
pi = (pi(-1) - lambda_p*c(-1) + lambda_p*alpha*k(-1) - (alpha + psi)*lambda_p*n(-1) + lambda_p*a(-1)) / beta;
c = c(-1) - theta_n*n(-1) + theta_n*phi_b*b(-1) - (phi_pi/sigma_tilde)*pi(-1) - (1/sigma_tilde)*pi + theta_n*n - theta_t*phi_b*b - theta_t*phi_g*(rho_g - 1)*g(-1);
n = ( (1-alpha)*n - gamma_c*c(-1) - (1-gamma_c_tilde)*k(-1) + ((1-gamma_c_tilde)*eta*phi_pi)*pi(-1) - (omega-(beta*gamma_c))*c + (omega+beta*(1-gamma_c_tilde-alpha))*k - ((1-gamma_c_tilde)*eta)*pi - (1-(beta*rho_g))*g(-1) + (1-(beta*rho_a))*a(-1) ) / (omega*(1+psi)+beta*(1-alpha));
b = (1 + rho)*(1 - phi_b)*b(-1) + (1 + rho)*(1 - phi_g)*g(-1);
g = rho_g*g(-1) + epsilon_g;
a = rho_a*a(-1) + epsilon_a;
end;
initval;
k = 9;
pi = 1;
c = .76;
n = .3;
b = .76;
a = 1;
g = 0;
end;
steady;
shocks;
var a = .007;
var g = .005;
end;
stoch_stimul;