real :: a = 0.0   ! initialize A to be zero