Rotation of Basic Shapes 2

xxxxxxxxxx

 

1

import pylab as pl; import matplotlib.patches as pt

2

import numpy as np; import random as rd

3

def rotate_xy(k,x,y):

4

    return np.array([[np.cos(i*np.pi/k)*x-np.sin(i*np.pi/k)*y,

5

                      np.sin(i*np.pi/k)*x+np.cos(i*np.pi/k)*y] 

6

                     for i in range(2*k)])

7

def circle_mandala(k,n,m):

8

    c=np.array([rotate_xy(n,.5+k*i,.5+k*i) for i in range(m)])

9

    r=[rd.randint(300,1800)/1000 for i in range(m)]

10

    l=.1+r[m-1]+c.max()

11

    col=[[rd.randint(100,900)/1000 

12

          for j in range(2)]+[1] for i in range(m)]

13

    fig,ax=pl.subplots(figsize=(6,6))

14

    for i in range(m):

15

        for j in range(2*n):

16

            ax.add_patch(pt.Circle(c[i][j],r[i],

17

                         alpha=.1,color=col[i]))

18

            ax.add_patch(pt.Circle(c[i][j],r[i],

19

                         fill=False,edgecolor=col[m-i-1]))

20

    pl.title(r'$k=%.3f; \ m=%d; \ n=%d$'%(k,m,n))

21

    pl.axis('off'); pl.xlim(-l,l); pl.ylim(-l,l)

22

    pl.tight_layout(); pl.show()

23

k=rd.randint(100,900)/1000; 

24

n,m=rd.randint(5,15),rd.randint(3,9)

25

circle_mandala(k,n,m)

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages