computation of H assuming point masse seems to work