www.pudn.com > cghost.rar > JMAXN0.C

 
  #include "stdio.h" 
  #include "jmaxn.c" 
  main() 
  { int js[2]; 
    double eps,x[3]; 
    eps=0.000001; x[0]=0.0; x[1]=0.0; 
    jmaxn(x,2,eps,10,js); 
    printf("\n"); 
    printf("js(0)=%3d\n",js[0]); 
    printf("x(0)=%e\n",x[0]); 
    printf("x(1)=%e\n",x[1]); 
    if (js[1]<0.0) printf("MIN:  "); 
    if (js[1]>0.0) printf("MAX:  "); 
    printf("z=%e\n",x[2]); 
    printf("\n"); 
  } 
 
  double jmaxnf(x,n,j) 
  int n,j; 
  double x[]; 
  { double y; 
    n=n; 
    switch(j) 
      { case 0: y=(x[0]-1.0)*(x[0]-10.0) 
                  +(x[1]+2.0)*(x[1]+2.0)+2.0; 
                break; 
        case 1: y=2.0*(x[0]-1.0); break; 
        case 2: y=2.0*(x[1]+2.0); break; 
        defaut: {} 
      } 
    return(y); 
  }