d=input('ingresa los datos de la siguiente forma: [x1 f(x1) f`(x1);...;xn f(xn) f´(xn)], si no hay derivada ingrese 0 ');
n=length(d);
c=inline('chop(x,5)','x')
x=c(d(:,1));
a=zeros(n-1,n-1);
for j=1:n-1;
    if j==1;
        for k=1:n-j;
            if x(k)==x(k+1);
                a(k,j)=c(d(k,3));
            else
                a(k,j)=c((c(d(k+1,2))-c(d(k,2)))/(x(k+1)-x(k)));
            end
        end
    else
        for k=1:n-j;
            a(k+j-1,j)=c((c(a(k+j-1,j-1))-c(a(k+j-2,j-1)))/(x(k+j)-x(k)));
        end
    end
end
a
%Algoritmo para evaluar el polinomio
P=[d(1,2) diag(a)'];
s=length(P);
% v=x(1):0.01:x(n);
% lv=length(v);
% for k=1:lv;
t=input('introduzca el x en el cual desea evaluar= ');
e=P(s);
for j=1:s-1;
    y=(t-x(s-j))*e;
    w=P(1,s-j)+y;
    e=w;
end
%y(k)=w;
w
%end
%plot(d(:,1),d(:,2),'o',v,y,'r');