matlab做数值分析Romberg积分预测近似值

%fprintf('使用Romberg积分，对于积分[048]上f(x)=sqrt(1+(cosx)^2)计算近似值\\n') %%

%确定R11R21R31R41R51，并用这些近似值预言积分值 %确定R22R33R44R55的值，并修改你的预言 %确定R61R62R63R64R65R66并修改你的预言 %确定R77R88R99R1010并作出最后的预言 %%

a=0;b=48; h=zeros(101); h(1)=b-a; n=10;

R=zeros(1010);

R(11)=h(1)/2*(sqrt(1+(cos(a))^2)+sqrt(1+(cos(b))^2)); h(2)=h(1)/2; for k=2:n s=0;

h(k)=(b-a)/2^(k-1); for i=1:2^(k-2)

x=a+(2*i-1)*h(k);

s=s+sqrt(1+(cos(x))^2); end

R(k1)=(R(k-11)+h(k-1)*s)/2; end

for j=2:n for k=j:n

R(kj)=R(kj-1)+(R(kj-1)-R(k-1j-1))/(4^(j-1)-1); end end

sprintf('对于积分[048]上f(x)=sqrt(1+(cosx)^2)计算积分预测值分别为：\\n\a\\\b\\\c\\\d\\\n%f\%f\%f\%f'R(51)R(55)R(66)R(1010))

更多相关数值分析实用代码：不动点迭代，Newton插值，Nevilles插值，矩阵LDL'分解，

LL'，Choleski分解，求解线性方程组，hermit插值，样条插值，scant迭代，false position法。有意者私聊