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