數(shù)值分析上機(jī)實習(xí)報告-資料下載頁
2025-03-23 08:39本頁面
【正文】 D(4*i2,j)。Ah3(i,j)=D(4*i1,j)。Ah4(i,j)=D(4*i,j)。endAxk(i)=xkk(4*i3)。endCh1=[Axk39。 Ah1]Ch2=[Axk39。 Ah2]Ch3=[Axk39。 Ah3]Ch4=[Axk39。 Ah4]Cheng=[Ch1(1,[2 3 4])。Ch2(1,[2 3 4])。 Ch3(1,[2 3 4])。Ch4(1,[2 3 4])選做題一：Romberg算法function [I,step] = Rg(f,a,b,eps)%UNTITLED4 Summary of this function goes here% Detailed explanation goes here if nargin==3 eps=1e4。 end M=1。 tol=10。 k=0。 T=zeros(1,1)。 h=ba。 T(1,1)=(h/2)*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))。 while toleps k=k+1。 h=h/2。 Q=0。 for i=1:M x=a+h*(2*i1)。 Q=Q+subs(sym(f),findsym(sym(f)),x)。 end T(k+1,1)=T(k,1)/2+h*Q。 M=2*M。 for j=1:k T(k+1,j+1)=T(k+1,j)+(T(k+1,j)T(k,j))/(4^j1)。 end tol=abs(T(k+1,j+1)T(k,j))。 end I=T(k+1,k+1)。 step=k。end選做題二：newton算法function [y] = newton(x)%UNTITLED5 Summary of this function goes here% Detailed explanation goes hereformat long。x1=xfunc2_1(x)/func2_1_1(x)。if(abs(x1)) delt=abs(x1x)。else delt=abs((x1x)/x1)。endwhile(delt1e6||abs(func2_1(x))1e6) x=x1。 x1=xfunc2_1(x)/func2_1_1(x)。 if(abs(x1)) delt=abs(x1x)。 else delt=abs((x1x)/x1)。 end if func2_1(x1)==0 break end endy=x。endfunction [y] = func2_1(x)%UNTITLED6 Summary of this function goes here% Detailed explanation goes herey=xlog(x)2。endfunction [y] = func2_1_1(x)%UNTITLED7 Summary of this function goes here% Detailed explanation goes herey=11/x。end
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1