function w = BernsteinWeights(d, j)% w = BernsteinWeights(d, j) returns a matrix of B-spline scaling function% weights for Bernstein polynomials of degree d, level j.w = eye(2^j + d);if d == 0  return;end;u = Knots(d, j);g = Greville(d, u);for i = 1:2^j-1  for r = 1:d    [u, g, w] = InsertKnot(d, u, g, w, i/2^j);  end;end;return;