Change subs().evalf() to evalf( subs={} )

This eliminates the need for declaring numeric variables rho_max,
u_max, u_star and then doing an additional substitution. It's also
recommended way to evalf.