Tarea 2
Tarea 2
Tarea 2
Problema 1
A continuacion se muestra el problema a tratar y posteriormente mostraremos el codigo de programacion
realizado en R studio.
Datos x0 = 1, x1 = 4, x2 = 6 y x3 = 5 se utilizaron para estimar ln(x) mediante el polinomio de interpolacion
de Newton de tercer grado.
Solucion
El polinomio de tercer grado es
1,7917591,386294
f [x2 , x1 ] = 64 = 0,2027326
1,609481,791759
f [x3 , x2 ] = 56 = 0,1823216
0,20273260,4620981
f [x2 , x1 , x0 ] = 61 = 0,05187311
0,18232160,2027326
f [x3 , x2 , x0 ] = 54 = 0,02041100
0,02041100(0,05187311)
f [x3 , x2 , x1 , x0 ] = 51 = 0,007865529
rm(list = ls())
graphics.off()
n<-length(x)
d<-matrix(0,n,n)
10 for(i in 1:n){
d[i,1]<-y[i]
}
for(j in 2:n){
for(i in j:n){
15 d[i,j]<-(d[i-1,j-1]-d[i,j-1])/(x[i-j+1]-x[i])
}
}
a<-diag(d)
20 print(a)
Tabla
Pagina 2 de 2