[sage-support] strange n()

2010-04-10 Thread bb
Just experimenting with Sage syntax I found something strange: sage: n(pi) 3.14159265358979 sage: n(pi,20) 3.1416 sage: n(pi,29) 3.1415927 sage: n(pi,59) 3.1415926535897932 sage: n(pi,0x59) 3.1415926535897932384626434 sage: pi.n(digits=17) 3.1415926535897932 sage: Is there any explanation?

Re: [sage-support] strange n()

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 1:02 AM, bb bblo...@arcor.de wrote: Is there any explanation? Could you be more specific in your question? Everything there looks normal to me. n(pi, 20) means to compute using 20 bits of precision. --Mike -- To post to this group, send email to

Re: [sage-support] strange n()

2010-04-10 Thread Alex Ghitza
On Sat, 10 Apr 2010 01:08:12 -0700, Mike Hansen mhan...@gmail.com wrote: On Sat, Apr 10, 2010 at 1:02 AM, bb bblo...@arcor.de wrote: Is there any explanation? Could you be more specific in your question? Everything there looks normal to me. n(pi, 20) means to compute using 20 bits of

Re: [sage-support] Re: Invoking Lisp from within Sage

2010-04-10 Thread Burcin Erocal
On Fri, 9 Apr 2010 22:51:30 -0700 (PDT) Alec Mihailovs alec.mihail...@gmail.com wrote: On Apr 10, 1:32 am, Adam Getchell adam.getch...@gmail.com wrote: Hi all, I realize this maybe a bit of an insane question, but I'm looking for a way to use ecl within sage besides: ./sage -ecl

Re: [sage-support] strange n()

2010-04-10 Thread bb
Mike Hansen schrieb: On Sat, Apr 10, 2010 at 1:02 AM, bb bblo...@arcor.de wrote: Is there any explanation? Could you be more specific in your question? Everything there looks normal to me. n(pi, 20) means to compute using 20 bits of precision. --Mike Ok, I see - the argument

[sage-support] Arrangements().list() nice format

2010-04-10 Thread bb
In deviation from the example given in the Sage Reference v4.3.5 » Combinatorics » http://www.sagemath.org/doc/reference/sage/combinat/permutation.html I get the elements in one line instead of one below the other. sage: set=['adam', 'eva', 'otto', 'emil'] sage: Arrangements(set,4).list()

[sage-support] Gap packages

2010-04-10 Thread Stochastix
I downloaded and installed sage 4.3.5 on my MacBook Pro. Everything went like a charm. I'm interested in using the simplicial homology GAP packages available at http://www.cis.udel.edu/~dumas/Homology/ It doesn't seem to be in the optional GAP packages available for SAGE. Can I install it

Re: [sage-support] Gap packages

2010-04-10 Thread David Joyner
On Sat, Apr 10, 2010 at 9:06 AM, Stochastix laurent.decreusef...@gmail.com wrote: I downloaded and installed sage 4.3.5 on my MacBook Pro. Everything went like a charm. I'm interested in using the simplicial homology GAP packages available at http://www.cis.udel.edu/~dumas/Homology/ It

[sage-support] Re: Gap packages

2010-04-10 Thread Simon King
Hi! On 10 Apr., 15:06, Stochastix laurent.decreusef...@gmail.com wrote: It doesn't seem to be in the optional GAP packages available for SAGE. Can I install it directly in SAGE or should I have a stand-alone installation of GAP ? Sage *has* a stand-alone installation of GAP. If you run sage

Re: [sage-support] Re: Problems to susbstitute a function

2010-04-10 Thread Eckhard Kosin
Hi Alec, Am Freitag, den 09.04.2010, 16:52 -0700 schrieb Alec Mihailovs: For example, var('t R_u c') map(function,('R_b', 'psi', 'm_z')) r = R_b(t)*sin(psi(t)) z = R_b(t)*cos(psi(t))+m_z(t) Dr, Dz = r.diff(t), z.diff(t) v=vector([Dr,Dz]) * vector([cos(psi(t)),sin(psi(t))])

[sage-support] Re: Problems to susbstitute a function

2010-04-10 Thread Alec Mihailovs
On Apr 10, 10:37 am, Eckhard Kosin e...@mathematik-service-kosin.de wrote: I tried to understand my mistake and I realised that I have difficulties to distinguish between symbolic expressions and functions and when to use either. The general rule is that it is better to use expressions if

Re: [sage-support] Re: Eigenvectors of a matrix

2010-04-10 Thread Mike Witt
On 04/09/2010 10:05:08 PM, Alec Mihailovs wrote: For M, you could do something like M=matrix([[1,1+i],[1-i,-1]])/sqrt(3) html.table(maxima(M).eigenvectors().sage()) or html.table([[r,(M-matrix(2,2,r)).right_kernel().basis_matrix()] for r in M.eigenvalues()]) And numerically the eigenvectors

