program factorial
      real fact 
      integer n
      print*,'Enter n of factorial'
      read*,n
      fact=1
      do i=1,n
         fact=fact*float(i) 
      enddo
      print*,'factorial of ',n,' is ',fact 
      end