double precision acc
      parameter(acc=1d-4)