欧拉方法python代码实现

最新推荐文章于 2024-10-08 12:37:10 发布

方天又画了个戟

最新推荐文章于 2024-10-08 12:37:10 发布

阅读量5.8k

点赞数 7

分类专栏：笔记分享文章标签：欧拉公式 python

本文链接：https://blog.csdn.net/The_Shiled/article/details/121367221

版权

如有不对，请批评指正！

1. 显式Euler和隐式Euler

在这里插入图片描述

2. 改进的Euler

在这里插入图片描述

3. 龙格-库塔（Runge-Kutta）

在这里插入图片描述

4. 举个栗子

在这里插入图片描述
代码仅为了说明计算过程，并未进行优化！

import math
import matplotlib.pyplot as plt

y_e = []  # 显式
y_i = []  # 隐式
y_pro = []  # 改进
y_p = []  # 改进
y_rk = []  # 龙格-库塔
y_c = []  # 解析解
x = []
h = 0.1  # 步长
num = 20


def fun(x0, y0):
    return -2 * x0 / y0 + y0


def Euler():
    print('{:>3s}{:>8s}{:>9s}{:>12s}{:>9s}{:>9s}'.format('x', 'y_e', 'y_i', 'y_pro', 'y_rk', 'y_c'))
    for i in range(num):
        # 显示Euler
        y_e