real :: pi = 4.*atan(1.)
real :: tau

include "inc2.f90"