BSpline.m
041b8992f8b7master

BSpline.m
View Options

			function [ result ] = BSpline( n, x, shift, width )
			% BSpline( n, x, shift, width ) evaluates the nth order B-spline at points
			% x. The spline is centered in shift an shrinked to support [shift - width,
			% shift + width].

			x = (x - shift) / width;

			result = zeros(1, length(x));

			for k=0:n+1
			lastTerm = x - k + 0.5 * (n + 1);
			result = result + (lastTerm >= 0) .* nchoosek(n+1, k) * ((-1)^k) .* lastTerm.^n;
			end

			result = result / factorial(n);

			end