DATA_SECTION
  init_int nobs;
  init_matrix data(1,nobs,1,2)
  vector stock(1,nobs);
  vector recruit(1,nobs);

INITIALIZATION_SECTION

PARAMETER_SECTION
  init_bounded_number alpha(0.,10.);
  init_bounded_vector stockdevs(1,nobs,-5.,5.,2);
  init_number beta;
  vector pred_recruit(1,nobs)
  vector true_stock(1,nobs)
  objective_function_value f;

PRELIMINARY_CALCS_SECTION
  // get the data out of the columns of the data matrix 
  stock=column(data,1);
  recruit=column(data,2);

PROCEDURE_SECTION
  true_stock=elem_prod(mfexp(stockdevs),stock);
  pred_recruit= alpha*elem_prod(true_stock,mfexp(-beta*(true_stock/true_stock(1)-1.)));
  f=robust_regression(log(recruit),log(pred_recruit),.7);
  f+=12.5*norm2(stockdevs);
  //f+=square((alpha-.5))

REPORT_SECTION
  report <<pred_recruit<<endl;
  report <<recruit<<endl;
  report <<true_stock<<endl;
  report <<stock<<endl;