home *** CD-ROM | disk | FTP | other *** search
- \ MandelPot
- \ Calculates the M-potential of a complex number
- \ stored in creal,cimag
-
- var creal, cimag, zreal, zimag, ztemp
- var iter, potential
-
- creal=0; cimag=1.01
- zreal=0; zimag=0; iter=0
-
- repeat
- ztemp=zreal
- zreal=zreal*zreal-zimag*zimag+creal
- zimag=2*ztemp*zimag+cimag
- iter=iter+1
- until zreal*zreal+zimag*zimag>4 OR iter>127
-
- if iter<=127
- potential=0.5*LOG(zreal*zreal+zimag*zimag)/(2^iter)
- else
- potential=0
- endif
-
- show potential
-