Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 9,783 of 10,432    |
|    antispam@math.uni.wroc.pl to clicliclic@freenet.de    |
|    Re: More teething help    |
|    25 Jan 18 17:17:41    |
      clicliclic@freenet.de wrote:       >       > Albert Rich schrieb:       > >       > > On Tuesday, January 23, 2018 at 4:40:12 AM UTC-10, anti...@m       th.uni.wroc.pl wrote:       > > >       > > > I do not see Albert's message, only this reply so here I am       > > > answering to Albert.       > > >       > > > AFAICS x/((x^3+9)*sqrt(x^3-1)) has no elementary antiderivative.       > > > I am surprised that you go that way. Examples that Martin gave       > > > are essentially diofantine phenomana depending on arithmetic       > > > of elliptic curves. ATM Rubi contains rules for handful       > > > of curves. But there are infinitely many essentially different       > > > elliptic curves and tables giving "simple" ones contain       > > > thousands of entries (google Cremona). Each such curve leads       > > > to its own family of pseudoelliptics. You may be able to       > > > find some patterns, but apparently number theorists believe       > > > that core behaviour is essentially random. So this task       > > > is like table of prime numbers, only entries have much       > > > more complicated nature.       > > >       > >       > > Sorry, I made a typo. I meant to type x/((x^3+8)*sqrt(x^3-1)), the       > > integrand of Martin's second pseudo-elliptic integral example.       > >       > > As far as the density of pseudo-elliptic integrals, Martin is better       > > equipped to address that question than I. However, I would bet that       > > only a finite number of rules are required to find elementary       > > antiderivatives for pseudo-elliptic integrands of the form       > >       > > x^m * (a+b*x^3)^n * (c+d*x^3)^(p/2)       > >       > > where m, n and p are integers.       > >       > > Note that rather than finding optimal antiderivatives for all       > > possible expressions, the more modest goal of Rubi is to find such       > > antiderivatives for all elements of well-defined classes, like the       > > form given above.       > >              _Generic_ c and d in (c+d*x^3)^(p/2) give you _single elliptic curve.       Let me recall theory: integrals as above can be reduced to sum       of three different kinds:              (c+d*x^3)^(1/2)       x*(c+d*x^3)^(1/2)       (x - a)^(-1)*(c+d*x^3)^(1/2)              Note that rules present in Rubi are doing this, so there is no       point at looking triples different than (m, n, p) = (-1, -1, -1).              First forms above (first and second kind) are no problem: single pattern       for each. Problem is with third kind. First, to have any chance       of elementary integral one have to subract apropriate multiple       of integral of first kind. Second, corrected single term is       elementary if (a, (c + d*a^3)^(1/2)) is a point of finite order       on elliptic curve. Third, modulo derivatives of elementary functions       multiple terms of form              (c + d*a_i)^(1/2)*(x - a_i)^(-1)*(c+d*x^3)^(1/2)              can be reduced to single such term. So if you split rational       part into partial fractions you are likely to loose elementary       integrability.              Generic integral doable by single logarithm is of form              f = D((a*log((p + q*del)/(p - q*del)), x)              where a is a constant, p and q are polynomials and del = sqrt(P)       where P is polynomial giving elliptic curve.              f can be written as R/del where R is rational function. Now,       the first condition is that denominator d of R is of form              d = p^2 - q^2*P              When deg(P) = 3 the simples possible case is deg(d) = 3 which       leads to constant q and linear p. So in this case checking       that d is of requested form is reasonably easy and gives you       one integrability condition. Once you know p and q you can       compute numerator of R up to multiplicative constant. In       particular for integrability your integrand must agree with the       computed form, so you get extra 3 integrbility conditions.       In principle one can give explicit formulas for all steps       and turn this into a rule or set of rules.              However, many practical examples are of nongeneric type,       where there is cancelation between numrator and denominator       of R. As I wrote for "1/(x -a)" term you need torsion       points on elliptic curve. Martin has posted formulas for       multiples a point on elliptic curve and for each k one can use       such formula to verify if "a" above corresponds to a torsion       point of order k. But to integrate you also need to find p and       q above and that is more work. Several examples include       multiple logaritms (pairs, triples, ...) and then you need       to determine residues and then check if each corresponding       points (divisors) is a tersion point.              Even for _very_ simple integrals corresponding rules are       likely to be incomplete: for curves defined over rationals       you may get points of order 12 which would imply rather       hairy rule. In practice square roots can appear for       various reasons and than instead of 12 you need to consider       higher bound.              > Thanks to the (still unpublished?) Masser-Zanier counterexample Waldek       > informed us about, we know that there are infinitely many (though       > increasingly complicated algebraic) u^2 for which       >       > INT(x/((x^2 - u^2)*SQRT(x^3 - x)), x)       >       > has a solution (also of increasing complexity) in terms of elementary       > functions. Consequently, via the usual Moebius transformation of the       > integration variable, there are also infinitely many v = (u^2 - 1)/       > (u^2 + 1) for which       >       > INT((y^2 - 1)/((y^2 + 2*v*y + 1)*SQRT(y^4 - 1)), y)       >       > has such an elementary solution. I too expect these patterns to be       > exceptions rather than the rule, the rule being that the number of       > pseudo-elliptic cases is finite and small. But does this pair already       > exhaust the store of patterns involving simple square-root radicands       > and giving rise to infinitely many pseudo-elliptic integrals?              Well, any torsion point leads to a pseudo-elliptic integral. Over       algebraically closed field there are infintely many torsion points       on any curve. So example with _single_ parameter may be exceptional,       put once you allow algebraic extentions there is see of strange       examples -- just two parameters family:              integrate((1/(x - a) - b)/sqrt(P), x)              has no elementary integral valid for continuous family of a and b       and infintely many integrable cases when P is a polynomial of       degree 3 without multiple factors.              --        Waldek Hebisch              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca