04/25/2008, 11:25 PM

I found something interesting with Sage recently, and I'd thought I'd share. It seems that Sage implements complex exponentiation differently, or incorrectly, I'm not quite sure, but this plot is different than the same one produced with Mathematica, so I'm wondering if we need to start worrying about how Sage implements (^) and if this might cause problems...

Running the Sage code:

produces:

which, as you can see, is equivalent to the one made in Mathematica, except that the side with negative real part is the same as the side with positive real part. Is there any reason why this should be?

Andrew Robbins

Running the Sage code:

Code:

`contour_plot(imag((x+y*i)**(1/(x+y*i))), (x, -1, 3), (y, -2, 2), `

fill=False, contours=0, plot_points=1000)

which, as you can see, is equivalent to the one made in Mathematica, except that the side with negative real part is the same as the side with positive real part. Is there any reason why this should be?

Andrew Robbins