There use to be a very nice gallery of examples of complex valued functions drawn using mpmath library available in Sage through SymPy, but it is too bad, the link is now broken. Maybe it can be restored by exporting the google code repository?
UPDATE: Thanks to Samuel Lelièvre and Fredrik Johansson, the gallery is back online: http://mpmath.org/gallery/
matplotlib graphs may be nice and come in handy. However, you have to be aware that they are not Sage "graphical objects" : you are on your own when it comes to using them.
Second point : the current 3D abilities of matplotlib are a bit limited. Especially, there is no, as far as I know, dynamic rendering, allowing scaling/rotation as allowed by the dynamic viewers used by Sage. What you get is a static image, similar to what Tachyon gives you.
Hi @Sébastien, the gallery is now back from brokenland, see http://mpmath.org/gallery/ -- it's a bit thanks to your answer, which prompted me to ask the mpmath list about the gallery, which made Fredrik Johansson revive it. Many thanks to Ask Sage, to @nebuckandazzer, to @Sébastien, to the mpmath list, and to Fredrik!
There is a nice complex_plot function which maps the argument of the function's value to the color of the representative point, and its modulus to tthe brightness. This may or may not be what you want.
Yon can also give a 3D representation using a colored 3D surface, mapping modulus to z and argument to color. The proble is, of course, what to do with poles ?
I have written this a couple times already, but not yet filed a ticket for that. If you're interested, I may post my code here, but don't expect miracles : it will be _s-l-o-w_...
No ticket yet. I'll keep you posted here... But I already have my load of late tickets...
The problem is that 3D surface coloring and the ability to let "holes" in them are founded on _undocumented_ abilities of our 3D infrastructure. Writing new code using them is risky if those features remain undocumented.
@Emmanuel Charpentier, open the ticket anyway, and another ticket for documenting the features.
Done : this is Trac#24335, depending on Trac#24331.
You should give more details on how is your function defined, the answer depends on this. Let me assume that you have a symbolic function like:
f(x,y) = cos(x+I*y)
You can see its norm as follows:
plot3d(lambda x,y: abs(f(x,y)),[-5,5],[-1,1])
You can see its argument as follows (note that you will have a discontinuity at -pi=pi):
plot3d(lambda x,y: arg(f(x,y)),[-5,5],[-1,1])
Asked: 2017-12-04 10:08:48 +0200
Seen: 851 times
Last updated: Dec 06 '17
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.