function f = Nijk_param (p, pi, pj, pk)
% 娭悢 N_{i,j,k}(p)
% p = [x;y]
    f = det([p-pj,p-pk])/det([pi-pj,pi-pk]);
end
