[sage-support] Re: matplotlib & animations

2013-01-20 Thread LFS
http://sage.math.canterbury.ac.nz/home/pub/246 (Sorry posted notebook and not published version.) Linda, I wasn't able to access your notebook. > -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-suppor

[sage-support] Re: matplotlib & animations

2013-01-20 Thread David Morawski
Thanks, John. That worked just fine for me. After posting, I realized PNG was my problem, but none of the other formats compatible with savefig are animate-able either, so I was at a loss. Linda, I wasn't able to access your notebook. It would be great if there were a help document/tutorial de

[sage-support] parametric_plot and xmin/xmax

2013-01-20 Thread Robert Jacobson
Setting xmin/xmax for parametric_plot doesn't seem to do anything, but ymin/ymax work as expected. What am I doing wrong? t = var('t') parametric_plot( (cos(t), sin(t)), (t, 0, 2*pi), xmin=-2, xmax=2, ymin=-2, ymax=2) -- You received this message because you are subscribed to the Google Groups

Re: [sage-support] Fast Fourier Transform

2013-01-20 Thread David Joyner
On Sun, Jan 20, 2013 at 1:49 PM, Jose Guzman wrote: > Would anyone we so kind to provide a minimal example with Sage? There are examples in fft.pyx and dft.py in the docstrings. > > > On 20 January 2013 19:48, David Joyner wrote: >> >> On Sun, Jan 20, 2013 at 12:29 PM, Santanu Sarkar >> wrot

Re: [sage-support] Fast Fourier Transform

2013-01-20 Thread Jose Guzman
Would anyone we so kind to provide a minimal example with Sage? On 20 January 2013 19:48, David Joyner wrote: > On Sun, Jan 20, 2013 at 12:29 PM, Santanu Sarkar > wrote: > > Dear all, > > We know using Fast Fourier Transform we can evaluate a polynomial f(x) > > of degree n at n points x_1, .

Re: [sage-support] Fast Fourier Transform

2013-01-20 Thread David Joyner
On Sun, Jan 20, 2013 at 12:29 PM, Santanu Sarkar wrote: > Dear all, > We know using Fast Fourier Transform we can evaluate a polynomial f(x) > of degree n at n points x_1, .., x_n in O(n) time. Is it implemented in > Sage? It is in the gsl module http://hg.sagemath.org/sage-main/file/9519a7bb2f

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Jeroen Demeyer
On 2013-01-20 16:38, Jeroen Demeyer wrote: > This looks like > http://trac.sagemath.org/sage_trac/ticket/13672 This was due to a bad call to PARI (a ring object was passed as variable name where PARI expects a string). Needs review: http://trac.sagemath.org/sage_trac/ticket/13672 -- You receive

[sage-support] Re: bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Dima Pasechnik
On 2013-01-20, Jeroen Demeyer wrote: > This looks like > http://trac.sagemath.org/sage_trac/ticket/13672 indeed: sage: K.=GF(5)[] sage: R.=K[] sage: S.=GF(5)[] sage: f=w^10+2*w^6+2*w^5+w+2 sage: S(f).factor() (y + 3)^6 * (y^4 + 2*y^3 + 4*y^2 + 3*y + 3) sage: (f-t).discriminant().factor() (4) * t^

[sage-support] Re: bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Volker Braun
On a related note, it is usually a bad idea to create polynomial rings over polynomial rings if you really want multivariate polynomials. If anything, the proper multivariate polynomials will be much faster: sage: R. = GF(5)[] sage: f = x^10+2*x^6+2*x^5+x+2-t sage: f = x^10+2*x^6+2*x^5+x+2

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Jeroen Demeyer
This looks like http://trac.sagemath.org/sage_trac/ticket/13672 -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegroups.com. To unsubscribe from this group, send email to sage-support+unsub

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Jeroen Demeyer
On 2013-01-20 16:25, Jeroen Demeyer wrote: > On 2013-01-20 15:56, Kannappan Sampath wrote: >> In Sage 5.5 and the 5.6.5c1, I get 4 as the output. > ...which is still wrong, the correct result would be 0. So we replaced > a wrong answer by a different wrong answer, progress! I meant to say, the di

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Jeroen Demeyer
On 2013-01-20 15:56, Kannappan Sampath wrote: > In Sage 5.5 and the 5.6.5c1, I get 4 as the output. ...which is still wrong, the correct result would be 0. So we replaced a wrong answer by a different wrong answer, progress! -- You received this message because you are subscribed to the Google G

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Peter Mueller
Sorry, I meant to say Sage version 5.4.1. In Sage 5.5 and the 5.6.5c1, I get 4 as the output. May be upgrade Sage? :) > I just updated to 5.5, with the same result as you. So the OLD bug was replaced by a NEW bug, because the discriminant of f-t must have the root 0, so it cannot be a nonzero c

Re: [sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Kannappan Sampath
Hello, sage: version() 'Sage Version 5.5, Release Date: 2012-12-22' On Sun, Jan 20, 2013 at 7:36 PM, Peter Mueller wrote: > I believe the following Sage code (version 4.5.1) exhibits a bug: > > sage: K.=GF(5)[] > sage: R.=K[] > sage: S.=GF(5)[] > sage: f=x^10+2*x^6+2*x^5+x+2 > sage: > sage: S(f)

[sage-support] bug in discriminant of polynomial over polynomial ring

2013-01-20 Thread Peter Mueller
I believe the following Sage code (version 4.5.1) exhibits a bug: sage: K.=GF(5)[] sage: R.=K[] sage: S.=GF(5)[] sage: f=x^10+2*x^6+2*x^5+x+2 sage: sage: S(f).factor() (y + 3)^6 * (y^4 + 2*y^3 + 4*y^2 + 3*y + 3) The code and its up to here correct results show that f is inseparable, so 0 should

[sage-support] Re: Sage versus Excel (spreadsheets in general): Are my arguments correct and complete?

2013-01-20 Thread Javier López Peña
On Sunday, January 20, 2013 6:27:16 AM UTC, Dima Pasechnik wrote: > For serious statistics work, one can just use R (which is distributed > with Sage, by the way). http://www.r-project.org/ > Or if you want to keep things pythonic just install the python data analysis library with easy_instal