FUNCTION SSINT(F,XS,XF,EPS,N,FS)
      EXTERNAL F
      COMMON /INTEG/ KOD
      DATA AEPS/5./
C     BÛ×ÈCËEHÈE ÈHTEÃPAËA 3-Mß METOÄAMÈ;
C     SPR  - ÏPßMOÓÃOËÜHÈKOB (S1)
C     STR  - TPAÏEÖÈÈ        (S2)
C     SSIM - CÈMÏCOHA        (S3)
C
C     F     - ÏOÄÈHTEÃPAËÜHAß ÔÓHKÖÈß
C     XS,XF - ÈHTEPBAË ÈHTEÃPÈPOBAHÈß
C     E     - TO×HOCTÜ 
C
      E=EPS*AEPS
      N=0
      AF=F(XF)
      S3=AF+FS
      S2=S3*.5
      FS=AF
      H=XF-XS
C
   50 S1=0
      H2=H*.5
      X=XS+H2
   51 FX=F(X)
      S1=S1+FX
      X=X+H
      IF(X-XF) 51,52,52
C
   52 SPR=S1*H
      STR=S2*H
      SSIM=(S3+4.*S1)*(H/6.)
      SINT=(SPR+STR)*.5
      IF(KOD.GE.7) WRITE(6,222)N,SPR,STR,SINT,SSIM
  222 FORMAT(' ***7 SSINT*** N SPR STR SSINT SSIM =',I5,4F12.5)
      IF(ABS(STR-SPR)-E) 53,53,54
   54 S2=S2+S1
      S3=S3+2.*S1
      H=H2
      N=N+1
      GOTO 50
C
   53 CONTINUE
      SINT=(SPR+STR)*.5
      IF(KOD.GE.6) WRITE(6,20)XS,XF,E,N,SPR,STR,SSIM,SINT
   20 FORMAT(' ***6 SSINT*** XS XF E N SPR STR SSIM SSINT ='
     *,2F10.2,F8.4,I5,4F10.3)
      SINT=SSIM
      SSINT=SINT
      RETURN
      END