function D=density2(y)
N=length(y);
for i=1:N
    h=y(i);
    if h<11000
        D(i)=101290*(1.000208275-0.2252846432e-4*h)^5.256/...
            (83007.3712-1.8696392*h);
    elseif (h>11000 & h<25000) 
        D(i)=.3629246369*exp(1.73-0.157e-3*h);
    elseif h>25000
        D(i)=2488/((.6550784855+0.1380424746e-4*h)^11.388...
            *(40875.6712+.8613592*h));
    else
        D(i)=0;
    end
end