double precision ferr, berr
double precision phtmpin(maxn), phtmpout(maxn)
double precision phxout(maxn)
double precision tmpvec1(maxn), tmpvec2(maxn)
double precision phrwork(3*maxn)