from pylab import *
#import pylab
from mpl_toolkits.mplot3d.axes3d import Axes3D
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import axes3d


def drange(start, stop, step):
    r = start
    while r < stop:
        yield r
        r += step


a=[]
b=[]
z1=[]
c=[]
d=[]
z2=[]




plt.ion()

fig = plt.figure(dpi=100)
ax = axes3d.Axes3D(fig)#fig.add_subplot(111, projection = '3d')

t = 0
dt = 1
tn= 50

for i in drange (t+dt, tn,dt):
    aa = sin(i)
    bb = cos(i)
    cc = aa*bb
    aa1 = 1.5 *sin(i)
    bb1 = 1.5*cos(i)
    cc1 = aa1*bb1
    a.append(aa)
    b.append(bb)
    z1.append(cc)
    c.append(aa1)
    d.append(bb1)
    z2.append(cc1)
    #clf()
    ax.plot(a,b,z1)#linestyle='None', marker='o', markersize = 5)
    ax.plot(a,b,z1, marker='o', linestyle='None')
    ax.plot(c,d,z2)
    ax.plot(c,d,z2, marker='o', linestyle='None')
    plt.draw()

#grid()
#show()
