How to get error out of astropy constants
While I'm not familiar with the package, a very quick look at the documentation suggests that you want

In [90]: c.M_sun.uncertainty


instead.  I've just checked and this appears to be correct.

&gt; python -c "from astropy import constants as c ; print c.M_sun.uncertainty"

5e+25