diff options
Diffstat (limited to 'solution_visualization.m')
| -rw-r--r-- | solution_visualization.m | 6 |
1 files changed, 2 insertions, 4 deletions
diff --git a/solution_visualization.m b/solution_visualization.m index 68b7a55..f5d69d4 100644 --- a/solution_visualization.m +++ b/solution_visualization.m @@ -5,15 +5,13 @@ q0=waste_r2q(w0,r0,lambda); x=linspace(x0,xf,1000); % we will calculate beam profile between x0 and xf printf('======== Forward propagation ======\n') q_forward=gbeam_propagation(x,q0,x0,optics); -w_forward=q2waste(q_forward, lambda); -r_forward=q2radius(q_forward,lambda); +[w_forward,r_forward]=q2wr(q_forward, lambda); printf('======== Backward propagation =====\n') q_backward0=waste_r2q(wf,rf,lambda); q_backward=gbeam_propagation(x,q_backward0,xf,optics); -w_backward=q2waste(q_backward, lambda); -r_backward=q2radius(q_backward,lambda); +[w_backward,r_backward]=q2wr(q_backward, lambda); printf('========= Plotting ================\n') plot (x,w_forward, '1;forward propagation;', x, w_backward, '2;backward propagation;', ... |