From 4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99 Mon Sep 17 00:00:00 2001 From: Matt Argao Date: Sun, 31 Mar 2013 22:35:49 -0400 Subject: Added descriptions to each function. --- gbeam_propagation_froward_only.m | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) (limited to 'gbeam_propagation_froward_only.m') diff --git a/gbeam_propagation_froward_only.m b/gbeam_propagation_froward_only.m index afc1316..2cdc1d8 100644 --- a/gbeam_propagation_froward_only.m +++ b/gbeam_propagation_froward_only.m @@ -1,4 +1,4 @@ -function q = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) +function [q] = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) % calculate the 'q' parameter of the Gaussian beam propagating through optical % 'optics_elements' only in the positive direction along 'x' axis at points 'x_pos' % takes the gaussian beam with initial q_in parameter at x_in @@ -8,6 +8,10 @@ function q = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) % % all x_pos must be to the right of x_in % x_pos must be monotonic! + + %Initialize q_lens (q at position of lens) + q_lens = zeros(1,size(optics_elements,2)); + if (any(x_pos < x_in)) error('all beam positions must be to the right of the x_in'); end @@ -45,8 +49,7 @@ function q = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) index = (x_last_calc <= x_pos); x=x_pos(index); q(index)= q_after_free_space(q_last_calc,x-x_last_calc); - - + end -- cgit v1.2.3