Practica 06

Apellidos y Nombres: Mallqui Rivera Paul Lenyn Ciclo: IV Fecha:

07/04/2021

Ejercicio 01
Programo un programa para obtener las diferencias divididas:

In[79]:= DD[a_,b_]:=Module{x,y,n},
x= a;
y=b;
n=Length[a];
d= ConstantArray[0,{n,n+1}];
d[[All,1]]=x;
d[[All,2]]=y;
Forj=3,j<=n+1,j++,

Fori=j-1,i<=n,i++,
(d[[i-1,j-1]]-d[[i,j-1]])
d[[i,j]]= ;
(d[[i-j+2,1]]-d[[i,1]])


;
dd = Take[d,All,{2,n+1}];

Diagonal[dd]


Ejercicio 02
La tabla siguiente se muestra la estimación de la concentración atmosférica media de Dióxido de
carbono en atmósfera terrestre, partes por millón por volumen. Encuentre el polinomio de interpo-
lador de Newton que interpole a los datos y úselo para estimar la concentración de CO2 en a) En el
año 1950; b) Para el año futuro del 2050.

In[27]:= data = {{1800,280},{1850,283},{1900,291},{2000,370}};

Como tenemos 4 puntos nuestro polinomio de Newton sera de grado 3 de la forma:

2 Practica_06.nb

In[11]:= P[x_]:= a1 + a2*(x-1800) + a3*(x-1800)(x-1850) + a4*(x-1800)(x-1850)(x-1900);

Lo resolvemos mediante métodos rápidos Solve

In[17]:= coeficientes=Solve[{P[1800]⩵280,P[1850]⩵283,P[1900]⩵291,P[2000]⩵370},{a1,a2,a3,a4}]

3 1 1
Out[17]= a1 → 280, a2 → , a3 → , a4 → 
50 1000 62 500

In[20]:= ( a1 a2 a3 a4 ) = {a1, a2, a3, a4} /. coeficientes

3 1 1
Out[20]= 280, , , 
50 1000 62 500

Tenemos ya el polinomio particular

In[22]:= P[x]//Simplify

3213 x 439 x2 x3
Out[22]= - 97 730 + - +
20 5000 62 500

Graficamos el polinomio y lo interpolamos con los puntos

In[34]:= d=ListPlot[data,PlotStyle→PointSize[0.02]];
p =Plot[P[x],{x,1700,2000},PlotRange→{0,400},AxesLabel→ {"x",y},PlotStyle→Red];
Show[p,d]

y
400

300

Out[36]= 200

100

x
1700 1750 1800 1850 1900 1950 2000

Con instrucciones rápidas

In[48]:= Fitdata,1,x,x2 ,x3 ,x,WorkingPrecision→∞

3213 x 439 x2 x3
Out[48]= - 97 730 + - +
20 5000 62 500

Con diferencias divididas en Mathematica

 Practica_06.nb 3

1800 280
1850 283
In[49]:= DD , 
1900 291
2000 370

3 1 1
Out[49]= {280},  ,  ,  
50 1000 62 500

Con diferencias divididas en Matlab

Como podemos comprobar con ambos software nos salen los mismos coeficientes de Newton por
lo que encontramos nuestro polinomio interpolador.
En el inciso a nos solicita estimar la concentración de CO2 en el año 1950 por lo que reemplazamos
los datos:

In[50]:= P[1950]

Out[50]= 316

a) En el año 1950 la concentración de Co2 fue de 316 ppm.

En el inciso a nos solicita estimar la concentración de CO2 en el año 2050 por lo que reemplazamos
los datos :

In[51]:= P[2050]

Out[51]= 465

b) En el año 2050 se estima una concentración de Co2 fue de 465 ppm.

In[53]:= ClearAll["Global`*"]
borra todo

Ejercicio 03
En siguiente tabla horas: 13 14 15 16 17 18; grados: 18 18 17 16 15 14 se muestra las temperaturas
que fueron medidas cada hora durante un lapso total de 5 horas en la ciudad de Huacho el día de
ayer.
a) Hallar el polinomio interpolador de Newton

