Esope source files

   1 : C SINI      SOURCE    CHAT      05/01/13    03:19:45     5004           
   2 : C SINI SOURCE                                                                   
   3 :       SUBROUTINE SINI(X,SINTX)                                                  
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   6 : C                                                                               
   7 : C     FONCTION SINUS INTEGRAL                                                   
   8 : C                                                                               
   9 : C                                                                               
  10 : -INC CCREEL                                                                     
  11 : C                                                                               
  12 :       SINTX=0                                                                   
  13 : C                                                                               
  14 :       SAGN=SIGN(1.D0,X)                                                         
  15 :       X=ABS(X)                                                                  
  16 : C                                                                               
  17 :       IF(X.LE.1.7D0)THEN                                                        
  18 :         SINTX=X-X**3/18+X**5/600                                                
  19 :       ELSE                                                                      
  20 :         IF(X.GE.6.5D0)THEN                                                      
  21 :           SINTX=XPI/2.D0-COS(X)/X-SIN(X)/X**2+2.D0*COS(X)/X**3                  
  22 : C                                                                               
  23 :         ELSE                                                                    
  24 :           NPAS=1000                                                             
  25 :           NPASM1=999                                                            
  26 :           XSN=X/NPAS                                                            
  27 :           SINTX=0.5D0*(SIN(X)+1.D0)/NPAS                                        
  28 :           DO 20K=1,NPASM1                                                       
  29 :              SINTX=SINTX+SIN(K*XSN)/K                                           
  30 :  20       CONTINUE                                                              
  31 :         ENDIF                                                                   
  32 :       ENDIF                                                                     
  33 : C                                                                               
  34 :       SINTX=SAGN*SINTX                                                          
  35 : C                                                                               
  36 :       RETURN                                                                    
  37 :       END                                                                       
  38 :  
  39 :  
  40 :