When I calculate the below, I get an unexpected output, not giving me the density, but rather a value M * ..
Clear("Global`*") ;
R = UnitConvert(Quantity(2.7828000000000004`*^7, "Meters"),
"JupiterEquatorialRadius");
i == 1.5135776449524154`; (*rad*)
M == Quantity(0.080 , "JupiterMass")/Sin(1.5135776449524154) ;
(Rho) == M/(4 (Pi) R^3) // N
The out put is then:
(Rho) == M (Quantity(1.34932, 1/("JupiterEquatorialRadius")^3))
Why is M not represented as a numerical value?