Plotting the airy disk pattern in MATLAB

Question

So the formula for the airy pattern is

with J being the Bessel function besselj(1,x) in MatLab. How would i plot this and approximate it for just the x values of 0-10?

Answer 1

x = [0:0.1:10]; % kasin(phi)
I0 = 1;
I = I0*(2*besselj(1,x)./x).^2
plot(x,I)