Random Coefficients & Colors 3

xxxxxxxxxx

 

1

import numpy as np,pylab as pl; pi=np.pi

2

n=24; t=np.arange(0,2*pi,.1^3*2*pi/n)

3

def randi(nmin,nmax): 

4

    return np.random.randint(nmin,nmax)

5

a,b,c=randi(5,11),randi(12,24),randi(25,81)

6

d,m=randi(216,256),randi(100,300)

7

pl.figure(figsize=(6,6)); ax=pl.gca()

8

ax.set_facecolor('black')

9

for i in range(n):

10

    f1=(a+.9*np.cos(b*t+2*pi*i/n))*(1+.1*np.cos(c*t+2*pi*i/n))

11

    f2=(1+.01*np.cos(d*t+2*pi*i/n))*(1+np.sin(t+2*pi*i/n))

12

    x=f1*f2*np.cos(t); y=f1*f2*np.sin(t)

13

    pl.scatter(x,y,s=3*.1^3,c=np.random.random(3))

14

st='a=%d; b=%d; c=%d; d=%d; m=%d'

15

pl.title(st%(a,b,c,d,m)); pl.grid(c='white',alpha=.3)

16

pl.tight_layout(); pl.show()

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

---

Messages