#######################################################################################################################################
#                                                                                                                                     #
#         P A R A D O X E   D E   B E R T R A N D  :                                                                                  #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirk/.BERT.11.1.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20050122112952).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z        setParam            _____Iterations1              200

$Z        SETParam            _____Iterations2              1000000

$Z        setParam            _____Seuil                    `calcul (((4*$pi) - (3*sqrt(3))) / (12*$pi))`

$Z        setParam            _____Nombre                   10

$Z        $DELETE             $xTV/RAPPORT.1

$Z        $xrk/ParadoxeB.01$X                                                                                                           \
$Z                            iterations=$_____Iterations1                                                                              \
$Z                            E=$_____ImageE                                                                                            \
$Z                            D=$_____ImageD                                                                                            \
$Z                            multicolore=VRAI                                                                                          \
$Z                            niveau_inferieur_N=$GRIS_1                                                                                \
$Z                            niveau_superieur_N=$GRIS_8                                                                                \
$Z                            editer=FAUX proportions=VRAI probabilite=FAUX                                                             \
$Z                            R=$xTV/BERTRAND.1                                                                                         \
$Z                                                                                                                      $formatI        \
$Z                                                                                                            > $xTV/RAPPORT.1

$Z        $DELETE             $xTV/RAPPORT.2

$Z        $xrk/ParadoxeB.01$X                                                                                                           \
$Z                            iterations=$_____Iterations2                                                                              \
$Z                            E=$_____ImageE                                                                                            \
$Z                            D=$_____ImageD                                                                                            \
$Z                            multicolore=VRAI                                                                                          \
$Z                            niveau_inferieur_N=$GRIS_1                                                                                \
$Z                            niveau_superieur_N=$GRIS_8                                                                                \
$Z                            editer=FAUX proportions=VRAI probabilite=FAUX                                                             \
$Z                            R=$xTV/BERTRAND.2                                                                                         \
$Z                                                                                                                      $formatI        \
$Z                                                                                                            > $xTV/RAPPORT.2

$Z        $xci/neutre$X                                                                                                                 \
$Z                            A=$xTV/BERTRAND.1                                                                                         \
$Z                            R=$xTV/BERTRAND$ROUGE                                                                                     \
$Z                                                                                                                      $formatI
$Z        $xci/convol.01$X                                                                                                              \
$Z                            A=$_____ImageD                                                                                            \
$Z                            points=81 circulaire=VRAI                                                                                 \
$Z                            R=$xTV/BERTRAND$VERTE                                                                                     \
$Z                                                                                                                      $formatI
$Z        $xci/init$X                                                                                                                   \
$Z                            niveau=$NOIR                                                                                              \
$Z                            R=$xTV/BERTRAND$BLEUE                                                                                     \
$Z                                                                                                                      $formatI


$Z        $DELETE             $xTV/RAPPORT.3

$Z        set                 CoMpLeMeNt_SeUiL=`$xcg/SOUS.01$X a=1 b=$_____Seuil`

$Z        $xrv/INCLcc.11$X                                                                                                              \
$Z                            ne=0                                                                                                      \
$Z                            fichierX=$xTV/RAPPORT.2                                                                                   \
$Z                            fichierA=$_____Seuil                                                                                      \
$Z                            fichierB=$CoMpLeMeNt_SeUiL                                                                                \
$Z                                                                                                            > $xTV/RAPPORT.3

$Z        set                 AiRe_E=`$xci/comptage$X A=$_____ImageE origine=$NOIR_PLANCHER extremite=$BLANC Prme=VRAI`
$Z        set                 AiRe_D=`$xci/comptage$X A=$_____ImageD origine=$NOIR_PLANCHER extremite=$BLANC Prme=VRAI`

$Z        set                 InTeR_2=`$CA $xTV/RAPPORT.2 | $WCl`
$Z        set                 InTeR_3=`$CA $xTV/RAPPORT.3 | $WCl`

$Z        echo                "AIRE(D)/AIRE(E) droites.................. = $AiRe_D/$AiRe_E (="`calculNS $AiRe_D/$AiRe_E`")"
$Z        echo                "Nombre de droites intersectant le Domaine = $InTeR_2"
$Z        echo                "Probabilite.............................. = $InTeR_3/$InTeR_2 (="`calculNS $InTeR_3/$InTeR_2`")"
$Z        echo                "Nombre de droites visualisees............ = "`$CA $xTV/RAPPORT.1 | $WCl | $AW ' { print $1 } '`

$Z        $DELETE             $xTV/HISTOGRAMME

$Z        @                   NoMbRe = 2 * $_____Nombre + 1
$Z                                      # Il est preferable d'utiliser un nombre impair d'elements pour l'histogramme...              #

$Z        $xrv/histogram.01$X                                                                                                           \
$Z                            ne=0                                                                                                      \
$Z                            fichier=$xTV/RAPPORT.2                                                                                    \
$Z                            nombre=$NoMbRe                                                                                            \
$Z                            forcer=VRAI minimum=0 maximum=1                                                                           \
$Z                                                                                                            > $xTV/HISTOGRAMME



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