[SCILAB] Newton

Mam do rozwiązania taki problem:

Napisz program obliczający metodą Newtona wartość funkcji zadanej w postaci tabeli:

x 0.0 0.2 0.4 0.6 0.8 1.0
f(x) 0.000000 0.587785 0.951057 0.951057 0.587785 0.000000

dla wartości x należacej do przedziału [0.0 , 1.0]. wartość X wprowadzana z klawiatury. wynik porównać z wartością f(x)=sin(%pi*x) oraz wynikiem użycia funkcji "interp".

Proszę o jakikolwiek kod. Nie mam pojęcia jak się za to zabrać, może ma ktoś jakiś schemat? Proszę o pomoc.

r=0.0:0.2:1.0
q=[0.000000; 0.587785; 0.951057; 0.951057; 0.587785; 0.000000]
x=input('podaj wrtość x= ');
while (x<0|x>1)
    disp('podaj nową wartość x ');
    x=input('x= ');
k=1:6
for j=2:6
    for i=
        a(k,1)=q(k);
        end
end