Given the local coordinates of a parented object, how do I get its global coordinates using the Python API?
The formula is

default_cube.matrix_world * sphere.matrix_basis


You can directly get the global coordinates (including constraints) of the sphere from 

sphere.matrix_world.translation


See Object.matrix_basis, mathutils.Matrix