改进欧拉公式求解常微分方程（c++）

#include<iostream>

#include<iomanip>

using namespace std;

int main()

{

      

      double x,y,h,temp,f;     

      x = 0;      //对x赋初值

      y = 1;      //对y赋初值

      h = 0.1;      //步长设置为0.1

      

      cout << setiosflags(ios::left);

      cout << setw(20) << "y的计算值";

      cout << setw(20) << "y的理论值";

      cout << setw(20) << "x的值";

      cout << setw(20) << "误差" << endl;

      cout << setw(20) << y;

      cout << setw(20) << y;

      cout << setw(20) << x;

      cout << setw(20) << 0 << endl;

      

      for (int i = 0; i < 10; i++)

      {

            temp=y;

            f = y - (2 * x / y);      //保存未改变的f（x,y）的值

            y = y + h*(y - (2 * x / y));      //得到预估值

            x += h;      //迭代得到最新的x值

            y=temp + h / 2 *(f + (y - (2 * x / y) ) );

            cout << setw(20) << y;      //输出y的新值

            cout << setw(20) << sqrtf(1 + 2 * x);      //计算y的理论值

            cout << setw(20) << x;      //输出x的新值

            cout << setw(20) << y - sqrtf(1 + 2 * x) << endl;      //计算误差

      }

      return 0;

}