aboutsummaryrefslogtreecommitdiff
path: root/solution_visualization.m
diff options
context:
space:
mode:
Diffstat (limited to 'solution_visualization.m')
-rw-r--r--solution_visualization.m14
1 files changed, 7 insertions, 7 deletions
diff --git a/solution_visualization.m b/solution_visualization.m
index 71e9724..cfce0ab 100644
--- a/solution_visualization.m
+++ b/solution_visualization.m
@@ -2,30 +2,30 @@ trial_solution;
q0=wr2q(w0,r0,lambda);
x=linspace(x0,xf,1000); % we will calculate beam profile between x0 and xf
-printf('======== Forward propagation ======\n')
+fprintf('======== Forward propagation ======\n')
q_forward=gbeam_propagation(x,q0,x0,optics);
[w_forward,r_forward]=q2wr(q_forward, lambda);
-printf('======== Backward propagation =====\n')
+fprintf('======== Backward propagation =====\n')
q_backward0=wr2q(wf,rf,lambda);
q_backward=gbeam_propagation(x,q_backward0,xf,optics);
[w_backward,r_backward]=q2wr(q_backward, lambda);
-printf('========= Plotting ================\n')
+fprintf('========= Plotting ================\n')
plot (x,w_forward, '1;forward propagation;', x, w_backward, '2;backward propagation;', ...
x,-w_forward, '1;;', x, -w_backward, '2;;')
-printf('======= final check =============\n')
-printf('following are desired values: \n')
+fprintf('======= final check =============\n')
+fprintf('following are desired values: \n')
w0
r0
x0
wf
rf
xf
-printf('======= after propagation ========\n')
-printf('values below should match wf and rf: \n')
+fprintf('======= after propagation ========\n')
+fprintf('values below should match wf and rf: \n')
[waste_at_the_end,radius_at_the_end] = q2wr(q_forward(end), lambda);
waste_at_the_end
radius_at_the_end