Análisis y Métodos Numéricos:

Derivación, Diferencias Finitas y EDP

Deisy Chaves
[email protected]
 Diferenciación Numérica

 La diferenciación numérica aproxima la derivada de una

función en un punto cuando solo se conocen los valores
de la función en algunos puntos

 La diferenciación numérica es útil en la solución de

ecuaciones diferenciales y los métodos de las
diferencias finitas
 Serie de Taylor

 La serie de Taylor trunca, después de un

termino de orden N, se expresa como:
2
hn−1
f (x) = f (a) + hf ′ (a) +
h ′′
2!
f (a) + ... +
(n −1)!
( )
f (n−1) (a) + O hn

donde N es el orden del polinomio que

aproxima la función y

( )
n
h
O hn = f n (ξ )
n!
 Aproximación de la Primera Derivada

 Veamos como se aproxima la primera

derivada de una función f(x) cuando se
conocen los valores de la función en dos
puntos x y x+h, usando el desarrollo de la serie
de Taylor truncada en N=1
2
+ h ′′
f (x + h) = f (x) + hf ′(x) f (ξ )
2

f (x + h) − f (x)
f ′(x) ≈
h
 Aproximación de la Primera Derivada

 Entonces el error de aproximación de la primera

derivada, cuando se usan los dos primeros
términos del desarrollo de series de Taylor, es

h2 2
E= f ′′(ξ ) ≤ h max f ′′(x)
2 2 a≤ x≤b

 Si conocemos los valores que toma la función

en los puntos x-h y x+h podemos obtener una
mejor aproximación de la derivada de la función
en x usando el desarrollo de la serie de Taylor
truncada en N=2
 Aproximación de la Primera Derivada
2 3
+h + h ′′′
f (x + h) = f (x) + hf ′(x) f ′′(h) f (ξ )
2! 3!
2 3
+h − h ′′′
f (x − h) = f (x) − hf ′(x) f ′′(h) f (η )
2! 3!
h 3 ′′′
f (x + h) − f (x − h) = 2h f ′(x) + ( f (ξ ) + f ′′′(η))
3!
La primera derivada se aproxima como:
f (x + h) − f (x − h)
f ′(x) ≈
2h
El error de la aproximación esta dado por
h 2 ′′′ φ 2
E= f ( ) ≤ h max f ′′′(x)
6 6 a≤ x≤b
 Aproximaciones de la Primera Derivada

 Con la notación usada en para el calculo de las

diferencias divididas tenemos que:
Formula de las diferencias
f (xi+1 ) − f (x) divididas hacia adelante
f ′(xi ) ≈
h
Formula de las diferencias
f (xi ) − f (x i−1 ) divididas hacia atrás
f ′(xi ) ≈
h

f (x i+1 ) − f (x i−1 ) Formula de las diferencias

f ′(xi ) ≈ divididas centradas hacia
2h
adelante
 Ejemplo

 Calcular la primera derivada en x=0.5 con

h=0.25, de la función
f (x) = −0.1x4 −0.15x3 − 0.5x2 − 0.25x +1.2
Solución:
xi−2 = 0 f (xi−2 ) = 1.2
xi−1 = 0.25 f (xi−1 ) = 1.10527344
xi = 0.5 f (xi ) = 0.934375
xi+1 = 0.75 f (xi+1 ) = 0.65214844
xi+2 = 1 f (xi+2 ) = 0.2
 Ejemplo

 Usando las siguientes expresiones tenemos

que:
0.65214844 − 0.934375
f ′(0.5) ≈ = −1.12890625
0.25
Hacia adelante

0.934375 −1.10527455
f ′(0.5) ≈ = −0.68359375
0.25
Hacia atrás

0.65214844 −1.10527344
f ′(0.5) ≈ = −0.90625
0.5 Centradas
 Ejemplo: Errores

 Calculemos el error de aproximación usando

diferencias divididas hacia adelante y hacia
atrás
f ′′(0) = −1 E=
h 2 ′′ ξ 2
f ( ) ≤ h max f ′′(x)
f ′′(0.25) = −1.24375 2 2 a≤ x≤b

