function Bx = exemplarCubicBSpline(x) in1 = find(x <= -2) ; Bx = zeros(1,length(in1)) ; in2 = find(x > -2 & x <= -1) ; valIn2 = (x(in2) + 2).^3 ; in3 = find(x > -1 & x <= 0) ; valIn3 = 1 + 3*(x(in3)+1) + 3*(x(in3)+1).^2 - 3*(x(in3)+1).^3 ; in4 = find(x > 0 & x <= 1) ; valIn4 = 1 + 3*(1-x(in4)) + 3*(1-x(in4)).^2 - 3*(1-x(in4)).^3 ; in5 = find(x > 1 & x <= 2) ; valIn5 = (2 - x(in5)).^3 ; in6 = find(x > 2) ; valIn6 = zeros(1,length(in6)) ; Bx = [Bx,valIn2,valIn3,valIn4,valIn5,valIn6] ;