program calcul_puissance

  implicit none

  ! Declarations
  real (kind=8) :: x         ! valeur de x (DONNEE)
  integer :: n               ! valeur de n (DONNEE)
  real (kind=8):: puissance  ! valeur de x^n (RESULTAT)
  integer :: i               ! variable INTERMEDIAIRE
  
  
  ! Instructions

  print*, "Donner la valeur de x (reelle):"
  read*, x
  print*, "Donner la valeur de n (entier strictement positif):"
  read*, n
  
  puissance = 1.
  
  if (n<0) then 
    do i=1,-n
      puissance = puissance * x
    end do
    print*,"x^n =",1./puissance
  else
    do i=1,n
      puissance = puissance * x
    end do
    print*,"x^n =",puissance
  end if

end program calcul_puissance