clc;
clear;
n=5;
a=-3.2;
b=2;

X=[-3.12; -2.09 ; -1.07 ; 0.61 ; 1.37 ; 1.91 ]
x_k=0.11;
m=13;
function y=f(X)
    y=sin(X);
endfunction

Y=f(X);
function w=omega(X,x_k,n)
    w=1;
    for i=0:n
        w=w*(x_k - X(i+1));
    end
endfunction

h=(b-a)/m;

function y=fp(X)
    y=-sin(X);
endfunction
for i=0:m
    W(i+1)=a+i*h
end

function w=msup(W,m)
    for i=0:m
        A(i+1)=fp(W(i+1));
    end
    w=max(A);
endfunction

wynik=msup(W,m)/factorial(n+1)*abs(omega(X,x_k,n));
disp(wynik);