#!/usr/bin/python

import numpy as np
import matplotlib.pyplot as plt 
import math

fig = plt.figure()


# initial values
Ax = 1.0
Ay = 1.0

omegax = 1
omegay = 2
phaseshift = 0.0*math.pi


tmax=10.0
h=0.02
n = int((tmax)/h)

t = np.linspace(0,tmax,num=n)
x = np.zeros(n)
y = np.zeros(n)


for i in range(0,len(t)):
    x[i] = Ax * math.sin(omegax * t[i] + phaseshift)
    y[i] = Ay * math.sin(omegay * t[i])


plt.plot(x,y)


plt.show()