Procedures | Cast3M

   1 : * PROBDENS  PROCEDUR  CHLOE     98/11/18    21:17:30     3360           
   2 : 'DEBPROC' PROBDENS MU*'FLOTTANT' SIGMA*'FLOTTANT' 
   3 :                    RB1*'FLOTTANT' B2*'FLOTTANT';
   4 : *|=====================================================================
   5 : *|
   6 : *| << OBJET >> :
   7 : *|
   8 : *| Procedure permettant de définir la densite de probabilite et
   9 : *| la fonction de repartition d'une variable aletoire connaissant
  10 : *| ses parametres statistiques (moyenne, ecart-type, symetrie,
  11 : *| aplatissement)
  12 : *|
  13 : *|
  14 : *| NB : Les évolutions résultats sont données entre
  15 : *|           Mu - 6 * Sigma et Mu + 6 * Sigma 
  16 : *|
  17 : *|
  18 : *| << ENTREE >> :
  19 : *| -- ARGUMENTS OBLIGATOIRES DANS TOUS LES CAS --
  20 : *|
  21 : *| MU = objet de type REEL designant la valeur moyenne
  22 : *|                    de la variable aleatoire.
  23 : *|
  24 : *| SIGMA = objet de type REEL designant l'ecart-type
  25 : *|                       de la variable aléatoire.
  26 : *|
  27 : *| RB1 = objet de type REEL designant la valeur du coefficient 
  28 : *|                    de symétrie de la variable aleatoire.
  29 : *|
  30 : *| B2 = objet de type REEL designant la valeur du coefficient
  31 : *|                    d'aplatissement de la variable aleatoire.
  32 : *| << SORTIE >> :
  33 : *|
  34 : *|
  35 : *| PDF = objet de type EVOLUTION contenant la densite de probabilite
  36 : *|                     de la variable aléatoire
  37 : *|
  38 : *| CDF = objet de type EVOLUTION contenant la fonction de repartition 
  39 : *|                     de la variable aléatoire
  40 : *|
  41 : *|=====================================================================
  42 : *                                                                               
  43 : * ______________________                                                        
  44 : *                                                                               
  45 : * DECLARATIONS GENERALES                                                        
  46 : * ______________________                                                        
  47 : *                                                                               
  48 : SI = 'MOT' 'SI' ; SINON = 'MOT' 'SINON' ; FINSI = 'MOT' 'FINSI' ;               
  49 : REPETER = 'MOT' 'REPETER'; QUITTER = 'MOT' 'QUITTER'; FIN = 'MOT' 'FIN';        
  50 : EXISTE = 'MOT' 'EXISTE' ; EGA = 'MOT' 'EGA' ;                                   
  51 : ERREUR = 'MOT' 'ERREUR' ;                                                       
  52 : ENTIER   = 'MOT' 'ENTIER  ' ;                                                   
  53 : FLOTTANT = 'MOT' 'FLOTTANT' ;                                                   
  54 : MOT      = 'MOT' 'MOT     ' ;                                                   
  55 : * _____                                                        
  56 : *                                                                               
  57 : * TESTS                                                        
  58 : * _____                                                        
  59 : *                                                                               
  60 : * _________________                                                        
  61 : *                                                                               
  62 : * DEBUT DES CALCULS                                                        
  63 : * _________________                                                        
  64 : *
  65 : X1 = PROG;
  66 : Y1 = PROG;
  67 : Y2 = PROG;
  68 : oldy = 0.;
  69 : *
  70 : NPAS = 100;
  71 : REPE BOUC NPAS;
  72 :   absx = (MU - (6 * SIGMA)) + ((&BOUC*(12*SIGMA))/NPAS);
  73 :   X1 = X1 ET (PROG absx);
  74 :   newy = PROB MU SIGMA RB1 B2 absx;
  75 :   Y1 = Y1 ET (PROG ((newy - oldy)/((12*SIGMA)/NPAS)));
  76 :   Y2 = Y2 ET (PROG newy);
  77 :   oldy = newy;
  78 : FIN BOUC;
  79 : CDF = evol 'MANU' 'Valeur' X1 'Repartition' Y2 ;
  80 : PDF = evol 'MANU' 'Valeur' X1 'Densite' Y1 ;
  81 : 'FINPROC' PDF CDF;
  82 :  
  83 :