profile on
  [X,y,z,Z,muu,itersinfo,meritfn]=theta(n,indsE,toler,crosstoler);
profile report
