Re: [sage-support] implicit_plot3d and solve

2011-07-12 Thread Steven Pollack
Thanks. I might work on some sort of script, I suppose, because maple has an "intersectplot" function, and that's definitely something I need for this project. How hard would it be to scan through the list of solutions to "solve([nodalSet1,nodalSet2], [x,y,z], dict_solution=True)" and use som

Re: [sage-support] Re: integrate() problem

2011-07-12 Thread William Stein
On Tue, Jul 12, 2011 at 10:37 PM, robin hankin wrote: > OK so it's a bug in maxima.  Quite a bad one, I'd say!  The answer > given by maxima is unequivocal, and wrong. > > What is the procedure for getting this fixed?  Do I write to the > maxima developers or is there a mechanism to upload such re

[sage-support] Re: psi function.

2011-07-12 Thread robin
Hi again [replying to self] I didn't make myself clear here. What I meant to ask was, I think, "Look, sage doesn't seem to know this fact about psi(). Does sage store a list of known facts about the psi function anywhere? Because, if it does, I would like to suggest that the following list of

[sage-support] Re: implicit_plot3d and solve

2011-07-12 Thread Steven Pollack
Thanks for the reply. In my particular instance, there are a lot of constants, and the problem looks a bit difficult to automate. Here are the specifics: var('x1,x2,x3,x4,k,x,y,z', domain=RR) # Definition of P = Im[(x1 + i x2)^k] P(k, x1, x2, x3, x4) = (x1^2 + x2^2)^(k/2)*sin(k*arctan2(x2,x1)) #

Re: [sage-support] Re: integrate() problem

2011-07-12 Thread Michael Orlitzky
On 07/12/2011 04:37 PM, robin hankin wrote: > OK so it's a bug in maxima. Quite a bad one, I'd say! The answer > given by maxima is unequivocal, and wrong. > > What is the procedure for getting this fixed? Do I write to the > maxima developers or is there a mechanism to upload such reports > fr

Re: [sage-support] psi function.

2011-07-12 Thread Michael Orlitzky
On 07/12/2011 06:11 PM, robin hankin wrote: > Hi. > > > The psi function is defined over the entire complex plane. > > But sage doesn't seem to know anything about its properties. > > For example, I happen to know that psi(1-z) == pi/tan(pi*z) + psi(z) > > So I would expect the following sage

[sage-support] ultraspherical (Gegenbauer polynomials)

2011-07-12 Thread Steven Pollack
Hi, I know that sage has ultraspherical(n,a,x) implemented, however if a is not a number, ultraspherical(n,a,x) returns the error: NameError: name 'a' not defined (even if I write a = var('a')). This, partly, flies in the face of the fact that the Gegenbauer polynomials are functions of a. Wors

[sage-support] Re: implicit_plot3d and solve

2011-07-12 Thread achrzesz
Example; intersection of an elipsoid and sphere: var('x,y,z') solve([x^2 +y^2+z^2 ==1,x^2+y^2+2*z^2 ==1],[x,y,z]) #[[x == r1, y == -sqrt(-r1^2 + 1), z == 0], [x == r2, y == sqrt(-r2^2 + 1), z == 0]] var('r1') p1=parametric_plot3d([r1,-sqrt(-r1^2 + 1),0], (-1,1),thickness=10,color='red') p2=paramet

[sage-support] psi function.

2011-07-12 Thread robin hankin
Hi. The psi function is defined over the entire complex plane. But sage doesn't seem to know anything about its properties. For example, I happen to know that psi(1-z) == pi/tan(pi*z) + psi(z) So I would expect the following sage command: sage: (psi(1-z)- psi(z)-pi/tan(pi*z)).full_simplify()

[sage-support] implicit_plot3d and solve

2011-07-12 Thread Steven Pollack
Hi everyone, I'm new to SAGE, so I'm sorry if this is an amateur question, however I've been trying to find the simplest way to plot the intersection of two surfaces. The impression I'm under is that I "should" be able to do this with implicit_plot3d and solve. More specifically, I have two functi

Re: [sage-support] Re: integrate() problem

2011-07-12 Thread robin hankin
OK so it's a bug in maxima. Quite a bad one, I'd say! The answer given by maxima is unequivocal, and wrong. What is the procedure for getting this fixed? Do I write to the maxima developers or is there a mechanism to upload such reports from the sage community? Robin On Tue, Jul 12, 2011 at 1

[sage-support] Re: Generating extension fields

2011-07-12 Thread Maarten Derickx
The way i'd expect it to work is not working as i'd expect. sage: F.=GF(25) sage: F2. Display all 114 possibilities? (y or n) sage: F2.=F.extension(x^12-a) sage: F[b^4][b^6] Univariate Quotient Polynomial Ring in a over Univariate Quotient Polynomial Ring in a over Finite Field in a of size 5^2 w

[sage-support] Re: Maybe not a sage problem?

2011-07-12 Thread Graham Enos
I have a similar problem on Mac OS X 10.6 with Sage 4.6. I'm not what's up either, but I've just been resorting to the notebook interface when I want 3D graphics. Any help others may offer would be appreciated. Best, Graham On Jul 11, 2:17 pm, Jotace wrote: > Hi all, > > I'm using sage 4.7 on ub

[sage-support] Re: list_plot does not plot numpy arrays any more

2011-07-12 Thread Jason Grout
On 7/12/11 6:15 AM, Stan Schymanski wrote: Dear all, since I upgraded from 4.6 to 4.7 (I also moved from a Mac to Ubuntu), I cannot use a 1-D numpy array in list_plot any more. Is this intended, or is it a bug? It breaks backwards compatibility and I don't see the rationale behind it. In sear

[sage-support] Implementation of Lenstra-Lenstra-Lovesz Algorithm

2011-07-12 Thread Santanu Sarkar
I have written following Lattice Reduction Algorithm. However, it does not work properly (does not matche with the LLL algorithm implemented in Sage). It will be great for me if any one check the program. # LLL Algorithm M=matrix(ZZ,4,4, [1,57,67,75, 3,4,98,98, 34,23,267,111, 6,134,125,68 ])

[sage-support] Re: list_plot does not plot numpy arrays any more

2011-07-12 Thread achrzesz
``list_plot`` takes either a single list of data, a list of tuples, or a dictionary and plots the corresponding points. sage: list_plot(data.tolist()) A Ch -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+u

[sage-support] list_plot does not plot numpy arrays any more

2011-07-12 Thread Stan Schymanski
Dear all, since I upgraded from 4.6 to 4.7 (I also moved from a Mac to Ubuntu), I cannot use a 1-D numpy array in list_plot any more. Is this intended, or is it a bug? It breaks backwards compatibility and I don't see the rationale behind it. Here is an example (just copy and paste into a workshe

[sage-support] Re: integrate() problem

2011-07-12 Thread achrzesz
It seems that Maxima has problem here but mpmath has not: sage: from mpmath import * sage: mp.pretty=True sage: quad(lambda x:(x^2)*exp(x)/(1+exp(x))^2,[-inf,+inf]) 3.28986813369645 sage: n(pi^2/3) 3.28986813369645 A Ch -- To post to this group, send email to sage-support@googlegroups.com To uns