Простая функция определения плотности воздуха

% Модель атмосферы в диапазоне от 0 до 100 км с погрешностью 1,5%
% В. А. ЯРОШЕВСКИЙ "АППРОКСИМАЦИЯ МОДЕЛИ СТАНДАРТНОЙ АТМОСФЕРЫ"
% Ученые записки ЦАГИ, т. XL, №3, 2009.
%
% density = rho(h)
% Atmospheric density
%   h [meters]
%   density [kg/m^3]
function density = rho(h)

a  = [-6.3759,-7.3012,-1.1817];
b  = [-0.4754,-0.0096,-0.0068,-0.0120,0.0042];
c  = [ 0.1803, 0.0872,-0.0153, 0.0145,0     ];


x = (h*0.001)/50.0-1; 
sa = a(1) + a(2)*x + a(3)*x*x;

%rho = sum(arrayfun(@(i) b(i)*cos(i*pi*x) + c(i)*sin(i*pi*x), 1:5));
rho = 0;
for i=1:5
    rho = rho + b(i)*cos(i*pi*x) + c(i)*sin(i*pi*x);
end

density = exp(sa + rho);

end