Как расширить в ряд Тейлора состав функций?

Question

Я хотел бы расширить в ряд Тейлора функцию типа: f(x+f(x)) вокруг x=a в случае, когда f(a)=0.

(%i1) atvalue(f(x),[x=a],0)$

Прямое исчисление дает:

(%i2) taylor(f(x+f(x)),x,a,2);

(%o2)/T/ f(a)+(at('diff(f(f(x)+x),x,1),x=a))*(x-a)+((at('diff(f(f(x)+x),x,2),x=a))*(x-a)^2)/2+...

Если я определю промежуточную функцию:

(%i3)define(tf(x),taylor(f(x),x,a,2))$

Затем развернутое в ряд Тейлора я получу:

(%i4) taylor(f(x+tf(x)),x,a,2);

(%o4) 0+...

Я ожидаю следующий результат: f(1+f'(a))f'(a)(x-a)+(x-a)^2 f''(a)[f'(a)+(1+f'(a))^2/2]+o(x-a)^2

Как я мог решить эту проблему?

Answer 1

Вы можете использовать gradef для упрощения записи.

gradef(f(x),  f1(x)) $
gradef(f1(x), f2(x)) $

atvalue(f(x), x = a,  0) $

e: f(x+f(x)) $
e: taylor(e, x, a, 2) $
e: expand(e, 0, 0)$ /* 'taylor' form to ordinar expression */
e: ev(e, nouns);    /* f(a) to 0 */

возвращает

          2                                          2
       (f1 (a) f2(a) + 3 f1(a) f2(a) + f2(a)) (x - a)
(%o7) -----------------------------------------------
                              2
                                                          2
                                                     + (f1 (a) + f1(a)) (x - a)