From 87191666b4659ddb5d13d9899048a9bc2d6a7a20 Mon Sep 17 00:00:00 2001
From: tTh <tth@konrad.lan>
Date: Sat, 11 Feb 2023 16:54:55 +0100
Subject: [PATCH] trigofest again

---
 GrafAnim/README.md     |  2 +-
 GrafAnim/trigofest.f90 | 12 +++++++-----
 2 files changed, 8 insertions(+), 6 deletions(-)

diff --git a/GrafAnim/README.md b/GrafAnim/README.md
index 20d2682..5c83cf7 100644
--- a/GrafAnim/README.md
+++ b/GrafAnim/README.md
@@ -11,7 +11,7 @@ un peu foireux sur les tracés de ligne...
 
 ## trigofest
 
-Distorsions de la courbe de Lissajous.
+Distorsions approximatives de la courbe de Lissajous.
 
 ## doubledice
 
diff --git a/GrafAnim/trigofest.f90 b/GrafAnim/trigofest.f90
index 5ae9f9c..731b714 100644
--- a/GrafAnim/trigofest.f90
+++ b/GrafAnim/trigofest.f90
@@ -19,15 +19,17 @@ program trigofest
   W = 512     ;       H = 342
   allocate(picz(W,H), stat=errcode)
 
-  blouber = 0.01
+  blouber = 0.1
   do loop=0, 359
     call spirale(picz, blouber, loop*9)
-    blouber = blouber * 1.0259
+    blouber = blouber + 0.3333
     write (filename, "(a, i5.5, a)") "F/spi/", loop, ".pgm"
     call spit_as_pgm_8(picz, trim(filename))
     print *, loop, blouber
   enddo
 
+  deallocate(picz)
+
   STOP ': WORLD FINISHED'
 
 contains        !------------------------------------------
@@ -102,13 +104,13 @@ subroutine spirale(pic, inirad, param)
     angle = real(foo) * 0.01745329252
 
     ! rx = radius * sin(angle) * 1.21
-    kx = sin(angle*2)
+    kx = 1.55 * sin(angle+(0.04*radius))
     rx = radius * kx
 
     ! ry = radius * cos(angle)
-    ky = cos(angle) + (0.5*cos(angle*4.0))
+    ky = cos(angle) + (0.5*cos(angle*6.0))
     ry = radius * ky
-    radius = radius + 0.025
+    radius = radius + 0.0245
     ix = int(rx) + W/2
     iy = int(ry) + H/2