|
|
A280572
|
|
E.g.f. satisfies: A(x - Integral A(x) dx) = x + Integral 4*A(x) dx.
|
|
10
|
|
|
1, 5, 40, 525, 10025, 253475, 8015725, 305359050, 13645726250, 701304298375, 40822454374125, 2658840618527250, 191861336190647375, 15213199343853357500, 1316408013706224687500, 123576861126283832953125, 12521371849855149886590625, 1363361618975383978443843750, 158900334287408210286438971875, 19755940413686794723417400000000, 2612146114817877629253999384562500, 366294181903982533559997504649828125
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
E.g.f. A(x) satisfies:
(1) A(x - Integral A(x) dx) = x + Integral 4*A(x) dx.
(2) A(x) = x + 5 * G( (A(x) + 4*x)/5 ), where G(x) = Integral A(x) dx.
(3) A(x) = -4*x + 5 * Series_Reversion(x - Integral A(x) dx).
(4) A( (A(x) + 4*x)/5 ) = (A'(x) - 1)/(A'(x) + 4).
(5) A'(x - Integral A(x) dx) = (1 + 4*A(x))/(1 - A(x)).
a(n) = Sum_{k=0..n-1} A277410(n,k) * 5^(n-k-1).
|
|
EXAMPLE
|
E.g.f.: A(x) = x + 5*x^2/2! + 40*x^3/3! + 525*x^4/4! + 10025*x^5/5! + 253475*x^6/6! + 8015725*x^7/7! + 305359050*x^8/8! + 13645726250*x^9/9! + 701304298375*x^10/10! + 40822454374125*x^11/11! + 2658840618527250*x^12/12! +...
Let G(x) = Integral A(x) dx, then A(x - G(x)) = x + 4*G(x) where
G(x) = x^2/2! + 5*x^3/3! + 40*x^4/4! + 525*x^5/5! + 10025*x^6/6! + 253475*x^7/7! + 8015725*x^8/8! + 305359050*x^9/9! + 13645726250*x^10/10! + 701304298375*x^11/11! + 40822454374125*x^12/12! +...
Also, A(x) = x + 5 * G( (A(x) + 4*x)/5 ).
RELATED SERIES.
We have (A(x) + 4*x)/5 = Series_Reversion( x - Integral A(x) dx ), where
(A(x) + 4*x)/5 = x + x^2/2! + 8*x^3/3! + 105*x^4/4! + 2005*x^5/5! + 50695*x^6/6! + 1603145*x^7/7! + 61071810*x^8/8! + 2729145250*x^9/9! + 140260859675*x^10/10! + 8164490874825*x^11/11! + 531768123705450*x^12/12! +...
Further, A( (A(x) + 4*x)/5 ) = (A'(x) - 1)/(A'(x) + 4), which begins
A( (A(x) + 4*x)/5 ) = x + 6*x^2/2! + 63*x^3/3! + 1045*x^4/4! + 24105*x^5/5! + 716195*x^6/6! + 26137820*x^7/7! + 1134457060*x^8/8! + 57203895725*x^9/9! + 3292221321425*x^10/10! + 213282348138700*x^11/11! + 15380885339509825*x^12/12! +...
|
|
MATHEMATICA
|
m = 23; A[_] = 0;
Do[A[x_] = -4 x + 5 InverseSeries[x - Integrate[A[x], x] + O[x]^m], {m}];
|
|
PROG
|
(PARI) /* A(x) = x + (p+q)*G((p*A(x) + q*x)/(p+q)) ; G(x) = Integral A(x) dx: */
{a(n, p=1, q=4) = my(A=x, G); for(i=1, n, G = intformal(A +x*O(x^n)); A = x + (p+q)*subst(G, x, (p*A + q*x)/(p+q)) +x*O(x^n)); n!*polcoeff(A, n)}
for(n=1, 30, print1(a(n, 1, 4), ", "))
(PARI) /* A(x - Integral p*A(x) dx) = x + Integral q*A(x) dx: */
{a(n, p=1, q=4) = my(A=[1], F=x); for(i=1, n, A=concat(A, 0); F=x*Ser(A); G=intformal(F); A[#A] = -polcoeff(subst(F, x, x - p*G) - q*G, #A) ); n!*A[n]}
for(n=1, 30, print1(a(n, 1, 4), ", "))
(PARI) /* Informal code to generate the first N terms: */
{N=20; p=1; q=4; A=x; for(i=1, N, G=intformal(A +x*O(x^N)); A = x + (p+q)*subst(G, x, (p*A + q*x)/(p+q))); Vec(serlaplace(A))}
|
|
CROSSREFS
|
Cf. A277410, A210949, A277403, A279843, A279844, A279845, A280570, A280571, A280573, A280574, A280575.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|