我的M文件如下: YXFUZ9a#e
function [G5]=fun(x); r;~2NxMF/
syms n1 n2 k1 k2 d; &3_S+.JO
[subwave,nsub]=textread('subout.txt','%f %f'); usFhcU
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); fxcCz 5
[p,q]=size(subwave); A4|L;z/A[h
for j=1:p; fRow@DI\
wavel(j)=subwave(j,1); DO5H(a
n(j)=n1+n2/((wavel(j))^2); S0zk<