pickover++

Fraktalism/fraktals.f90

@@ -52,9 +52,10 @@ subroutine simple_julia(pic, cx, cy, maxiter)
 
 end subroutine simple_julia
 !-----------------------------------------------------
-!
-!
-!
+! 
+!  d'après les pages 91/92 du livre de Roger T Stevens
+!       "Fractal programming in C"
+! 
 subroutine pickover_0(pic, count)
     implicit none
     integer, intent(inout), dimension (:,:)  ::  pic
@@ -62,12 +63,13 @@ subroutine pickover_0(pic, count)
 
     double precision      :: xa, ya, za, xb, yb, zb
     double precision      :: ka, kb, kc, kd
-    integer               :: i, px, py
+    integer               :: i, w, h, px, py
 
     ka = 2.24 ; kb = 0.43 ; kc = -0.65 ; kd = -2.43
-
     xa = 0.00 ; ya = 0.00 ; za = 0.0
 
+    w  = ubound(pic, 1)
+    h = ubound(pic, 2)
 
     do i=1, count
 
@@ -75,8 +77,13 @@ subroutine pickover_0(pic, count)
       yb = za*sin(kc*xa) - cos(kd*ya)
       zb = sin(xa)
 
-      print *, i, xb, yb, zb
+      px = (xb * (w/4.05)) + (w / 2)
+      py = (yb * (h/4.05)) + (h / 2)
 
+      pic(px, py) = 200              ! WARNING COREDUMP
+
+      print *, xb, yb, zb
+      
       xa = xb ; ya = yb ; za = zb 
 
     enddo