- curvpp(0) * curvpp(0) ) / dist;     h0d(0, 0) = (0 - curvpp(0) * curvpp(0) ) / dist;        vmat(0,0) = planee0(0);   vmat(0,0) = planee0(0);   vmat(0,0) = planee0(0);   vmat(0,0) = planee0(0);   vmat(0,0) = planee0(0);   vmat(0,0) = planee0(0);   