parameter (nmax=165000) double precision x(nmax) end