�� : ��
SymPy

�� (>0), �� (>=0).
�� .
import sympy as sp
import matplotlib as plt

v0, alpha, t, s, h = sp.symbols('v0 alpha t s h')

v_0 = int(input())
alp = int(input())
alp = (alp%360)*sp.pi/180
h_0 = int(input())

eqt = sp.Eq(s/(v0*sp.cos(alpha)), t)
eqh = sp.Eq(h, h_0 + v0*t*sp.sin(alpha) - 0.5*9.8*t**2)

sol = sp.solve(eqh.subs([(h, 0), (v0, v_0), (alpha, alp)]))
ss = sp.solve(eqt.subs([(t, sol[-1]), (v0, v_0), (alpha, alp)]))

exprt = eqh.subs([(t, s/(v0*sp.cos(alpha))), (v0, v_0), (alpha, alp)])
exprt

sp.plot(h_0 + v_0*(s/(v_0*sp.cos(alp)))*sp.sin(alp) - 0.5*9.8*(s/(v_0*sp.cos(alp)))**2, (s, 0, ss[0]), showcase=False)
��