%本程序與caculate_thermalense.m程序配合使用，本程序是第一個程序
clear;
wl=1.033e-3;          % wl為Yb:GdCOB的波長，單位為mm
n=1.74;              % Yb:GdCOB的折射率

r=50;                % r為輸出鏡的曲率半徑
l1=1.5;                % l1為晶體離入射全反鏡的距離
lc=2;                %晶體長2mm
l3=54.5;
syms f;

T1=[1,l1
    0,1];
T2=[1,lc/n/2
    0,1];
T3=[1,0
    -1/f,1];
T4=[1,lc/n/2
    0,1];
T5=[1,l3
    0,1];
T6=[1,0
    -2/r,1];

M=T3*T4*T5*T6*T5*T4*T3*T2*T1*T1*T2;
ADhalf=(M(1,1)+M(2,2))/2
w=sqrt(wl/pi*abs(M(1,2))/sqrt(-1*M(1,2)*M(2,1)-((M(1,1)-M(2,2))/2)^2))