%Analyze bandwidth
wl=[400:0.5:750];
A=W1(wl);
I1=find(A==max(A));
wl_max=wl(I1)
I=find(A>=max(A)/2);
FWHM=length(I)/2