改进的欧拉法计算常微分方程初值问题

改进的欧拉法计算常微分方程初值问题

标签： 计算方法实验

---

/*
    本实验用改进的欧拉公式计算常微分方程初值问题
    1) y' = y / x - y * y / x, 1 <= x <= 1.5
    2) y(1) = 0.5
    取步长h = 0.1(即n = 5).
*/
#include <stdio.h>
#include <math.h>

double f(double x, double y){  //微分方程右端函数
    return (y / x - y * y / x);
}

void Euler(double x0, double xn, double y0, int n){  //改进的欧拉公式