Revision history [back]

I think the problem is that you're trying to mix numpy and matplotlib and Sage in ways that are not supported. The documentation for plt.streamplot says:

Parameters
**********

x, y : 1D/2D arrays
   Evenly spaced strictly increasing arrays to make a grid.  If 2D,
   all rows of *x* must be equal and all columns of *y* must be equal;
   i.e., they must be as if generated by "np.meshgrid(x_1d, y_1d)".

u, v : 2D arrays
   *x* and *y*-velocities. The number of rows and columns must match
   the length of *y* and *x*, respectively.

It looks like you're trying to pass grid in the place of both x and y, but grid is a function, not a pair of arrays of values. You also define x and y using np.meshgrid, but the later assignment eqn(x,y) = K*Q / (x**2 + y**2)^(1/2) redefines x and y to be symbolic variables. I suggest instead:

xx,yy = np.meshgrid (gridx, gridy) # don't redefine x and y

and

G = [grad(x,y) for (x,y) in zip(gridx,gridy)]
Gx = [a[0] for a in G]
Gy = [a[1] for a in G]

and finally

plt.streamplot (Gx, Gy, xx,yy, color='green')

This doesn't work, either, but the error message is different. Maybe you need to do more conversion to get xx and Gx in the same format (convert both to floats? convert Gx to a numpy array?).