高斯光束matlab仿真

%% Intensity distribution of Gaussian with paraxial approximation 
 %近轴高斯光束分布　
clc 
clear 
close all 
 
w0=40e-3;                  % waist size 束腰半径
lambda=30e-3;              % wavelength 　波长
k0=2*pi/lambda;            % wavenumber in free space　波数 
I0=1; 
zR=pi*w0^2/lambda;         %  瑞利长度
r=linspace（-200200401）*1e-3;  
z=linspace（-400400801）*1e-3; 
wz=w0*sqrt（1+（z/zR）.^2）;    % 光斑半径 
Rz=z+zR.^2./z;              % 光波波前的曲率半径
phiz=atan（z./zR）;           %Gouy phase shift 
I=zeros（length（r）length（z））; 
E=zeros（length（r）length（z））; 
 
for ii=1:le