File: julia.ode

package info (click to toggle)
xppaut 5.85-3.3
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 11,568 kB
  • ctags: 7,729
  • sloc: ansic: 80,804; cpp: 965; makefile: 271
file content (10 lines) | stat: -rw-r--r-- 224 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
 
# julia set for z -> z^2+c
par cx=0,cy=0
r=sqrt(sqrt((x-cx)^2+(y-cy)^2))
th=atan2(y-cy,x-cx)/2
s=sign(ran(1)-.5)
x'=s*r*cos(th)
y'=s*r*sin(th)
@ total=2000, meth=disc,transient=100,xlo=-2,ylo=-2,xhi=2,yhi=2,xp=x
@ yp=y
done