Computing numerical value of formula

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"), 
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?