scilab 求微分_数值计算:微分方程的数值解法的C++程序

微分方程的各种算法的C++实现，是学数值计算方法的实验，有欧拉法，二阶龙格库塔法和四阶龙格库塔法。 程序使用了我自己写的库(已上传)。

Linux下编译方法：g++ df.cpp -o df -lm

附图为把程序对y'=cos(x)在[0,10]区间内插值运行结果复制到scilab里画的图，实线是原始函数图像，点实线为欧拉法的解，*号为RK2的解，虚线为RK4结果。

#include #include #include #include #include using namespace std;

//微分方程 df = f(x,y);

double df(double x, double y)

{

return cos(x); // df/dx = cos(x); 解为f(x) = sin(x) + C

}

//欧拉法解微分方程，分块数量为n, 解的区间为[a,b], 解向量为(x,y),方程初值为(x0,y0)

bool Eluer(vector& x, vector& y, const int & n,

const double & a, const double & b,

const double & x0, const double & y0)

{

double h = (b-a)/n;

x.clear();

y.clear();

x.resize(n+1);

y.resize(n+1);

x[0] = x0;

y[