function [mwa,densa]=mwdens(x)
global nc mw dens
densa=0;
mwa=0;
for j=1:nc
    mwa=x(j)*mw(j)+mwa;
    densa=x(j)*dens(j)+densa;
end