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

#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");
}

转载于:https://www.cnblogs.com/neoLin/p/10999688.html

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

一阶欧拉近似matlab,MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程.doc
MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程姓名:樊元君学号:2012200902 日期:2012.11.06 一.实验目的掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题 ...
数值分析C++实现用四阶龙格-库塔(Runge-Kutta)方法求解常微分方程初值问题
问题:用四阶龙格-库塔(Runge-Kutta)方法求解常微分方程初值问题. 算法描述算法的程序框图: 具体算法: (1)读取a,b,n,f (2)计算步长h = (b-a)/n,x0=a,y0=f ...
python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法
我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...
四阶龙格库塔（Runge-Kutta）求解微分方程-多种编程语言
前期是分享了matlab下面实现四阶龙格库塔(Runge-Kutta)求解微分方程,这期分享一下C++.C.Java.Python下面的四阶龙格库塔(Runge-Kutta)求解微分方程. 前文传送门 ...
通过swagger下载的文件乱码解决方法，求解
通过swagger下载的文件乱码解决方法,求解参考文章: (1)通过swagger下载的文件乱码解决方法,求解 (2)https://www.cnblogs.com/shuiqian/p/10568 ...
一阶常微分方程的数值解法（二阶显式、隐式 Adams 公式及 Milne 方法）
一阶常微分方程的数值解法这里我们介绍二阶显式.隐式 Adams 公式及 Milne 方法求解方程. 题目: 对初值问题 u ′ = u − t 2 , 0 ≤ t ≤ 1 , u ( 0 ) = 0 ...
经典龙格-库塔法（四阶龙格-库塔法）求解求一阶常微分方程相应的特解的Python程序
基本原理例题代码 #四阶龙格-库塔法 #求一阶常微分方程,相应的特解 #x变量的区间 a = 0 b = 1 #已知条件 X = [0] Y = [1] h = 0.2 #设置步长 n = (b- ...
蒙特卡罗方法近似求解定积分-Go语言
蒙特卡罗方法近似求解定积分-Go语言蒙特卡罗方法蒙特卡洛方法是一种用大数统计频率来模拟概率的方法.也就是说,当样本越多,越有可能获得及接近最优解. 定积分的几何意义定积分的几何意义实际上就是求函 ...
Java的开方和次方Math.sqrt()方法和Math.pow()方法以及求解有多少种连续正整数之和为N的算法详解
1. Java的开方和次方运算 java.lang.Math.sqrt(double a) 返回一个数字a的正平方根,返回结果是double型 java.lang.Math.pow(double a, ...

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

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

最新文章

热门文章