init = [
#model.shortL0
#  1;
#model.shortL1
#  1;
#model.shortLength
  -10;
#model.openC0
  1;
#model.openC1
  0;
#model.openC2
#  0;
#model.openLength
#  0;
#model.loadR
#  50;
#model.loadR1 
#  2e-9;
#model.loadL0=
#  1;
#model.loadL1=
#  -1e-24;
#model.loadLength=
  1; 
  ]
  
  
 t_init = 20;
 t_min = 1;
 mu_t = 1.002;
 iters_fixed_t = 10;
 max_rand_step = [0.1; 0.1; 0.1; 0.1; 0.1];
 
  opt= optimset ("max_rand_step", max_rand_step, "t_init", t_init, "T_min", t_min, "mu_t", mu_t,"iters_fixed_T", iters_fixed_t); ##, "trace_steps", true)
  
  f = @ residue;
  [p, residuals, cvg, outp] = nonlin_residmin (f, init, opt) ;
  p