h=0.2
x=[0:h:1]
f=-3*x.^3+x.^2+2*x+2
y=zeros(6,6)
for k=1:6
    y(k,1)=f(k)
end

m=6
for i=2:6
    m=m-1
    for j=1:m
        y(j,i)=y(j+1,i-1)-y(j,i-1)
    end
end

d1=0
for n=1:5
    d1=d1+((-1).^(n-1))*y(1,n+1)*(1/n)
end
disp('Elso derivalt a 0 helyen:')
d1=d1/0.2