Universidad Tecnológica Nacional

Facultad Regional Resistencia

Matemática Superior Aplicada

Trabajo Práctico N°1:
Números Complejos, Interpolación, Derivación e
Integración Numérica

INTEGRANTES:
Aguirre, Claudia Verónica
Avalos Llano, Ivana Natali
Gomez, Daiana Soledad
Lopez Croce, Tomas
Pibernus, Ariana Agustina
PROFESORES:
Milena Balbi
Nápoles Valdés

FECHA: 2/11/22

AÑO 2022
1-NUMEROS COMPLEJOS
A)- HALLAR LA RELACIÒN QUE DEBEN VERIFICAR LOS COEFICIENTES A Y B REALES PARA QUE
LAS RAICES DE LA ECUACIÓN 𝑨 ∗ 𝒁𝟐 + 𝑩 ∗ 𝒁 = 𝟏 − 𝒊
𝑍1 = 𝑥1 + 𝑦1 𝑖 ; 𝑍2 = 𝑥2 + 𝑦2 𝑖

𝐴 ∗ (𝑥1 + 𝑦1 𝑖)2 + 𝐵 ∗ (𝑥2 + 𝑦2 𝑖) = 1 − 𝑖

𝐴 ∗ (𝑥1 2 + 2𝑥1 𝑦1 𝑖 + 𝑦1 2 𝑖 2 ) + 𝐵𝑥2 + 𝐵𝑦2 𝑖 = 1 − 𝑖

𝐴𝑥1 2 + 𝐴2𝑥1 𝑦1 𝑖 − 𝐴𝑦1 2 + 𝐵𝑥2 + 𝐵𝑦2 𝑖 = 1 − 𝑖

(1) 𝐴𝑥1 2 − 𝐴𝑦1 2 + 𝐵𝑥2 = 1

(2) 𝐴2𝑥1 𝑦1 𝑖 + 𝐵𝑦2 𝑖 = −𝑖
De (1) despejo B:

𝐴𝑥1 2 − 𝐴𝑦1 2 + 𝐵𝑥2 = 1

𝐴(𝑥1 2 − 𝑦1 2 ) + 𝐵𝑥2 = 1

1 − 𝐴(𝑥1 2 − 𝑦1 2 )
𝐵=
𝑥2
Luego reemplazo B en (2):

1 − 𝐴(𝑥1 2 − 𝑦1 2 )
𝐴2𝑥1 𝑦1 𝑖 + ( )𝑦2 𝑖 = −𝑖
𝑥2

𝑦2 𝑖 𝐴(𝑥1 2 − 𝑦1 2 ) ∗ 𝑦2 𝑖
𝐴2𝑥1 𝑦1 𝑖 + − = −𝑖
𝑥2 𝑥2

𝑦2 (𝑥1 2 − 𝑦1 2 )𝑖 𝑦2 𝑖
𝐴 (2𝑥1 𝑦1 𝑖 − ) = (−𝑖 − )
𝑥2 𝑥2
𝒚 𝒊
(−𝒊 − 𝒙𝟐 )
𝟐
𝑨= 𝟐 − 𝒚 𝟐 )𝒊
𝒚 (𝒙
(𝟐𝒙𝟏 𝒚𝟏 𝒊 − 𝟐 𝟏 𝒙 𝟏 )
𝟐

Finalmente remplazando A en B:
𝒚𝟐 𝒊
(−𝒊 − )
𝒙𝟐
𝟏−[ ] (𝒙𝟏 𝟐 − 𝒚𝟏 𝟐 )
𝒚𝟐 (𝒙𝟏 𝟐 − 𝒚𝟏 𝟐 )𝒊
(𝟐𝒙𝟏 𝒚𝟏 𝒊 − )
𝒙𝟐
𝑩=
𝒙𝟐
B)-HALLA EL LUGAR GEOMÉTRICO QUE SATISFACE:
i)- |𝑍 + 1| > |𝑍 − 𝑖| + 1
|𝑥 + 𝑦𝑖 + 1| > |𝑥 + 𝑦𝑖 − 1| + 1

