home *** CD-ROM | disk | FTP | other *** search
/ Fractal Frenzy 2 / Fractal_Frenzy_Volume_II_Walnut_Creek_August_1995.iso / pc / formulas / mandne.frm < prev    next >
Encoding:
Text File  |  1993-04-10  |  1.3 KB  |  40 lines
  1. M-SetInNewton(XAXIS) {; use float=yes
  2.                       ; jon horner 100112,1700, 12 feb 93
  3.   z = 0,  c = pixel,  cminusone = c-1:
  4.   oldz = z,  nm = 3*c-2*z*cminusone, dn = 3*(3*z*z+cminusone),
  5.   z = nm/dn+2*z/3,   |(z-oldz)|>=|0.01|
  6.  }
  7.  
  8. F'M-SetInNewtonA(XAXIS) {; use float=yes
  9.                       ; jon horner 100112,1700, 12 feb 93
  10.   z = 0,  c = fn1(pixel),  cminusone = c-1:
  11.   oldz = z, nm = p1*c-2*z*cminusone, dn = p1*(3*z*z+cminusone),
  12.   z = nm/dn+2*z/p1, |(z-oldz)|>=|0.01|
  13.  }
  14.  
  15. ; #: 14384 S4/Fractal Sources
  16. ;    16-Feb-93 15:55:32
  17. ;Sb: #14099-Newtons Mandelbrot
  18. ;Fm: Chuck Ebbert 76306,1226
  19. ;To: Jon Horner 100112,1700 (X)
  20.  
  21. ;Jon, your formula led to these:
  22.  
  23. F'M-SetInNewtonB(XAXIS) { ; use float=yes, periodicity=no
  24.   ; set p1 >= 3, 1e-30 < p2 < .01
  25.   z=0, c=fn1(pixel), cm1=c-1, cm1x2=cm1*2, twoop1=2/p1, p1xc=c*real(p1):
  26.    oldz = z,
  27.    z= (p1xc - z*cm1x2 )/( (sqr(z)*3 + cm1 ) * real(p1) ) + z*real(twoop1),
  28.     |z - oldz| >= p2
  29.  }
  30.  
  31.  F'M-SetInNewtonC(XAXIS) { ; same as F'M-SetInNewtonB except for bailout
  32.   ; use float=yes, periodicity=no
  33.   ; (3 <= p1 <= ?) and (1e-30 < p2 < .01)
  34.   z=0, c=fn1(pixel), cm1=c-1, cm1x2=cm1*2, twoop1=2/p1, p1xc=c*real(p1):
  35.    z = (p1xc - z*cm1x2 )/( (sqr(z)*3 + cm1 ) * real(p1) ) + z*real(twoop1),
  36.     abs(|z| - real(lastsqr) ) >= p2
  37.  }
  38.  
  39.  
  40.