[sage-support] Re: Some bugs and some wishes

2012-02-20 Thread Rob Beezer
On Feb 16, 1:43 am, Manuel Kauers wrote: > 7. Nullspace for matrices over finite fields is unreasonably slow > > sage: M = MatrixSpace(GF(2^31-1), 1000, 1001).random_element(); > sage: %time M.right_kernel(); > CPU times: user 165.71 s, sys: 0.01 s, total: 165.73 s > Wall time: 166.20 s > > Math

[sage-support] Re: Add new/custom integral code to sage

2012-02-20 Thread kcrisman
On Feb 20, 2:52 pm, Burcin Erocal wrote: > > On Feb 8, 10:38 pm, ObsessiveMathsFreak > > wrote: > > > I have a certain integration result which Sage is currently unaware > > > of. I need a way to make sage aware of it in some fashion, via > > > substitution or anything else. > > > > For example

Re: [sage-support] Re: Add new/custom integral code to sage

2012-02-20 Thread Burcin Erocal
> On Feb 8, 10:38 pm, ObsessiveMathsFreak > wrote: > > I have a certain integration result which Sage is currently unaware > > of. I need a way to make sage aware of it in some fashion, via > > substitution or anything else. > > > > For example, sage currently cannot perform the following integral

[sage-support] Re: Simplification Issue Implicates Canonical Form

2012-02-20 Thread Nils Bruin
On Feb 18, 5:24 pm, Mark Rahner wrote: > converts 1/sqrt(5) to 1/5*sqrt(5) so I suspect that this issue can be > traced to the GiNaC canonical form. Yes, it does so for a very good reason: By simplifying expressions this way, you're sure to recognize equal expressions. Compare 2/ sqrt(6) and sqrt

Re: [sage-support] Keeping coefficients as floats when using solve

2012-02-20 Thread Michael Orlitzky
On 02/20/12 11:35, Burcin Erocal wrote: > > > Or you can do this: > > sage: t = -2/3*x + 4/3 > sage: t._convert(RR) > -0.667*x + 1.33 > Where were you a few weeks ago? =) Would anyone be opposed to making this a visible method? -- To post to this group, send email to

[sage-support] Re: polygon graph is too small (wrong scale)

2012-02-20 Thread Goebbe
Sorry Jason, it seems that I missed to define my Normal distribution. Here is the complete second example: NV01 = RealDistribution('gaussian',1) signi = 0.1 lu = NV01.cum_distribution_function_inv(signi/2) lo = NV01.distribution_function(lu) ru = NV01.cum_distribution_function_inv(1-(signi/2)) ro

[sage-support] Re: Coefficients of multivariate polynomials along with degrees

2012-02-20 Thread john_perry_usm
On Feb 19, 6:38 am, ObsessiveMathsFreak wrote: > Basically, I want to restrict the variables of a multivariate > polynomial to a certain set of symbols. > > B.variables() should return > [x,y] > instead of > [a,x,y] > > How can this be done. Maybe I misunderstand you, but this approach works fine

Re: [sage-support] Re: gcd's of numbers mod N

2012-02-20 Thread William Stein
On Mon, Feb 20, 2012 at 10:43 AM, Ken Ribet wrote: > I'll add here is that sage, when asked to find gcd(Mod(5,6),5), will > just echo back 5.  Of course it then gives the same answer when asked > to find gcd(Mod(11,6),5).  In number theory courses, students are told > that these quantities don't m

[sage-support] Re: Simplification Issue Implicates Canonical Form

2012-02-20 Thread kcrisman
On Feb 18, 8:24 pm, Mark Rahner wrote: > I appreciate that background info.  I hadn't tried invoking maxima > because I read somewhere that simplify() used maxima.  I must've been > reading outdated material.  As you stated, maxima does the correct > thing.  Because Sage can invoke maxima, perha

[sage-support] Re: Keeping coefficients as floats when using solve

2012-02-20 Thread kcrisman
On Feb 20, 10:39 am, Michael Orlitzky wrote: > On 02/19/12 19:58, Mike wrote: > > > When I run: > > > x,y=var('x,y', domain=RR) > > solve(2.0*x+3.0*y==4.0, y) > > > I get > > > [y == -2/3*x + 4/3] > > > but I would like to get > > > [y == -0.666*x + 1.3] > > > How can I d

