
  params.sigma = 0.1;  % standard deviation of shock
  params.nK = 2000;
  params.nZ = 11;
  [r2,sol2] = ha1(params);

  
