[sage-support] Re: Differential algebra support

2009-07-30 Thread Martin Rubey
William asked me to forward his reply... (One remark: William always developed for Axiom. In Sage, the variant of Axiom usually provided is FriCAS. To the best of my knowledge, all libraries developed for Axiom is provided by FriCAS as well.) "William Sit" writes: > Dear Martin: > > I just n

[sage-support] Re: [sage-edu] Re: Black body graph is sage with interact

2009-07-30 Thread David Joyner
On Wed, Jul 29, 2009 at 11:27 PM, Minh Nguyen wrote: > > Hi Offray, > > On Thu, Jul 30, 2009 at 8:23 AM, Offray Vladimir Luna > Cárdenas wrote: >> >> Hi, >> >> Two students are interested in making the graph of the black body[1][2] >> with Sage and making it interactive in the web. [2] contains, a

[sage-support] 2d plot

2009-07-30 Thread surfer
Hi, could someone help me to plot the function def f(x,y): if x<=y: return x*(1-y) elif y<=x: return y*(1-x) for 0<=x,y<=1 in 2D? Thanks in advance. --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com T

[sage-support] Re: 2d plot

2009-07-30 Thread Harald Schilly
sage: density_plot(f, (0,1), (0,1)) looks like an adequate visualization in 2d. h On Jul 30, 5:04 pm, surfer wrote: > Hi, > > could someone help me to plot the function > > def f(x,y): >     if x<=y: >        return x*(1-y) >     elif y<=x: >        return y*(1-x) > > for 0<=x,y<=1 in 2D? > >

[sage-support] Re: 2d plot

2009-07-30 Thread surfer
> sage: density_plot(f, (0,1), (0,1)) > looks like an adequate visualization in 2d. thanks it looks pretty, but i would like to have a curve. --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, s

[sage-support] Finding the containing elliptic curve

2009-07-30 Thread VictorMiller
Suppose that q is a prime power, and I have an elliptic curve E over GF (q) (say created by E = EllipticCurve(coefficient_list)) and P,Q = E.gens() How can I find E just given P (say if I pass just P to a function)? If I say print P.parent() I get something like Abelian group of point on Ell

[sage-support] Re: Finding the containing elliptic curve

2009-07-30 Thread William Stein
On Thu, Jul 30, 2009 at 8:51 AM, VictorMiller wrote: > > Suppose that q is a prime power, and I have an elliptic curve E over GF > (q) > (say created by E = EllipticCurve(coefficient_list)) > > and P,Q = E.gens() > > How can I find E just given P (say if I pass just P to a function)? > > If I say

[sage-support] Y axis non at y=0.

2009-07-30 Thread Laurent
Hello everyone The following small example draws the function in such a way that it "suggest" that f is negative between 0.8 and 1. var('x') f(x)=x+1/sqrt(x) plot(f,(0,2)) Indeed, the X axis is drawn at y=2. Since nothing warn that behavior, it is quite disturbing. Is it willing ? I think tha

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread kcrisman
On Jul 30, 1:05 pm, Laurent wrote: > Hello everyone > > The following small example draws the function in such a way that it > "suggest" that f is negative between 0.8 and 1. > > var('x') > f(x)=x+1/sqrt(x) > plot(f,(0,2)) > > Indeed, the X axis is drawn at y=2. > Since nothing warn that behavi

[sage-support] Unsubscribe from Sage-Support

2009-07-30 Thread Craig Pellegrini
Hi, How can I unsubscribe from sage-support? -- Craig Pellegrini Department of Physics University of Illinois at Chicago (312) 996-1323 --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send

[sage-support] Re: Unsubscribe from Sage-Support

2009-07-30 Thread David Joyner
Read the bottom of the email you just posted (or any other email posted to the group). It says: "To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com" On Thu, Jul 30, 2009 at 3:11 PM, Craig Pellegrini wrote: > > Hi, > > How can I unsubscribe from sage-support?

[sage-support] Re: Unsubscribe from Sage-Support

2009-07-30 Thread William Stein
On Thu, Jul 30, 2009 at 12:11 PM, Craig Pellegrini wrote: > > Hi, > > How can I unsubscribe from sage-support? You should also be able to edit your membership here: http://groups.google.com/group/sage-support E.g., you could stay subscribed (so you can post), but elect to read all messages

[sage-support] Re: Why is R running?

2009-07-30 Thread Simon King
Hi! On Jul 28, 3:15 am, William Stein wrote: ... > > in _expect_expression. So, perhapsRisn't involved, but it is just a > > misleading warning message? > > It is a bug in the warning message indeed, in line 820. > > > Do you think I should post a ticket, changing it into "Control-C > > pressed.

