polinomio

function x=polinomio(a,b,xb)

n=length(a);

m=zeros(n);

p=n;

for i=1:n

    for j=1:n

        m(i,j)= (a(i)^(p-1));

        p=p-1;

    end

    p=n;

end

x=inv(m)*b';

p=n;

r=0;

for i=1:n

    r=r+x(i)*(xb^(p-1));

    p=p-1;

end

disp('los coef del polinomio')

disp(x)

disp('el resultado es de p(x)=')

disp(r);

por ejemplo si a=[0 1 2 3] y b=[0 1 4 9] y xb=2.6
el resul est

los coef del polinomio
   -0.0000
    1.0000
   -0.0000
         0
el resultado es de p(xb)=
    6.7600