function y = sscumc(x,param)

hrs    = x(1);
emp    = x(2);
c1     = x(3);
k      = x(4);

alpha  = param(1);
beta   = param(2);
delta  = param(3);
gamma  = param(4);
lam    = param(5);
omega  = param(6);
pa     = param(7);
pb     = param(8);
psi    = param(9);
sig    = param(10);
theta  = param(11);
ig     = param(12);
bg     = param(13);
army   = param(14);

z      = 1; 
ome1   = 1-omega;
thet1  = 1-theta;
gamm1  = 1-gamma;
hwgt   = (1-gamma)*(1-omega);
cwgt   = gamma*(1-omega);
cwgt1  = cwgt-1;
hwgt1  = hwgt-1;

zkk    =  z*lam*k^theta;
c0     =  c1*(1-hrs)^(hwgt/cwgt1);
%dpn   =  pa+2*pb*emp;
d      =  (emp-pa)*pb;
dpn    =  3.d0*d*d*(1+sign(d))*pb;
y(1,1) =  gamm1/gamma*emp^theta*c1-zkk*(1-hrs);
y(2,1) =  c1^cwgt*(1-hrs)^hwgt/ome1-c0^cwgt/ome1-dpn- ...
          gamma*c1^cwgt1*(1-hrs)^hwgt*(c1-c0-zkk*thet1*emp^(-theta)*hrs);
y(3,1) =  emp*c1+(1-emp)*c0+delta*k+bg-zkk*hrs*emp^thet1;
y(4,1) =  beta*(1-delta+theta*zkk/k*emp^thet1*hrs)-1;

