romb0.c

徐士良写的《数值分析与算法》

  #include "stdio.h"
  #include "romb.c"
  main()
  { double a,b,eps,t,rombf();
    a=0.0; b=1.0; eps=0.000001;
    t=romb(a,b,eps,rombf);
    printf("\n");
    printf("t=%e\n",t);
    printf("\n");
  }

  double rombf(x)
  double x;
  { double y;
    y=x/(4.0+x*x);
    return(y);
  }