script transition from cartesian to cylindric issues

Hi

I wrote a code for beam propagation in a medium but for some reason the cartesian vs the cylindric is different.

I think it's somewhere in the fourier transformations but I'm not sure.

Notice in the polar that I only use 'r' since my beam is assumed to be symmetrical so to count the 'theta' I just multiply by 2pi

this is the line of how i look at the beam in the radial:

plot(abs(E_r_t(:,length(T)/2)).^2)

and it looks bad already after first propagation in air. so i assume i'm making something bad. can look at line 135 in the radial for example and insert the code above.

for examle instead of 'fft2' i'm doing fft(E_r_f,[],1) for the radial since it's a different dimention.

Attaching the two codes.