ts.c

来自「两种曲线拟合的算法」· C语言 代码 · 共 49 行

#include <stdio.h>
#include <math.h>
#define f(x) (cos(x)+sin(x)*sin(x))
tixing(float h,float n)
{int i;
 double s,a;
 a=0;
 s=f(a)+f(a+(n-1)*h);
 for(i=1;i<n-1;i++)
	 s=s+2*f(a+i*h);
 
 s=s*h/2;
 printf("梯形公式结果为：T%f=%lf",n-1,s);
}

simpson(float h,float n)
{int i;
 double s,a,b,x;
 a=0;b=1;
 s=f(a)-f(b);
 x=a;
 for(i=1;i<n;i++)
 {
	 x=x+h/2;
	 s=s+4*f(x);
	 x=x+h/2;
	 s=s+2*f(x);
 }
 s=s*h/6;
 printf("辛浦生公式结果为：s%f=%lf",n-1,s);
}

main()
{float h,n;
 
 printf("请输入接点个数：");
 scanf("%f",&n);
 printf("\n");
 h=1/(n-1);
 tixing(h,n);
 printf("\n");
 simpson(h,n);
 printf("\n");
}