FORTRAN PROGRAM 2 Trapezoidal Rule

Q write a program using fortram to find area using trapezoidal rule
Ans  
         program trap
implicit none
         REAL :: H,Y1,YLAST,y,A,S
       
         INTEGER :: i,N
       
         WRITE(*,*) "ENTER NUMBER OF TERMS"
         READ(*,*) N
         WRITE(*,*) "ENTER STEP SIZE"
         READ(*,*) H
         WRITE(*,*) "ENTER FIRST TERM"
         READ(*,*) Y1
         WRITE(*,*) "ENTER LAST TERM"
         READ(*,*) YLAST
         S=Y1+YLAST
         DO i=1,N-1
         WRITE(*,*) "ENTER L",i,"TH TERM"
         READ(*,*) y
         S=S+2*y
   
         END DO
       
         A=(H*S)/2
         WRITE(*,*) "AREA IS"
         WRITE(*,*) A
         READ(*,*)
         end program trap