Интерполяционный многочлен, построенный на узлах Чебышёва

Код на Matlab.

Начало кода:

a = -2;

b = 4;

function y = func(x)

    y = -exp(x)+x.^2;

endfunction

for n = 2:2:16

    k1(n/2-1) = ...

    k2(n/2-1) = ...

    % страшный индекс обусловлен тем, что он равен 1,2,3,...

end