欧拉法求解常微分方程（c++）

#include<iostream>

#include<iomanip>

using namespace std;

int main()

{

double x, y, h;      //，x为对应的每一步x的值,其中y为对应的每一步y的值

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++)

{

y = y + h*(y - (2 * x / y));      //迭代得到y的新值

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

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

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

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

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

}

return 0;

}

运行结果：