#######################################################################################################################################
# #
# 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 E T O I L E ) : #
# #
# #
# Author of '$xiirk/.SYRA.A1.2.$U' : #
# #
# Jean-Francois Colonna (LACTAMME, 20130121100229). #
# #
#######################################################################################################################################
$Z setParam _____U0 7
$Z setParam _____TransformationX neutre$X
$Z # Parametre introduit le 20130125131220, 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 20130124083406 (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.1
$Z FilSTmpE FSyracuse
$Z set Nparts=`$CA $FCoOrD.RHO.1 | $WCl`
$Z $xci/valeurs_inte$X p=1 d=$Nparts \
$Z vD=0 vA=0 \
$Z lineaire=VRAI \
$Z > $FCoOrD.RHO.2
$Z $PAST $FCoOrD.RHO.2 $FCoOrD.RHO.1 | \
$Z $R "$K_TAB" "$K_NL" \
$Z > $FCoOrD.RHO
$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 $xrv/neutre$X \
$Z ne=0 \
$Z fichier== \
$Z exemplaires=2 \
$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 $xrv/neutre$X \
$Z ne=0 \
$Z fichier== \
$Z exemplaires=2 \
$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`
#20150415100248____:$Z $xci/valeurs_inte$X p=1 d=$Nparts \ #
#20150415100248____:$Z vD=$NOIR_PLANCHER vA=$BLANC \ #
#20150415100248____:$Z lineaire=VRAI \ #
#20150415100248____:$Z entier=VRAI \ #
#20150415100248____:$Z > $FCoOrD.NIVEAU #
$Z set ExEmPlAiReS=2
$Z @ NpartsS2 = $Nparts / $ExEmPlAiReS
$Z $xci/valeurs_inte$X p=1 d=$NpartsS2 \
$Z vD=$NOIR_PLANCHER vA=$BLANC \
$Z lineaire=VRAI \
$Z entier=VRAI | \
$Z $xrv/neutre$X ne=0 \
$Z fichier== \
$Z exemplaires=$ExEmPlAiReS \
$Z > $FCoOrD$NIVEAU