作业实验5.docxVIP

实验5.1 复化求积公式计算定积分(1)用复化梯形公式、复化辛普森公式利用复化高斯-勒让德公式计算，误差限在1/2e-7，计算各个方法的应用步长。以第一题为例输入n1=trapen(f,a,b,wucha)n2=simpn(f,a,b,wucha)n3=Gausslegn(f,a,b,wucha);%%n1为复化梯形公式需要的步长%%n2为复化辛普森公式需要的步长%%n3为复化高斯勒让德公式需要的步长。输出n1=9310n2=85n3=29 (2)分别用复化梯形公式、复化辛普森公式利用复化高斯-勒让德公式计算输入Tn=trape(f,a,b,n1)Sn=simp(f,a,b,n2)Gn=Gaussleg(f,a,b,n3)输出Tn=0.693147187950964Sn=0.693147183740232Gn超出计算机运算范围(3)将计算结果与精确解比较，并比较各种算法的计算量输入T1=log(2)-TnS1=log(2)-Sn输出T1=-7.3910e-09S1=-3.1803e-09梯形公式计算时间为5.931496s辛普森公式计算时间为0.470143s同理可得其他三题的步长估计、积分计算值和计算量整理计算结果得表如下：步长统计复化梯形复化辛普森复化高斯勒让德题19310852932457144分值统计复化梯形复化辛普森复化高斯勒让德题