SIMPSON 3/8TH RULE

!        SIMPSON 3/8 RULE
         WRITE(*,*)'GIVE THE LIMTS'
         READ(*,*)TL,TU
         WRITE(*,*)'GIVE THE NO OF INTERVAL'
         READ(*,*)N
         H=(TU-TL)/N
         SUM=F(TU)+F(TL)+3*F(TU-H)+3*F(TU-2*H)
         DO 10 I=1,N-5,3
         SUM=SUM+3*F(TL+I*H)+3*F(TL+(I+1)*H)+2*F(TL+(I+2)*H)
  10     CONTINUE
         VALUE=3*(SUM*H)/8
         WRITE(*,20)VALUE
  20     FORMAT(5X,'THE VALUE OF INTEGRAL IS==',F8.3)
         STOP
         END
       
       
         FUNCTION F(X)
         FN=X**3+2*X**2-3*X+4
         FD=X**2+5*X-6
         F=FN/FD
         RETURN
         END