[sage-support] Re: 2d plot

2009-07-30 Thread David Joyner
I don't understand the question. Do you want a contour plot? An implicit plot? On Thu, Jul 30, 2009 at 11:30 AM, surfer wrote: >> sage: density_plot(f, (0,1), (0,1)) > >> looks like an adequate visualization in 2d. > > thanks it looks pretty, but i would like to have a curve. > > > > > > --~--~

[sage-support] Re: 2d plot

2009-07-30 Thread surfer
sorry, I would like to obtain the curve for fixed x in (0,1). (or for fixed y in (0,1)). --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread Laurent
> > and that sometimes the graph might be "busy" in any particular > location one chose. Some users might also find the labeling of axes > distracting. Anyway, please let us know what your ideas would be. > * For me, it would be satisfactory to draw the X axis in grey instead of black when

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread Robert Bradshaw
On Jul 30, 2009, at 2:01 PM, Laurent wrote: > > >> >> and that sometimes the graph might be "busy" in any particular >> location one chose. Some users might also find the labeling of axes >> distracting. Anyway, please let us know what your ideas would be. >> > * For me, it would be satisfactor

[sage-support] The speed of random_element in finite fields versus Magma

2009-07-30 Thread VictorMiller
I just did a test of SAGE versus Magma on the same computer. I had a finite field GF(19991^2), and timed generating a random element in SAGE and in Magma. I found, much to my surprise, that Magma was a factor of 7 times faster. Does anyone know what method they use? Victor --~--~-~--~-

[sage-support] Re: 2d plot

2009-07-30 Thread David Joyner
I think this is with y = 1/4: sage: h = Piecewise([[(0,1/4),3*x/4],[(1/4,1),(1-x)/4]]) sage: h.plot() On Thu, Jul 30, 2009 at 4:45 PM, surfer wrote: > sorry, I would like to obtain the curve for fixed x in (0,1). (or for fixed > y in (0,1)). > > > > --~--~-~--~~~---

[sage-support] Re: The speed of random_element in finite fields versus Magma

2009-07-30 Thread Robert Bradshaw
On Jul 30, 2009, at 2:50 PM, VictorMiller wrote: > I just did a test of SAGE versus Magma on the same computer. > > I had a finite field GF(19991^2), and timed generating a random > element in SAGE and in Magma. > I found, much to my surprise, that Magma was a factor of 7 times > faster. Does an

[sage-support] Re: The speed of random_element in finite fields versus Magma

2009-07-30 Thread VictorMiller
Thanks, I found that it even gets much worse with a bigger finite field. For example if p = next_prime(100) and I try the same experiment in GF(p^2), Magma is over 100 times faster! Victor On Jul 30, 6:07 pm, Robert Bradshaw wrote: > On Jul 30, 2009, at 2:50 PM, VictorMiller wrote: > > >

[sage-support] Re: The speed of random_element in finite fields versus Magma

2009-07-30 Thread Robert Bradshaw
I'm not surprised. These fields should be written to use NTL as a backend directly (and probably using the same ideas as polynomial template, and there's enough finite fields that they should probably be put into their own directory while all the re-factoring is going on...) but no one's f

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread kcrisman
On Jul 30, 5:17 pm, Robert Bradshaw wrote: > On Jul 30, 2009, at 2:01 PM, Laurent wrote: > > > > >> and that sometimes the graph might be "busy" in any particular > >> location one chose. Some users might also find the labeling of axes > >> distracting. Anyway, please let us know what your id

[sage-support] Re: Is it possible to replace the list built-in type in Sage?

2009-07-30 Thread Offray Vladimir Luna Cárdenas
Hi, Bill Page escribió: >> [ x ] No, I can read the above just fine. It is crystal clear. >> > > ... but of course unnecessarily verbose. In my opinion a more common > notation in Sage: > > sage: x=2*vector(range(10))+vector(10*[3]) > sage: list_plot(map(lambda a:[cos(a),sin(a)],x/max(x)))

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread Jason Grout
Laurent wrote: > >> and that sometimes the graph might be "busy" in any particular >> location one chose. Some users might also find the labeling of axes >> distracting. Anyway, please let us know what your ideas would be. >> > * For me, it would be satisfactory to draw the X axis in grey in

[sage-support] Re: 2d plot

2009-07-30 Thread MaxTheMouse
On Jul 30, 10:45 pm, surfer wrote: > sorry, I would like to obtain the curve for fixed x in (0,1). (or for fixed > y in (0,1)). A general solution is to return a function of only one axis. For example, for a constant y you can do the following: def slice_y(y): def g(x): return f