下面是使用Python求解在区间[0,0.51]上初值问题的数值解的示例代码,望采纳:
import numpy as np
def diff_eq(x, y):
return x - 2 * y
def solve_ivp(x0, y0, x_end, h):
# 计算步数
n = int((x_end - x0) / h)
# 初始化 x 和 y 数组
x = np.zeros(n + 1)
y = np.zeros(n + 1)
x[0] = x0
y[0] = y0
# 使用 Euler 方法求解
for i in range(n):
x[i + 1] = x[i] + h
y[i + 1] = y[i] + h * diff_eq(x[i], y[i])
return x, y
# 求解初值问题
x, y = solve_ivp(0, 1, 0.51, 0.1)
# 输出结果
for i in range(len(x)):
print(f'x = {x[i]:.2f}, y = {y[i]:.2f}')
你熟悉什么编程语言呢?