MATLAB: Numerical value for sum of bessel functions

for loopMATLABsumvectorization

I would like to find as numerical value the sum of bessel functions A0 .
Is it possible to speed it ?
x= linspace(-2,2,100); [X,XX]=meshgrid(x,x)
% X matrix,
syms k
A0 = symsum(((-1)^k).*((X).^(2*k)).*besselj(2*k,X), k, 0, 10);

Best Answer

k = 0:10;
x = 3;
A0 = sum((-1).^k.*x.^(2*k).*besselj(2*k,x))
If x is a vector, use
A0 = arrayfun(@(x)sum((-1).^k.*x.^(2*k).*besselj(2*k,x)),x)
Related Question