使用C++，用四阶Runge-Kutta的方法来求解一阶常微分方程

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了使用C++，用四阶Runge-Kutta的方法来求解一阶常微分方程，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含679字，纯文字阅读大概需要1分钟。

内容图文

#include <iostream>
using namespace std;

/*
dy/dx = y - 2x/y, 0< x <= 1
步长h = 0.2
*/
const double h = 0.2;

//待求解的一节常微分方程，dy/dx = f_x
double f_x(double x, double y)
{
    return (y - 2 * x / y);
}

int main() {
    double k[5];//其实k[0]没有用到，为了看起来更简洁
    double y[6];
    y[0] = 1;
    double x[6] = { 0.0,0.2, 0.4, 0.6, 0.8, 1.0 };
    //计算出的是y[i+1]所以到 i = 5就行了。
    for (int i = 0; i < 5; ++i) {
        k[1] = f_x(x[i], y[i]);
        k[2] = f_x(x[i] + h / 2, y[i] + h / 2 * k[1]);
        k[3] = f_x(x[i] + h / 2, y[i] + h / 2 * k[2]);
        k[4] = f_x(x[i] + h, y[i] + h * k[3]);
        y[i + 1] = y[i] + h / 6*(k[1] + 2 * k[2] + 2 * k[3] + k[4]);
    }
    for (int i = 0; i < 6; ++i) {
        cout << y[i] << endl;
    }
    system("pause");
}

内容总结

以上是互联网集市为您收集整理的使用C++，用四阶Runge-Kutta的方法来求解一阶常微分方程全部内容，希望文章能够帮你解决使用C++，用四阶Runge-Kutta的方法来求解一阶常微分方程所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端