MATLAB: Different order Bessel function in each column of an array

bessel functionbesseljMATLABvectorization

Is there a way of calculating the array J, below, without using a loop?

X = 0:0.1:20;
J = zeros(5,201);
for i = 0:4
    J(i+1,:) = besselj(i,X);
end

From the documentation, "If nu and Z are arrays of the same size, the result is also that size." Make arrays for nu and Z using meshgrid.

x = 0:0.1:20;
y = min(12,x);
plot(x,y,'-+')

In a more general case, use logical indexing, e.g.:

x = 0:0.1:20;
y = x;
y(y>12) = 12;

Why not just use FZERO?

 zroot = fzero(@(z) besselj(1, z), 3)
 besselj(1, zroot) % Should be small