density of mixture
function rho_mix input Real parts_mass[:]; input Real rho_parts[size(parts_mass, 1)]; output Real rho; end rho_mix;