function BP = derbes(B)
% BP = derbes(B)
% BP = derivative of bessel function B

N = length(B);

D = [-B(2),B(1:(N-2))];

BP = (D - B(2:N))*0.5;