#######################################################################################################################################
#                                                                                                                                     #
#         E T U D E   D E   L A   C O N J E C T U R E   D E   S Y R A C U S E                                                         #
#         E N   C O O R D O N N E E S   P O L A I R E S   ( E N   S P I R A L E )  :                                                  #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirk/.SYRA.I1.2.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20211026112009).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  setParam            _____U0                       7

$Z                  setParam            _____Parite0                  -0.02
$Z                  setParam            _____Parite1                  +0.02

$Z                  setParam            _____TPhi                     0

$Z                  SETParam            _____FRayon                   0
$Z                                      # Parametre introduit le 20211028165549...                                                    #

$Z                  FilSTmpB            FCoOrD
$Z                                      # ATTENTION : les 'FilSTmpE's ont lieu ailleurs...                                            #

$Z                  FilSTmpB            FSyracuse

$Z                  $xci/valeurs_Syra$X                                                                                                 \
$Z                                      p=1 d=100000                                                                                    \
$Z                                      u0=$_____U0                                                                                     \
$Z                                      editer_pariteparite=VRAI                                                                        \
$Z                                      stop_421=VRAI                                                                             |     \
$Z                  $SE                                                                                                                 \
$Z                                      -e 's/^0$'"/$_____Parite0/"                                                                     \
$Z                                      -e 's/^1$'"/$_____Parite1/"                                                                     \
$Z                                                                                                                                |     \
$Z                  $AC                                                                                                                 \
$Z                                                                              >> $FSyracuse

$Z                  $xrv/ADD2.11$X      ne=0                                                                                            \
$Z                                      fichier1=$FSyracuse                                                                             \
$Z                                      fichier2=$_____TPhi                                                                             \
$Z                                                                              >  $FCoOrD.PHI

$Z                  set                 Nparts=`$CA $FSyracuse | $WCl`

$Z                  FilSTmpE            FSyracuse

$Z                  set                 RhO_0=0.00

$Z                  if                  ($_____FRayon == 0) then
$Z                                      set                 RhO_N=1.00
$Z                  else
$Z                                      set                 RhO_N=`calcul $_____FRayon*$Nparts`
$Z                  endif

$Z                  $xci/valeurs_inte$X p=1 d=$Nparts                                                                                   \
$Z                                      vD=$RhO_0 vA=$RhO_N                                                                             \
$Z                                      lineaire=VRAI                                                                                   \
$Z                                                                              >  $FCoOrD.RHO
$Z                                      # Ainsi, tous les 'U(0)' seront sur un cercle de rayon 1 alors que les sequences {4,2,1}      #
$Z                                      # seront au centre...                                                                         #

$Z                  $xci/valeurs_inte$X p=1 d=$Nparts                                                                                   \
$Z                                      vD=$pis2 vA=$pis2                                                                               \
$Z                                      lineaire=VRAI                                                                                   \
$Z                                                                              >  $FCoOrD.THETA
$Z                                      # Avec ces valeurs, on prepare donc :                                                         #
$Z                                      #                                                                                             #
$Z                                      #                   X = rho.cos(phi).sin(theta)                                               #
$Z                                      #                   Y = rho.sin(phi).sin(theta)                                               #
$Z                                      #                   Z = rho.cos(theta)                                                        #
$Z                                      #                                                                                             #
$Z                                      # soit, puisque 'theta' est egal a pi/2 :                                                     #
$Z                                      #                                                                                             #
$Z                                      #                   X = rho.cos(phi)                                                          #
$Z                                      #                   Y = rho.sin(phi)                                                          #
$Z                                      #                                                                                             #
$Z                                      # nous placant ainsi en coordonnees polaires...                                               #
$Z                                      #                                                                                             #
$Z                                      #                                                                                             #
$Z                                      # Mais ATTENTION : on notera le 20211026155108 que {RHO,PHI[,THETA]} definissent des          #
$Z                                      # coordonnees spheriques dans l'espace "absolu" {X,Y,Z} et non pas une veritable marche       #
$Z                                      # pseudo-aleatoire ou en chaque point on tournerait de 'PHI' en se deplacant de 'RHO'...      #

$Z                  $xci/valeurs_inte$X p=1 d=$Nparts                                                                                   \
$Z                                      vD=$BLANC vA=$NOIR_PLANCHER                                                                     \
$Z                                      lineaire=VRAI                                                                                   \
$Z                                      entier=VRAI                                                                                     \
$Z                                                                              >  $FCoOrD$NIVEAU



Copyright © Jean-François Colonna, 2021-2021.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2021-2021.