[sage-devel] Re: strange ECL RunTimeError when integrating

2011-12-14 Thread achrzesz
On 14 Gru, 04:36, Dan Drake wrote: > On Tue, 13 Dec 2011 at 06:38PM -0800, kcrisman wrote: > > > Proposed bugfix: symbolic integration should reject floats. > > > Yikes!   So we wouldn't want > > > sage: integrate(1.0*x,x,0,2) > > 2 > > > to work any more?  (Pace the error messages below.) > > >

[sage-devel] Re: strange ECL RunTimeError when integrating

2011-12-13 Thread achrzesz
On Dec 13, 8:27 pm, Nils Bruin wrote: > On Dec 13, 11:07 am, kcrisman wrote: > > > Yes, we try to turn it off as much as possible - > > seehttp://hg.sagemath.org/sage-main/file/9e29a3d84c48/sage/interfaces/ma... > > and the keepfloat:true.  What an unhelpful error message, then! > > I think th

[sage-devel] Re: strange ECL RunTimeError when integrating

2011-12-13 Thread achrzesz
On Dec 13, 3:48 pm, kcrisman wrote: > On Dec 12, 10:11 pm, Dan Drake wrote: > > > > > I'm doing some integrals: > > > sage: a, b, t = var('a b t') > > sage: f(a,b,t) = sin(t)^2/(a + b*cos(t))^2 > > sage: integrate(f(3/2,1,t), (t,0,2*pi)) > > -2/5*(sqrt(5) - 3)*pi*sqrt(5) > > > Okay, that's fine

[sage-devel] Re: solve an algebraic expression

2011-05-15 Thread achrzesz
Maybe the following workaround attempt will be helpful? var('x,mp,me') eq=mp-sqrt(x^2+me^2)==x assume(mp>0) solve(eq^2,x) [x == 1/2*(me^2 - mp^2)/mp, x == -1/2*(me^2 - mp^2)/mp] On 15 Maj, 01:34, Dox wrote: > Hi again list, > > This time I come with a problem of solving an algebraic expression

[sage-devel] Re: Maxima plot_list example is broken

2011-04-18 Thread achrzesz
With precision(1) also works Sage 4.6.2 64bit Fedora13 Dell Vostro1720 On 18 Kwi, 20:17, achrzesz wrote: > After replacing precision(1) by precision(53) it works for me > > On 18 Kwi, 17:36, Jeroen Demeyer wrote: > > > Executing the following "not tested" doct

[sage-devel] Re: Maxima plot_list example is broken

2011-04-18 Thread achrzesz
After replacing precision(1) by precision(53) it works for me On 18 Kwi, 17:36, Jeroen Demeyer wrote: > Executing the following "not tested" doctest from > sage/interfaces/maxima.py gives trouble, anybody has a clue what is > going wrong? > I doesn't work neither on sage-4.6.2 nor on sage-4.7.alp

[sage-devel] Re: Series expansion.

2011-02-18 Thread achrzesz
Continuing my example: sage: maxima(f).powerseries('r',infinity) -4*r*'sum((2^(2*i3-1)-1)*2^(2*i3-1)*bern(2*i3)*r^(2*i3-1)/ (2*i3)!,i3,0,inf) (NOT CHECKED) A Ch -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+uns

[sage-devel] Re: Series expansion.

2011-02-18 Thread achrzesz
sage: var('r'); sage: f=2*r/sinh(2*r) sage: f.taylor(r,0,5) 14/45*r^4 - 2/3*r^2 + 1 sage: maxima(f).powerseries('r',0) -4*r*'sum((2^(2*i2-1)-1)*2^(2*i2-1)*bern(2*i2)*r^(2*i2-1)/ (2*i2)!,i2,0,inf) Andrzej Chrzeszczyk -- To post to this group, send an email to sage-devel@googlegroups.com To unsubs

[sage-devel] Re: Definite integral on 2005 Putnam Math competition

2011-02-16 Thread achrzesz
sage: import scipy.integrate sage: scipy.integrate.quad(lambda x:log(1+x)/(x^2+1),0,1) (0.27219826128795022, 3.0220077694481673e-15) A Ch -- 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...@googlegroup

[sage-devel] Re: Definite integral on 2005 Putnam Math competition

2011-02-16 Thread achrzesz
On 16 Lut, 18:54, Marshall Buck wrote: > In the Putnam solutions by Manjul Bhargava, Kiran Kedlaya, et al, it is > observed that the definite integral > >  integrate(log(1+x)/(x^2+1),(x,0,1)) > > can be evaluated successfully by mathematica.  This is true, and the answer > is  log(2)*pi/8 > sag

[sage-devel] Problems with solve (or n() )

2011-02-10 Thread achrzesz
sage: s=solve(3*x^3-9*x+10==0,x,solution_dict=True) sage: [n(t[x]) for t in s] [1.06780542232902 - 1.84949324407141*I, # 0.0277635108030695 + 1.24902476648341*I, # WRONG! -1.09556893313209 + 0.600468477588001*I]# sage: s=solve(3*x^3-9*x +10==0,x,solution_dict=True,to_poly_solve='forc

[sage-devel] SeLinux alert

2009-11-20 Thread achrzesz
Hello Can I ignore the SeLinux alert which I obtain from time to time starting sage 4.2.1 compiled from sources in Fedora 12 on QuadCore 2.4GHz, 4GB RAM Podsumowanie: (Summary:) SELinux is preventing /home/andy/Pobrane/sage-4.2.1/local/bin/python from making the program stack executable. Szczegó

[sage-devel] Re: Fwd: sage-devel "exact" numerical integration

2008-03-05 Thread achrzesz
In the present example all data are exact and the results is OK In the previous example you have mixed symbolic and numerical data (0.1) and that was risky Andrzej On 5 Mar, 23:59, David Harvey <[EMAIL PROTECTED]> wrote: > Begin forwarded message: > > > > > From: Andrzej Chrzęszczyk <[EMAIL PROTE