Yesterday I somehow managed to write this nice little program...
sync on : sync rate 0 : hide mouse
do
ox=320+sin(sin(b#)*sin(b#)*a#)*sin(sin(b#)*sin(b#)*a#-r#)*250
oy=240+cos(sin(b#)*sin(b#)*a#)*sin(r#-sin(b#)*sin(b#)*a#)*150
cls
`ink rgb(255,255,255),0 : print screen fps()
for i=1 to 50
for p=1 to 40
x=ox+sin(sin(b#)*a#)*i+sin(wrapvalue(cos(b#)*(a#+p*11.5-i*2)))*(60+sin(r#)*i*5+50*sin(wrapvalue(r#+p*4+i*4)))
y=oy+sin(sin(b#)*a#)*i+cos(wrapvalue(cos(b#)*(a#+p*11.5-i*2)))*(60+sin(r#)*i*5+50*sin(wrapvalue(r#+p*4+i*4)))
ink rgb(155+sin(wrapvalue(a#+p*7+i*2))*100,155+sin(wrapvalue(a#+r#+p*5+i*2))*100,155+sin(wrapvalue(r#+p*4+i*2))*100),0
dot x,y
next p
next i
sync
a#=wrapvalue(a#+0.7)
r#=wrapvalue(r#+1.1)
b#=wrapvalue(b#+0.35)
if mousemovey()0 or mousemovex()0 then exit
loop