implicit real*8 (a-h,o-z)
      parameter(maxvar=100)
      double precision p(maxvar+1,maxvar),pb(maxvar),y(maxvar+1)