Yes, that's it. I think my way of expressing a constant was much nicer than the one in the bug report, though. ;)
Thanks for pointing this out, Carl! Maybe I don't need to plot this after all! Stan Rob Beezer wrote: > This might be it. > > http://trac.sagemath.org/sage_trac/ticket/4384 > > On Mar 2, 7:47 am, Carl Witty <carl.wi...@gmail.com> wrote: > >> On Mon, Mar 2, 2009 at 7:35 AM, mabshoff <mabsh...@googlemail.com> wrote: >> >>> On Mar 2, 6:15 am, Stan Schymanski <schym...@gmail.com> wrote: >>> >>>> Dear all, >>>> >>>> I get a ZeroDivisionError: float division if I try to plot the >>>> following function: >>>> >>>> var('kab') >>>> fun = kab*((2*kab + 1)/(600*kab + 200) - (2*kab + 1)/(400*(kab + 1) + >>>> 200*kab))*((400*(kab + 1) + 200*kab)/(2*kab + 1) - (600*kab + 200)/ >>>> (2*kab + 1)) + 400*((2*kab + 1)/(400*(kab + 1) + 200*kab) - 1/400) + >>>> 200*((2*kab + 1)/(600*kab + 200) - 1/300) >>>> plot(fun,0.01,1) >>>> >>>> None of the denominators can turn 0 for positive kab, so I suppose >>>> that this is a bug. >>>> >>> It is a problem in the axis code for some reason. >>> >> In particular: >> sage: fun.simplify_full() >> 1/3 >> it's a bug in the axis code when you try to plot a constant function. >> >> Carl >> > > > -- ________________________________________ Stan Schymanski Scientist Max Planck Institute for Biogeochemistry Postfach 10 01 64 D-07701 Jena Phone: +49.3641.576264 Fax: +49.3641.577274 WWW: http://www.bgc-jena.mpg.de/~sschym Biospheric Theory and Modelling Group http://www.bgc-jena.mpg.de/bgc-theory/ _________________________________________ --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
