Hi:
In preparing a patch to fix the bugs in special.py reported a
few days ago, I ran across the following odd feature:
sage: f = lambda x: hypergeometric_U(1,1,x)
sage: f(0.1)
2.0146425447084515
sage: f(0.2)
1.4933487469322395
sage: f(0.9)
0.63994922663929976
sage: P = plot(f,0.1,0.9)
hangs. Computes values just fine (and quickly). Hangs when it
tries to plot.
Can anyone explain this?
Please ignore the following issues for now:
1. The values above are different than the ones you will get.
(I fixed a bug - now hypergeometric_U(1,1,0.9) returns the same value
as gp.eval("hyperu(1,1,0.9)"), for example.)
2. The plot of this fuction differs from that obtained via Mma
http://functions.wolfram.com/webMathematica/FunctionPlotting.jsp?name=Hypergeometric1F1
I'll try to figure out how these two are related but that isn't the question.
- David Joyner
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/
-~----------~----~----~----~------~----~------~--~---
