In [18]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib






    




/home/uri/anaconda3/lib/python3.5/site-packages/IPython/core/magics/pylab.py:161: UserWarning: pylab import has clobbered these variables: ['square']
`%matplotlib` prevents importing * from pylab and numpy
  "\n`%matplotlib` prevents importing * from pylab and numpy"

In [19]:

    

triang = array([0,2*pi/3,4*pi/3,0])
square = array([0,pi/2,pi,3*pi/2,0])
square += pi/4
rad = 1
triangxs = rad*sin(triang)
triangys = rad*cos(triang)
squarexs = rad*sin(square)
squareys = rad*cos(square)
plot(triangxs,triangys,'C1')

plt.axes().set_aspect('equal', 'datalim')
ylim(-1,1)
axis('off')
savefig('triang.pdf')


    

In [20]:

    

plot(squarexs,squareys,'C0')

plt.axes().set_aspect('equal', 'datalim')
ylim(-1,1)
axis('off')
savefig('square.pdf')


    

In [21]:

    

plot(triangxs,triangys,'C1')
plot(squarexs,squareys,'C0')

plt.axes().set_aspect('equal', 'datalim')
ylim(-1,1)
axis('off')
savefig('twoshapes.pdf')


    

In [22]:

    

plot(squarexs,squareys,'C0')
#boundaries = [arange(x,x+pi/3,pi/60) for x in square]
#midpoints = []
#for x in boundaries:
#    midpoints.append(x[-1]+pi/12)
#boundaries = concatenate(boundaries)
#boundariesxs = rad*sin(boundaries)
#boundariesys = rad*cos(boundaries)

for i in arange(0,pi/3,pi/60):
    square+=pi/60
    squarexs = rad*sin(square)
    squareys = rad*cos(square)
    plot(squarexs,squareys,'C0')
    
triang-=pi/12
triangxs = rad*sin(triang)
triangys = rad*cos(triang)
plot(triangxs,triangys,'C1')
#plot(boundariesxs,boundariesys,'k')

plt.axes().set_aspect('equal', 'datalim')
ylim(-1,1)
axis('off')
savefig('twoshapesslack.pdf')


    

In [ ]: