Procedures | Cast3M

   1 : * @STAT     PROCEDUR  AM        93/11/16    21:15:14     978            
   2 : *-----------------------------------------------------------            
   3 : ******              PROCEDURE @STAT                   ******
   4 : *-----------------------------------------------------------            
   5 : *                                                                       
   6 : *    CETTE PROCEDURE A ETE MISE GRACIEUSEMENT                           
   7 : *   A DISPOSITION DE LA COMMUNAUTE  CASTEM2000                          
   8 : * PAR  DELERUYELLE Fr.  (SOCOTEC-INDUSTRIE à l'IPSN/DES)
   9 : *                                                                       
  10 : *-----------------------------------------------------------            
  11 : *           OPERATIONS STATISTIQUES COURANTES
  12 : *----------------------------------------------------------
  13 : 'DEBP' @STAT LX*'LISTREEL' LY/'LISTREEL' ;
  14 : XI = 0. ; XI2 = 0. ; NX = 'DIME' LX ;
  15 : 'SI' ('EXISTE' LY) ;
  16 : * DEUX LISTES
  17 :  YI = 0. ; YI2 = 0. ; XIYI = 0. ; NY = 'DIME' LY ;
  18 :  'SI' (NX 'EGA' NY) ;                                                    
  19 :   K = 0 ;                                                                
  20 :   'REPE' BOU1 NX ;                                                       
  21 :    K = K + 1 ;                                                           
  22 :    XK = 'EXTR' LX K ; YK = 'EXTR' LY K ;
  23 :    XI  = XI    +  XK ;
  24 :    XI2 = XI2   + (XK * XK) ;
  25 :    YI  = YI    +  YK ;
  26 :    YI2 = YI2   + (YK * YK) ;
  27 :    XIYI = XIYI + (XK * YK) ;
  28 :   'FIN' BOU1 ;                                                           
  29 :   XM = XI / NX ; YM = YI / NY ;
  30 :   VX =    (XI2 / NX) - (XM * XM) ;
  31 :   VY =    (YI2 / NY) - (YM * YM) ;
  32 :   COXY = (XIYI / NX) - (XM * YM) ;
  33 :   EX = VX ** 0.5 ; EY = VY ** 0.5 ;                                     
  34 :   A = COXY / VX ; B = YM - (A * XM) ; R = COXY / (EX * EY) ;
  35 :  'SINON' ;                                                               
  36 :    'MESS' '    ' ;                                                       
  37 :    'MESS' '*** STAT ==> Les listes lx et ly n ont pas meme dimension' ;  
  38 :    'MESS' '             Dim(lx)=' ('DIME' LX) 'Dim(ly)=' ('DIME' LY) ;   
  39 :    'MESS' '    ' ;                                                       
  40 :  'FINSI' ;                                                               
  41 :  'RESP' XM EX YM EY A B R ;
  42 : 'SINON' ;
  43 : * UNE SEULE LISTE
  44 :   K = 0 ;                                                                
  45 :   'REPE' BOU2 NX ;
  46 :    K = K + 1 ;                                                           
  47 :    XK = 'EXTR' LX K ;
  48 :    XI  = XI    +  XK ;
  49 :    XI2 = XI2   + (XK * XK) ;
  50 :   'FIN' BOU2 ;
  51 :   XM = XI / NX ;
  52 :   VX =    (XI2 / NX) - (XM * XM) ;
  53 :   EX = VX ** 0.5 ;
  54 :  'RESP' XM EX ;
  55 : 'FINSI' ;
  56 : 'FINP' ;