how to make these vector with polar coordinate r into a matrix with coordinate phi?

Hello, i have a simple vector in coordinate r (polar) without the coordinate phi (polar), i would like to make this a matrix that will show a circles from the middle. how to do this?

lambda=520*10^-9;
k=2*pi/lambda;
D=2.5*10-2;
R=D/2;
f=7.5*10^-2;
NA=R/f;
limit=0;
r=[-1e-6:0.00001*1e-6:-limit limit:0.00001*1e-6:1e-6];
arg=k*NA*r;
I=2*besselj(1,arg)./arg;
I=I.^2;
plot(r,I)
ylim([0 10*10^-4])
xlim([-0.5*10^-7 0.5*10^-7])
grid on
ylabel('I(r)')
xlabel('r [m]')

[dpb]

Not following what the expected plot should look like...can you post a picture that's similar (even a hand sketch would likely do)?