Получить n-член пос-ти по производящей ф-ии [Matlab/Wolfram]

[GAA]

maximkarimov в сообщении #1474277 писал(а):

Ммм.. а как нужно изменить вторую строчку, чтобы не весь вектор, а только один коэффициент получить?

R2013b:

Используется синтаксис Matlab M

>> syms x
>> C = coeffs(taylor(x/(1 - x - x^2), x, 0, 'Order', 46)); C(45)
ans = 1134903170

maximkarimov в сообщении #1474277 писал(а):

P.S. Пробовал применить синтаксис MuPAD, который Вы показали выше, но внутри Matlaba он "не работает".

Как было приведено в предыдущей теме в R2013b в командном окне (или в m-файле) можно вызвать функцию MuPAD при помощи feval(symengine, …)

Используется синтаксис Matlab M

>> P = feval(symengine,'taylor', x/(1 - x - x^2), x, 0, 45);
P = Series::Puiseux::create(1, 1, 46, [1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, 610, 987, 1597, 2584, 4181, 6765, 10946, 17711, 28657, 46368, 75025, 121393, 196418, 317811, 514229, 832040, 1346269, 2178309, 3524578, 5702887, 9227465, 14930352, 24157817, 39088169, 63245986, 102334155, 165580141, 267914296, 433494437, 701408733, 1134903170], x, 0, Undirected)

В данном случае в этом мне не видно смысла.

[maximkarimov]

Если:

Используется синтаксис Matlab M

syms x
P = feval(symengine,'taylor', x/(1 - x - x^2), x, 0, 45);
 

то получаем не одно значение, а весь вектор.

Ну да ладно, поскольку производительность в данном конкретном случае для меня не важна, вышел из положения так:

Используется синтаксис Matlab M

syms x
a=coeffs(taylor(x/(1 - x - x^2), x, 0, 'Order', 45));
a=a(size(a,2));
 

Спасибо!