Repository navigation
Fractional powers of a ± a sin(y) are integrated by the half angle at which they are squares - #1405
Merged
Merged
Conversation
… which they are squares (A + C sin(y)^2)/((c - c sin(y))^(3/2) sqrt(a + a sin(y))) and its kin ran past any budget: the substitution search has no radical of a trigonometric function it can rationalise, and the closed rule for (a + b sin)^n with a^2 = b^2 reads one such power alone. 1 + sin(y) is 2 sin(u)^2 and 1 - sin(y) is 2 cos(u)^2 for u = y/2 + pi/4, and 1 ± cos(y) the same with u = y/2, so each such power is a whole power of a sine or cosine of u times its sign -- a constant between the zeros, written in front -- and what is left is rational in sin(u) and cos(u), which the rules answer. Exact for any a, since 1 + sin(y) is not negative and (a q)^p = a^p q^p for the principal powers whatever a is; a^2 = b^2 is decided, not assumed; a lone positive power goes to the closed rule first, which answers it in one form across the zeros. Before the substitution search, which spent twenty seconds on the radicals. Rubi's family 4 sample: 261 to 277 of 422, sixteen rows and none lost, six of them timeouts before. The 1774-problem suite is unchanged at 1706 with 0 wrong. Part of #718. Co-Authored-By: Claude Opus 5 (1M context) <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura
Member
Author
|
The Ubuntu leg failed twice on |
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.
The largest remaining pocket of Rubi's family 4 was fractional powers of
a ± a sin(y)anda ± a cos(y)beside other things —(A + C sin^2)/((c - c sin)^(3/2) sqrt(a + a sin)),(A + B sin)/((a + a sin)(c - c sin)^(5/2)),sec^2/(a + a cos)^(5/2)— every one a timeout or a decline: the substitution search has no radical of a trigonometric function it can rationalise, andSolveAHalfPowerOfOnePlusASinereads one such power alone, to a positive exponent, witha^2 = b^2decided by evaluation (sosqrt(a + a sin(x))with a symbol was not even that).The rule.
1 + sin(y)is2 sin(u)^2and1 - sin(y)is2 cos(u)^2foru = y/2 + pi/4;1 ± cos(y)the same withu = y/2. So(a + a sin(y))^pisa^p 2^p |sin(u)|^(2p)=a^p 2^p sgn(sin(u)) sin(u)^(2p)for a half-integerp, the sign a constant between the zeros of the sine that comes out in front of the integral (top level only, as every rule that writes a sign for a function), and withsin(y) = 2 sin(u)^2 - 1,cos(y) = 2 sin(u) cos(u)anddx = 2 du/fwhat is left is rational insin(u)andcos(u), which the chain answers. Exact for anya:1 + sin(y)is not negative, so(a q)^p = a^p q^pfor the principal powers whateverais — the test pinsa = -3as well asa = 2.a^2 = b^2is decided by simplification, not assumed; a lone positive power still goes to the closed rule first, which answers it in one form across the zeros. Before the substitution search, which spent twenty seconds on the radicals before declining.Measured. Rubi's family 4 sample (422, six per file): 261 → 277, sixteen rows and none lost, six of them timeouts before —
(A + C sin^2)/((c - c sin)^(3/2) sqrt(a + a sin))at 498 ms,(a + a sin)^(7/2)and1/(a + a sin)^(5/2)(which the closed rule could not decide for a symbolica),sin^2 (a + a sin)^(3/2)(c - c sin)at 30 ms where it had run 24 s, and the cosine and secant kin. The 1,774-problem independent suite: 1706/1774, 0 wrong, 0 timeouts, verdicts identical to master, 89.6 s → 85.0 s. All four suites green (11,517);HalfAngleSquareTest(6 rows × two values ofa) with its.editorconfigsection; BREAKING-CHANGES entry.Not this rule's: a second radical beside them (
sqrt(g sin) sqrt(a + a sin)), a symbolic exponentm, and(a + a sec)^(m/2)shapes, which want the secant read first.Part of #718.
🤖 Generated with Claude Code
https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura