Repository navigation
Integrate a rational function of x/sqrt(P) for an even quartic P, high powers of two linears, and nine more algebraic shapes - #2023
Merged
Conversation
…h powers of two linears, and nine more algebraic shapes Eleven shapes of Rubi's family 1, each answered where the rules declined or ran out of time: - a rational function beside a half-odd power of P = a + b x^2 + c x^4, where it is one of t = x/sqrt(P), since dt/dx = (a - c x^4)/P^(3/2); - (c x)^t, t not whole, as x^t beside (c x)^t x^(-t), constant wherever defined; - two roots whose radicands' quotient is a constant times an even power of x, as one root up to a factor whose square is constant; - a polynomial over high whole powers of two linears, under t = L1/L2; - an algebraic function of 1/x alone, under x = 1/u; - x^m (a + b x^n)^p, m and n whole and p <= -2, by the reduction to p = -1; - a power of a sum of powers of x beside a symbolic power of a binomial, multiplied out into binomial differentials; - a difference of radicands that is one power of x, written as one; - a rational function beside an odd root of a linear past the square, split first where a linear below meets the root at an irrational value; - an odd integrand beside a high power of a quadratic left to u = x^2, and a quartic in one shifted linear to u = d + k x. Chain against master f7f2a7a: pocket of 1,079 rows of family 1, 1,013 -> 1,038 solved, 0 wrong; samples unchanged but family 1's 40-a-file, 1,375 -> 1,379; the 25 rows that moved, alone, 0 -> 25 solved, each in 18-772 ms; no row's allocation grew past 1.3x. x^2/((1 - a x)^16 (1 + a x)^11), which master answers in 20 s and 8.5 GB as 16,292 characters (exactly right; the probe's double evaluation calls it wrong), is 1,380 in 0.12 s. Part of #233. Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Part of #233.
Eleven shapes of Rubi's family 1, each answered where the rules declined or ran out of time. One test class per shape.
P = a + b x^2 + c x^4, where it is a rational function oft = x/sqrt(P). Sincedt/dx = (a - c x^4)/P^(3/2),(a - c x^4)/((a d + a e x^2 + c d x^4) sqrt(P))isdt/(d + (a e - b d) t^2). Whether an integrand qualifies is read iny = x^2moduloy^2 - w y + a/c; an elliptic integral does not pass.(c x)^t,tnot whole, writtenx^tbeside(c x)^t x^(-t), which is constant wherever it is defined.t = L1/L2.1/xalone, underx = 1/u.x^m (a + b x^n)^pfor wholemandnandp <= -2, by the reduction that raisespto-1.u = x^2.u = d + k x.Chain against master f7f2a7a, with the harness at its current revision on both arms:
Re-measured with the probe on f7f2a7a and on this branch:
(a - c*x^4)/((a*d + a*k*x^2 + c*d*x^4)*sqrt(a + b*x^2 + c*x^4))(g + h*x)*sqrt(a + b*x + c*x^2)/(a*d + b*d*x + c*d*x^2)^(3/2)x^2/((1 - a*x)^16*(1 + a*x)^11)1/((a + b/x)^(5/2)*(c + d/x)^2)(c*x)^(-1 - 2*n - n*p)*(a + b*x^n)^p1/(x^2*(a + b*x^3)^5)(d + k*x)^4/(a + b*(d + k*x)^2 + c*(d + k*x)^4)1/(x*(1 - x^2)^(1/3)*(3 + x^2))The probe calls master's answer to the third row wrong. It is exactly right: the difference is 0 at 50, 300 and 2,000 digits at two points. The probe's double evaluation of a 16,292-character expression is what fails. One answer gets longer, the last row.
One shape here,
(c x)^tasx^ttimes a constant, is close to #2012's. The rows that moved are the ones #2012 leaves.The BREAKING-CHANGES.md entry's "Now" column is re-measured on this build, and its comparison is against f7f2a7a.
🤖 Generated with Claude Code
https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura