I really like this example from it: def newton(f, z, precision=0.001) : while abs(f(x=z)) >= precision: z = z - f(x=z) / diff(f)(x=z) return z
complex_plot(lambda z : newton(f, z), (-1,1), (-1,1)) I've wanted to be able to plot basins of attraction for my students easily, but didn't have time to think about the best way to do it. If you think that this solution for getting a complex basin of attraction plot is close to optimal (perhaps should create g = diff(f) so we don't keep calling Pynac?), please open a ticket and put a version of code you like for this. We could call it basins_of_attraction(f,z,prec,xrange,yrange)... or maybe not. What do you think? - kcrisman On Dec 1, 9:43 am, Sébastien Labbé <sla...@gmail.com> wrote: > Hi Sage-Devel, > > Last Monday, Franco Saliola and myself gave a talk on Sage at Montreal > Python 17. I think between 60 and 70 persons came. Franco presented > the first part with slides in English. I did a demo using the Notebook > in French. Slides are here : > > http://www.thales.math.uqam.ca/~labbes/pdf/2010-11-python-montreal-sa...http://www.thales.math.uqam.ca/~labbes/pdf/2010-11-python-montreal-sa...http://www.thales.math.uqam.ca/~labbes/pdf/2010-11-python-montreal-sa... > > We received many positive feedbacks, people really enjoyed it. > > They gave us the book *Beginning Python Visualization: Creating Visual > Transformation Scripts* which is really neat! Organizers even proposed > us to organize a coding sprint on Sage tickets with them. So if you > have any ticket a Python developer can fix, just tell me, I will > construct a list of ticket to fix for them. They have already done it > by the past on other open source projects. This could happen during > the spring. > > Sébastien > LaCIM, Montréal -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
