Пока я разбираюсь с индексами, возможно вам будет проще понять по листингу проверки равенства в "Maple 12". Тут
![$x=5, y=4$ $x=5, y=4$](https://dxdy-01.korotkov.co.uk/f/8/6/6/866bd865e34e5a182ceb2df9032866bb82.png)
, таким образом, результат равен
![$625$ $625$](https://dxdy-04.korotkov.co.uk/f/b/c/b/bcb1a31cbdd1601588434a73561d889f82.png)
.
Код:
restart;
with(combinat);
n := 5;
s1 := array(1 .. n, 1 .. n);
for i to n do
for j to n do
s1[i, j] := stirling1(i, j)/factorial(i-1)
end do
end do;
s1 := evalm(s1);
sm := convert(s1, Matrix);
C := array(1 .. n, 1 .. n);
for i to n do
for j to n do
C[i, j] := (-1)^(i+j)*binomial(i-1, j-1)
end do
end do;
C := evalm(C);
P := convert(C, Matrix);
with(LinearAlgebra);
S[1] := Transpose(sm);
f := proc (j) options operator, arrow; j^y end proc;
y := 4;
k := Vector(5, f);
Y := MatrixVectorMultiply(P, k);
F := proc (j) options operator, arrow; x^(j-1) end proc;
x := 5;
k := Vector[row](5, F);
X := VectorMatrixMultiply(k, S[1]);
Multiply(X, Y);
Можно просто вставить код в Maple и поиграть со значениями
![$x$ $x$](https://dxdy-04.korotkov.co.uk/f/3/3/2/332cc365a4987aacce0ead01b8bdcc0b82.png)
и
![$y$ $y$](https://dxdy-02.korotkov.co.uk/f/d/e/c/deceeaf6940a8c7a5a02373728002b0f82.png)
.