From 12948ff5436c6c385d239572f5eb20ed19f93392 Mon Sep 17 00:00:00 2001 From: Matt Argao Date: Wed, 20 Feb 2013 14:15:44 -0500 Subject: Displays lenses in beam propogation graph. --- gbeam_propagation.m | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) (limited to 'gbeam_propagation.m') diff --git a/gbeam_propagation.m b/gbeam_propagation.m index d81f70e..eeb6a3d 100644 --- a/gbeam_propagation.m +++ b/gbeam_propagation.m @@ -1,4 +1,4 @@ -function q = gbeam_propagation(x_pos, q_in, x_in, optics_elements) +function [q, q_lens] = gbeam_propagation(x_pos, q_in, x_in, optics_elements) % calculate the 'q' parameter of the Gaussian beam propagating through optical % 'optics_elements' array along 'x' axis at points 'x_pos' % takes the gaussian beam with initial q_in parameter at x_in @@ -8,7 +8,7 @@ function q = gbeam_propagation(x_pos, q_in, x_in, optics_elements) if any(x_pos >= x_in) % Forward propagation to the right of x_in - q(x_pos >= x_in) = gbeam_propagation_froward_only(x_pos(x_pos>=x_in), q_in, x_in, optics_elements); + [q(x_pos >= x_in), q_lens] = gbeam_propagation_froward_only(x_pos(x_pos>=x_in), q_in, x_in, optics_elements); end if any(x_pos < x_in) @@ -45,8 +45,10 @@ function q = gbeam_propagation(x_pos, q_in, x_in, optics_elements) % final assignment of the backwards propagating beam % which we need to flip back q(x_pos