In [ ]:
#masses in kg
m = np.zeros([9,1])
m[0] = 0.33*10**24 # merc
m[1] = 4.87*10**24 # venus
m[2] = 5.97*10**24 # earth
m[3] = 0.642*10**24 # mars
m[4] = 1898*10**24 # jup
m[5] = 568*10**24 # sat
m[6] = 86.8*10**24 # uran
m[7] = 102*10**24 # nep
m[8] = 0.0166*10**24 # eris