Skip to content

Commit a1f6426

Browse files
Split a rational function of a sine over powers of linears in it by the residues at their roots (#2001)
cot(c + d x)^6/(a + b sin(c + d x))^3 is (1 - s^2)^3/(s^6 (a + b s)^3) in s = sin(x), over powers of linears in s: the sine, 1 +- s from an even power of the cosine, and a + b s. It is split by the residues at their roots, each block a Taylor series at a root, into constants over powers of one linear that the rules for csc^n, 1/(1 +- sin)^n and 1/(a + b sin)^n answer. Under the half-angle tangent, where master took it, every linear but the sine is a quadratic, and the cube of a symbolic one ran out of time. Only where the half-angle is that slow. And a repeated quadratic with Gaussian rational roots is split by the residues at one root, those at the other their conjugates: sin(c + d x)^6/(a + b tan(c + d x))^2. Rubi pocket (1,817 rows of family 4): 1769 -> 1783 solved; the moves re-run alone gain 15 rows of 4.1.1.2, 4.1.1.3, 4.1.2.1, 4.1.2.2 and 4.3.1.3. Part of #233. Claude-Session: https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura Co-authored-by: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
1 parent efd72ac commit a1f6426

5 files changed

Lines changed: 610 additions & 17 deletions

File tree

‎BREAKING-CHANGES.md‎

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -763,6 +763,33 @@ unevaluated are now answered ([#718](https://github.com/asc-community/AngouriMat
763763
| `"e^(4*atanh(a*x))*(k-a*k*x)^p".ToEntity().Integrate("x")` | `integral(...)` | 434 characters |
764764
| `"e^atanh(x)*x*(1+x)^2".ToEntity().Integrate("x")` | `integral(...)` | 1,110 characters |
765765

766+
### A rational function of a sine over powers of linears in it is split by the residues at their roots
767+
768+
**Answers where there were none.** `cot(c + d x)^6/(a + b sin(c + d x))^3` was declined, and so were eight more of
769+
Rubi's 4.1.1.3, 4.1.2.1, 4.1.2.2 and 4.3.1.3 ([#233](https://github.com/asc-community/AngouriMath/issues/233)):
770+
- In `s = sin(x)` that integrand is `(1 - s^2)^3/(s^6 (a + b s)^3)`, over powers of linears in `s`: the sine, `1 ± s`
771+
from an even power of the cosine, and `a + b s`. It is split by the residues at their roots, each block a Taylor
772+
series at a root, and every piece is a constant over a power of one linear, which the rules for `csc(x)^n`,
773+
`1/(1 ± sin(x))^n` and `1/(a + b sin(x))^n` answer. Under the half-angle tangent, where the unreleased master took
774+
them, every linear but the sine is a quadratic, and the cube of a symbolic one ran out of time. Only where the
775+
half-angle is that slow: a linear with a symbolic root to the third power or more, and the half-angle's denominator
776+
past the eleventh degree. Not for the cosine, whose half-angle answer is the shorter.
777+
- A repeated quadratic with Gaussian rational roots `r = ρ + σ i` and `r̄` is split by the residues at `r`, one
778+
division of series per order, with those at `r̄` their conjugates. `sin(c + d x)^6/(a + b tan(c + d x))^2` is
779+
`t^6/((1 + t^2)^4 (a + b t)^2)` in `t = tan(c + d x)`, and was declined.
780+
781+
`csc(c + d x)^2 sec(c + d x)^2/(a + b sin(c + d x))^3` was declined by 2.5.0 too, and the unreleased master answered
782+
it in eleven seconds and 104,032 characters, which are a third of a second and 7,926 now.
783+
784+
| Input | Was (2.5.0) | Now |
785+
|---|---|---|
786+
| `"tan(c + d*x)^4/(a + b*sin(c + d*x))^3".ToEntity().Integrate("x")` | `integral(...)` | 9,096 characters |
787+
| `"cot(c + d*x)^6/(a + b*sin(c + d*x))^3".ToEntity().Integrate("x")` | `integral(...)` | 9,316 characters |
788+
| `"1/((a + b*sin(k + f*x))^3*(c + d*sin(k + f*x))^3)".ToEntity().Integrate("x")` | `integral(...)` | 15,245 characters |
789+
| `"sin(c + d*x)^6/(a + b*tan(c + d*x))^2".ToEntity().Integrate("x")` | `integral(...)` | 1,818 characters |
790+
| `"sin(c + d*x)^4/(a + b*tan(c + d*x))^4".ToEntity().Integrate("x")` | `integral(...)` | 2,586 characters |
791+
| `"csc(c + d*x)^2*sec(c + d*x)^2/(a + b*sin(c + d*x))^3".ToEntity().Integrate("x")` | `integral(...)` | 7,926 characters, where the unreleased master gave 104,032 |
792+
766793
### A sine's binomial cube below the bar is split over the cube's two factors
767794

768795
**Answers where there were none.** `csc(c + d x)/(a + b sin(c + d x)^3)` was declined, and so were the second to fifth

0 commit comments

Comments
 (0)