Во-первых, я вижу ошибку, когда вы пишете a(n+1)
, для n
переходя от 0
до 6
. Когда n
равно 6
, у вас будет a(7)
, которого не существует. a
имеет только 6 элементов.
Тогда у вас есть другая проблема, потому что вы индексируете a
с помощью символа, и символическая индексация не разрешена (см. этот пост ).
Действительно ли n
должен быть символом, и вам действительно нужно использовать symsum
? Если нет, вы можете попробовать:
a = [1,2,3,40,51,61];
n = 0:5;
syms u
S1(u) = sum( sym( (a(n+1).*(-u).^n)./factorial(n) ) )
Возвращает:
S1(u) =
- (61*u^5)/120 + (17*u^4)/8 - (20*u^3)/3 + (3*u^2)/2 - 2*u + 1