sigma_N=zeros(1,50)
wave_N=zeros(1,50)
for j=2:50
[a1,a2,a3,err]=diff_backscatter(355,250,j,-1,'nitrogen');
sigma_N(j)=a1;
wave_N(j)=a3;
end
sigma_O=zeros(1,50)
wave_O=zeros(1,50)
for j=2:50
[a1,a2,a3,err]=diff_backscatter(355,250,j,-1,'oxygen');
sigma_O(j)=a1;
wave_O(j)=a3;
end
sN=sigma_N;
sO=sigma_O;
figure
bar(wave_N(2:end),.78*sN(2:end),.1,'r')
hold on
bar(wave_O(2:end),0.21*sO(2:end),.1,'g')

xlim([353 355])
for j=2:length(wave_N)
    if (wave_N(j)>=353 && wave_N(j)<=355)
    text(wave_N(j),1.05*.78*sN(j),num2str(j),'Color','red','FontSize',10)
    end
    if (wave_O(j)>=353 && wave_O(j)<=355 && sO(j)~=0)
    text(wave_O(j),1.05*.21*sO(j),num2str(j),'Color','green','FontSize',10)
    end
end
