function [uret, gret, pret] = InsertKnot(d, u, g, p, unew)

% [uret, gret, pret] = InsertKnot(d, u, g, p, unew) inserts a new knot at
% unew for B-spline scaling functions of degree d, thereby modifying knot
% vector u, Greville abscissas g, and synthesis matrix p.

uret = sort([u unew]);
gret = Greville(d, uret);
pret = PolyEval(g, p, gret);
return;