首先,很抱歉,如果这是一个菜鸟问题,我是编码新手。我正在尝试用Euler方法解决这个方程:
期望我们可以像这样解决它(来自示例):
def ODE_Euler(ic,h,endpoint,f):
n = int((1/0.1) + 1)
x = np.linspace(start=0,stop=1,num=n)
y = np.zeros(n)
y[0] = ic # initial condition
for i in range(10):
y[i+1] = y[i] + h*f(x[i],y[i])
return x,y
def dydx(x,y):
return #the equation
x,y = ODE_Euler(1.0,0.1,1,dydx)
import pandas as pd
result = pd.DataFrame({
"x":x,"y":y
})
result
有人可以帮忙吗?