f ′′(0.5) = −1.6
f ′′(0.75) = −2.18125 2
h ′′ ) = 0.05
E= f (xm
f ′′(1) = −3.1 2
h2
E= max
f ′′(x) = 0.096875
2 a≤ x≤b
 Ejemplo: Errores

 Calculemos el error de aproximación usando

diferencias divididas centradas
f ′′′(0) = −0.9 h2
f ′′′(φ )
2
E= h max f ′′′(x)
6 ≤ 6 a≤ x≤b

f ′′′(0.25) = −1.5
2
h f ′′′(x m ) = 0.021875
f ′′′(0.5) = −2.1 E=
6
f ′′′(0.75) = −2.7
h2
E= max
f ′′′(1) = −3.3 f ′′′(x) = 0.034375
6 a≤ x≤b
 Resumiendo

f ′(0.5)adelante = −1.1289 Eaprox = 0.05 Ereal = 0.2164

f ′(0.5)atras = −0.68359 Eaprox = 0.05 Ereal = 0.2289
f ′(0.5)centrada = −0.906 Eaprox = 0.02812 Ereal = 0.00625

 El error real se calcula usando el valor

verdadero de la primera derivada evaluada en el
punto 0.5
f (x) = −0.1x4 −0.15x3 − 0.5x2 − 0.25x +1.2

f ′(0.5)=−0.9125
 Aproximación de la Segunda Derivada

 Veamos como se aproxima la segunda

derivada de una función f(x) cuando se
conocen los valores de la función en los
puntos x-h, x y x+h, usando el desarrollo de la
serie de Taylor truncada en N=3
h2 h3 h4
f (x + h) = f (x) + hf ′(x) + f ′′(x) + f ′′′(x) + f (4) (ξ )
2! 3! 4!

h2 h3 h4
f (x − h) = f (x) − hf ′(x) + f ′′ (x) − f ′′′(x) + f (4)
(η)
2! 3! 4!
 Aproximación de la Segunda Derivada

h2 h3 h4
f (x + h) = f (x) + hf ′(x) + f ′′(x) + f ′′′(x) + f (4) (ξ )
2! 3! 4!
+h 2 −h 3 +h 4
f (x − h) = f (x) − hf ′(x) f ′′(x) f ′′′(x) f (4) (η )
2! 3! 4!

f (x + h) + f (x − h) = 2 f (x) +
h2
2!
f ′′(x) +
2h 4 (4)
4!
(
f (ξ ) − f (4) (η))
La segunda derivada se aproxima con
f (x + h) − 2 f (x) + f (x − h)
f ′′ (x) ≈
h2
El error de la aproximación esta dado por

E= h 2
h 2
f (4) (θ ) ≤ max (4)
f (x)
12 12 a≤ x≤b
 Aproximación de la Segunda Derivada

 Con la notación usada en para el calculo de las

diferencias divididas tenemos que:
Formula de las diferencias
′′(x f (x i+2 ) − 2 f (x i+1 ) + f (x i ) divididas hacia adelante
f i) ≈
h2
Formula de las diferencias
′′(x f (xi ) − 2 f (x i−1 ) + f (x i−2 ) divididas hacia atrás
f i) ≈
h2
Formula de las diferencias
′′(x
f (x i+1 ) − 2 f (x i ) + f (x i−1 ) divididas centradas hacia
f i) ≈ adelante
h2
 Ejercicio

 Calcular la segunda derivada en x=0.5, con

h=0.25 de la función
f (x) = −0.1x4 −0.15x3 − 0.5x2 − 0.25x +1.2

Recordemos que:
xi−2 = 0 f (xi−2 ) = 1.2 f ``(0.5)adelante = -2.3625
xi−1 = 0.25 f (xi−1 ) = 1.103516 f ``(0.5)atras = -0.304687
xi = 0.5 f (xi ) = 0.925 f ``(0.5)centrada =−1.76

xi+1 = 0.75 f (xi+1 ) = 0.6363281

xi+2 = 1 f (xi+2 ) = 0.2
Derivadas de Orden Superior Hacia Adelante
Derivadas de Orden Superior Hacia Atrás
Derivadas de Orden Superior Centradas
 Derivadas Parciales

 La derivada parcial de una función f(x, y) se logra

aplicando una derivada ordinaria a la función con
respecto a una de las variables independientes,
considerando a la otra constante

 Ejemplo: Calcular la primera derivada con respecto

a x, en x=0.5,y=1 con h=0.25, de la función
f (x,y) = −x4 y3-0.15x3y − 0.5x2 − 0.25x +1.2

𝑑𝑑𝑑𝑑 𝑓𝑓 𝑥𝑥𝑖𝑖+1 , 𝑦𝑦𝑖𝑖 − 𝑓𝑓(𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑖𝑖 )

𝑥𝑥 , 𝑦𝑦 =
𝑑𝑑𝑑𝑑 𝑖𝑖 𝑖𝑖 ℎ
 Derivadas Parciales

 La derivada parcial de una función f(x, y) se logra

aplicando una derivada ordinaria a la función con
respecto a una de las variables independientes,
considerando a la otra constante

 Ejemplo: Calcular la primera derivada con respecto

a x, en x=0.5,y=1 con h=0.25, de la función
f (x,y) = −x4 y3-0.15x3y − 0.5x2 − 0.25x +1.2
𝑑𝑑𝑑𝑑 𝑓𝑓 𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑖𝑖+1 − 𝑓𝑓(𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑖𝑖 )
𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑖𝑖 =
𝑑𝑑𝑑𝑑 ℎ

𝑑𝑑𝑑𝑑 𝑓𝑓 0.5,1.25 − 𝑓𝑓(0.5,1)

𝑥𝑥 , 𝑦𝑦 =
𝑑𝑑𝑑𝑑 𝑖𝑖 𝑖𝑖 0.25
 Ecuaciones Diferenciales

 Toda ecuación que establece la dependencia de

una variable respecto a otra u otras mediante
derivadas es una ecuación diferencial
 Ecuaciones Diferenciales Ordinarias (EDO)

 Las ecuaciones diferenciales ordinarias están

compuestas por la cantidad que se está derivando,
conocida como variable dependiente, la cantidad
con respecto a la cual se está derivando, conocida
como variable independiente
 Ejemplo: ecuación de razón, usada para calcular
la velocidad ∂v c
=g− v
∂t m
donde g es la constate gravitacional, m es la masa
y c es el coeficiente de arrastre
 Ecuaciones Diferenciales Ordinarias (EDO)

 Las ecuaciones diferenciales ordinarias están

compuestas por la cantidad que se está derivando,
conocida como variable dependiente, la cantidad
con respecto a la cual se está derivando, conocida
como variable independiente
 Ejemplo: ecuación de razón, usada para calcular
la velocidad ∂v c
=g− v
∂t m
donde g es la constate gravitacional, m es la masa
y c es el coeficiente de arrastre
 v es la variable dependiente
 t es la variable independiente
 Ecuaciones diferenciales Parciales (EDP)

 Las ecuaciones diferenciales parciales

involucra dos o mas variables independientes
 Ecuaciones diferenciales Parciales (EDP)

 Las ecuaciones diferenciales parciales

involucra dos o mas variables independientes

 Ejemplo: la siguiente ecuación describe la

dependencia de r respecto de x, y y z

∂r ∂r ∂r
+ − 2 xy = 1
∂x ∂y ∂z
 Ecuaciones diferenciales Parciales (EDP)

 El orden de una ecuación diferencial parcial es

el orden de la derivada mayor que aparece en
ella
 Ejemplo: EDP de orden 2
𝜕𝜕2 u 𝜕𝜕2 u 𝜕𝜕𝜕
A 2 + 𝐵𝐵 + 𝐶𝐶 2 + 𝐹𝐹 = 0
𝜕𝜕x 𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝑦𝑦

Clasificación:
o Si B2 – 4AC < 0, la ecuación se llama elíptica
o Si B2 – 4AC = 0, la ecuación se llama parabólica
o Si B2 – 4AC > 0, la ecuación se llama hiperbólica
 EPD Eliptica

 Ecuación de Laplace utilizada para modelar

problema de la transmisión de calor en una
placa

𝜕𝜕 2 𝜕𝜕 2 𝑓𝑓 = 0 : Laplace
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 𝑓𝑓(𝑥𝑥, 𝑦𝑦) �
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦 𝑓𝑓 ≠ 0: Poisson
 Ejercicio (EDP)

Se tiene una delgada placa de acero rectangular

de 10 cm x 20 cm. Uno de los bordes de 10 cm se
mantiene a 100⁰c y los otros tres bordes a 0⁰c ,
estas temperaturas son producto de una fuente
de calor colocada debajo de ella.
 Ejercicio (EDP)

Se tiene una delgada placa de acero rectangular

de 10 cm x 20 cm. Uno de los bordes de 10 cm se
mantiene a 100⁰c y los otros tres bordes a 0⁰c ,
estas temperaturas son producto de una fuente
de calor colocada debajo de ella.

𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦
Ejercicio (EDP)

𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦
Ejercicio (EDP)

𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦
Ejercicio (EDP)

𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦
 Ejercicio (EDP)

Se tiene una delgada placa de acero rectangular de 10 cm x

20 cm. Uno de los bordes de 10 cm se mantiene a 100⁰c y
los otros tres bordes a 0⁰c, estas temperaturas son producto
de una fuente de calor colocada debajo de ella. Use ∆x =5 y
∆y=5

Cuales son las condiciones

de Frontera?
 Ejercicio (EDP)

Se tiene una delgada placa de acero rectangular de 10 cm x

20 cm. Uno de los bordes de 10 cm se mantiene a 100⁰c y
los otros tres bordes a 0⁰c, estas temperaturas son producto
de una fuente de calor colocada debajo de ella. Use ∆x =5 y
∆y=5

Cuales son las condiciones

de Frontera?
 Ejercicio (EDP)

Se tiene una delgada placa de acero rectangular de 10 cm x

20 cm. Uno de los bordes de 10 cm se mantiene a 100⁰c y
los otros tres bordes a 0⁰c, estas temperaturas son producto
de una fuente de calor colocada debajo de ella. Use ∆x =5 y
∆y=5

Cuales son las condiciones

de Frontera?
Ejercicio (EDP)

𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

P1 (x2,y2)

P2(x3,y2)

P3(x4,y2)
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

P1 (x2,y2) 1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0

P2(x3,y2)

P3(x4,y2)
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

P1 (x2,y2) 1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0

1/25[0 +U3,2+ 0 + 0 -4U2,2]=0

P2(x3,y2)

P3(x4,y2)
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

P1 (x2,y2) 1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0

1/25[0 +U3,2+ 0 + 0 -4U2,2]=0

P2(x3,y2) 1/25[U2,2 +U4,2+U3,1+U3,3 -4U3,2]=0

1/25[U2,2 +U4,2+ 0 + 0 -4U3,2]=0
1/25[U3,2 +U5,2+U4,1+U4,3 -4U4,2]=0
P3(x4,y2)
1/25[U3,2 +U5,2+ 0 + 0 -4U4,2]=0
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

Reemplazando Wi=U(Pi)
1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0
P1 (x2,y2)
1/25[ 0 +W2 +0 +0 -4W1 ]=0

P2(x3,y2)

P3(x4,y2)
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

Reemplazando Wi=U(Pi)
1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0
P1 (x2,y2)
1/25[ 0 +W2 +0 +0 -4W1 ]=0
1/25[U2,2 +U4,2+U3,1+U3,3 -4U3,2]=0
P2(x3,y2)
1/25[ W1 + W3+ 0 + 0 -4W2 ]=0

P3(x4,y2) 1/25[U3,2 +U5,2+U4,1+U4,3 -4U4,2]=0

1/25[ W2 +100+ 0 + 0 -4W3 ]=0
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

Reemplazando Wi=U(Pi)
1/25[U1,2 +U3,2+U2,1+U2,3 -4U2,2]=0
P1 (x2,y2)
1/25[ 0 +W2 +0 +0 -4W1 ]=0
1/25[U2,2 +U4,2+U3,1+U3,3 -4U3,2]=0
P2(x3,y2)
1/25[ W1 + W3+ 0 + 0 -4W2 ]=0

P3(x4,y2) 1/25[U3,2 +U5,2+U4,1+U4,3 -4U4,2]=0

1/25[ W2 +100+ 0 + 0 -4W3 ]=0
Ejercicio (EDP)
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦

Sustituyendo ∆x =5 y ∆y=5

Describa el cambio en la temperatura en los puntos

P1 (x2,y2); P2(x3,y2); P3(x4,y2)

Reemplazando Wi=U(Pi)

Reordenando y escribiendo como un sistema de

ecuaciones
1/25[-4W1 + W2 + 0 + 0 + 0 ]=0
1/25[ W1 -4W2 + W3+ 0 + 0 ]=0
1/25[ 0 + W2 -4W3+ 0 +100]=0
Ejercicio (EDP)
Describa el cambio en la temperatura en los puntos
P1 (x2,y2); P2(x3,y2); P3(x4,y2)

Reemplazando Wi=U(Pi)

Reordenando y escribiendo como un sistema de

ecuaciones
1/25[-4W1 + W2 + 0 + 0 + 0 ]=0
1/25[ W1 -4W2 + W3+ 0 + 0 ]=0
1/25[ 0 + W2 -4W3+ 0 +100]=0

1/25[-4W1 + W2 + 0 ]=0
1/25[ W1 -4W2 + W3 ]=0
1/25[ 0 + W2 -4W3 ]=-100/25

−4/25 1/25 0 𝑤𝑤1 0

1/25 −4/25 1/25 𝑤𝑤2 = 0
0 1/25 −4/25 𝑤𝑤3 −4
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

eje de la tubería
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

y
𝜕𝜕 2 𝜕𝜕 2
𝑈𝑈(𝑥𝑥, 𝑦𝑦) + 2 𝑈𝑈(𝑥𝑥, 𝑦𝑦) = 0
𝜕𝜕𝑥𝑥 2 𝜕𝜕𝑦𝑦
x

eje de la tubería
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

30 30 30
CONDICIONES DE FRONTERA

30 30 30
20

10 15

0 0 0 0 0 0
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

30 30 30 Puntos/Nodos
NODOS a calcular
A CALCULAR
CONDICIONES DE FRONTERA

1 2 30 30 30
20

3 4 5 6
10 15

0 0 0 0 0 0
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

30 30 30 NODOS A CALCULAR
CONDICIONES DE FRONTERA

1 2 30 30 30
20

3 4 5 6
10 15

0 0 0 0 0 0

U0,1 + U2,1 + U1,0 + U1,2 -4U1,1 =0

Reemplazando Wi=U(Pi) 20 + W2 + 30 + W3 -4W1 =0
 Ejercicio (EDP)
Calcule la temperaturas en la siguiente tubería.
Use ∆x =1 y ∆y=1

W1

30 + W2 + W3 +20 - 4W1 = 0 W2
30 + 30 + W4 + W1 - 4W2 = 0 W3
W1 + W4 + 0 + 10 - 4W3 = 0
W4
W2 + W5 + 0 + W3 - 4W4 = 0
30 + W6 + 0 + W4 - 4W5 = 0 W5
30 + 15 + 0 + W5 - 4W6 = 0 W6
 Lecturas Complementarias

 Numerical Analysis Mathematics of Scientific Computing

Capítulo 9.1

 Computational Problems for Physics

Capítulo: 4.6 (Proyecto: Simulación de fluido incompresible utilizando
ecuaciones de Navier-Stokes)