数值分析实验数值积分和数值微分.doc

黑河学院数学系
实验报告
课程名称: 数值分析
实验项目: 数值积分和数值微分
实验室名称: 数学系实验室
实验台号: 02
指导教师: 王国贤
学生学号: 2011562003
学生姓名: 柳雪梦
实验日期:__2013___年__11___月___28__日
:
要求学生掌握牛顿-科特斯求积公式、复化辛甫生求积法、高斯型求积分方法,初步掌握用Matlab、C语言等软件求解问题。
:
应用自适应复化辛甫生求积法实习、应用软件编写程序,求给定问题的数值积分,并与真实值作比较,进行验证,同时达到熟练掌握和应用的程度。
:
用复化的辛甫生公式计算积分 M=4
给定下列表格值:
X
50
55
60
65
Y




复化辛甫生公式
利用C语言编程,程序源代码如下:
#include <>
#include <>
float f(float x)
{
return(sqrt(x));
}
void main()
{
float a,b,h,x[10],sum=0,sum1=0,sum2=0,sum3=0;
int i,m;
printf("input a:");
scanf("%f",&a);
printf("input b:");
scanf("%f",&b);
printf("input m:");
scanf("%d",&m);
h=(b-a)/m;x[0]=a;x[m]=b;
for(i=1;i<m;i++)
x[i]=a+h*i;
sum1=f(x[0])+f(x[m]);
for(i=1;i<=m-1;i++,i++)
sum2=sum2+f(x[i]);
for(i=2;i<=m-2;i++,i++)
sum3=sum3+f(x[i]);
sum=(sum1+4*sum2+2*sum3)*h/3;
printf("the result is %f",sum);
getch();
}
运行后,得到如下结果:
input a:0
input b:9
input m:4
the result is