《计算方法》 李晓红等 第八章 常微分方程数值解 例题 代码

#include "stdafx.h"

// 常微分方程数值解

 

// dy/dx = f(x,y), 初值条件 y(x0) = y0

// 解存在唯一性条件:

// f(x,y) 关于y 在区域D={(x,y)| a<=x<=b, y belong to R}上

// 满足利普希茨条件 |f(x,y1)-f(x,y2)| < L |y1-y2|, L>0

 

static double fval(double x, double y)

{

    return ((x*y)/2);

}

 

void EulerMethod()

{

    printf("欧拉法:\n");

    int n=10;

    double a=0,b=1;

    double