V[T_]:=999.87-.06426T+.0085043T^2-.0000679T^3
V'[T]
Solve[V'[T]==0,T]

Plot[V[T],{T,0,30}]
(* where the volume is a minimum, the density will be a maximum (since rho=mass/volume)
*)