lambda = 1E-6; % m

w0 = 0.35E-3; % m

zR = pi*w0^2/lambda;



z = linspace(0.1*zR,6*zR);

R = z.*(1+(zR./z).^2);

plot(z,R);

xlabel('Distance from waist, z [m]');

ylabel('Wavefront curvature, R(z) [m]');

grid on;

title('Gaussian beam wavefront curvature, with \lambda=1{\mu}m & w_0 = 0.35mm');