fongostevДа! Я про эту ссылку говорил. Вот мой вариант кода
from matplotlib import pyplot as plt import numpy as np from scipy.integrate import odeint G=6.67*10**-8 M=2*10**33. def motion(x_vec, t): x, y = x_vec # tuple unpacking return [y, G*M/y**2.] y0 = [10**13., 10**13.] t_output = np.arange(0, 10**7, 10**3) y_result = odeint(motion, y0, t_output) plt.plot(y_result[:,0],y_result[:,1]) plt.show()
OlmerТам идет распаковка набора аргументов. Грубо говоря, там происходит следующее:
a, b, c = (a, b, c)
fongostevВот здесь? “x, y = x_vec ” Немного не догоняю…
x, y = x_vec
OlmerДа, именно. Просто в интерпретаторе введи
t = (1,2) a, b = t
fongostevНо x_vec, что такое? Что это за вектор?
Olmerэто твой набор независимых переменных, в твоем случае х и х с точкой.
fongostevА ты можешь записать, как будет выглядеть тело функции, или у меня правильно написано?
return [y, G*M/x**2]
fongostevОн рисует, что-то похожее. Но если поменять знак на минут перед G*M, то получается страшное решение.
Ты просто неаккуратно расписал систему. Замени строку и попробуй. Должно быть так: