import numpy as np
pi=np.pi
ts=12 ; U=5 ; g=9.8 ; f=2 ; h=10
tp=np.linspace(0,ts,500)
alfa0=pi/3 ; A=0.1
f1=U*(np.cos(alfa0)+A*np.cos(2*pi*f*tp))*(ts-tp)
f2=h+(U*(np.sin(pi/2)+0.1*np.sin(2*pi*f*tp)))*(ts-tp)-g*(ts-tp)**2/2
import matplotlib.pylab as plt
plot=plt.plot
plot(f1,f2,'o') ; plt.show();plt.axis([-2,10,0,13])
plt.xlabel('x') ; plt.ylabel('y') ;
plt.savefig('.')