用龙格库塔法解微分方程组.doc

《用龙格库塔法解微分方程组.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#includeusingnamespacestd;intmain(){doubley1[100];doubley2[100];doubley3[100];intt=0;doubleh=0.1;doublek[3][4];intn;y1[0]=0;y2[0]=0;y3[0]=0;cout<

2、0]=-800*y1[n-1]-80*y2[n-1]-24*y3[n-1]+1;k[0][1]=y2[n-1]+h/2*k[1][0];k[1][1]=y3[n-1]+h/2*k[2][0];k[2][1]=-800*(y1[n-1]+h/2*k[0][0])-80*(y2[n-1]+h/2*k[1][0])-24*(y3[n-1]+h/2*k[2][0])+1;k[0][2]=y2[n-1]+h/2*k[1][1];k[1][2]=y3[n-1]+h/2*k[2][1];k[2][2]=-800

3、*(y1[n-1]+h/2*k[0][1])-80*(y2[n-1]+h/2*k[1][1])-24*(y3[n-1]+h/2*k[2][1])+1;k[0][3]=y2[n-1]+h*k[1][2];k[1][3]=y3[n-1]+h*k[2][2];k[2][3]=-800*(y1[n-1]+h*k[0][2])-80*(y2[n-1]+h*k[1][2])-24*(y3[n-1]+h*k[2][2])+1;y1[n]=y1[n-1]+h/6*(k[0][0]+2*k[0][1]+2*k[0]

4、[2]+k[0][3]);y2[n]=y2[n-1]+h/6*(k[1][0]+2*k[1][1]+2*k[1][2]+k[1][3]);y3[n]=y3[n-1]+h/6*(k[2][0]+2*k[2][1]+2*k[2][2]+k[2][3]);cout<