Rewrote Symbol.coeff() and added __array__ attribute so binary operations with a numpy array work both ways.