source("abcd.m")
source("cavity.m")

lambda=.532e-6
R1=.01
R2=.01

cavity_abcd = abcd_lens(R1/2) * abcd_free_space(.007) * abcd_free_space(.01/1.8) * abcd_free_space(0.0036)

global A = cavity_abcd(1,1)
global B = cavity_abcd(1,2)
global C = cavity_abcd(2,1)
global D = cavity_abcd(2,2)

q_at_mirror=find_q_at_the_mirror(1/R1+3i)
w_at_mirror=q2waste(q_at_mirror,lambda)
R_at_mirror=q2radius(q_at_mirror,lambda)