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