for i = 0:4
    answer = bspliner (3,0,i);
    answer2 = bsplinerd (3,0,i);
    fprintf ('For i = % d, B (3,0,i) = % f and (dy/dx)B (3,0,i) = % f \
\n', i, answer,answer2);
end
 
function blah = bspliner (k,i,x) 
% fprintf ('k = % d, i = % d, x = % d \n',k,i,x);
 
if (k == 0) 
    if (i<= x && x < (i+1)) 
         blah = 1;
    else 
         blah = 0;
    end
else
 
    blah = ((x - i)./k)*bspliner (k-1,i,x) \
+(((i+k+1)-x)./((i+k+1)-(i+1)))*bspliner (k-1,i+1,x);
end
 
function blah = bsplinerd (k,i,x) 
% fprintf ('k = % d, i = % d, x = % d \n',k,i,x);
 
if (k == 0) 
    if (i<= x && x < (i+1)) 
         blah = 1;
    else 
         blah = 0;
    end
 
else
 blah = (k./k)*bsplinerd (k-1,i,x) - (k./k)*bsplinerd (k-1,i+1,x);
end