As stated by pavel above you need to write out the bloch vector [u,v,w], here's a python function that returns it from the density matrix
def uvw_representation(rho):
w=rho[0,0]-rho[1,1]
u=rho[0,1]+rho[1,0]
v=-1j*(rho[1,0]-rho[0,1])
return [u.real,v.real,w.real]
Remember that this is only meant for two level systems, that is 2x2 density matrices. For more complicated systems the bloch sphere is not adequate
some other methods for those systems are explained
here