Metodo de interpolacion de newton.

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;

...