数值分析各种求积公式的应用(2)

数值分析第二次作业

班级：双控 姓名：王雪萍 学号：201102208

options={'积分下限a','积分上限b' ,'插入点相关的值n'}; topic='seting'; lines=1;

def={'0','1','1000'};

h=inputdlg(options,topic,lines,def); a=eval(h{1});%积分下限 b=eval(h{2});%积分上限 n=eval(h{3});%子区间个数

%********************************************

f='func';%用f来调用被积函数func h=(b-a)/(n); t=0;

for k=1:1:n-1 x=a+h*(k); t=t + feval(f,x);

end

format long

t=h*( feval(f,a) +feval(f,b)+2*t)/2;%s是梯形规则的总计 end

%定义被积函数func function y=func(x)

y=4*7782.5*sqrt(1-(972.5/7782.5)^2*sin(x).^2); end

运行上面的程序得周长为：4.870743851190014e+004。

由上述计算结果可知编程计算与用积分函数计算的结果是一样的。与精确值之间的误差也在允许的范围内。

终上所述，卫星轨道的周长约为48707km.