SUBROUTINE my_square (x)
IMPLICIT DOUBLE PRECISION (a-h, o-z)
INTENT (INOUT) :: x
x = x * x
RETURN
END