Re: [sage-support] strange n()

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 3:37 AM, bb bblo...@arcor.de wrote: sage: n(sqrt(2.), digits=40) 1.414213562373095145474621858738828450441 sage: n(sqrt(2), digits=40) 1.414213562373095048801688724209698078570 sage: If you wanted this to be more like Maxima, the appropriate thing to do would some

[sage-support] loading maxima

2010-04-10 Thread Mathieu Roux
Hello everybody, I have just installed Sage Version 4.3.5, Release Date: 2010-03-28 on my laptop MacBook with MacOS X 10.5.8, by compiling the source code. There was no problem during the compilation. But now i can't use maxima. For exemple, if i enter f = maxima.eval('ij_entry[i,j] := i/j')

Re: [sage-support] loading maxima

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 11:09 AM, Mathieu Roux mthr...@gmail.com wrote: RuntimeError: Unable to start maxima because the command 'maxima- noreadline --userdir=/Users/mathieuroux/.sage/maxima -p / Applications/sage-4.3.5/local/bin/sage-maxima.lisp' failed. What do you get if you try running the

[sage-support] Re: loading maxima

2010-04-10 Thread Mathieu Roux
Mike thank you for your quick answer. The answer of my system is: Mathieu:~ mathieuroux$ sage -maxima --userdir=/Users/ mathieuroux/.sage/maxima -p /Applications/sage-4.3.5/local/bin/sage- maxima.lisp /Applications/sage-4.3.5/local/bin/sage-sage: line 325: maxima: command not found see you

[sage-support] Re: Arrangements().list() nice format

2010-04-10 Thread Alec Mihailovs
For example, sage: set=['adam', 'eva', 'otto', 'emil'] sage: a=Arrangements(set,4) sage: for x in range(0,a.cardinality(),2): print a[x],a[x+1] : ['adam', 'eva', 'otto', 'emil'] ['adam', 'eva', 'emil', 'otto'] ['adam', 'otto', 'eva', 'emil'] ['adam', 'otto', 'emil', 'eva'] ['adam', 'emil',

[sage-support] Re: Gap packages

2010-04-10 Thread Alec Mihailovs
On Apr 10, 9:06 am, Stochastix laurent.decreusef...@gmail.com wrote: I downloaded and installed sage 4.3.5 on my MacBook Pro. Everything went like a charm. I'm interested in using the simplicial homology GAP packages available at http://www.cis.udel.edu/~dumas/Homology/ It doesn't seem to

[sage-support] Re: Arrangements().list() nice format

2010-04-10 Thread Alec Mihailovs
Another way, sage: set=['adam', 'eva', 'otto', 'emil'] sage: a=Arrangements(set,4).list() sage: for i in zip(a[0:-1:2],a[1:-1:2]): i : (['adam', 'eva', 'otto', 'emil'], ['adam', 'eva', 'emil', 'otto']) (['adam', 'otto', 'eva', 'emil'], ['adam', 'otto', 'emil', 'eva']) (['adam', 'emil', 'eva',

[sage-support] Re: Invoking Lisp from within Sage

2010-04-10 Thread Robert Dodier
On Apr 9, 11:32 pm, Adam Getchell adam.getch...@gmail.com wrote: We've got some rather neat causal dynamical triangulation (2d quantum gravity) code running in Lisp. The lisp environment lacking certain facilities, I thought it would be interesting to find a way to run it within Sage and take

[sage-support] points on elliptic curves mod N

2010-04-10 Thread Kenneth A. Ribet
Hi, I'd like to present Lenstra's elliptic curve factoring method to a class. This means that I'd like to define an elliptic curve over Integers(N), where N is composite, and then add points on that curve in sage. I may be doing something stupid, but I'm getting a NotImplementedError with

[sage-support] Re: points on elliptic curves mod N

2010-04-10 Thread Ken Ribet
The key trick is: sage: R = Integers(2010) sage: # Make Sage think that R is a field I'm stunned. Thanks! Ken -- 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 For more options,

[sage-support] Re: Problems to susbstitute a function

2010-04-10 Thread Alec Mihailovs
On Apr 10, 10:37 am, Eckhard Kosin e...@mathematik-service-kosin.de wrote: sage: sin.diff() --- AttributeError ... There is an interesting way to fix that, sin(x)=sin(x) sin.diff() x |-- cos(x) Alec --

[sage-support] Re: points on elliptic curves mod N

2010-04-10 Thread Alec Mihailovs
On Apr 10, 8:27 pm, Robert Bradshaw rober...@math.washington.edu wrote: This particular point however seems to have order 3 on both E(GF(67))   and E(GF(89)). E=EllipticCurve([1,R(1)]) seems to be working, 7*E([0,1]) Traceback (click to the left of this block for traceback) ...