[sage-support] Re: gcd's of numbers mod N

2012-02-20 Thread Ken Ribet
I'll add here is that sage, when asked to find gcd(Mod(5,6),5), will just echo back 5. Of course it then gives the same answer when asked to find gcd(Mod(11,6),5). In number theory courses, students are told that these quantities don't make any sense. Ken -- To post to this group, send email t

Re: [sage-support] Keeping coefficients as floats when using solve

2012-02-20 Thread Burcin Erocal
On Mon, 20 Feb 2012 10:39:11 -0500 Michael Orlitzky wrote: > On 02/19/12 19:58, Mike wrote: > > When I run: > > > > x,y=var('x,y', domain=RR) > > solve(2.0*x+3.0*y==4.0, y) > > > > I get > > > > [y == -2/3*x + 4/3] > > > > but I would like to get > > > > [y == -0.666*x + 1.33

Re: [sage-support] Keeping coefficients as floats when using solve

2012-02-20 Thread Michael Orlitzky
On 02/20/12 10:39, Michael Orlitzky wrote: > For a workaround, someone recently showed me this. You would call > `symbolic_approx` on your result. Whoops, you'll need this, too: from sage.symbolic.expression_conversions import Converter -- To post to this group, send email to sage-support@go

Re: [sage-support] Keeping coefficients as floats when using solve

2012-02-20 Thread Michael Orlitzky
On 02/19/12 19:58, Mike wrote: > When I run: > > x,y=var('x,y', domain=RR) > solve(2.0*x+3.0*y==4.0, y) > > I get > > [y == -2/3*x + 4/3] > > but I would like to get > > [y == -0.666*x + 1.3] > > How can I do this? > Wild guess: the float coefficients are coerced to

[sage-support] Re: polygon graph is too small (wrong scale)

2012-02-20 Thread Jason Grout
On 2/18/12 7:30 PM, Goebbe wrote: Using sage 4.8 in the notebook environment I have problems to scale polygons: For example, the following code generates a graphs that seems to be to small (not high enough): polygon2d([[0,0],[0,1],[-1,0]], rgbcolor=(1,0,1),xmin=-4, xmax=4, ymin=0, ymax=1) Use

[sage-support] Re: polygon graph is too small (wrong scale)

2012-02-20 Thread Goebbe
Thanks a lot for your answers. The figsize option did not work. The problem is that the graph is not tall enough (has not enough height). figsize just scales the whole image but keeps the wrong aspect ratio. Resizing the y-axis did not do the trick neither. Look at my first example, I already spec

[sage-support] Re: Simplification Issue Implicates Canonical Form

2012-02-20 Thread Dox
Hi Mark, I agree with you, and by experience I could say that SAGE could be a bit frustrating to learn at first... but is worth! I'm not sure if it'd help you, but notice that there are two different routines of simplification, say, simplify() and full_simpllify(). Try using the second one with y

Re: [sage-support] Flipping an equation/inequality

2012-02-20 Thread Ivan Andrus
On Feb 20, 2012, at 9:05 AM, ancienthart wrote: > Guys, is there a method or function to flip an equation or inequality? > > I'd like to turn 2*x + 3*y < 8 into 8 > 2*x + 3*y automatically. > > Joal Heagney If you just mean to reverse the order in which it is printed without changing the mathe

Re: [sage-support] gcd's of numbers mod N

2012-02-20 Thread John Cremona
On 20 February 2012 05:01, William Stein wrote: > On Sun, Feb 19, 2012 at 8:25 PM, D. S. McNeil wrote: >>> sage seems to think that the gcd of 6 and (-2 mod 6) is -2 mod 6, which it >>> converts to 4.  A mathematician would say that the gcd is 2. >>> Is this a bug, or does sage have a higher pur

[sage-support] Flipping an equation/inequality

2012-02-20 Thread ancienthart
Guys, is there a method or function to flip an equation or inequality? I'd like to turn 2*x + 3*y < 8 into 8 > 2*x + 3*y automatically. Joal Heagney -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@goog