高斯-勒让德积公式.docVIP

实 习 论 文 题 目 高斯勒让德积分公式 专 业 信息与计算科学 班 级 计算092 学 号 3090811065 学 生 周吉瑞 指导教师 秦新强 2011 年 高斯勒让德积分公式 专 业： 信息与计算科学 学 生： 周吉瑞 指导老师： 秦新强 摘要 关于数值积分公式，除了用误差来分析其精度以外，还可以用代数精度来判断其代数精度的高低，已知n+1点Newton-Cotes型积分公式，当n为奇数时，其代数精度为n，当n为偶数时，其代数精度达到n+1。 n+1点的Newton-Cotes型积分公式属于插值积分型积分公式，一般地，若对随机选取的n+1个节点作插值型积分公式也仅有n次代数精度，但是，如果求积节点选取适当，就有可能提高数值积分的代数精度，高斯型积分公式就可以实现这一目标。 关 键 词：数值积分，代数精度，高斯型积分公式 一、目的意义 构造Gaoss型求积公式除需要求出正交多项式外，还需要求出正交多项式的零点和求积系数，当时，这些工作均很困难，因此给出高斯-勒让德积分公式的零点和系数。 二、公式 高斯-勒让德积分公式：； 零点及系数： n x A 1 0 2 2 0.5773502692 1 3 0.77459666920 0 0.5555555556 0.8888888889 4 0.8611363116 0.3399810436 0.3478548451 0.6521451549 5 0.9061798459 0.53845931010 0 0.2369268851 0.4786286705 0.5688888889 6 0.9324695142 0.6612093865 0.2386191816 0.1713244924 0.3607615730 0.4679139346 7 0.9491079123 0.7415311856 0.4058451514 0 0.1294849662 0.2797053915 0.3818300505 0.4179591834 8 0.9602898565 0.7966664774 0.5255324099 0.1834346425 0.1012285363 0.2223810345 0.3137066459 0.3626837834 三、算法流程 Step1:输入所用的点数n； Step2:对i=1，2，···，n循环执行步3； Step3:I= I+ ; Step4:输出I；结束。 四、算法程序 #includestdio.h #includemath.h double Leg(double x) { double z; z=8/(4+(1+x)*(1+x)); return z; } void main() { double x[9],A[9],I=0; int i,n; printf(请输入点数n:); scanf(%d,n); switch(n) { case 1: x[1]=0,A[1]=2;break; case 2: x[1]=0.5773502692,x[2]=-0.5773502692,A[1]=A[2]=1;break; case 3: x[1]=0.77459666920,x[2]=0,x[3]=-0.77459666920,A[1]=A[3]=0.5555555556,A[2]=0.8888888889;break; case 4: x[1]=0.8611363116,x[2]=0.3399810436,x[3]=-0.8611363116,x[4]=-0.3399810436; A[1]=A[3]=0.3478548451,A[2]=A[4]=0.6521451549;break; case 5: x[1]=0.9061798459,x[2]=0.53845931010,x[3]=0,x[4]=-0.9061798459,x[5]=-0.53845931010; A[1]=A[4]=0.2369268851,A[3]=0.5688888889,A[2]=A[5]=0.4786286705;break; case 6: x[1]=0.9324695142,