peyo написа:In [49]: print(latex(integrate(F_0*(t/(m_0-k*t)) - g*t,(t,0,t))))
$h(t)=- \frac{F_{0} t}{k} + \frac{F_{0} m_{0} \log{\left(- m_{0} \right)}}{k^{2}} - \frac{F_{0} m_{0} \log{\left(k t - m_{0} \right)}}{k^{2}} - \frac{g t^{2}}{2}$
Странно наистина, но това трябва да дава реални положителни числа от 0 до някой момент.
Да проверим дали наистина това дава очаквана графика!?
Нека ракетата да тежи 100 кг. $F_0$ трябва да е по-голяма от 100*9.81 = 980 N, да вземем 1050 N, а константата $k$ да е 1, така ракетата ще се самоизяде за 100 секунди. Да проверим графиката от 0 до 200 секунди:
In [72]: h = integrate(F_0*(t/(m_0-k*t)) - g*t,(t,0,t))
In [77]: print(latex(h.subs({m_0: 100, F_0: 1050, k: 1, g: 9.81})))
$- 4.905 t^{2} - 1050 t - 105000 \log{\left(t - 100 \right)} + 105000 \log{\left(100 \right)} + 105000 i \pi$
In [73]: plot(h.subs({m_0: 100, F_0: 1050, k: 1, g: 9.81}), (t,0,200))
Out[73]: <sympy.plotting.plot.Plot at 0x18036ac87d0>

- Figure_rocket.png (19.58 KiB) Прегледано 1671 пъти
Изглежда сякаш много добре. Първо тази формула дава реални числа въпреки, че има имагинерен компонент. В началото височината расте много бавно, но с намалаването на тежестта расте много бързо. След 100 секунди няма графика, защото масата на ракетата е паднала под нулата. В 100 секунди графиката върви към безкрайност, защото масата клони към 0.