iloop=1;
g=inline('0.32*exp(-((353.35-x).^2)/(.1802)^2)');
XL=150;
XH=350;
dX=(XH-XL)/200;
X=linspace(150,350,200);
for iloop=1:length(X)
    NN2=0.;
    for nline=11:20
       [a1,a2,a3,err]=diff_backscatter(354.7,X(iloop),nline,-1,'nitrogen');
       NN2=NN2+0.78*g(a3)*a1;
    end
    NO2=0.
    for oline=15:27
        [a1,a2,a3,err]=diff_backscatter(354.7,X(iloop),oline,-1,'oxygen');
        NO2=NO2+0.21*g(a3)*a1;
    end
    Q_NN2(iloop)=NN2;
    Q_NO2(iloop)=NO2;
    Q_comb(iloop)=NN2+NO2;
end

figure
plot(X,(Q_comb).^-1)
hold on
%plot(X,1/Q_NN2,'r')
%plot(X,1/Q_NO2,'g')
legend
