var dd dnd lambdad lambdand i kd knd bd yd ynd;

predetermined_variables kd knd bd;

varexo e;

varexo_det z;

parameters beta delta omeg gam theta;
    beta = 0.8;
    delta = 0.03;
    omega = 3;
    gam = 1;
    theta = 0.7;
    r = 0.04;

model;
yd = gam*z*kd^(alpha);
ynd = gam*z*knd^(alpha);

kd(+1) = (1-delta)*k + i;
dd =