數(shù)值分析上機(jī)實(shí)驗(yàn)報(bào)告-資料下載頁(yè)
2025-08-03 23:50本頁(yè)面
【正文】 j]=A[j][i]。} for(k=0。k9。k++) {y[k]=0。 for(j=0。j9。j++) y[k]+=A[k][j]*u[j]。 y[k]/=a。} kr=0。 for(k=0。k9。k++) kr+=y[k]*u[k]。 kr/=2*a。 for(k=0。k9。k++)q[k]=y[k]kr*u[k]。 for(k=0。k9。k++) {for(j=0。j9。j++) A[k][j]=u[k]*q[j]+u[j]*q[k]。} }}void expunction(double A[9][9],double b[9],double x[9]){int i,j,k。 double B[9][10]。 double z[3]。 double t1=0,t2=0,t3=0。 for(i=0。i8。i++) {if(A[i+1][i]A[i][i]) {for(j=i,k=0。ji+3。j++,k++) z[k]=A[i][j]。A[i][j]=A[i+1][j]。A[i+1][j]=z[k]。 t1=b[i]。b[i]=b[i+1]。b[i+1]=t1。} t2=A[i+1][i]。 for(j=i。ji+3。j++) A[i+1][j]=A[i+1][j]A[i][j]*t2/A[i][i]。 b[i+1]=b[i+1]b[i]*t2/A[i][i]。} x[8]=b[8]/A[8][8]。 for(i=7。i=0。i) {for(j=i+1。j9。j++) t3=t3+A[i][j]*x[j]。 x[i]=(b[i]t3)/A[i][i]。 t3=0。}} 運(yùn)行結(jié)果 MATLAB上機(jī)程序unction [x]=mgauss2(A,b,flag)if nargin3,flag=0。endn=length(b)。for k=1:(n1) [ap,p]=max(abs(A(k:n,k)))。 p=p+k1。 if pk A([k p],:)=A([p k],:)。 b([k p],:)=b([p k],:)。 end m=A(k+1:n,k)/A(k,k)。 A(k+1:n,k+1:n)=A(k+1:n,k+1:n)m*A(k,k+1:n)。 b(k+1:n)=b(k+1:n)m*b(k)。 A(k+1:n,k)=zeros(nk,1)。 if flag~=0,Ab=[A,b],endendx=zeros(n,1)。x(n)=b(n)/A(n,n)。for k=n1:1:1 x(k)=(b(k)A(k,k+1:n)*x(k+1:n))/A(k,k)。endformat longA=[,,。,,。,,。,,。,,。,,。,,。,,。,,]。b=[,,]39。x=mgauss2(A,b)。x=x39。：。，數(shù)組中的下標(biāo)從零開始。，編程時(shí)必須注意。1. 若不給自己設(shè)限，則人生中就沒(méi)有限制你發(fā)揮的藩籬。2. 若不是心寬似海，哪有人生風(fēng)平浪靜。在紛雜的塵世里，為自己留下一片純靜的心靈空間，不管是潮起潮落，也不管是陰晴圓缺，你都可以免去浮躁，義無(wú)反顧，勇往直前，輕松自如地走好人生路上的每一步3. 花一些時(shí)間，總會(huì)看清一些事。用一些事情，總會(huì)看清一些人。有時(shí)候覺(jué)得自己像個(gè)神經(jīng)病。既糾結(jié)了自己，又打擾了別人。努力過(guò)后，才知道許多事情，堅(jiān)持堅(jiān)持，就過(guò)來(lái)了。4. 歲月是無(wú)情的，假如你丟給它的是一片空白，它還給你的也是一片空白。歲月是有情的，假如你奉獻(xiàn)給她的是一些色彩，它奉獻(xiàn)給你的也是一些色彩。你必須努力，當(dāng)有一天驀然回首時(shí)，你的回憶里才會(huì)多一些色彩斑斕，少一些蒼白無(wú)力。只有你自己才能把歲月描畫成一幅難以忘懷的人生畫卷。學(xué)習(xí)參考
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1