Metodo de interpolacion de newton.
Enviado por Jorge Urueta • 30 de Noviembre de 2016 • Apuntes • 290 Palabras (2 Páginas) • 163 Visitas
clear,clc;
fprintf(' Programa echo por jorge urueta & jesus garcia \t\n');
pause(2);
clear,clc;
fprintf('-------------------------------------- \t\n');
fprintf(' Algoritmo de Interpolacion de newton \t\n');
fprintf('-------------------------------------- \t\n');
pause(2);
fprintf ('\n');
n=input ('ingrese el numero de parejas: ');
for p=1:n
fprintf ('\n');
fprintf ('ingrese el Valor X en la posicion %d: ',p);
x=input ('');
parx(p)=x;
fprintf ('ingrese el Valor Y en la posicion %d: ',p);
y=input ('');
pary(p)=y;
end
%acava de llenar los vectores con los valores
fprintf ('\n');
val=input('ingrese el valor a Interpolar: ');
clc;
%hallaremos las Bi
d=zeros(n);
d(:,1)=pary';
fprintf('las bi del polinomio son: \n');
pause (1);
fprintf ('\n');
for k=2:n
for j=1:n+1-k
d(j,k)=(d(j+1,k-1)-d(j,k-1))/(parx(j+k-1)-parx(j));
end
end
disp(d);
%mostraremos el resultado de newtom
fprintf ('\n');
disp('El resultado del polinomio de newton es: ')
pause (1);
aux=0;
for i=1:n
if (i==1)
aux=aux+d(1,i);
else
f=1;
for h=1:(i-1)
f=f*(val-parx(h));
end
aux=(aux + (d(1,i)*f));
end
end
fprintf ('\n');
disp(aux);
%polinomio de newton
fprintf ('\n');
disp('El Polinomio de newton generado es: ')
pause (1);
p=num2str(d(1,1));
xx=parx*-1;
for q=2:n
signo= '';
if d(i,q) >= 0
signo= '+';
end
xt = '';
for v=1:q-1
signo2= '';
if xx(i) >= 0;
signo2= '+';
end
xt=strcat(xt,'*(X',signo2,num2str(xx(v)),')');
end
p=strcat(p,signo,num2str(d(1,q)),xt);
end
fprintf ('\n');
fprintf ('Pn = %s\n',p);
fprintf ('\n');
...