Added descriptions to each function.v3.0

1 files changed, 5 insertions, 33 deletions

diff --git a/self_gbeam_propagation.m b/self_gbeam_propagation.m
index 10bcf05..824be5a 100644
--- a/self_gbeam_propagation.m
+++ b/self_gbeam_propagation.m
@@ -1,11 +1,11 @@
 function [ w, w_pos ] = self_gbeam_propagation( w0, x_lens, f, x0, lambda )
-%SELF_GBEAM_PROPAGATION Summary of this function goes here
-%   Detailed explanation goes here
+%Beam propagation based on Sidney A. Self's "Focusing of spherical Gaussian
+%beams". Applied Optics, Vol. 22, Issue 5, pp. 658-661 (1983)
 
 n_lens = 3;
 
-w = w0;
-w_pos = x0;
+w = w0; %Initial waist
+w_pos = x0; %Initial waist position
 
 for i = 1:n_lens
     x_from_next_lens = x_lens(i) - w_pos(i);
@@ -14,32 +14,4 @@ for i = 1:n_lens
     w_pos(i+1) = x_lens(i) + s_new;
 end
 
-% x_pos= x_lens;
-% x(1)=w_pos(1);
-% waist_prop(1)=w(1);
-% 
-% x_temp=linspace(x0 ,x_pos(1));
-% waist_prop=[waist_prop, w(1)*sqrt(1+((x_temp-w_pos(1))/(pi*w(1)^2/lambda)).^2)];
-% x=[x x_temp];
-% 
-% for i=1:n_lens
-%     x_temp=linspace(x_pos(i),x_pos(i+1));     
-%     waist_prop=[waist_prop, w(i)*sqrt(1+((x_temp-w_pos(i))/(pi*w(i)^2/lambda)).^2)];  
-%     x=[x x_temp];
-% end
-% 
-% x_temp=linspace(x_pos(end),x_pos(end)+abs(f(end)));  
-% waist_prop=[waist_prop, w(end)*sqrt(1+((x_temp-w_pos(end))/(pi*w(end)^2/lambda)).^2)];
-% x=[x x_temp];
-% 
-% 
-% figure(1)
-% plot(x,waist_prop,'r',x,-waist_prop,'r'); 
-% 
-
-
-
-
-
-end
-
+end
\ No newline at end of file