常微分方程数值解法（2）

 

#include <stdio.h>

#include <cmath>

#define N 20//离散点个数取20，对应步长0.1

#define low -1.0//定义域x区间下限

#define high  1.0//定义域x区间上限

#define h (high-low)/N//步长

double q(double x)//定义函数q(x)

{

      return 1+pow(x,2);

}

double gamma(double x)//定义函数gamma(x)

{

      return 0;

}

int main()

{

      double x[N+1];//自变量

      double y[N+1];//y

<