Update Transmission_FP_f.m

Transmission_FP_f.m

@@ -1,5 +1,13 @@
 function[T,R]=Transmission_FP_f(lambda,Gain,L,n1,n2,n3)
 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% "Optoelectronic", Cambridge Books Online
+% Prof. Emmanuel Rosencher,
+% Complement to Chapter 9
+% 9.D Fabry–Perot cavities and Bragg reflectors, page 434
+% http://dx.doi.org/10.1017/CBO9780511754647.028
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
 c=2.99792458e8;             %% speed of light [m/s]
 k=2*pi/lambda;              %% wave vector    [m-1]
 w=2*pi*c/lambda;            %% pulsation      [s-1]
@@ -35,4 +43,4 @@
 end
 
 R = squeeze( abs( M(1,2,:)./M(2,2,:) ).^2 );
-T = squeeze( n3/n1./abs(M(2,2,:)).^2 );
+T = squeeze( n3/n1./abs(M(2,2,:)).^2 );