0
Этот код в порядке?Какой будет код python для метода runge kutta second?
def rKN(x, fx, n, hs):
k1 = []
k2 = []
k3 = []
k4 = []
xk = []
for i in range(n):
k1.append(fx[i](x)*hs)
for i in range(n):
xk.append(x[i] + k1[i]*0.5)
for i in range(n):
k2.append(fx[i](xk)*hs)
for i in range(n):
xk[i] = x[i] + k2[i]*0.5
for i in range(n):
k3.append(fx[i](xk)*hs)
for i in range(n):
xk[i] = x[i] + k3[i]
for i in range(n):
k4.append(fx[i](xk)*hs)
for i in range(n):
x[i] = x[i] + (k1[i] + 2*(k2[i] + k3[i]) + k4[i])/6
return x
ли работа? Откуда вы знаете? –
См. Https://stackoverflow.com/questions/35071393/runge-kutta-code-not-converging-with-builtin-method (или любую из связанных ссылок на боковой панели) для более компактного кода RK4. – LutzL