Re: [sage-support] Matrices and Cyclotomics

2012-04-26 Thread Keshav Kini
On Friday, April 27, 2012 6:38:50 AM UTC+8, William Stein wrote: > > E.g., in sage-5.0-pre15: > > sage: sage: s=CyclotomicField(24,'s').gen() > sage: sage: (8*s^6-1)^10 > -1098715216*s^6 - 372960063 > sage: sage: xb=matrix(1,1,[8*s^6-1]) > sage: sage: xb^10 > [-1098715216*s^6 - 372960063] >

[sage-support] Solution of the system of differential equations

2012-04-26 Thread Doaa El-Sakout
Hi, I have a system of differential equations and sage give me the solution as follows, The system: x10,x20,x30,t,f1,f2,c11,c22,c12,c23,k1,k2,x1,x2,x3=var('x10,x20,x30,t,f1,f2,c11,c22,c12,c23,k1,k2,x1,x2,x3') v1=f1*(x1-k1*x2)/(1+c11*x1+c12*x2) v2=f2*(x2-k2*x3)/(1+c22*x2+c23*x3) x1= function('x1',t

Re: [sage-support] Matrices and Cyclotomics

2012-04-26 Thread William Stein
On Wed, Apr 25, 2012 at 8:45 AM, Graham Gerrard wrote: > Finding occasional inconsistencies when using matrices with cyclotomic > entries, though works well most of the time... > > sage: s=CyclotomicField(24,'s').gen() > sage: (8*s^6-1)^10 > -1098715216*s^6 - 372960063 > sage: xb=matrix(1,1,[8*s^6

Re: [sage-support] Matrices and Cyclotomics

2012-04-26 Thread JamesHDavenport
The fact that the discrepancy is 46273*46153 (both primes) makes me suspect that there's a factor of 2 missing in the CRT bounds, to allow for $\pm$. But I don't have the source here to check. On Thursday, 26 April 2012 11:19:49 UTC+1, Alastair Irving wrote: > > On 25/04/2012 16:45, Graham Gerra

Re: [sage-support] Re: Solving a set of quadratic inequalities (continuous)

2012-04-26 Thread Volker Braun
I thought the variables are real in the OP, non? Then you can split n linear inequalitiies abs(sum(p_i))>1 into 2^n cases sum(p_i))>1 or sum(p_i))<-1 On Wednesday, April 25, 2012 2:56:45 AM UTC-4, Nathann Cohen wrote: > > Hello !!! > > > Do you mean to say that you have complex num

Re: [sage-support] Re: Solving a set of quadratic inequalities (continuous)

2012-04-26 Thread Dima Pasechnik
On Thursday, 26 April 2012 23:32:01 UTC+8, Nathann Cohen wrote: > > Hell !!! > > > it should not be hard to deal with the case when you don't have > inequalities > > like |p_j-p_k|>=D, but only |p_j-p_k|<=C. > > This would make your problem convex, etc. > > Indeed, but in this ca

Re: [sage-support] Re: Solving a set of quadratic inequalities (continuous)

2012-04-26 Thread Nathann Cohen
Hell !!! > it should not be hard to deal with the case when you don't have inequalities > like |p_j-p_k|>=D, but only |p_j-p_k|<=C. > This would make your problem convex, etc. Indeed, but in this case all my problems would have a trivial solution ==> all variables equal to zero :-) N

Re: [sage-support] Re: Solving a set of quadratic inequalities (continuous)

2012-04-26 Thread Dima Pasechnik
On Wednesday, 25 April 2012 14:56:45 UTC+8, Nathann Cohen wrote: > > Hello !!! > > > Do you mean to say that you have complex numbers p_j and your > inequalities > > are of the form > > |p_j-p_k|<=C and |p_j-p_k|>=D, and that you also have > > some equations on Re(p_j) and Im(p_j) ?

[sage-support] Re: About SAGE

2012-04-26 Thread William Stein
On Apr 26, 2012 1:07 AM, "Luiz Roberto Meier" wrote: > Dear Stein, > > I have missed (!) the date of the next release of the SAGE package. I use > it in LiveCD and installed in my Fedora 16 (Sage 4.8). Do you know when > they will release the next version? I'm sure that I saw that in the website

Re: [sage-support] Matrices and Cyclotomics

2012-04-26 Thread Alastair Irving
On 25/04/2012 16:45, Graham Gerrard wrote: Finding occasional inconsistencies when using matrices with cyclotomic entries, though works well most of the time... sage: s=CyclotomicField(24,'s').gen() sage: (8*s^6-1)^10 -1098715216*s^6 - 372960063 sage: xb=matrix(1,1,[8*s^6-1]) sage: xb^10 [103692