real, parameter :: pi = 3.141593