龙格-库塔法解常微分方程（c++）

用龙格库塔法计算

 

#include <iostream>

#include<iomanip>

#include <cmath>

using namespace std;

int main()

{

      double a = 1, b = 3;      //a,b表示[a,b]求解区间