[sage-support] Re: Numerical integration and parametic curves

tions" that avoid this, but those "special functions" have usually been created to avoid this kind of problems. Even trig functions can be seen as "special functions" curtaining e^i*x, but their geometric interpretation predates largely (by about two millenia) their analyt

[sage-support] Re: solve(-(1/2*sqrt((4*w+1)+1))*t+w==0,w)

("w,t") (w, t) sage: solve(-(1/2*sqrt((4*w+1)+1))*t+w==0,w) [w == 1/2*t*sqrt(4*w + 2)] sage: maxima.version() '5.33.0' What are you talking about, guys ? HTH, -- Emmanuel Charpentier Le jeudi 16 octobre 2014 22:49:48 UTC+2, kcrisman a écrit : > > > sage: solve(-(1/2*sq

Re: [sage-support] Re: solve(-(1/2*sqrt((4*w+1)+1))*t+w==0,w)

Le vendredi 17 octobre 2014 16:37:55 UTC+2, vdelecroix a écrit : > > 2014-10-17 10:09 UTC, Emmanuel Charpentier >: > > Ahem ! > > > > On one machine : > > > > sage: sage.version.version > > '6.4.beta4' > > sage: var("

Re: [sage-support] Re: solve(-(1/2*sqrt((4*w+1)+1))*t+w==0,w)

ma) to recognize the existence of "biradical" terms and act accordingly ? HTH, -- Emmanuel Charpentier Le vendredi 17 octobre 2014 23:00:28 UTC+2, Emmanuel Charpentier a écrit : > > > > Le vendredi 17 octobre 2014 16:37:55 UTC+2, vdelecroix a écrit : >>

[sage-support] Is that Sage vs Maxima inconsistency known ?

e reason. Similarly, setting Maxima's domain to "real" (via maxima.domain("real")) does not allow to solve the problem. One notes that, conversively, setting domain to complex in Sage's Maxima does not create a similar problem. It simply returns unsimplified expres

[sage-support] Noticed in sage 6.4rc1 : sage_mode broken

I just noticed that sage 6.4rc1 breaks sage_mode : launching sage with M-x sage never returns. Emacs is unresponsive in all of its buffers. : you have to kill emacs to get out of this mess. Now compiling rc2 to test this a bit further. HTH, -- Emmanuel Charpentier -- You received this

[sage-support] Re: Noticed in sage 6.4rc1 : sage_mode broken

A bit later : Same problem with 6.4rc2. Hints for debugging welcome... BTW : still no typeset mode possible in the new notebook, as far as I can tell. I understand that this was to be expected. HTH, -- Emmanuel Charpentier Le mercredi 12 novembre 2014 21:36:44 UTC+1, Emmanuel Charpentier a

Re: [sage-support] Re: Noticed in sage 6.4rc1 : sage_mode broken

e a look at this (tempus adjuvante...). A quick try was unsuccessfull, but I might have misinterpreted the instructions... -- Emmanuel Charpentier > -Ivan > > On Nov 12, 2014, at 2:20 PM, Emmanuel Charpentier > wrote: > > A bit later : > > Same problem with 6.4rc2. Hints

[sage-support] Re: Why solve(5^( x -1) == (0.04)^(2*x), x) returns empty set?

Why 0.04 ? Th notebook says : S=(5^( x -1) == (0.04)^(2*x)).subs({0.04:1/25}).log().solve(x) ; S [x == log(5)/(2*log(25) + log(5))] bool(S[0].rhs()==1/5) True (The last step is easily done by mental computation ; this is only a check.). HTH, -- Emmanuel Charpentier Le dimanche 16 novembre

[sage-support] Re: Why solve(5^( x -1) == (0.04)^(2*x), x) returns empty set?

by hand to get an expression easier to handle. HTH, -- Emmanuel Charpentier On Sunday, November 16, 2014 12:54:20 PM UTC-6, RRogers wrote: > > Apparently the default solver doesn't do logarithms. > For the default try: > solve(log(5^( x -1)) == log((0.04)^(2*x)), x) > >

[sage-support] Re: Why solve(5^( x -1) == (0.04)^(2*x), x) returns empty set?

Le mardi 18 novembre 2014 17:10:42 UTC+1, Chris Seberino a écrit : > > Emmanuel > > Any way to make Sage act like it can't find the solution (emit question > back to user) INSTEAD of emitting the empty set? > > "I can't find the solution" and "There is no solution" are NOT the same > thing? >

[sage-support] Re: Why solve(5^( x -1) == (0.04)^(2*x), x) returns empty set?

;true", "false" or "unknown") and uses it, so it's at least conceptually doable. And useful ! But with very deep consequences. Shouldn't this be discussed on sage-devel ? HTH, -- Emmanuel Charpentier > > Jakob > > > > Am Mittwoch

[sage-support] Re: Solving sistems of two equations

())^2, solve(eq0,psi_d)), map(lambda t:(sin(t).trig_expand())^2, solve(eq1,psi_d))) [cos(psi_d)^2 + sin(psi_d)^2 == l^2*phi^2/V_f^2 + v^2/V_f^2] leading you to the condition l^2*phi^2/V_f^2 + v^2/V_f^2==1. Beware : since we squared the previous results, this eqial

