remz0.c

来自「数值分析算法描述与习题解答」· C语言 代码 · 共 25 行

  #include "math.h"
  #include "stdio.h"
  #include "remz.c"
  main()
  { int i;
    double a,b,eps,p[5],remzf();
    a=-1.0; b=1.0; eps=1.0e-10;
    remz(a,b,p,4,eps,remzf);
    printf("\n");
    for (i=0; i<=3; i++)
      printf("p(%2d)=%e\n",i,p[i]);
    printf("\n");
    printf("MAX(p-f)=%e\n",p[4]);
    printf("\n");
  }

  double remzf(x)
  double x;
  { double y;
    y=exp(x);
    return(y);
  }