#######################################################################################################################################
# #
# 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.C1.2.$U' : #
# #
# Jean-Francois Colonna (LACTAMME, 20130121175907). #
# #
#######################################################################################################################################
$Z setParam _____U0 7
$Z setParam _____TransformationX neutre$X
$Z # Parametre introduit le 20130125131230, en notant que : #
$Z # #
$Z # PUIX.01$X exposant=0.02 #
$Z # #
$Z # permet de prendre mieux en compte la "dynamique" du processus... #
$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 stop_421=VRAI | \
$Z $xrv/$_____TransformationX \
$Z ne=0 \
$Z fichier== | \
$Z $xrv/normalise.01$X \
$Z ne=0 \
$Z fichier== \
$Z >> $FSyracuse
$Z # Le 20130124083225 (avec retard...) fut introduit "p=1 d=100000" afin d'atteindre {4,2,1}... #
$Z $xrv/neutre$X ne=0 \
$Z fichier=$FSyracuse \
$Z > $FCoOrD.RHO
$Z FilSTmpE FSyracuse
$Z set Nparts=`$CA $FCoOrD.RHO | $WCl`
$Z @ NpartsP1 = $Nparts + 1
$Z $xci/valeurs_inte$X p=1 d=$NpartsP1 \
$Z vD=$pis2 vA=$pis2 \
$Z lineaire=VRAI | \
$Z $HEA --lines=-1 \
$Z > $FCoOrD.THETA
$Z $xci/valeurs_inte$X p=1 d=$NpartsP1 \
$Z vD=0 vA=$dpi \
$Z lineaire=VRAI | \
$Z $HEA --lines=-1 \
$Z > $FCoOrD.PHI
$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 set Nparts=`$CA $FCoOrD.RHO | $WCl`
$Z $xci/valeurs_inte$X p=1 d=$Nparts \
$Z vD=$NOIR_PLANCHER vA=$BLANC \
$Z lineaire=VRAI \
$Z entier=VRAI \
$Z > $FCoOrD$NIVEAU