[sage-support] Re: Solving sistems of two equations

Wups : I forgot to specify that I looked for *real* solutions. If we work in complexes, the "obvious" conditions do not hold... -- Emmanuel Charpentier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from

[sage-support] Possible bug in Maxima interface

ppens that maxima.domain() exists, and the doc asserts that it's default is ... "real". So why did I get the "complex" behaviour ? It seems that trying to set Maxima's working value of domain via this interface is ... ineffective ... Is that a known problem ? HTH,

[sage-support] Re: buggy integral

I can't reproduce your problem (using the server at www.sagenb.org). I ask : var('x, p, q') assume(p,'integer'); assume(p>0) assume(q,'integer'); assume(q>0) fun = exp(-x^2) * x^(2*p) * x^q A = integral(fun, (x,-oo,oo)) B = integral(fun, (x,-oo,0)) + integral(fun, (x,0,oo)) then I ask for A : 1/

[sage-support] Symbolic functions (?) problems

I am in the process of learning Sage, coming from Maxima (and Mathematica, which I do not like much...). Cut'n'pastes from a notebook running on sagenb.org version() version() ==> 'Sage Version 5.4, Release Date: 2012-11-09' var('t,a,b,d') ## beta density dbeta(t,a,b)=t^(a-1)*(1-t)^(b-1)/bet

[sage-support] Re: Symbolic functions (?) problems

Dear Sir, Thank you for your prompt advice. Some comments below : Le vendredi 21 décembre 2012 15:03:45 UTC+1, kcrisman a écrit : > > >> I am in the process of learning Sage, coming from Maxima (and >> Mathematica, which I do not like much...). >> Cut'n'pastes from a notebook running on sagenb.o

[sage-support] Re: Symbolic functions (?) problems

Dear list, Le vendredi 21 décembre 2012 21:22:31 UTC+1, kcrisman a écrit : > > >> That's a substantial difference, IMHO. >> >> If you do a+b, then Python calls a.__add__(b). So, Python being object >> oriented, you can easily overload the a.__add__ method. Sage has the >> class sage.structure.e

Re: [sage-support] Re: Symbolic functions (?) problems

Le vendredi 21 décembre 2012 20:31:42 UTC+1, KnS a écrit : > > Emmanuel wrote: "Please let me know how to comment a ticket, and I will > report this." > > I think the usual procedure is to request a TRAC account (the details of > how to request are outlined in the SAGE TRAC homepage). > > HTH,

[sage-support] Re: Symbolic functions (?) problems

Dear Simon, Thank you for that prompt answer. I have a couple of comments below : Le samedi 22 décembre 2012 13:18:32 UTC+1, Simon King a écrit : > > Hi Emmanuel, > > On 2012-12-22, Emmanuel Charpentier > > wrote: > > So, if I follow you, Sage's add was des

[sage-support] Mathematica interface has changed in V9 ?

Note to (potential) users of the sage interface to Mathematica : something seems to have changed in Mathematica version 9 interface with "the rest of the world". Setup(s) : Debian wheezy with self-compiled sage v 5.4.1 then v 5.5, Mathematica Linux 64 bits V8 then V9. (1) sage v 5.4 <--> Mathem

Re: [sage-support] Call sage from R

I am somewhat skeptical about getting such an R package in CRAN : the dependency on Sagemath is probably a bit heavy for its platforms... and heavily platform dependent (we have serious implementation differences between Linux Mac and Windows versions). A few remarks below : Le lundi 5 avril 2

[sage-support] Trouble understanding `sage -optional` answer

Dear list Trying to list the installed optional packages, I innocently typed : charpent@zen-book-flip:~$ sage -optional | grep -v not_instal The answer left me stimyed : /usr/local/sage-9/local/lib/python3.9/site-packages/sage/misc/package.py:114: UserWarning: failed to fetch the version of

[sage-support] Re: Symbolic Fourier transform in sagemath.

This can be computed “by hand” using (one of) the textbook definition(s) : sage: var("omega, s") (omega, s) sage: integrate(sin(x^2)*e^(-I*s*x), x, -oo, oo) 1/2*sqrt(2)*sqrt(pi)*cos(1/4*s^2) - 1/2*sqrt(2)*sqrt(pi)*sin(1/4*s^2) Both sympy and giac have implementations of this transform : sage:

[sage-support] Re: memory leak in trivial calculation?