In[54]:= data = {{13,18},{14,18},{15,17},{16,16},{17,15},{18,14}};

Como tenemos 6 puntos nuestro polinomio de Newton sera de grado 5 de la forma:

4 Practica_06.nb

In[55]:= P[x_]:= a1 + a2*(x-13) + a3*(x-13)(x-14) + a4*(x-13)(x-14)(x-15) + a5*(x-13)(x-14)(x-15)(x-

Lo resolvemos por métodos rápidos Solve

In[58]:= coeficientes=Solve[{P[13]⩵18,P[14]⩵18,P[15]⩵17,P[16]⩵16,P[17]⩵15,P[18]⩵14},{a1,a2,a3,a4,a5

In[59]:= ( a1 a2 a3 a4 a5 a6 ) = {a1, a2, a3, a4, a5, a6} /. coeficientes

1 1 1 1
Out[59]= 18, 0, - , , - , 
2 6 24 120

Tenemos ya el polinomio particular

In[62]:= P[x]//Simplify

1
Out[62]= - 1 024 320 + 323 724 x - 40 720 x2 + 2555 x3 - 80 x4 + x5 
120

Graficamos el polinomio y lo interpolamos con los puntos

In[66]:= d=ListPlot[data,PlotStyle→PointSize[0.02]];
p =Plot[P[x],{x,10,20},PlotRange→{0,20},AxesLabel→ {"Horas","Grados"},PlotStyle→Red];
Show[p,d]

Grados
20

15

Out[68]= 10

Horas
10 12 14 16 18 20

b) Lo mismo pero con diferencias divididas luego utilizar el programa

Creamos la matriz que contendrá todas las operaciones como contamos con 6 puntos la matriz
sera de 6x7
13 18     
14 18 0    
-1
15 17 -1 2
  
A= 16 16 -1 0 1
 
6
1
17 15 -1 0 0 - 24

1
18 14 -1 0 0 0 120

Realizamos las operaciones para ir llenando la matriz

18 - 18
A23 = =0
13 - 14
 Practica_06.nb 5

18 - 17 0+1 -1
A33 = = -1 A34 = =
14 - 15 13 - 15 2

17 - 16 -1 + 1 -1 / 2 - 0 1
A43 = = -1 A44 = = 0 A45 = =
15 - 16 14 - 16 13 - 16 6

A53 =
16 - 15 -1 + 1 0-0 1/6-0 1
= -1 A54 = = 0 A55 = = 0 A56 = = -
16 - 17 15 - 17 14 - 17 13 - 17 24

15 - 14 -1 + 1
A63 = = -1 A64 = =
17 - 18 16 - 18
0-0 0-0 -1 / 24 - 0 1
0 A65 = = 0 A66 = = 0 A67 = =
15 - 18 14 - 18 13 - 18 120

Por lo tanto el Polinomio interpolador de Newton de grado 5 es el siguiente :

-1 1 1
P (x) = 18 + (x - 13) (x - 14) + (x - 13) (x - 14) (x - 15) + - (x - 13)
2 6 24
1
(x - 14) (x - 15) (x - 16) + (x - 13) (x - 14) (x - 15) (x - 16) (x - 17);
120
Utilizando el programa de diferencias divididas en Mathematica
13 18
14 18
15 17
In[80]:= DD , 
16 16
17 15
18 14

1 1 1 1
Out[80]= {18}, {0}, - ,  , - ,  
2 6 24 120

Utilizando el programa de diferencias divididas en Matlab

c) Resolver con instrucciones rápidas

In[82]:= Fitdata,1,x,x2 ,x3 ,x4 ,x5 ,x,WorkingPrecision→∞//Simplify

1
Out[82]= - 1 024 320 + 323 724 x - 40 720 x2 + 2555 x3 - 80 x4 + x5 
120
6 Practica_06.nb

In[83]:= ClearAll["Global`*"]
borra todo

Ejercicio 04
Si el número de artículos producidos por semana es x = [0, 2, 4, 6] medidos en miles, la función
costo de un fabricante es C = [2, 46/15, 38/15, 4/5] medidos en miles de soles.
a) Hallar función Costo por interpolación de Newton y diferencias divididas.

da