source("abcd.m")
########################################## 
lambda= 1.064E-6 
Ltot= 0.55
r0= 1e100
w0= 0.00675
f1=0.1526
lns1.abcd=abcd_lens( f1) 
lns1.x= 0.2
f2=0.019
lns2.abcd=abcd_lens( f2) 
lns2.x= lns1.x+ f1 + f2
lns3.abcd=abcd_lens( 0.0629) 
lns3.x= 0.500
# befrore we calculate final waste we need to take in account propagation through the thick mirror
# since we try to match to the waste right in front of front coupling mirror
r_at_mirror=.01
w_at_mirror=83.729e-6
q_at_mirror=waste_r2q(w_at_mirror, r_at_mirror, lambda)
n_fused_cilica=1.44968
R_of_mirror=.01
mirror_thickness=6.35e-3
f_of_mirrr_formed_lense=1/((n_fused_cilica-1)*(1/(-R_of_mirror)))
abcd_of_mirrr_formed_lense=abcd_lens(f_of_mirrr_formed_lense)
abcd_of_propagation_through_mirror=abcd_free_space(mirror_thickness/n_fused_cilica)
abcd_of_thick_mirror=abcd_of_propagation_through_mirror*abcd_of_mirrr_formed_lense
# q right in front of front coupling mirror of the cavity 
q_infront_of_mirror=q_after_element(q_at_mirror, abcd_of_thick_mirror)
wf=q2waste(q_infront_of_mirror,lambda)
# we need minus for rf since we will solve injection to the cavity problem
rf= - q2radius(q_infront_of_mirror,lambda)
##########################################