diff options
| -rw-r--r-- | solution_visualization.m | 14 |
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 |