Nice one… Indeed: sage: f.expand()(y=i) -(3735/1394*I + 405/1394)*z^3 - (606/697*I - 942/697)*z^2 - (8681/6970*I + 15973/6970)*z + 1 sage: f.partial_fraction(y)(y=i) -(3735/1394*I + 405/1394)*z^3 - (606/697*I - 942/697)*z^2 - (8681/6970*I + 15973/6970)*z + 1 sage: f.simplify_full()(y=i) -(373

Re: [sage-support] Re: Symbolic Fourier transform in sagemath.

Le vendredi 28 mai 2021 à 19:01:38 UTC+2, dim...@gmail.com a écrit : > On Fri, May 28, 2021 at 5:38 PM Hongyi Zhao wrote: > > > > > > > > On Friday, May 28, 2021 at 8:19:07 PM UTC+8 Emmanuel Charpentier wrote: > >> > >> This can be c

Re: [sage-support] Re: Symbolic Fourier transform in sagemath.

Le samedi 29 mai 2021 à 01:16:21 UTC+2, hongy...@gmail.com a écrit : > On Saturday, May 29, 2021 at 1:01:38 AM UTC+8 dim...@gmail.com wrote: > >> On Fri, May 28, 2021 at 5:38 PM Hongyi Zhao wrote: >> > >> > >> > >> > On Friday, May 28, 202

[sage-support] QQbar operations extremely slow ; is this normal ?

As the first part of a demonstration on eigensystems, I was surprised to see that computing QQbar polynomial roots was way slower that computing the roots of the same polynomial expressed as a symbolic expression, or solveing it : # Relative timings of QQbar.roots and solve from time import t

[sage-support] Re: QQbar operations extremely slow ; is this normal ?

] [0 1/4 1/2*I*sqrt(23) + 1/2] HTH, ​ Le vendredi 11 juin 2021 à 09:51:46 UTC+2, Emmanuel Charpentier a écrit : > As the first part of a demonstration on eigensystems, I was surprised to > see that computing QQbar polynomial roots was way slower that computing > the roo

[sage-support] Re: QQbar operations extremely slow ; is this normal ?

