Interpolacion de newton

INSTITUTO POLITECNICO NACIONAL
UNIDAD INTERDICIPLINARIA DE INGENIERIA CAMPUS GUANAJUATO



TALLER DE METODOS NUMERICOS


METODO INTERPOLADOR DE NEWTON



OLMOSORTEGA ISAAC ANTONIO
TREJO RODRIGUEZ JUAN BERNARDO
HERNANDEZ SUARES ADRIAN
ORETGA ESCOBEDO BRIAN DANIEL


INGENIERIA EN SISTEMAS AUTOMOTRICES 2SV4INTRODUCCION:
El método de interpolación de Newton es un poco más complicado que el de LaGrange, pero esa complicación en los cálculos nos lleva a una solución más precisa.
Como todo método tieneun desarrollo teórico para llegar a la ecuación general:    

.
OBJETIVO:
Implementar en Scilab el algoritmo para hallar el polinomio interpolador de Newton, utilizarlo en un conjunto deDatos para hacer estimaciones y hacer comparaciones con otros modelos.
DESCRIPCION DEL METODO
Paso 1: Construye en Scilab una función que implemente el algoritmo del polinomio interpolador de NewtonPaso 2: Realizar un diagrama de dispersión con la tabla que se nos dio en la práctica.
Paso 3: Utilizando el método de Newton, halla una estimación para los siguientes cilindrares: 3,
3.5, 4, 4.5,5, 5.5, 6, 6.5, 7, 7.5 y 8 litros.
Paso 4: añadir la información del carro que se agrega en el paso 4 de la practica y volver a introducir los datos en el programa.
Paso 5: comparar los datosobtenidos del paso 4 con los anteriores y hacer una grafica con los datos obtenidos.

Líneas de código
//Funcion

function N=pnewton(X, Y)
n=length(X);
dfd=zeros(n,n);
dfd(:,1)=Y(:);for j=2:n
for i=j:n
dfd(i,j)=(dfd(i,j-1)-dfd(i-1,j-1))/(X(i)-X(i-j+1))
end
end
N=dfd(1,1);
for i=2:n
N=N+dfd(i,i)*poly(X(1:(i-1)),"x");
end
function y=f(x)y=horner(N,x);
endfunction
x=min(X):(max(X)-min(X))/100:max(X);
plot(X,Y,'*');
fplot2d(x,f);
endfunction

//Datos

x=[1 3 5 7 13];
y=[800 2310 3090 3940 4755];

//Funcionamiento...