function beta=back_RR(temp)
ntemp=length(temp);
for n=1:ntemp
g=inline('exp(-((353.35-x).^2)/(.1802)^2)/0.3194');
    NN2=0.;
    gg=0.;
    for nline=6:30
       [a1,a2,a3,err]=diff_backscatter(354.7,temp(n),nline,-1,'nitrogen');
       gg=gg+g(a3);
       NN2=NN2+g(a3)*a1;
    end
    NN2=NN2/gg;
    NO2=0.;
    gg=0.;
    for oline=6:30
        [a1,a2,a3,err]=diff_backscatter(354.7,temp(n),oline,-1,'oxygen');
        gg=gg+g(a3);
        NO2=NO2+(0.21/.78)*g(a3)*a1;
    end
    NO2=NO2/gg;
    Q_NN2=NN2;
    Q_NO2=NO2;
    Q_comb=NN2+NO2;
    beta(n)=Q_comb;
end
end