/* Funkcni hodnota a gradient 2-dim. Rosenbrockovy funkce. */

double f(long int ndim, double x[], double grad[])
{
  double r1,r2;
  
  r1= x[1]-x[0]*x[0];
  r2= 1-x[0];
  grad[0]= -400*x[0]*r1 -2*r2;
  grad[1]= 200*r1;
  return(100*r1*r1 +r2*r2);
}