miércoles, 12 de octubre de 2011

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

No hay comentarios:

Publicar un comentario