以光工書里20頁左右部分的例子，將其輸入到matlab中進行光線追跡，獲得近軸計算的系統焦距。如下為程序： 9`!#5i)VU8  
rLcXo%w  
clc;clear all;close all; XlJ+:st  
Pz-=Eq  
x5,|kJ9S  
%%透鏡組結構(透鏡曲率半徑、厚度、折射率) ^?T,>ZI  
r=[26.67 189.67 -49.66 25.47 72.11 -35.00]; JB HnJm  
d=[5.20 7.95 1.6 6.7 2.8 0]; 'h 7n}  
n=[1 1.6140 1 1.6475 1 1.6140 1]; \8`^QgV`@  
w^/"j_p@  
uR82},r$m  
%%鏡面數 aWb5w  
L=length(r); _RX*Ps=  
P> ~Lx  
6H:EBj54?  
%%入射角、孔徑角、位置 e
_Q(
l'f  
i=zeros(1,L); j{P3o<l&`  
i_=zeros(1,L); cQ/T:E7$`  
u=zeros(1,L); &XSe&1  
u_=zeros(1,L); #4!6pMW(&7  
l=zeros(1,L); XU$\.g p-  
l_=zeros(1,L); e{d$OzT) V  
Y+"Gx;F>  
V,:~FufM^  
%%物方參數 5@~5RNrq2  
l(1)=-inf;                            %%物體位于無窮遠 =eG:Scoug?  
u(1)=0;                               %%平行光入射，孔徑角為零 !_h<w