我的M文件如下： [63;8l
}  
function [G5]=fun(x); x)5v8kgf  
syms n1 n2 k1 k2 d; mf'1.{  
[subwave,nsub]=textread('subout.txt','%f %f'); X*q C:]e  
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); ]YzAcB.R  
[p,q]=size(subwave); dEL"(e#0s4  
for j=1:p; .NC}TFN|  
    wavel(j)=subwave(j,1); jU~%5R  
    n(j)=n1+n2/((wavel(j))^2); []Ea0jYu  
    k(j)=k1*exp(k2/wavel(j)); 8PS:yBkA|  
    N(j)=n(j)-i*k(j); ?R}