改进后的欧拉方法.cpp

来自「数值分析（计算方法）：欧拉方法、高斯－赛德尔迭代法、拉格朗日及分段线性插值、龙贝」· C++ 代码 · 共 24 行

#include "iostream"
#include "iomanip"
using namespace std;

void main()
{double x0=0,y0=1,x1,y1,yp,yc,h;
int i,n;
cout<<"输入n:";
cin>>n;
cout<<"输入h:";
cin>>h;
i=1;
while(i<=n)
{x1=x0+h;
yp=y0+h*(y0-2*x0/y0);
yc=y0+h*(yp-2*x1/yp);
y1=(yp+yc)/2;
cout<<setprecision(8)<<x1<<' ';
cout<<setprecision(8)<<y1<<endl;
x0=x1;
y0=y1;
i++;
}
}