Re: [sage-devel] Computing Resultant of Polynomials and Equation Solving

2023-08-04 Thread Emmanuel Briand
You want to eliminate variables in a system of polynomials equations. Gröbner bases are more efficient than resultants for this task. In sagemath, see especially the "elimination_ideal" method at https://doc.sagemath.org/html/en/reference/polynomial_rings/sage/rings/polynomial/multi_polynomial_idea

Re: [sage-devel] Isomorphic multi digraphs have distinct number of spanning trees, bug?

2023-06-18 Thread Emmanuel Briand
Have you taken into account the following? ignature: G1.spanning_trees_count(root_vertex=None)Docstring: Return the number of spanning trees in a graph. In the case of a digraph, counts the number of spanning out-trees rooted in "root_vertex". Default is to set first vertex as root

Re: [sage-devel] Discrepancy between integral and numerical_integral

2023-05-01 Thread Emmanuel Briand
This is not an answer but I have tried to reproduce this bug with a simpler example and I obtained this: u(x) = sqrt((1-tan(x)^2)) f(x) = pi/2-arccos(u(x)); r(x) = f(x)*cos(x) integral(r(x), (x, 0, pi/4)) 1.11072073453959 But numerical_integral(r(x), 0, pi/4) (0.785398166410623, 6.11034898000

Re: [sage-devel] RealField isn't doing it right

2023-04-17 Thread Emmanuel Briand
Interesting. The two inputs have different types type(RR(11/10)) type(RR(1.1)) Then from the documention at https://doc.sagemath.org/html/en/reference/rings_numerical/sage/rings/real_mpfr.html *class *sage.rings.real_mpfr.RealLiteral Bases: RealNumber

Re: [sage-devel] Possibly incorrect result in groebner basis over the integers

2023-03-10 Thread Emmanuel Briand
Not a bug, but a (nice) feature: Gröbner bases for rings of polynomials with integers coefficients, documented in: https://doc.sagemath.org/html/en/reference/polynomial_rings/sage/rings/polynomial/multi_polynomial_ideal.html El vie, 10 mar 2023 a las 11:57, Georgi Guninski () escribió: > I thin

Re: [sage-devel] Final call for GitHub user names

2023-01-19 Thread Emmanuel Briand
My trac username is: ebriand and my GitHub username is: EmmanuelJeanBriand Emmanuel Briand El jue, 19 ene 2023 a las 3:14, Kwankyu Lee () escribió: > Dear readers, > > Attached unmapped_users.pdf is the list of usernames encountered in the > preparation of the migration from Tr

[sage-devel] Re: Bug with a very specific integral computed with giac

2022-01-30 Thread Emmanuel Briand
duce().exponentialize().factor() >>>> -1/2*I*(Ei(x - I)*e^(2*I) - Ei(x + I))*e^(-I) >>>> sage: integrate(e^x/(x^2+1), x, algorithm="fricas") >>>> -1/2*I*(Ei(x - I)*e^(2*I) - Ei(x + I))*e^(-I) >>>> >>>> Mathematica gi

[sage-devel] Bug with a very specific integral computed with giac

2022-01-28 Thread Emmanuel Briand
s of other functions are ok. The problem does not show up when avoiding giac: sage: integrate(e^x/(x^2+1), x, -pi, pi, algorithm='maxima') integrate(e^x/(x^2 + 1), x, -pi, pi) Asking for the giac version gives no clue: sage: giac.version() "Done" Should I open a tick