√(𝑥 + 1)2 + (𝑦)2 > √(𝑥 − 1)2 + (𝑦)2 + 1

(𝑥 + 1)2 + (𝑦)2 > (√(𝑥 − 1)2 + (𝑦)2 + 1 )2

𝑥 2 + 2𝑥 + 1 + 𝑦 2 > 1 + (𝑥 − 1)2 + (𝑦)2 + 2√(𝑥 − 1)2 + (𝑦)2

𝟒𝒙 − 𝟏 > 2√(𝑥 − 1)2 + (𝑦)2

Elevo al cuadrado
(4𝑥 − 1)2 > 4((𝑥 − 1)2 + (𝑦)2 )
16𝑥 2 − 8𝑥 + 1 > 4(𝑥 2 − 2𝑥 + 1 + 𝑦 2 )
16𝑥 2 − 8𝑥 + 1 > 4𝑥 2 − 8𝑥 + 4 + 4𝑦 2
16𝑥 2 − 4𝑥 2 − 8𝑥 + 8𝑥 + 1 − 4 − 4𝑦 2 > 0
12𝑥 2 − 3 − 4𝑦 2 > 0
12𝑥 2 − 4𝑦 2 > 3
𝑥2 𝑦2
− > 3/3
3/12 3/4
𝑥 2 3𝑦 2
− >1
4 4
Gràfico:
ii)- 0 < |𝑧 − 2 + 𝑖| < 4
0 < |𝑥 + 𝑦𝑖 − 2 + 𝑖| < 4

0 < √(𝑥 − 2)2 + (𝑦 + 1)2 < 4

0 < (𝑥 − 2)2 + (𝑦 + 1)2 < 42
Gràfico:

C)- SEA LA FUNCIÒN 𝑾 = 𝒇(𝒛)/𝒘 = 𝒛𝟐 + 𝒛 − 𝟏. Si x=0¿QÙE PUNTOS CORRESPONDEN

SEGÙN W?

𝑾 = 𝒇(𝒛)/𝒘 = 𝒛𝟐 + 𝒛 − 𝟏; 𝒔𝒊 𝑿 = 𝟎
𝒘 = 𝒛𝟐 + 𝒛 − 𝟏
𝒘 = (𝒙 + 𝒚𝒊)𝟐 + (𝒙 + 𝒚𝒊) − 𝟏
𝒘 = 𝒙𝟐 + 𝟐. 𝒙. 𝒚𝒊 − 𝒚 + 𝒙 + 𝒚𝒊 − 𝟏
𝑤 = (𝒙𝟐 + 𝒙 − 𝒚 − 𝟏) + 𝒊(𝟐𝒙. 𝒚 + 𝒚)
𝒖 = 𝒙𝟐 + 𝒙 − 𝒚 − 𝟏
𝒗 = 𝟐𝒙. 𝒚 + 𝒚
Si x=0
𝒖(𝟎, 𝒚) = −𝒚 − 𝟏 = 𝒖 (𝟏)
𝒗(𝟎, 𝒚) = +𝒚 = 𝒗(𝟐)
De (2)= +𝒚 = 𝒗
Sustituye en (1)
𝒖 = −𝒗 − 𝟏
Función en el Plano XY

Funciòn en el Plano W
2- INTERPOLACIÓN
A) DEFINE INTERPOLACIÓN DE POLINOMIOS Y SUS USOS E IMPORTANCIA. ¿SE DEFINEN
OTROS TIPOS DE “INTERPOLACIÓN FUNCIONAL”? NOMBRA CUÁLES Y EJEMPLIFICA
RELACIONANDO CON ALGUNAS DE LAS ACTIVIDADES REALIZADAS EN CLASE.
Interpolación de Polinomios:
En análisis numérico, la interpolación polinómica es una técnica de interpolación de un conjunto
de datos o de una función por un polinomio. Es decir, dado cierto número de puntos obtenidos
por muestreo o a partir de un experimento se pretende encontrar un polinomio que pase por
todos los puntos.

Definición
Dada una función f de la cual se conocen sus valores en un número finito de abscisas
𝑋0 , 𝑋1 , 𝑋2 , . . . , 𝑋𝑛 , se llama interpolación polinómica al proceso de hallar un polinomio 𝑃𝑚 (𝑋) de
grado menor o igual a m, cumpliendo. 𝑃𝑚 (𝑋𝑘 ) = 𝑓(𝑋𝑘 ), ∀𝑘 = 0,1, . . . , 𝑚.
A este polinomio se le llama Polinomio interpolador de grado m de la función f.
Se dispone de varios métodos generales de interpolación polinómica que permiten aproximar
una función por un polinomio de grado m. El primero de estos polinomios es el método de las
diferencias divididas de Newton. Otro de los métodos es la interpolación de Lagrange, y por
último, la interpolación de Hermite.

Método de las diferencias divididas de Newton

Sea 𝑓𝑛 una variable de n elementos y sea 𝑋𝑛 otra variable de n elementos los cuales
corresponden, por parejas, a la imagen u ordenada y abscisa de los datos que se quieran
interpolar, respectivamente, tales que:
𝑓(𝑋𝑘 ), ∀𝑘 = 0,1, . . . , 𝑛
Este método es muy algorítmico y resulta sumamente cómodo en determinados casos, sobre
todo cuando se quiere calcular un polinomio interpolador de grado elevado.
El polinomio de grado 𝑛 − 1 resultante tendrá la forma
𝑛−1

∑ 𝑎𝑗 𝑔𝑗 (𝑋)
𝑗=0
definiendo 𝑔𝑗 (𝑋) como
𝑗−1
𝑔𝑗 (𝑋) = ∏𝑖=0 (𝑥 − 𝑥𝑖 )
y definiendo 𝑎𝑗 como
𝑎0 = 𝑓(𝑋0 ), 𝑎1 = 𝑓(𝑋0 , 𝑋1 ), . . . , 𝑎𝑗 = 𝑓 (𝑋0 , 𝑋1 , . . . , 𝑋𝑗−1 , 𝑋𝑗 )

Los coeficientes 𝑎𝑗 son las llamadas diferencias divididas.

Una vez se hayan realizado todos los cálculos, nótese que hay (muchas) más diferencias divididas
que coeficientes 𝑎𝑗 . El cálculo de todos los términos intermedios debe realizarse simplemente
porque son necesarios para poder formar todos los términos finales. Sin embargo, los términos
usados en la construcción del polinomio interpolador son todos aquellos que involucren a 𝑋𝑛 .
Estos coeficientes se calculan mediante los datos que se conocen de la función 𝑓.
𝑓 (𝑋0 , 𝑋1 , . . . , 𝑋𝑗−1 , 𝑋𝑗 queda definido, como:

Se muestra ahora una tabla mnemotécnica con las diferencias divididas de una cierta función
dada para construir un polinomio interpolador de grado 2:

Ejemplo:
Calcular la tabla de diferencias divididas finitas con los siguientes datos:

X -3 -2 0 4

Y 5 8 4 2

Y usar la información en la tabla, para construir el polinomio de interpolación de Newton.

8−5 −2−3 −1,75−−0,25
-3 5 =3 = −1,666 = −0.3
−2−−3 0−−3 3−−2
9−6 −6−1
-1 6 =1 = −1,75
2−−1 3−−1
3−9
2 9 = −6
3−2
3 3

Por lo tanto el polinomio de interpolación de Newton nos queda:

𝑓(𝑥) = 5 + 3(𝑥 + 3) − 1,666(𝑥 + 3)(𝑥 + 2) − 0,20238(𝑥 + 3)(𝑥 + 2)(𝑥)

𝑓(𝑥) = −0,20238 𝑥 3 − 2,6779 𝑥 2 − 6,54428 𝑥 + 4,004

Interpolación de Lagrange:
Sea 𝑓 la función a interpolar, sean 𝑋0 , 𝑋1 , 𝑋2 , . . . , 𝑋𝑛 las abscisas conocidas de 𝑓 y sean 𝑓0 ,𝑓1,...𝑓𝑛 ,
los valores que toma la función en esas abscisas, el polinomio interpolador de grado 𝑛 de
Lagrange es un polinomio de la forma:

∑ 𝑓𝑗 𝑙𝑗 (𝑋), 𝑛 ≤ 𝑚
𝑗=0

donde 𝑙𝑗 (𝑋) son los llamados polinomios de Lagrange, que se calculan de este modo:
𝑋 − 𝑋𝑖 (𝑋 − 𝑋0 )(𝑋 − 𝑋1 ). . . (𝑋 − 𝑋𝑗−1 )(𝑋 − 𝑋𝑗+1 ). . . (𝑋 − 𝑋𝑛 )
𝑙𝑗 (𝑋) = ∏ =
𝑋𝑗 − 𝑋𝑖 (𝑋𝑗 − 𝑋0 ) (𝑋𝑗 − 𝑋1 ). . . (𝑋𝑗 − 𝑋𝑗−1 )(𝑋𝑗 − 𝑋𝑗+1 ). . . (𝑋𝑗 − 𝑋𝑛 )
𝑖≠𝑗

Los coeficientes 𝑙𝑗 (𝑋) están bien definidos y son siempre distintos de cero.
ejemplo:

Calcular el polinomio de Lagrange usando los siguientes datos:

X 1 3 5 47

Y -2 1 2 -3

𝒇(𝒙) = 𝒚𝟎 𝒍𝟎 (𝒙) + 𝒚𝟏 𝒍𝟏 (𝒙) + 𝒚𝟐 𝒍𝟐 (𝒙) + 𝒚𝟑 𝒍𝟑 (𝒙)

𝒇(𝒙) = −𝟐𝒍𝟎 (𝒙) + 𝒍𝟏 (𝒙) + 𝟐𝒍𝟐 (𝒙) − 𝟑𝒍𝟑 (𝒙)
donde:
(𝒙 − 𝟑)(𝒙 − 𝟓)(𝒙 − 𝟕) (𝒙 − 𝟑)(𝒙 − 𝟓)(𝒙 − 𝟕) 𝟏
𝒍𝟎 (𝒙) = = = (−𝒙𝟑 + 𝟏𝟓𝒙𝟐 − 𝟏𝟕𝒙 + 𝟏𝟎𝟓)
(−𝟐)(−𝟒)(−𝟔) (−𝟒𝟖) 𝟒𝟖
(𝒙 − 𝟏)(𝒙 − 𝟓)(𝒙 − 𝟕) (𝒙 − 𝟏)(𝒙 − 𝟓)(𝒙 − 𝟕) 𝟏 𝟑
𝒍𝟏 (𝒙) = = = (𝒙 − 𝟏𝟑𝒙𝟐 + 𝟒𝟕 𝒙 − 𝟑𝟓)
(𝟐)(−𝟐)(−𝟒) (𝟏𝟔) 𝟏𝟔
(𝒙 − 𝟏)(𝒙 − 𝟑)(𝒙 − 𝟕) (𝒙 − 𝟏)(𝒙 − 𝟑)(𝒙 − 𝟕) 𝟏 𝟑
𝒍𝟐 (𝒙) = = = (𝒙 − 𝟏𝟏𝒙𝟐 + 𝟑𝟏𝒙 − 𝟐𝟏)
(𝟒)(𝟐)(−𝟐) (𝟏𝟔) 𝟏𝟔
(𝒙 − 𝟏)(𝒙 − 𝟑)(𝒙 − 𝟓) (𝒙 − 𝟏)(𝒙 − 𝟑)(𝒙 − 𝟓) 𝟏 𝟑
𝒍𝟑 (𝒙) = = = (𝒙 − 𝟗𝒙𝟐 + 𝟐𝟑𝒙 − 𝟏𝟓)
(𝟔)(𝟒)(𝟐) (𝟒𝟖) 𝟒𝟖
𝒇(𝒙) = −𝟐𝒍𝟎 (𝒙) + 𝒍𝟏 (𝒙) + 𝟐𝒍𝟐 (𝒙) − 𝟑𝒍𝟑 (𝒙) =
𝒙𝟑 𝟓𝒙𝟐 𝟏𝟕𝒙 𝟑𝟓 𝒙𝟑 𝟏𝟑𝒙𝟐 𝟒𝟕𝒙 𝟑𝟓 𝒙𝟑 𝟏𝟏𝒙𝟐 𝟑𝟏𝒙 𝟐𝟏
𝒇(𝒙) = − + − + − + − + − + −
𝟐𝟒 𝟖 𝟐𝟒 𝟖 𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔
𝒙𝟑 𝟑𝒙𝟐 𝟐𝟑𝒙 𝟓
+ − + −
𝟒𝟖 𝟏𝟔 𝟒𝟖 𝟏𝟔
 𝟑𝒙𝟑 𝟑𝟕𝒙𝟐 𝟗𝟕𝒙 𝟏𝟑𝟏
𝒇(𝒙) = − + −
𝟏𝟔 𝟏𝟔 𝟏𝟔 𝟏𝟔
B) DEFINE UN POLINOMIO QUE DESCRIBA UN ESBOZO DEL SIGUIENTE DIBUJO
Esbozo de Imagen en con polinomios realizados en GeoGebra:

Función Polinómica Puntos

comprendidos:
A(x) = -0.01x⁵ + 0.15x⁴ - 1.24x³ + 4.8x² - 7.5x + 11.6 Si(2 ≤ x ≤ 7)
B(x) = -0.16x³ + 3.6x² - 26.53x + 71.54 Si(7 ≤ x ≤ 8.73)
C(x) = -12.13x² + 211.29x - 913.31 Si(8.73 ≤ x ≤ 8.99)
D(x) = 1.89x³ - 53.72x² + 507.87x - 1589.38 Si(8.99 ≤ x ≤ 9.59)
E(x) = 0.13x² - 1.33x + 5.38 Si(5 ≤ x ≤ 9.59)
F(x) = 0.04x⁴ - 0.65x³ + 4.05x² - 11.33x + 14.2 Si(1.33 ≤ x ≤ 5)
G(x) = -14.11x² + 48.58x - 34.73 Si(1.33 ≤ x ≤ 1.68)
H(x) = -0.87x² - 0.14x + 9.75 Si(1.36 ≤ x ≤ 1.68)
I(x) = -0.11x² + 0.45x + 7.55 Si(1.36 ≤ x ≤ 2)
J: (x - 2.94)² + (y - 7.56)² = 0.09
K(x) = -0.13x⁴ + 1.84x³ - 9.59x² + 22.21x - 10.84 Si(1.68 ≤ x ≤ 5.51)
L(x) = 0.07x² - 0.71x + 5.33 Si(5.85 ≤ x ≤ 9.59)
M(x) = -0.15x³ + 2.76x² - 16.89x + 40.4 Si(5.51 ≤ x ≤ 7.38)
N(x) = 6.17x³ - 130.18x² + 915.22x - 2139.67 Si(6.89 ≤ x ≤ 7.38)
Ñ(x) = 4.67x² - 61.06x + 203.83) Si(6.49 ≤ x ≤ 6.89)
O(x) = 2.39x³ - 42.04x² + 247.01x - 481.12 Si(5.85 ≤ x ≤ 6.49)
P(x) = -0.13x³ + 2.72x² - 18.29x + 43.28 Si(4.27 ≤ x ≤ 5.85)
Q(x) = ) 0.96x³ - 15.08x² + 78.85x - 131.98) Si(4.27 ≤ x ≤ 5.51
3- ELIJE UNA SITUACIÓN EXTRAÍDA DE LOS LABORATORIOS DE QUÍMICA PARA REALIZAR
UNA APLICACIÓN DE LOS TEMAS DESARROLLADOS: NÚMEROS COMPLEJOS-
INTERPOLACIÓN- DERIVACIÓN E INTEGRACIÓN NUMÉRICA.

Integración Numérica:
Se sacó un problema de aplicación en el laboratorio de Química, más precisamente de la
cátedra de Fisicoquímica, en la cual necesitamos calcular las fugacidades de diferentes gases
de manera práctica y considerando la dependencia de los parámetros.

Un poco de teoría:
La fugacidad es una propiedad termodinámica auxiliar, en la que indica la presión parcial efectiva
de un gas real, que actúa como factor de corregir y transformar una ecuación teórica ideal en
una situación práctica real.

La fugacidad de un gas depende de la presión, temperatura, volumen molar (Vi) o factores de

compresibilidad (Z). Para evaluar las fugacidades en el laboratorio de un gas determinado, se
realizan mediciones de presiones para una sola temperatura a 0°C.
La ecuación general en para una tabla de datos experimental será:
𝑃
𝑉𝑖  −  𝑉𝑖 ⋅
𝑙𝑛 𝑙𝑛  𝑓𝑖  =   𝑙𝑛 𝑙𝑛  𝑃  +   ∫ 𝑑𝑃   =
0 𝑅𝑇
𝑃 𝑉𝑖 −𝑉𝑖⋅
𝑓𝑖  = 𝑒 [𝑙𝑛𝑙𝑛  𝑃 +  ∫0 ( 𝑅𝑇 )𝑑𝑃]

𝑃 𝑍 − 1
ó  𝑓𝑖 = 𝑒 [𝑙𝑛𝑙𝑛  𝑃 +∫0 ( 𝑃
)𝑑𝑃]

Para distintos valores de presión obtendremos distintos valores de fugacidad. Podemos obtener
datos a partir de valores de volumen molar o con factores de compresibilidad (Z). Para ello
necesitamos evaluar la integral definida realizada por los valores de tabla y que satisfagan las
ecuaciones anteriores.

Ejemplo:
A partir de los siguientes datos del factor de compresibilidad para CO2 a 150°C, calcule la
fugacidad del CO2.
P (bar) Z (Z-1) / P
10 0,985 -0,0015
40 0,942 -0,00145
100 0,869 -0,00131
300 0,765 -0,000793
500 0,91 -0,00018
Realizamos los gráficos de (Z-1) / P con respecto a P:

Regla de Simpson:
𝑏
𝛥𝑥
∫ 𝑓(𝑥). 𝑑𝑥   = (𝑦𝑜  +  4. 𝑦1  + 𝑦2)
𝑎 3
𝑃 (𝑍−1)
Para poder aproximar el área debajo de la curva: ∫0 𝑃
. 𝑑𝑃

(𝑏 − 𝑎) (500 − 10)
𝛥𝑥 = = = 245
2 2
Con la ecuación polinómica de aproximación P(L) (Polinomio de Lagrange) calculo y1

𝑦𝑜  = −1,5.10−3 ;  𝑦1 = −2,2806.10−3 ;  𝑦2 = −1,8.10−4

𝑏
𝛥𝑥 245
∫ 𝑓(𝑥). 𝑑𝑥 = (𝑦𝑜 + 4. 𝑦1 + 𝑦2)  =   [−1,5.10−3 + 4 ⋅ (−2,2806.10−3 ) + (−1,8.10−4 )] =
𝑎 3 3
𝑃 (𝑍
− 1)
∫ . 𝑑𝑃  = 0,8822
0 𝑃

Ahora calculamos la fugacidad de CO2 de 10 a 500 bar y 150°C que nos pide el problema:
𝑃 𝑍−1
[𝑙𝑛𝑙𝑛  𝑃 +  ∫0 ( 𝑃 ).𝑑𝑃]
𝑓𝑖 = 𝑒 = 𝑒 [𝑙𝑛𝑙𝑛  500 + 0,8822]
= 1208,10
Método de interpolación de Lagrange:
Se busca el polinomio de aproximación por el método de Lagrange.

P (bar) Z (Z-1) / P
10 0,985 -0,0015
40 0,942 -0,00145
100 0,869 -0,00131
300 0,765 -0,000793
500 0,91 -0,00018

(𝑥 − 40) ⋅ (𝑥 − 100) ⋅ (𝑥 − 300) ⋅ (𝑥 − 500) 𝑥 4 − 940𝑥 3 + 266000𝑥 2 − 24200000𝑥  + 600000000

𝑙𝑜 = =
(10 − 40) ⋅ (10 − 100) ⋅ (10 − 300) ⋅ (10 − 500) 383670000

(𝑥 − 10) ⋅ (𝑥 − 40) ⋅ (𝑥 − 300) ⋅ (𝑥 − 500) 𝑥 4 − 850𝑥 3 + 190400𝑥 2 − 7820000𝑥  + 60000000

𝑙2 = =
(100 − 10) ⋅ (100 − 40) ⋅ (100 − 300) ⋅ (100 − 500) 432000000

(𝑥 − 10) ⋅ (𝑥 − 40) ⋅ (𝑥 − 100) ⋅ (𝑥 − 500) 𝑥 4 − 650𝑥 3 + 80400𝑥 2 − 2740000𝑥  + 20000000

𝑙3 = =
(300 − 10) ⋅ (300 − 40) ⋅ (300 − 100) ⋅ (300 − 500) −3016000000

(𝑥 − 10) ⋅ (𝑥 − 40) ⋅ (𝑥 − 100) ⋅ (𝑥 − 300) 𝑥 4 − 450𝑥 3 + 50400𝑥 2 − 1660000𝑥  + 12000000

𝑙4 = =
(500 − 10) ⋅ (500 − 40) ⋅ (500 − 100). (500 − 300) 18032000.103

Hago el polinomio de lagrange:

𝑃(𝑙) = 𝑦1. 𝑙1 + 𝑦2 ⋅ 𝑙2 + 𝑦3 ⋅ 𝑙3 + 𝑦4. 𝑙4 = 𝑃(𝑙1) + 𝑃(𝑙2) + 𝑃(𝑙3) + 𝑃(𝑙4)

𝑃(𝑙1) = −3,77929.10−12 𝑥 4 + 3,55253.10−9 𝑥 3 − 1,00529.10−6 𝑥 2 + 0,0000914588𝑥 − 0,00226757

𝑃(𝑙2) = −3,03241.10−13 𝑥 4 + 2,57755.10−10 𝑥 3 − 5,7737.10−8 𝑥 2 + 0,0000237134𝑥 − 0,0000181944

𝑃(𝑙3) = 2,62931.10−13 𝑥 4 − 1,70905.10−10 𝑥 3 + 2,11397.10−8 𝑥 2 − 7,20431.10−7 𝑥 + 5,25862.10−6

𝑃(𝑙4) = −9,98225.10−15 𝑥 4 + 4,49201.10−12 𝑥 3 − 5,03106.10−10 𝑥 2 + 1,65705.10−8 𝑥  −  1,19787.10−7

𝑃(𝑙) = 𝑃(𝑙1) + 𝑃(𝑙2) + 𝑃(𝑙3) + 𝑃(𝑙4)

𝑃(𝑙) = −3,82958225 ⋅ 10−12 𝑥 4 + 3,6438 ⋅ 10−9 𝑥 3 − 1,042366 ⋅ 10−6 𝑥 2 + 1,14468 ⋅ 10−4 𝑥 − 2,280625.10−3

4- ECUACIONES DIFERENCIALES ELEGIR UNA SITUACIÓN PROBLEMÁTICA QUE SE RESUELVA

CON UNA ECUACIÓN DIFERENCIAL, DE CUALQUIER ORDEN, LUEGO RESOLVERLA CON
“MÉTODO CLÁSICO” (COMO SE RESOLVIERON EN AMII). VERIFICAR POR TRANSFORMADA DE
LAPLACE.

Situación problemática de ecuaciones diferenciales:

La pendiente de una curva en cualquier punto(x,y) vale x + 2y. Determinar la ecuación de dicha
curva si pasa por el origen.
𝑑𝑦
Datos: 𝑑𝑥 = 𝑥 + 2𝑦
 𝑑𝑦
1. EDO en forma estándar: 𝑦′ + 𝑃(𝑥) + 𝑄(𝑥) a 𝑑𝑥 − 2𝑦 = 𝑥
2. Obtengo el factor integrante: 𝑉(𝑥) = 𝑒 − ∫ 𝑃(𝑥).𝑑𝑥 𝑉(𝑥) = 𝑒 − ∫ −2.𝑑𝑥 = 𝑒 2𝑥
𝑑𝑢 𝑑𝑢
3. Reemplazo: ∗ 𝑉(𝑥) = 𝑄(𝑥) a ∗ 𝑒 2𝑥 = 𝑥 a ∫ 𝑑𝑢 = ∫ 𝑒 −2𝑥 𝑥𝑑𝑥
𝑑𝑥 𝑑𝑥
𝑒 −2𝑥
4. Integro por partes: 𝑢 = 𝑣 ; 𝑑𝑢 = 𝑑𝑥 ; 𝑑𝑣 = 𝑒 −2𝑥 𝑑𝑥; 𝑣 = − 2

∫ 𝑢 . 𝑑𝑣 = 𝑢. 𝑣 − ∫ 𝑣 . 𝑑𝑢

𝑥. 𝑒 −2𝑥 −𝑒 −2𝑥
∫ 𝑥. 𝑒 −2𝑥 𝑑𝑥 = − −∫ 𝑑𝑥
2 2

resuelvo la integral por sustitución:

𝑑𝑢 1
𝑢 = −2𝑥; =− → ∫ 𝑒 𝑢 𝑑𝑢 = 𝑒 𝑢 = 𝑒 −2𝑥
𝑑𝑥 2
𝑥 1
reemplazo en la original: ∫ 𝑥. 𝑒 −2𝑥 𝑑𝑥 = − 2 𝑒 −2𝑥 − 4 𝑒 −2𝑥 + 𝐶

−𝑥 1 𝑥 1
5. expreso 𝑦 = 𝑈(𝑥) ∗ 𝑉(𝑥) = ( 2 𝑒 −2𝑥 − 4 𝑒 −2𝑥 + 𝑐) ∗ 𝑒 2𝑥 = − 2 − 4 + 𝑐. 𝑒 2𝑥
1 1
6. Despejo constante: si f (0) = 0 a 𝑦(0) = − 4 + 𝑐 𝑐=4
𝑥 1 1
7. 𝑦 = − 2 − 4 + 4 . 𝑒 2𝑥

Verificación por Transformada de Laplace:

𝑑𝑦
= 𝑥 + 2𝑦 → 𝑓′(𝑡) − 2𝑓(𝑡) = 𝑡 𝑐𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 𝑖𝑛𝑖𝑐𝑖𝑎𝑙 𝑓(0) = 0
𝑑𝑥

1. Transformo:
1 1 1
𝑠𝐹(𝑠) − 2𝐹(𝑠) = 𝑠2 → 𝐹(𝑠)(𝑠 − 2) = 𝑠2 𝐹(s)=𝑠2(𝑠−2)

2. Aplico fracciones parciales:

1 𝐴 𝐵 𝐶 −1 1 1
𝐹(𝑠) = = 2+ + = + −
𝑠 2 (𝑠 − 2) 𝑠 𝑠 (𝑠 − 2) 2𝑠 2 4(𝑠 − 2) 4𝑠

3. Inversa de Laplace:
−1 1 1 −1 1 1
𝑓(𝑡) = 2𝑠2 + 4(𝑠−2) − 4𝑠 = ∗ 𝑡 + 4 ∗ 𝑒 2𝑡 − 4 ∗ 1
2
4. Verifico:
1 1 1 1 1
𝑓′(𝑡) − 2𝑓(𝑡) = 𝑡 → − 2 + 2 𝑒 2𝑡 − 2(− 2 𝑡 + 4 𝑒 2𝑡 − 4) = 𝑡

1 1 1 1
− + 𝑒 2𝑡 + 𝑡 − 𝑒 2𝑡 + = 𝑡 𝑡=𝑡
2 2 2 2