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