预校算法解微分方程

1. #include <iostream.h>
2. #include <math.h>
3. double function(double x,double y)
4. {
5.     double z;
6.     z=-2*x*y*y;
7.     return z;
8. }
9. main()
10. {
11.     double h,a,b,t,t1;
12.     double n;
13.     cout<<"请输入步长h:"<<endl;
14.     cin>>h;
15.     cout<<"请输入取值区间左端点:"<<endl;
16.     cin>>a;
17.     cout<<"请输入取值区间右端点:"<<endl;
18.     cin>>b;
19.     n=(b-a)/h;
22.     
23.     double *x = new double[n+2];// 构造动态数组
24.     double *y1 = new double[n+2];
25.     double *y = new double[n+2];
26.     cout<<"请输入y0 的值"<<endl;
27.     cin>>y[0];
28.     for (int i=0;i<=n;i++)
29.     {
30.         x[i]=a+i*h;
31.     }
32.     for( int k=0;k<=n;k++)
33.     {
34.         t=function(x[k],y[k]);
35.         y1[k+1]=y[k]+h*t;
36.         t1=function(x[k+1],y1[k+1]);
37.         y[k+1]=y[k]+h/2*(t+t1);
38.     }
40.         for (i=0;i<=n;i++)
41.     {
42.         cout<<"y["<<i<<"]="<<y[i]<<endl;
43.     }
46.     
47.     delete [] x;// 释放数组
48.     delete [] y1;
49.     delete [] y;
50.     return 0;
51. }

 

 .