alpha=0.3;
beta=0.99;
delta=.025;
gamma=1.1;
rho=0.8;
Ve=1.;
lg=0;


[gk, ga, ge, gsig, gkk, gaa, gee, gsigsig, gka, gke, gksig, gae, ...
        gasig, gesig, kss] = second(alpha,beta,delta,gamma,rho,Ve,lg);