eps,fopt,x[3],xm[3],g[3];

 double eps,fopt,x[3],xm[3],g[3];
 int itmax = 50, n = 3,it;
 double f();
 void df();
 eps = 1e-7;