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

2010-04-10 Thread Alec Mihailovs
On Apr 10, 8:27 pm, Robert Bradshaw 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) ... ZeroDivisionError: Inverse of 5092 does

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

2010-04-10 Thread Alec Mihailovs
On Apr 10, 10:37 am, Eckhard Kosin wrote: > > sage: sin.diff() > --- > AttributeError ... There is an interesting way to fix that, sin(x)=sin(x) sin.diff() x |--> cos(x) Alec -- To post to this group, send

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

2010-04-10 Thread Robert Bradshaw
On Apr 10, 2010, at 5:08 PM, Kenneth A. Ribet wrote: 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 somet

[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,

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

2010-04-10 Thread William Stein
On Sat, Apr 10, 2010 at 5:08 PM, Kenneth A. Ribet wrote: > 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 doin

[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 the

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

2010-04-10 Thread Robert Dodier
On Apr 9, 11:32 pm, Adam Getchell 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 advantage of all th

[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: Gap packages

2010-04-10 Thread Alec Mihailovs
On Apr 10, 9:06 am, Stochastix 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 be in the optional GAP pa

[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', 'ev

[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 Mat

Re: [sage-support] loading maxima

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 11:09 AM, Mathieu Roux 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 following

[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') th

Re: [sage-support] strange n()

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 3:37 AM, bb 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 something like:

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 c

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

2010-04-10 Thread Alec Mihailovs
On Apr 10, 10:37 am, Eckhard Kosin 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 possible instead of functions. Th

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))]) > w=v.s

[sage-support] Re: Gap packages

2010-04-10 Thread Simon King
Hi! On 10 Apr., 15:06, Stochastix 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 -gap" then this installation

Re: [sage-support] Gap packages

2010-04-10 Thread David Joyner
On Sat, Apr 10, 2010 at 9:06 AM, Stochastix 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 be in the op

[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 directl

[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() [['ada

Re: [sage-support] strange n()

2010-04-10 Thread bb
Mike Hansen schrieb: On Sat, Apr 10, 2010 at 1:02 AM, bb 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 defines the number

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 wrote: > > On Apr 10, 1:32 am, Adam Getchell 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 > > For example, > > lisp("(def x

Re: [sage-support] strange n()

2010-04-10 Thread Alex Ghitza
On Sat, 10 Apr 2010 01:08:12 -0700, Mike Hansen wrote: > On Sat, Apr 10, 2010 at 1:02 AM, bb 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. ... which you can fi

Re: [sage-support] strange n()

2010-04-10 Thread Mike Hansen
On Sat, Apr 10, 2010 at 1:02 AM, bb 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 sage-support@googlegroups.com To unsu

[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? Rega