write a program to find the factorial of a given number (using do command)

implicit double precision (a-h,o-z)
    iprod=1
    write(*,*) 'n'
    read(*,*) n
      do i=1,n
    iprod=iprod*i
    enddo
    write(*,*) iprod
   
      stop
      end