double precision x_up,x_low
      common /xbound/  x_up(max_paropt),x_low(max_paropt)