幸卜生.cpp

这是一个数值分析算法的源程序,是用C来编写的,绝对是源码,而且也比较有用,因为只要是理工科的大学生大多要学数值分析,所以此程序有用

#include <stdio.h>
#include <math.h>
#define m 4
#define f(x) sqrt(x)

void xing()
{   int i;
	float a=0.00,b=9.00, x[5],h=(a+b)/m;
double t1[5],t;
for(i=0;i<5;i++)
{
	x[i]=a+i*h;
t1[i]=f(x[i]);
printf("%12.8f\n",t1[i]);
}
t=h*(f(x[0])+f(x[4])+4*(f(x[1])+f(x[3]))+2*f(x[2]))/3;
printf("%12.9f\n",t);
}
void tixing()
{ int i;
	float a=0.00,b=9.00, x[5],h=(a+b)/m;
double t1[5],t;
for(i=0;i<5;i++)
{
	x[i]=a+i*h;
t1[i]=f(x[i]);
printf("%12.8f\n",t1[i]);
}
t=h*(f(x[0])/2+f(x[4])/2+f(x[1])+f(x[3])+f(x[2]));
printf("%12.9f\n",t);
}
void main()
{ int a;
printf("\n");
printf("---------------------选择菜单----------------------\n");
printf("*******  1  ********复化梯形公式*****\n");
printf("*******  2  *******复化辛卜生公式****\n");
printf("input your choice:\n");

	scanf("%d",&a);
	switch(a)
	{case 1:  tixing(); break;
     case 2:  xing(); break;}
printf("input your choice again:\n");
		scanf("%d",&a);
	switch(a)
	{case 1:  tixing();break;
	case 2:  xing();break; }
printf("用积分方法知真实值为18\n");
printf("通过两者的对比易知复化辛卜生公式较为精确\n");
}