call HESSUPD( nn, nupdate, dd, jmin, jmax, xdiff,
subroutine HESSUPD( nn, mupd, dd, jmin, jmax, xdiff, lphprint )