SIMPSON'S 1/3RD RULE

!        SIMPSON 1/3 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)+4*F(TU-H)
         DO 10 I=1,N-3,2
         SUM=SUM+4*F(TL+I*H)+2*F(TL+(I+1)*H)
  10     CONTINUE
         VALUE=(SUM*H)/3
         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