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