h = 1.41421 [h^2-2=0], i = I > [i^2+1=0]} > > These representations are pretty similar to the SR solutions, just > displayed more compactly. You can use them for numerics: > > julia> AcbField(1000)(CP(x1)) > [+/- 4.69e-648] + [+/- 4.25e-648]*im > > julia> AcbField(100

Re: [sage-support] Re: QQbar operations extremely slow ; is this normal ?

mathematics with Sagemath*, which is more and more in order… What do you think ? ​ Le vendredi 11 juin 2021 à 21:27:05 UTC+2, vdelecroix a écrit : > Le 11/06/2021 à 20:22, Dima Pasechnik a écrit : > > On Fri, 11 Jun 2021, 18:25 Emmanuel Charpentier, < > > emanuel.c...@gmail.com> wrot

[sage-support] Possible bugs in SR

The same ask.sagemath question may have revealed two different bugs in symbolics handling. Input interpretation. Raw input, with spaces, indents and newlines : f(x) = (3/174465461165747500*pi*(-175*

[sage-support] Re: Possible bugs in SR

cision. I still think that something is wroing in SR handling, and I'll try to exhibit it in an understandable way. Stay tuned... Thanks a lot ! > On Wednesday, 7 July 2021 at 13:46:33 UTC-7 Emmanuel Charpentier wrote: > >> The same ask.sagemath question >> <https

[sage-support] Re: Possible bugs in SR

7;s handling of expressions is buggy in this case. Questions : * Does this deserve filing a ticket ? * Is so, how to file it efficiently ? Le jeudi 8 juillet 2021 à 19:39:48 UTC+2, Nils Bruin a écrit : > On Thursday, 8 July 2021 at 09:49:18 UTC-7 Emmanuel Charpentier wrote: > >> De

[sage-support] Re: Question about solving polynomial system over RR in sage

solve is a function working on symbolic expressions. Polynomials on RR as defined in your example, are *not* symbolic expressions. The set of operations applicable to them is different. And, BTW, symbolic expressions do have a *method* solve (not to be confused with the *function* solve), but

[sage-support] Re: UnicodeEncodeError: 'utf-8' codec can't encode characters in position 26-27: surrogates not allowed

I can't reproduce your problem on 9.5.beta0 installed from source on Debian testing running on core i7 + 16 GB RAM. - What version of Sage do you use ? - On which platform ? - How did-you install-it ? This information is decessary to try to understand and diagnose your problem. Le same

[sage-support] Re: UnicodeEncodeError: 'utf-8' codec can't encode characters in position 26-27: surrogates not allowed

Le dimanche 12 septembre 2021 à 11:02:07 UTC+2, Emmanuel Charpentier a écrit : > I can't reproduce your problem on 9.5.beta0 installed from source on > Debian testing running on core i7 + 16 GB RAM. > >- What version of Sage do you use ? >- On which platform ? >

[sage-support] Re: square_root_mod_prime runs indefinitely

I find the same thing. BTW, there is no square root of 12 modulo 17 : sage: R1.=Zmod(17)[] sage: %time (t^2-12).roots() CPU times: user 396 µs, sys: 2 µs, total: 398 µs Wall time: 407 µs [] Direct computation seems faster in this case, either in Zmod(17) sage: %time [t for t in Zmod(17) if t

Re: [sage-support] import

Je ne peux rien pour votre compte. En revanche, ceci devrait vous éclairer : sage: Ex = x > 3 sage: Ex.lhs() x sage: Ex.rhs() 3 sage: Ex.operator() # Kekcékçà ? sage: import_statements(Ex.operator()) from _operator import gt # Et d'ailleurs : sage: Ex.operands() [x, 3] À la vôtre ! ​ Le lundi

[sage-support] Re: Inverse of a Matrix in a Polynomial Quotient Ring

implementation of inverses in your ring in this version... -- Emmanuel Charpentier Le mardi 5 octobre 2021 à 21:52:25 UTC+2, Samanta a écrit : > Hi supporters, > I am using SageMath (version 8.9) in my Ubuntu 18.04 LTS and during the > calculation of inverse of a matrix in the quotient ring

[sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

A routine upgrade of Debian testing (where a few Python binaries were upgraded, broke a formerly functional (command line) Sage. It now crashes at startup : charpent@p-202-021:~$ sage ┌┐ │ SageMath version 9.5.beta2, Release D

Re: [sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

Le mercredi 10 novembre 2021 à 10:41:59 UTC+1, dim...@gmail.com a écrit : > it might be that you have to rebuild all the cython/python packages of > Sage. > That amounts to rebuilding from scratch, no ? > On Wed, 10 Nov 2021, 09:35 Emmanuel Charpentier, > wrote: > &g

Re: [sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

Le mercredi 10 novembre 2021 à 10:49:28 UTC+1, dim…@gmail.com a écrit : > > On Wed, 10 Nov 2021, 09:45 Emmanuel Charpentier, > wrote: > >> >> >> Le mercredi 10 novembre 2021 à 10:41:59 UTC+1, dim...@gmail.com a écrit : >> >>> it might be t

Re: [sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

org/ticket/32852#ticket>…, whose priority I’m uncertain about. Obvious workaround : ./configure --with-system-python3=no && make HTH, ​ Le mercredi 10 novembre 2021 à 12:31:06 UTC+1, Jan Groenewald a écrit : > Maybe sage -f ipython will fix it? > > > On Wed, 10 Nov 2021 at 1

Re: [sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

As reported in the ticket : charpent@zen-book-flip:~$ python --version Python 3.9.8 HTH, ​ Le mercredi 10 novembre 2021 à 19:35:20 UTC+1, dim...@gmail.com a écrit : > On Wed, Nov 10, 2021 at 5:41 PM Emmanuel Charpentier < > emanuel.c...@gmail.com> wrote:

Re: [sage-support] After (routine) upgrade of Python in Debian testing, command-line Sage crashes at startup

FWIW, I just positively reviewed the Mathias Koeppe's solution to the ticket. HTH, Le mercredi 10 novembre 2021 à 20:06:54 UTC+1, Emmanuel Charpentier a écrit : > As reported in the ticket : > > charpent@zen-book-flip:~$ python --version > Python 3.9.8 > > HTH, > ​

[sage-support] Re: subs

FWIW : sage: var("x, t") (x, t) sage: maxima_calculus.ratsubst(t,5^x, ex).sage() t^2 - 7*t + 4 sage: maxima_calculus.lratsubst([5^x==t], ex).sage() t^2 - 7*t + 4 See also this ask.sagemath.org question

[sage-support] Re: subs

Also : sage: import sympy sage: ex._sympy_().subs({(5^x)._sympy_():t._sympy_()})._sage_() t^2 - 7*t + 4 HTH, ​ Le mercredi 24 novembre 2021 à 17:20:10 UTC+1, Emmanuel Charpentier a écrit : > FWIW : > > sage: var("x, t") > (x, t) > sage: maxima_calculus.ratsubst(t,5^

[sage-support] Re: why does this not work? i am trying to create a pattern 1 2 3 4 5 6 which will repeat 6 times. i got this method from my lecturer but it is not working

Sage (i. e. Python) solution : list(range(1,9))*6 HTH, ​ Le samedi 13 novembre 2021 à 17:09:54 UTC+1, abdulra...@gmail.com a écrit : > -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails f

Re: [sage-support] I'm confused with symbolic fractions

Alternatives : sage: var("a, b, c") (a, b, c) sage: f=1/(a-b)+2/(b-c)+3/(c-a) sage: g=f*(a-b)*(b-c)*(c-a) sage: f 1/(a - b) - 3/(a - c) + 2/(b - c) sage: f.factor() (2*a^2 - 4*a*b + 3*b^2 - 2*b*c + c^2)/((a - b)*(a - c)*(b - c)) sage: g -(a - b)*(a - c)*(b - c)*(1/(a - b) - 3/(a - c) + 2/(b - c)

[sage-support] Re: Constrained optimization with strange result.

Variables of the form z_ are *integer* variables created by Maxima, which attempts to give you *also* the complex roots, if any, thus ignoring the assumptions on x, y and l. Note that : sage: solve(FOC[0], x) --- [ Sn

Re: [sage-support] Re: Constrained optimization with strange result.

tion assuming l>0 on the 3 conditions but it changes > nothing. > > > > - Mail d’origine - > De: Emmanuel Charpentier > À: sage-support > Envoyé: Mon, 29 Nov 2021 11:03:37 +0100 (CET) > Objet: [sage-support] Re: Constrained optimization with strange result. >

Re: [sage-support] Re: Constrained optimization with strange result.

all cases, only for Real > solutions. > Beware : There Might Be Tygers ! (Hint : I left you an exercise... ;-) > > - Mail d’origine - > De: Emmanuel Charpentier > À: sage-support > Envoyé: Tue, 30 Nov 2021 14:54:31 +0100 (CET) > Objet: Re: [sage-support] Re

[sage-support] A problem with the Sage solvers

ask.sagemat.org question demonstrating a problem common to all free equation solvers : solve $$ \begin{align *}-a{1}^{3} a{2} + a{1} a{2}^{2} \ -3 \, a{1}^{2} a{2} b{1} + 2 \, a{1} a{2} b{2} + a{2}

[sage-support] Re: A problem with the Sage solvers

ent in Sage 9.5.beta7). Sorry for the noise… ​ Le dimanche 5 décembre 2021 à 20:55:32 UTC+1, Emmanuel Charpentier a écrit : > ask.sagemat.org question > <https://ask.sagemath.org/question/59063/weird-c-values-from-solving-system-of-equations/> > > demonstrating a problem com

[sage-support] Transformations/functions of equalities (and possibly inequalities)

Sage can “distribute” many operations on equalities operands, such as : sage: var("a, b") (a, b) sage: (a==b)+3 a + 3 == b + 3 sage: 3*(a==b) 3*a == 3*b sage: (a==b)^3 a^3 == b^3 But not common functions : sage: log(a==b) log(a == b) sage: sin(a==b) sin(a == b) In both cases above, “distribut

Re: [sage-support] another "how to simplify" question

In Sage, this can be written wrong.maxima_methods().trigrat().expand(). HTH, ​ Le jeudi 9 décembre 2021 à 10:37:11 UTC+1, Daniel Volinski a écrit : > Hi All, > > In Maxima (embedded in SageMath) you can use: > > expand(trigrat(integrate(integrate(sin(x^2),x,y,1),y,0,1))); > > in order to get ex

[sage-support] Re: Sage crash report

ectory of the same version of the package will, by default, be deleted. Set the environment variable SAGE_KEEP_BUILT_SPKGS=yes to prevent this. make[1]: *** [Makefile:39 : all-start] Erreur 1 make[1] : on quitte le répertoire « /usr/local/sage-9 » make: *** [Makefile:13 : all] Erreur 2 ``` Suggestion

Re: [sage-support] Re: Sage crash report

r sure. > > On Fri, Dec 10, 2021 at 10:36 AM Emmanuel Charpentier > wrote: > > > > > > > > Le vendredi 10 décembre 2021 à 01:04:33 UTC+1, Matthias Koeppe a écrit : > >> > >> "make sagelib-clean" will fix this > > > > > > N

Re: [sage-support] Re: Sage crash report

anything when I stop-and-close a Jupyter notebook. Thoughts ? ​ Le vendredi 10 décembre 2021 à 16:39:53 UTC+1, Emmanuel Charpentier a écrit : > make giac-clean && make did it ! > > (*Alla cinese*) 10^4 thanks ! > ​ > Le vendredi 10 décembre 2021 à 11:56:00 UTC+1, d

[sage-support] Re: Deprecation warning when calling certain symbolic commands on Ubuntu 21.10 but not on 20.04 (SageMath 9.4)

Smells of an interface problem... Are you using a notebook, the command line or another interface (e. g. sage_shell_mode in emacs) ? If notebook, which browser do you use ? Le lundi 13 décembre 2021 à 20:31:48 UTC+1, Marcus Aichmayr a écrit : > Hi, > > I'm using SageMath 9.4 on both Ubuntu 20.0

[sage-support] Re: Deprecation warning when calling certain symbolic commands on Ubuntu 21.10 but not on 20.04 (SageMath 9.4)

Aichmayr a écrit : > It happens in command line and in jupyter notebook (firefox). > > On Monday, December 13, 2021 at 11:36:24 PM UTC+1 Emmanuel Charpentier > wrote: > >> Smells of an interface problem... >> >> Are you using a notebook, the command line or another i

[sage-support] Re: Bug in integrals

Le mercredi 15 décembre 2021 à 20:43:07 UTC+1, juanlui...@gmail.com a écrit : > See this example: > > f(x)=(x+sin(3*x))*exp(-3*x*I) > g(x)=f(x).expand() > integral(f(x)-g(x),(x,0,2*pi)) > > The answer is I*pi, but it should be 0. > Huh ? f has no poles ; therefore, the value of the integrate be

[sage-support] Re: Bug in integrals

A nice one, indeed. Here, Sage seems to use Maxima’s integrator : sage: table([[u,(f(x)-g(x)).integrate(x,algorithm=u)] for u in ["maxima", "sympy", "giac", "fricas", "mathematica_free"]], header_row=["Algorithm", "Indefinite integral"]) Algorithm Indefinite integral +--

[sage-support] Sage Crash report (again...)

Same circumstances as in the previous case : on Debian testing running on core i7 + 16 GB RAM, after upgrading Sage 9.5.beta7 compiled from a git tree  to 9.5.beta8, Sage crashes at startup. Crash report enclosed I'll try the same steps as in the last occurrence... -- Emmanuel Charpe

[sage-support] Re: Sage Crash report (again...)

n installation problem occurring on a previously-working system. Thank you for your attention. ​ Le samedi 18 décembre 2021 à 10:44:07 UTC+1, Emmanuel Charpentier a écrit : > Same circumstances as in the previous case : on Debian testing running on > core i7 + 16 GB RAM, after upgrading Sage

[sage-support] Re: Problem iniciating sage

You are trying to run Sage on a virtual machine. This Linix virtual machine requires you to log in. I note that this virtual machine runs Linux 2.6.32 on a 32-bit virtual machine, which is downright paleontologic... Are you trying to run a pre-packaged "appliance" ? Those weren't updated for a

Re: [sage-support] Re: Snowman

Homework ? Le vendredi 14 janvier 2022 à 19:41:05 UTC+1, iva.po...@gmail.com a écrit : > Sorry, here is the picture how snowman has to look like..I started drawing > with a Sphere, but I can't get it like in the picture. > > -- You received this message because you are subscribed to the Googl

Re: [sage-support] Problem in running SageMath application

Le samedi 15 janvier 2022 à 10:52:52 UTC+1, slelievre a écrit : > 2022-01-15 09:05:10 UTC, Dima Pasechnik: > > > > what is "Gilbert Source"? > > I would bet some auto-respelling software thought > that was a nice way to reinterpret "GitHub releases" > to make the conversation interesting. --Samue

[sage-support] Eigen spaces of algebraic matrices broken ?

Setup : Sage 9.5.rc1 running in Debian testing on core i7 + 16 GB RAM. def test(Size=2, Ring=QQ): from time import time as stime with seed(0): M = matrix(Ring, Size, Size, lambda u, v:Ring.random_element()) t0 = stime() SL = M.eigenspaces_left(algebraic_multiplicity=True)

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

rings/qqbar.py in pari_field(self) 3134 if self._pari_field is None: 3135 pari_pol = self._field.pari_polynomial("y") -> 3136 self._pari_field = pari_pol.nfinit(1) 3137 return self._pari_field 3138 cypari2/auto_gen.pxi in cypari2.ge

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

+1, Emmanuel Charpentier a écrit : > The promlem seems tolie wit (my use of) polynomial rings to compute the > eigenvectors. Manually, after executing > > dims = M.dimensions() > if dims[0] != dims[1]: raise DomainError("Not a square matrix !") > dim = dims[0] > BR

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

rday, 29 January 2022 at 13:51:14 UTC-8 Emmanuel Charpentier wrote: > >> /usr/local/sage-9/local/var/lib/sage/venv-python3.9/lib/python3.9/site-packages/sage/rings/qqbar.py >> >> in pari_field(self) >> >>>3134 if self._pari_field is N

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

, line 23641, in cypari2.gen.Gen_base.nfinit KeyboardInterrupt: 1.000?*v0 - 35.57125011095806?*v1 - 492.8896998473554? Similarly, V00.coefficient(v0).is_one() “never returns”. The current aritmetic on algebraics is therefore problematic for this kind of problems. Suggestions ? ​ Le

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

See lines #7-8 of the "Minimal example" : Le mercredi 2 février 2022 à 11:57:41 UTC+1, alan_thoma...@yahoo.co.uk a écrit : > > What is v0? When I run the above it isn't defined. > On Monday, January 31, 2022 at 11:19:49 PM UTC Emmanuel Charpentier wrote: > >&

[sage-support] Re: init_printing from sympy is no longer working in SageCell

WorksForMe(TM) in Sage 9.5, in both Jupyter and Jupyterlab... HTH, Le mercredi 2 février 2022 à 18:42:48 UTC+1, dsfitz...@gmail.com a écrit : > I'm teaching a linear algebra course where we use the Sympy Python package > for a lot of the computations. This includes a PreTeXt textbook where there

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

Le mercredi 2 février 2022 à 22:15:00 UTC+1, Nils Bruin a écrit : On Monday, 31 January 2022 at 15:19:49 UTC-8 Emmanuel Charpentier wrote: > >> As advertised, an atempt at a minimal (non-)working example : >> >> # Reproducible minimal example >> with seed(0): M =

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

, February 3, 2022 at 6:44:47 AM UTC Emmanuel Charpentier wrote: > >> Le mercredi 2 février 2022 à 22:15:00 UTC+1, Nils Bruin a écrit : >> >> On Monday, 31 January 2022 at 15:19:49 UTC-8 Emmanuel Charpentier wrote: >>> >>>> As advertised, an atempt at a mi

[sage-support] Re: Eigen spaces of algebraic matrices broken ?

That was my initial complaint... ;-) Le samedi 5 février 2022 à 18:05:42 UTC+1, alan_thoma...@yahoo.co.uk a écrit : > M.eigenvalues() never returns. > On Saturday, February 5, 2022 at 11:48:47 AM UTC Emmanuel Charpentier > wrote: > >> What exactly fails in the example ? &

Re: [sage-support] Re: Eigen spaces of algebraic matrices broken ?

Nins and me got *different* random matrices : his was composed of integers in (-2..2), and Sage could compute its eigenspaces ; mine was : sage: with seed(0): M = matrix(AA, 3, 3, lambda u,v: AA.random_element()) sage: M.apply_map(lambda u:u.radical_expression()) [ -sqrt(2) - 1

[sage-support] Inter-versions reproducibility problem of random elements.

Seen in this thread : On Sagecell : print(sage.version.version) with seed(0): M = matrix(AA, 3, 3, lambda u,v: AA.random_element()) M.apply_map(lambda u:u.radical_expression()) prints 9.4 [-2 2 -2] [-2

[sage-support] Re: Inter-versions reproducibility problem of random elements.

Also : coud you report the results on as many platforms and/or versions as possible ? Le dimanche 6 février 2022 à 20:23:39 UTC+1, Emmanuel Charpentier a écrit : > Seen in this thread > <https://groups.google.com/g/sage-support/c/IvjMhqryRQs> : > > On Sagecell <https://

[sage-support] Re: mathematica_free error

I have had similar problems. Wolfram may have changed (again !) something in their output format... Since I know zilch about HTML mysteries and miseries, I can't offer anything but my warmest condolences... Le lundi 14 février 2022 à 17:17:20 UTC+1, rodrigos...@gmail.com a écrit : > Good aftern

Re: [sage-support] nonlinear equation system

FWIW, executing : reset() # Don't scratch Sage's predefined identifiers, for sanity's sake... Vars= var('A B EE F II J RR T') eq1 = A*EE-B^2-B*F+EE^2==1 eq4 = A*II-B*J+II^2+RR^2==-1/2 eq5 = A*RR-B*T+2*RR*II==0 eq6 = B*II-EE*J+II*J+RR*T==0 eq8 = -B*RR+EE*T-RR*J-II*T==0 eq9 = EE*II-F*J+J^2+T

Re: [sage-support] nonlinear equation system

+ T^2 == 1/2, -(EE*RR) + F*T - 2*J*T == 0, II^2 - J^2 - RR^2 + T^2 == -1} sage: mathematica("Vars = {%s}"%", ".join([u._mathematica_init_() for u in Vars])) {A, B, EE, F, II, J, RR, T} sage: mathematica("Reduce[Sys, Vars]") False HTH, ​ Le jeudi 3 mars 2022 à 15:0

Re: [sage-support] nonlinear equation system

Or (unexpectedly much simpler) : ``` sage: mathematica.Reduce(Sys, Vars) False ``` Le jeudi 3 mars 2022 à 15:09:27 UTC+1, Emmanuel Charpentier a écrit : > And, BTW : > > sage: mathematica("Sys = {%s}"%", ".join([u._mathematica_init_() for u in > Sys])) > {

Re: [sage-support] Re: 2 questions on var

even more simpler :-) : sage: V=var("v", n=8) sage: V (v0, v1, v2, v3, v4, v5, v6, v7) sage: v2 v2 sage: V[2] v2 *“Who could ask for anything more ?”* ​ Le mardi 8 mars 2022 à 08:20:07 UTC+1, slelievre a écrit : > Even more practical, I find, is to name the tuple of indexed variables: > ``

[sage-support] Re: search_def(), search_src() are not working in Sage

How did you install Sage ? And on what platform ? Le mercredi 23 mars 2022 à 06:37:20 UTC+1, adarsh.k...@gmail.com a écrit : > Hello everyone, > I was going through some of the Sage functions and I wanted to look it up > in the source file definitions. For this I tried search_def() and then > s

Re: [sage-support] Solving logarithmic equations

BTW : sage: solve(1+2*log(x+1, 4)==2*log(x,2), x, to_poly_solve="force") [x == -sqrt(3) + 1, x == sqrt(3) + 1] HTH, ​ Le samedi 9 avril 2022 à 11:34:10 UTC+2, wdjo...@gmail.com a écrit : > On Sat, Apr 9, 2022 at 5:18 AM Paolo Robillos > wrote: > > > > Hi, > > > > I am trying to solve the fol

[sage-support] Re: \ZZ not defined

This means that Mathjax, which is used by Jupyter to interpret LaTeX expressions, has no definition of a `\ZZ` LaTeX macro. On 9.6.rc0, I get `\Bold𝑍/6\Bold𝑍`, which is not interpreted by `mathjax` either. I don't know how to work around this with Mathjax. Would you mind filing a ticket ? Le l

Re: [sage-support] Unexpected result for cos(a)^2 + sin(a)^2, question about variables

Also : sage: b.simplify_trig() 1 Using specific simplifications in a specific order is often the key to get interesting results that the brute-force simplify_full cannot. Such simplifications are : sage: import re sage: print(", ".join([v[0] for v in list(map(lambda u:re.findall(".*simplify

[sage-support] Re: Solving system of Linear Equatin over finite fields.

Your example has several problems : 1) You don’t define your polynomial indeterminates ; you should Rx.inject_variables(). 2) The syntax you use to substitute values in f is questionable… 3) f(v) is a polynomial in x0..x9 over GF(7), *not* a symbolic expression. Therefore f(v)==1 is *not* a

[sage-support] Sagetex version ?

I do not understand my installation of sagetex. As far as I understand it, Sage will use the files accessible via $SAGE_ROOT/venv/share/texmf/tex/latex/sagetex, and that’s what I symlinked to my LaTeX setup (in /usr/local/share/tex/latex/ in a Debian testing installation). However, this vers

[sage-support] Somethog in rotten is the state of Maxima...

... and all is not fresh in Sympy's realm either. Full details in [this ask.sagematjh.org question](https://ask.sagemath.org/question/64344/solving-a-system-of-linear-equations-with-complex-numbers-yields-false-solution/?answer=64363#post-id-64363). TL;DR : a linear system, perfectly solved by S

[sage-support] Re: Somethog in rotten is the state of Maxima...

It turns out that the original author has also posted this problem in sage-devel <https://groups.google.com/g/sage-devel/c/6jIKV1hPoCQ>, wher we should continue… ​ Le dimanche 9 octobre 2022 à 18:27:54 UTC+2, Nils Bruin a écrit : > On Sunday, 9 October 2022 at 08:45:36 UTC-7

[sage-support] Re: Plotting heat maps of scalar fields in SAGE

That can be done in Sage in a variety of ways . Here’s one : var("x, y") L = 3# Plotted function f =lambda x,y:cos(x)-2*y# Coloring# Colormap cm=colormaps["RdBu"]# We have to scale the colormap :# Range of values : possible shortcut via some analytical obviousnesses :# cmin = f(-3, 3).n()# cma

[sage-support] Re: Somethog in rotten is the state of Maxima...

s your execution *especially* how to convert the result back to Sage ? Martin > > On Sunday, 9 October 2022 at 18:43:03 UTC+2 Emmanuel Charpentier wrote: > >> It turns out that the original author has also posted this problem in >> sage-devel <https://groups.google.com/g/s

[sage-support] Re: Somethog in rotten is the state of Maxima...

age: Chk0f True You are right…. ​ Le jeudi 3 novembre 2022 à 16:51:31 UTC+1, Emmanuel Charpentier a écrit : > Le jeudi 3 novembre 2022 à 12:59:48 UTC+1, axio…@yahoo.de a écrit : > > @Emmanuel, why are you saying that FriCAS returns the same as sage? >> > I tried solve(Sys0, IVars

[sage-support] Re: Plotting the solution returned by desolve

Try this : sage: var("x, t") ## t will be used later (x, t) sage: y=function("y") ## Note : no default argument sage: sol=desolve(diff(y(x),x)==(x*y(x)^2 - cos(x)*sin(x) )/(y(x)*(1 - x^2)) ,y(x), ics=[0, 2]) ; sol ## Note : specify y argument 1/2*(x^2 - 1)*y(x)^2 + 1/2*cos(x)^2 == (-3/2) sage:

<    1   2   3   4   5   >