Tarea 2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

J.Reynaldo Garza M. CICFCM Aula 2 (Dr.

Omar Ibarra 4:00 pm): 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

f3 (x) = b0 + b1 (xx0 ) + b2 (xx0 )(xx1 ) + b3 (xx0 )(xx1 )(xx2 ) (1)

Las primeras diferencias divididas del problema son


1,3862940
f [x1 , x0 ] = 40 = 0,4620981

1,7917591,386294
f [x2 , x1 ] = 64 = 0,2027326

1,609481,791759
f [x3 , x2 ] = 56 = 0,1823216

Las segundas diferencias divididas son

0,20273260,4620981
f [x2 , x1 , x0 ] = 61 = 0,05187311

0,18232160,2027326
f [x3 , x2 , x0 ] = 54 = 0,02041100

La tercera diferencia dividida es

0,02041100(0,05187311)
f [x3 , x2 , x1 , x0 ] = 51 = 0,007865529

Los resultados de f [x1 , x0 ], f [x2 , x1 , x0 ] y f [x3 , x2 , x1 , x0 ] representan los coeficientes b1, b2 y b3 de la


ecuacion (1), respectivamente. Junto con b0 = f (x0 ) = 0,0, la ecuacion (1) es

f3 (x) = 0 + 0,4620981(x1)0,05187311(x1)(x4) + 0,007865529(x1)(x4)(x6) (2)

Problema 1 continua en la siguiente hoja. . . Pagina 1 de 2


J.Reynaldo Garza M. CICFCM Aula 2 (Dr. Omar Ibarra 4:00 pm): Tarea #2 Problema 1 (continua)

A continuacion se muestra el codigo de programcion realizado en R studio.

rm(list = ls())
graphics.off()

x<-c(1,4,6,5) #Intervalo en el cual la funcion sera evaluada


5 y<-log(x) #Funcion que sera aproximada

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)

Estos son los resultados obtenidos:


Coeficientes

Tabla

Pagina 2 de 2

También podría gustarte