我的M文件如下： +hKPOFa'  
function [G5]=fun(x); }HA2ce\  
syms n1 n2 k1 k2 d; 1,;qXMhK`;  
[subwave,nsub]=textread('subout.txt','%f %f'); v^lm8/}NO  
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); 9q0,K" x)  
[p,q]=size(subwave); ;hfG${l;  
for j=1:p; hF=V
?\  
    wavel(j)=subwave(j,1); 1!v >I"]  
    n(j)=n1+n2/((wavel(j))^2);