复化辛普生.c

来自「数值分析的几个C实现,比较小.」· C语言 代码 · 共 18 行

#include<stdio.h>
#include<math.h>
main()
{int k,n=6;float f();
float a=0,b=1,x=0,s=0,h;
/*scanf("%f%f%d\n",&a,&b,&n);*/
h=(b-a)/n;s=f(a)-f(b);x=a;
for(k=1;k<=n;k++)
{ x+=h/2;s+=4*f(x);
  x+=h/2;s+=2*f(x);
  }
 s*=h/6;
 printf("%.6f\n",s);getch();
 }
float f(float m)
 { return(m*sin(m));
 }