2 qq 23935107 qq_23935107 于 2016.02.26 12:02 提问

C++龙格库塔法步长问题

关于用龙哥库塔法求微分方程的问题:
在程序里每次改变步长h,求得结果都会不一样比如取步长为0.2和0.5,计算结果就不一样,不知道是什么原因,请大神帮帮我,谢谢,我把那个微分方程以及代码附上:

微分方程:dy/dx=y*y;0<=x<=1; y(0)=1;
代码:
double f(double x, double y)定义函数f,用来求微分方程
{
return y*y;
}
///////////////////////////
void CQweDlg::OnButton1()
{
// TODO: Add your control notification handler code here
double xx[100][100],yy[100][100],k[100][100];
double h;步长
int i;
char ch1[10];
GetDlgItem(IDC_EDIT1)->GetWindowText(ch1,10);
h=atof(ch1);

   long double T;
   T=1/h;循环次数
   xx[1][0]=0;
   yy[1][0]=1;
 for(i=0; i<=T;i++)
  {
    k[1][1]=f(xx[1][i],yy[1][i]);
        k[1][2]=f(xx[1][i]+h/2,yy[1][i]+h*k[1][1]/2);
        k[1][3]=f(xx[1][i]+h/2,yy[1][i]+h*k[1][2]/2);
        k[1][4]=f(xx[1][i]+h,yy[1][i]+h*k[1][3]);
        yy[1][i+1]=yy[1][i]+h*(k[1][1]+2*k[1][2]+2*k[1][3]+k[1][4])/6;
        xx[1][i+1]=xx[1][i]+h;
        CString str1,str2;
        str1.Format("%.20f",yy[1][1000]); 
        str2.Format("%.2f",h);
        GetDlgItem(IDC_EDIT2)->SetWindowText(str1);
}

}

1个回答

caozhy
caozhy   Ds   Rxr 2016.03.07 05:48
Csdn user default icon
上传中...
上传图片
插入图片
准确详细的回答,更有利于被提问者采纳,从而获得C币。复制、灌水、广告等回答会被删除,是时候展现真正的技术了!