研究生數(shù)值分析課后題上機編程部分答案-資料下載頁
2025-06-23 04:46本頁面
【正文】 for(j=1。(j*)=。j++) for(i=1。i=3。i++) /*本解法未使用0行，故從1開始取值*/{if(i=1) {k1=1。k2=1。k3=1。k4=1。 y[1][j]=y[1][j1]+h/6*(k1+2*k2+2*k3+k4)。} /*計算第一行數(shù)值*/ if(i=2) {k1=y[3][j1]。 k2=y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])。k3=y[3][j1]+*h*(10001000*(y[2][j1]+*h*y[3][j1])100*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])))。k4=y[3][j1]+h*(10001000*(y[2][j1]+*h*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])))100*(y[3][j1]+*h*(10001000*(y[2][j1]+*h*y[3][j1])100*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])))))。 y[2][j]=y[2][j1]+h/6*(k1+2*k2+2*k3+k4)。} /*計算第二行數(shù)值*/ if(i=3) {k1=10001000*y[2][j1]100*y[3][j1]。 k2=10001000*(y[2][j1]+*h*y[3][j1])100*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1]))。k3=10001000*(y[2][j1]+*h*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])))100*(y[3][j1]+*h*(10001000*(y[2][j1]+*h*y[3][j1])100*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1]))))。k4=10001000*(y[2][j1]+h*(y[3][j1]+*h*(10001000*(y[2][j1]+*h*y[3][j1])100*(y[3][j1]+*h*(10001000*y[2][j1]100*y[3][j1])))))100*(y[3][j1]+h*k3)。y[3][j]=y[3][j1]+h/6*(k1+2*k2+2*k3+k4)。} /*計算第三行數(shù)值*/} printf(y1(),y2(),y3()值分別為：%f,%f,%f\n,y[1][50],y[2][50],y[3][50])。 printf(y1(),y2(),y3()值分別為：%f,%f,%f\n,y[1][90],y[2][90],y[3][90])。 printf(y1(),y2(),y3()值分別為：%f,%f,%f\n,y[1][170],y[2][170],y[3][170])。 printf(y1(),y2(),y3()值分別為：%f,%f,%f\n,y[1][200],y[2][200],y[3][200])。 }（誤差分析、上機出現(xiàn)情況等）（1）本題的編寫過程相對前一題較為順利，過程中只出現(xiàn)一個較大問題，原程序中我在for循環(huán)體內(nèi)使用的是if(){} else if(){} else (){} 的結(jié)構(gòu)，無錯，但運行不出結(jié)果，改用if (){} if (){} if (){} 后出現(xiàn)結(jié)果，至今仍未找清楚原因。（2）另一個小問題是，在計算第三行數(shù)據(jù)時，由于要使用前一列二行的數(shù)據(jù)，剛開始直接將復(fù)制進來，稍后發(fā)現(xiàn)這時的K值是二行中的數(shù)據(jù)，放在第三行中將出錯，由于發(fā)現(xiàn)及時，未對后面程序編寫造成影響。（故后續(xù)將K值都由具體的算式替代，避免出錯）14
化學(xué)相關(guān)推薦
鄂ICP備17016276號-1