數(shù)值分析上機(jī)報(bào)告-資料下載頁
2025-08-04 00:46本頁面
【正文】 ,x)*pow(x,)*(5*x+7)*sin(x*x)。 return(y)。}三、計(jì)算結(jié)果如下：The result: y=四、問題討論：本程序可通過一個(gè)一維數(shù)組將進(jìn)行存儲(chǔ)，但考慮到對(duì)程序的可讀性，故將存儲(chǔ)在一個(gè)二維數(shù)組里。第六題：用定步長四階RungeKutta法求解h= 5，打印。一、分析如下：二、 程序如下：includeincludedefine N 3float y[N]={0,0,0},k[N][N+1]。float yd[N+1][N],a[N+1]={,}。/*以上是預(yù)處理和定義全局變量*/main(){int i,l,j=0,n=1。 float h=。do{for(l=0。lN。l++)for(i=0。i=N。i++)/*計(jì)算每一個(gè)k值*/k[l][i]=0。k[0][0]=h。k[1][0]=h*y[2]。k[2][0]=h*(10001000*y[1]100*y[2])。k[0][1]=h。k[1][1]=h*(y[2]+(1/)*k[0][0])。k[2][1]=h*(10001000*(y[1]+(1/)*k[1][0])100*(y[2]+(1/)*k[1][0]))。k[0][2]=h。k[1][2]=h*(y[2]+(1/)*k[1][0]+k[1][1])。k[2][2]=h*(10001000*(y[1]+(1/)*k[2][0]+k[2][1])100*(y[2]+(1/)*k[2][0]+k[2][1]))。k[0][3]=h。k[1][3]=h*(y[2]+k[1][0]k[1][1]+k[1][2])。k[2][3]=h*(10001000*(y[1]+k[2][0]k[2][1]+k[2][2])100*(y[2]+k[2][0]k[2][1]+k[2][2]))。for(i=0。iN。i++) y[i]=y[i]+(k[i][0]+3*k[i][1]+3*k[i][2]+k[i][3])/8。if(n==50||n==90||n==170||n==200)/*得到所需的數(shù)據(jù)*/{for(i=0。iN。i++) yd[j][i]=y[i]。 j++。}n++。}while(n=200)。printf(The result of sixth question is:\n)。for(j=0。jN+1。j++) {for(i=0。iN。i++) printf(y%d(%)=%.7f,i+1,a[j],yd[j][i])。printf(\b\n)。}}三、 計(jì)算結(jié)果如下：y1()= y2()= y3()=y1()= y2()= y3()=y1()= y2()= y3()=y1()= y2()= y3()=y1()= y2()= y3()=四、問題討論：在設(shè)計(jì)本程序初期，曾將再次計(jì)算時(shí)的初始值賦值給忽略了，造成只有第一個(gè)計(jì)算結(jié)果是正確的，而不能繼續(xù)進(jìn)行其它數(shù)據(jù)的計(jì)算，后經(jīng)修改，可實(shí)現(xiàn)對(duì)其它數(shù)據(jù)的連續(xù)計(jì)算。在設(shè)計(jì)過程中，由于將浮點(diǎn)數(shù)賦值給了整數(shù)，造成了計(jì)算誤差，后經(jīng)修改，使計(jì)算正確進(jìn)行。本程序可以通過修改主程序所調(diào)用的函數(shù)中的表達(dá)式來實(shí)現(xiàn)對(duì)其它函數(shù)的任意初值條件求微分計(jì)算。寧可累死在路上，也不能閑死在家里！寧可去碰壁，也不能面壁。是狼就要練好牙，是羊就要練好腿。什么是奮斗？奮斗就是每天很難，可一年一年卻越來越容易。不奮斗就是每天都很容易，可一年一年越來越難。能干的人，不在情緒上計(jì)較，只在做事上認(rèn)真；無能的人！不在做事上認(rèn)真，只在情緒上計(jì)較。拼一個(gè)春夏秋冬！贏一個(gè)無悔人生！早安！—————獻(xiàn)給所有努力的人.學(xué)習(xí)好
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1