Hello,
I'm trying to write a function which is calculating a B-Spline Function. Given is a knot vector U=[u₁,...u_m]. What i want to calculate know is the basis function of degree 0 which is defined by N_i,0(u)=1 if u_i <= u < u_i+1 and 0 otherwise. I tried to do this with two procs which are leading to the same problem:
procVec1 := proc (U)
m := nops(U);
N := [seq(0, i = 1 .. m-1)];
for i to m-1 do