M = 10000;
nu = linspace(-300,300,M); % frequency domain [GHz]
e0 = 769; % spacing in [micron]
R = 0.85; % per plate reflectivity
Dfp = 64; % diameter of clear aperture in [mm]
A = 0.001; % per plate absorption 1-R-T
defect = 25; % defect-limited finesse
n0 = 1.0003; % index of refraction of spacer material
lambda = 0.372; % wavelength in [micron]
etendue = pi*(0.001/2)^2*2*pi*(1-cos(asin(0.22))); % system etendue in [sr-m^2] (example: 1000 micron core, 0.22 NA)
source = 1; % 1/e half-width of source spectrum in [GHz]
tol = 1e-9; % model tolerance
T = etalon(nu,e0,R,Dfp,A,defect,n0,lambda,etendue,source,tol);
plot(nu,T,'k')
xlabel('Frequency offset (GHz)')
ylabel('Transmission')
title('Example model')
text(15,.7,sprintf('e0 = %d micron \nR = %1.2f \nD_{fpi} = %d mm \nEtendue = %1.1e sr-m^2',[e0 R Dfp etendue]),'fontsize',8)
