s=.095; s1=s/4; lamda=1270:1:1340;
a=1295;
b=1322;
d1=s1.*lamda.*(1-(a./lamda).^4));
d2=s1.*lamda.*(1-(b./lamda).^4));
[M11,I11]=max(d1)
[M21,I21]=max(d2)
[m11,i11]=min(d1)
[m21,i21]=min(d2)