Animacja równania rzutu ukośnego - Python

Mam do napisania program w Pythonie. Treść zadania jest następująca: zaanimuj rozwiązanie równania rzutu ukośnego w polu grawitacyjnym, w obszarze ograniczonym płaszczyznami prostopadłymi do płaszczyzny rzutu: jedną poziomą i dwiema ukośnymi pod kątem 60 stopni do poziomu. Jeśli chodzi o sam wykres wiem jakich pakietów i metod użyć, gorzej z poskładaniem wszystkiego w działającą całość, natomiast w kwestii ograniczenia obszaru płaszczyznami totalnie nie mam pojęcia jak powinno to wyglądać.
Na tę chwilę wymóżdżyłam coś takiego:

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation

x0,y0=0,0
vx0, vy0=10,10
g=9.81
t = 2*vy0/g
#y=odeint (polewekt, y0-punkt pocz trajektorii, t)

def polewekt(X,t):
    vx0,omega = X
    return [x0+vx0*t, y0+vy0*t-g*t^2/2 ]

t=np.arange(0,30,.01)    
rozw=odeint(polewekt,[y0+vy0*t-g*t^2/2],t)

print rozw
plt.plot(t,rozw)
plt.show()