Rotation of Basic Shapes

xxxxxxxxxx

 

1

def rotate_xy(k,x,y):

2

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

3

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

4

            for i in range(2*k)]

5

def circle_mandala(k,n,m):

6

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

7

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

8

    col=[[randint(300,900)/1000 

9

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

10

    c=sum([sum([circle([x,y],r[i],facecolor=col[i],

11

                       fill=True,alpha=.05)+\

12

                circle([x,y],r[i],color=col[m-i-1],thickness=.7) 

13

                for [x,y] in centers[i]]) for i in range(m)])

14

    c.show(axes=false,figsize=9,

15

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

16

k=randint(100,900)/1000; n=randint(5,15); m=randint(3,9)

17

circle_mandala(k,n,m)

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

---

Messages