文档介绍：实验三实验报告
贵州师范大学数学与计算机科学学院学生实验报告
课程名称： 数值分析 班级：数学（2）班 实验日期： 2013 年 10 月 11 日
学 n
2：对 k=1,2,...,n 执行 x=x+h/2;Sn=Sn+4f(x);x=x+h/2;Sn=Sn+2f(x).
h
3：Sn= Sn
6
4：输出 Sn.
题 2
（１）、算法描述：用变步长梯形法;
实验步骤：求积分 ，允许误差为ε 。
被积函数 f（x）,a,b,ε .
复合梯形积分值 T2n.
h<=b-a.
T1=h/2(f(a)+f(b)).
反复执行步 4—步 10。
S<=0;x<=a+h/2.
反复执行步 6—步 7.
S<=S+f（x）;x<=x+h.
若 x≥b,则退出本层循环。
T2<=T1/2+h/2S
e<=|T2-T1|;h<=h/2;T1<=T2.
若 e≤ε 则退出循环。
T2n<=T2.
输出 T2n.
六、调试过程及实验结果、
题 1
（1）调试过程：
a=[,,,,,,,,];
>> h=10;
>> T=(h/2)*(a(1)+a(6)+2*(a(2)+a(3)+a(4)+a(5)))
T =
+003
T=(h/2)*(a(1)+a(9)+2*(a(2)+a(3)+a(4)+a(5)+a(6)+a(7)+a(8)))
T =
+003
实验结果： 50s 时速度为 ;80s 时速度为
（2）调试过程：
x=0:10:80;
>> y=[,,,,,,,,];
>> cx=5:10:75;
>> cy=Lagrange(x,y,5,cx)
2 / 5实验三实验报告
cy =
>> h=10;
>>Sn=(h/6)*(y(1)+4*(cy(1)+cy(2)+cy(3)+cy(4)+cy(5))+2*(y(2)+y(3)+y(4)+y(5))+y(6))
Sn =
+003
>>Sn=(h/6)*(