円周率 (pi)
単精度用
real, parameter :: pi=3.14159265358979323846
倍精度用
real(8), parameter :: pi=3.14159265358979323846d0