Интерполяция методом Ньютона (Scilab)

[CraZZZy-GameRRR]

Здравствуйте, товарищи, нужна ваша помощь.
Задание: произвести интерполяцию по этой формуле:


Сделал так:

clc; clear;

X = [-8 -2 3 11];
Y = [2 7 9 4];
plot2d(X, Y, -1, leg="Точки");

// расчёт по формуле Ньютона
Y5 = [];
n = length(X);

C = [];
C(1) = Y(1);
for i = 1:n
    C(i+1) = Y(i+1) - C(1);
    for j = 1:i-1
        s = C(j+1);
        for k = 0:i-2
            s = s * (X(i+1) - X(k+1));
        end
        C(i+1) = C(i+1) - s;
    end
    
    for j = 0:i-1
        C(i+1) = C(i+1) / (X(i+1) - X(j+1));
    end
end

for x = -15:15
    p = C(1);
    for i = 1:n
        s = C(i+1);
        for j = 0:n-1
            s = s * (x - X(j+1));
        end
        p = p + s;
    end
    
    Y5 = [Y5, p];
end

plot2d(-15:15, Y5, style=[color("magenta")], leg="Метод Ньютона");

Компилятор находит ошибку:

     C(i+1) = Y(i+1) - C(1);
                    !--error 21 
Неправильный индекс.
at line      14 of exec file called by :    
exec('C:\Users\User\Documents\Scilab\lab2\newton.sce', -1)