double precision g0(4)
      data g0/-1d0,-1d0,-1d0,+1d0/
      save g0