An example from bug 1901044 (from Barton):
(%i1) integrate((5*x-2)^199,x);
(%o1) (5*x-2)^200/1000
(%i2) integrate((5*x-2)^3,x);
(%o2) (125*x^4)/4-50*x^3+30*x^2-8*x
Note the difference in the results.
In diffdiv in src/sin.lisp, there is a check to see if the exponent is
between 0 and 6. If so, the integrand is expanded and integrated
again, which gives the second result above. Otherwise, it uses the
general formula for (a*x+b)^n, which is the first answer above.
Any opinions on whether this is the right thing to do? And if so, why
6, and not some other value? Or should we just use the general
formula always?
Ray