lambda = 1500E-9; % m

R = 5E-2; % m



N = 1000;

L = linspace(0,R,N);

x = L/R;



w0 = sqrt(lambda*R/pi)*(x.*(1-x)).^(1/4);

w1 = sqrt(lambda*R/pi)*(x./(1-x)).^(1/4);

plot(x,w0*1e3); hold on;

plot(x,w1*1e3,'r');

xlabel('Plano-spheric resonator geometric ratio, L/R');

ylabel('Gaussian mode spot sizes [mm]');

title('Problem 1c: Beam radii of plano-spheric resonator, \lambda=1500nm, R=5cm');

grid on;

legend('Planar mirror spot size, w_0',...

       'Curved mirror spot size, w_1');