left out computeH.m