7nmtc.c

来自「常用算法的C程序。主要内容包括多项式的计算、复数运算、随机数的产生、矩阵运算、矩」· C语言 代码 · 共 29 行

  #include "stdlib.h"
  #include "math.h"
  #include "3rnd1.c"
  void nmtc(x,n,b,m,eps,f)
  int n,m;
  double x[],b,eps,(*f)();
  { int k,i;
    double a,r,*y,z,z1;
    y=malloc(n*sizeof(double));
    a=b; k=1; r=1.0; z=(*f)(x,n);
    while (a>=eps)
      { for (i=0; i<=n-1; i++)
          y[i]=-a+2.0*a*rnd1(&r)+x[i];
        z1=(*f)(y,n);
        k=k+1;
        if (z1>=z)
          { if (k>m) { k=1; a=a/2.0; }}
        else
          { k=1; 
            for (i=0; i<=n-1; i++) x[i]=y[i];
            z=z1;
            if (z<eps)  return;
          }
      }
    return;
  }