Esope source files

   1 : C EMME      SOURCE    CHAT      05/01/12    23:38:33     5004           
   2 : c***********************************************************************        
   3 : c subroutine per il calcolo dello scarto del coefficiente di Weibull   *        
   4 : c***********************************************************************        
   5 :       subroutine emme(n,rml,sig,sum1,diff)                                      
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 :       dimension sig(n)                                                          
   9 :       sum=0.                                                                    
  10 :       sum1=0.                                                                   
  11 :       sum2=0.                                                                   
  12 :       do 10 i=1,n                                                               
  13 :         sum = sum  + ( (sig(i)**rml) * dlog(sig(i)) )  
  14 :         sum1= sum1 + sig(i)**rml                                                     
  15 :         sum2= sum2 + dlog(sig(i))                                                    
  16 :   10  continue                                                                  
  17 :       rn=float(n)                                                               
  18 :       secmem=rn/(rn*sum/sum1-sum2)                                              
  19 :       diff=secmem-rml                                                           
  20 :       return                                                                    
  21 :       end                                                                       
  22 :