subroutine square(x,answer)
        integer, intent(in) :: x
        integer, intent(out) :: answer
        answer = x * x
      end