2.

- Resolver la ecuación de Laplace para la placa cuadrada de 100 cm de lado mostrada en

la figura, donde los lados izquierdo y derecho se mantienen a 0 y los lados inferior y superior
están dados por las funciones mostradas

Para resolver la ecuación de Laplace en una placa cuadrada con condiciones de contorno dadas,
podemos utilizar el método de separación de variables. En este caso, la ecuación de Laplace es:

𝛻²𝑢 = 0
Donde u es la función de temperatura en la placa y ∇² es el operador laplaciano.

Dado que los lados izquierdo y derecho de la placa se mantienen a 0, podemos establecer las
siguientes condiciones de contorno:

𝑢(0, 𝑦) = 0
𝑢(100, 𝑦) = 0
Además, los lados superior e inferior de la placa están dados por las funciones:

𝑢(𝑥, 100) = 𝑥²(100 − 𝑥)

𝑢(𝑥, 0) = 𝑥(100 − 𝑥)
Vamos a buscar una solución de la forma u(x, y) = X(x)Y(y) y sustituimos en la ecuación de Laplace:

𝑋′′(𝑥)𝑌(𝑦) + 𝑋(𝑥)𝑌′′(𝑦) = 0
Dividimos ambos lados por X(x)Y(y) para obtener:

(𝑋′′(𝑥) / 𝑋(𝑥)) + (𝑌′′(𝑦) / 𝑌(𝑦) ) = 0

De aquí

𝑋′′(𝑥) / 𝑋(𝑥) = −𝜆²

𝑌 ′′ (𝑦)/𝑌(𝑦) = 𝜆²

Donde λ es una constante de separación. Resolviendo esta ecuación diferencial, obtenemos las
soluciones para X(x):
 𝑋(𝑥) = 𝐴 𝑠𝑖𝑛(𝜆𝑥) + 𝐵 𝑐𝑜𝑠(𝜆𝑥)

Usando las condiciones de contorno X(0) = 0 y X(100) = 0, encontramos que B = 0 y λ = nπ/100,

donde n es un número entero distinto de cero.

Entonces, nuestras soluciones para X(x) son:

𝑋ₙ(𝑥) = 𝐴ₙ 𝑠𝑖𝑛(𝑛𝜋𝑥/100)
Ahora, consideramos la ecuación diferencial para Y(y):

𝑌’’(𝑦) / 𝑌(𝑦) = 𝜆²
𝑛𝜋
Donde utilizando los valores de 𝜆 = 100 se tiene
𝑛𝜋 𝑛𝜋
𝑌(𝑦) = 𝐶𝑒 100𝑦 + 𝐷𝑒 −100𝑦
Aplicando la condición del problema I
𝑛𝜋 𝑛𝜋
𝑌(0) = 𝐶𝑒 1000 + 𝐷𝑒 −1000 = 0
De donde teniendo la relación de C y D obtenemos

𝐷 = −𝐶
Por lo cual, factorizando C, obtenemos:
𝑛𝜋 𝑛𝜋 𝑛𝜋
𝑌(𝑦) = 𝐶 [𝑒 100𝑦 − 𝑒 −100𝑦 ] = 𝐶𝑠𝑒𝑛ℎ ( 𝑦)
100
De tal forma que la solución para el problema I queda, utilizando el principio de superposición,
como
∞
𝑛𝜋
𝑢𝐼 (𝑥, 𝑦) = ∑ 𝑏𝑛 𝑠𝑒𝑛ℎ ( 𝑦) 𝑠𝑖𝑛(𝑛𝜋𝑥/100)
100
𝑛=1

Donde, aplicando la última condición de frontera

𝑢(𝑥, 100) = 𝑥²(100 − 𝑥)

Se tiene
∞
𝑛𝜋
𝑢𝐼 (𝑥, 100) = 𝑥²(100 − 𝑥) = ∑ 𝑏𝑛 𝑠𝑒𝑛ℎ ( 100) 𝑠𝑖𝑛(𝑛𝜋𝑥/100)
100
𝑛=1

La cual es la serie de Fourier seno de la función

𝑥²(100 − 𝑥)
En donde
𝑛𝜋
𝑏𝑛 𝑠𝑒𝑛ℎ (100 100) son los coeficientes de la serie

Por lo cual estos coeficientes vienen dados por

100
2 𝑛𝜋
𝑏𝑛 𝑠𝑒𝑛ℎ(𝑛𝜋) = ∫ 𝑥²(100 − 𝑥)𝑠𝑖𝑛( 𝑥)𝑑𝑥
100 0 100
De donde:
100
2 𝑛𝜋
𝑏𝑛 = ∫ 𝑥²(100 − 𝑥)𝑠𝑖𝑛( 𝑥)𝑑𝑥
100𝑠𝑒𝑛ℎ(𝑛𝜋) 0 100

Resolvemos la integral separando términos:

100 100 100
2 𝑛𝜋 𝑛𝜋 𝑛𝜋
= [∫ 𝑥 2 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 + 100 ∫ 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 − ∫ 𝑥𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 ]
100𝑠𝑒𝑛ℎ(𝑛𝜋) 0 100 0 100 0 100

I II III

I
Resolvemos la primer integral por bloques:

(Polinomio) Derivamos Integramos (función trigonométrica)

𝑛𝜋
𝑥2 (+) 𝑠𝑒𝑛 (100 𝑥)

100 𝑛𝜋
2𝑥 (−) −( 𝑛𝜋 )𝑐𝑜𝑠 (100 𝑥)

100 𝑛𝜋
2 (+) −( 𝑛𝜋 )2 𝑠𝑒𝑛 (100 𝑥)
100 𝑛𝜋
0 ( 𝑛𝜋 )3 𝑐𝑜𝑠 (100 𝑥)

Dando como resultado de la primer integral y evaluando:

100 𝑛𝜋 100 2 𝑛𝜋 100 3 𝑛𝜋
100
I=−𝑥 2 ( 𝑛𝜋 ) cos (100 𝑥) + 2𝑥 ( 𝑛𝜋 ) 𝑠𝑒𝑛 (100 𝑥) + 2 ( 𝑛𝜋 ) cos (100 𝑥)
0
100 100 2 100 3
=−1002 ( 𝑛𝜋 ) cos(𝑛𝜋) + 2(100) ( 𝑛𝜋 ) 𝑠𝑒𝑛 (𝑛𝜋) + 2 ( 𝑛𝜋 ) cos(𝑛𝜋)

100 100 2 100 3

=−10000 ( 𝑛𝜋 ) cos(𝑛𝜋) + 200 ( 𝑛𝜋 ) 𝑠𝑒𝑛 (𝑛𝜋) + 2 ( 𝑛𝜋 ) cos(𝑛𝜋)

cos(𝑛𝜋) =(-1)n

II.

Resolvemos la segunda integral por separación de variables:

100
𝑛𝜋 100 𝑛𝜋
∫ 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 = − ( ) cos ( 𝑥)
0 100 𝑛𝜋 100
 𝑛𝜋
𝑢= 𝑥
100
𝑛𝜋
𝑑𝑢 = 𝑑𝑥
100
100
𝑑𝑥 = 𝑑𝑢
𝑛𝜋

Dando como resultado de la segunda integral y evaluando:

100
100 𝑛𝜋
II=− ( 𝑛𝜋 ) cos (100 𝑥)
0
100
=− ( 𝑛𝜋 ) cos(𝑛𝜋)

100
=− ( 𝑛𝜋 ) cos(𝑛𝜋)

cos(𝑛𝜋) =(-1)n

III.

Resolvemos la segunda integral de igual forma por bloques:

(Polinomio) Derivamos Integramos (función trigonométrica)

𝑛𝜋
𝑥 (+) 𝑠𝑒𝑛 (100 𝑥)

100 𝑛𝜋
1 (−) −( 𝑛𝜋 )𝑐𝑜𝑠 (100 𝑥)

100 𝑛𝜋
0 −( 𝑛𝜋 )2 𝑠𝑒𝑛 (100 𝑥)

Dando como resultado se la segunda integral y evaluando:

100 𝑛𝜋 100 2 𝑛𝜋 100

= − [− 𝑥𝑐𝑜𝑠 ( 𝑥) + ( ) 𝑠𝑖𝑛 ( 𝑥)]
𝑛𝜋 100 𝑛𝜋 100 0

100 100 2
= (100)𝑐𝑜𝑠(𝑛𝜋) − ( ) 𝑠𝑖𝑛(𝑛𝜋)
𝑛𝜋 𝑛𝜋
Simplificando términos de cada integral obtenemos:

I.
2 1,000,000 2,000,000
= [( ) (−1)𝑛 + ( ) (−1)𝑛 ]
100𝑠𝑒𝑛ℎ(𝑛𝜋) 𝑛𝜋 𝑛3 𝜋 3
 2 1,000,000 2,000,000
= {[( )+( )] (−1)𝑛 }
100𝑠𝑒𝑛ℎ(𝑛𝜋) 𝑛𝜋 𝑛3 𝜋 3

II.
2 100
=
100𝑠𝑒𝑛ℎ(𝑛𝜋)
[− (
𝑛𝜋
) (−1)𝑛 ]

III.
2 10,000
= {[ ] (−1)𝑛 }
100𝑠𝑒𝑛ℎ(𝑛𝜋) 𝑛𝜋
Juntando los tres términos resultantes de cada integral obtenemos:

2 1,000,000 100 10,000 2,000,000

[( )−( )+( )+( )] (−1)𝑛
100𝑠𝑒𝑛ℎ(𝑛𝜋) 𝑛𝜋 𝑛𝜋 𝑛𝜋 𝑛3 𝜋 3
1 1009900 2,000,000
= [( )+( )] (−1)𝑛
50𝑠𝑒𝑛ℎ 𝑛𝜋
( ) 𝑛𝜋 𝑛3 𝜋3

Entonces como solución obtenemos la siguiente:

1 1009900 2,000,000 𝑛𝜋 𝑛𝜋
𝑢(𝑥, 𝑦) = ∑∞
𝑛=1 [50𝑠𝑒𝑛ℎ(𝑛𝜋) [( )+( )] (−1)𝑛 ] 𝑠𝑒𝑛 (100 𝑥) 𝑠𝑒𝑛ℎ (100 𝑦)
𝑛𝜋 𝑛3 𝜋 3

Para el caso II al momento de calcular g(x) donde ocupamos las siguientes condiciones de frontera:

𝑢(x,100) = 𝑥(100 − 𝑥)

𝑢(0,y) = 0

𝑢(100,y) = 0
𝑢(x,b) = 0
Para resolver el problema utilizamos separación de variables:
Vamos a buscar una solución de la forma u(x, y) = X(x)Y(y) y sustituimos en la ecuación de Laplace:

𝑋′′(𝑥)𝑌(𝑦) + 𝑋(𝑥)𝑌′′(𝑦) = 0

Dividimos ambos lados por X(x)Y(y) para obtener:

(𝑋′′(𝑥) / 𝑋(𝑥)) + (𝑌′′(𝑦) / 𝑌(𝑦) ) = 0

𝑋 ′′ (𝑥) 𝑌 ′′ (𝑥) 𝜆=0
( )=− = −𝜆²
𝑋(𝑥) 𝑌(𝑥)

Se obtienen dos casos donde:

 i) 𝑋 ′′ (𝑥) + 𝜆²𝑋(𝑥)) = 0 𝑚2 + 𝜆² = 0
ii) 𝑌 ′′ (𝑦) − 𝜆²𝑌(𝑦)) = 0 𝑚2 − 𝜆² = 0

Entonces, nuestras soluciones para X(x) son:

𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛 (𝜆𝑥)

Ahora, consideramos la ecuación diferencial para Y(y):

𝑌(𝑦) = 𝐶3𝑒 𝜆𝑦 + 𝐶4𝑒 −𝜆𝑦

Aplicando condiciones de frontera:

𝑢(x,b) = 0
𝑌(𝑏) = 0 = 𝐶3𝑒 𝜆𝑦 + 𝐶4𝑒 −𝜆𝑦
𝑒 𝜆𝑦
𝐶4 = −𝐶3 = −𝐶3𝑒 2𝜆𝑦
𝑒 −𝜆𝑦
Sustituyendo en Y(y) tenemos que:

𝑌(𝑦) = 𝐶3[𝑒 𝜆𝑦 − 𝑒 2𝜆𝑦 𝑒 −𝜆𝑦 ]

Usando las condiciones de contorno X(0) = 0 y X(100) = 0:

Entonces, nuestras soluciones para X(0) son:

𝑋(0) = 𝐶2 𝑠𝑖𝑛(𝜆𝑥)
Entonces, nuestras soluciones para X(100) son:

𝑋(100) = 𝐶2 𝑠𝑖𝑛(𝜆100) = 0
𝑠𝑖𝑛(𝜆100) = 0
𝜆100 = 𝑛𝜋
𝑛𝜋
𝜆 = 100

De tal forma que la solución para el problema II queda, utilizando el principio de superposición,
como
∞
𝑛𝜋 𝑛𝜋 𝑛𝜋 𝑛𝜋
𝑢(𝑥, 𝑦) = ∑ 𝑏𝑛 [𝑒 100𝑦 − 𝑒 2100𝑏 𝑒 −100𝑦 ] 𝑠𝑖𝑛( 𝑥)
100
𝑛=1

Donde, aplicando la última condición de frontera

𝑢(𝑥, 100) = 𝑥(100 − 𝑥)

Se tiene
∞
𝑛𝜋 𝑛𝜋
𝑢𝐼𝐼 (𝑥, 100) = 𝑥(100 − 𝑥) = ∑ 𝑏𝑛 [1 − 𝑒 2100𝑏 ] 𝑠𝑖𝑛( 𝑥)
100
𝑛=1

Despejando bn:
 100
𝑛𝜋 2 𝑛𝜋
𝐵𝑛 = 𝑏𝑛 [1 − 𝑒 2100𝑏 ] = ∫ 𝑥(100 − 𝑥)𝑠𝑖𝑛( 𝑥)𝑑𝑥
100 0 100
100
2 𝑛𝜋
𝑏𝑛 = 𝑛𝜋 ∫ 𝑥(100 − 𝑥)𝑠𝑖𝑛( 𝑥)𝑑𝑥
100 [1 − 𝑒 2100𝑏 ] 0 100

Resolvemos la integral separando términos:

100 100 100
2 𝑛𝜋 𝑛𝜋 𝑛𝜋
= 𝑛𝜋 [∫ 𝑥𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 + 100 ∫ 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 − ∫ 𝑥𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 ]
2 𝑏 100 100 100
100 [1 − 𝑒 100 ] 0 0 0

I II III

I.

Resolvemos la primer integral por bloques:

(Polinomio) Derivamos Integramos (función trigonométrica)

𝑛𝜋
𝑥 (+) 𝑠𝑒𝑛 ( 𝑥)
100

100 𝑛𝜋
1 (−) −( 𝑛𝜋 )𝑐𝑜𝑠 (100 𝑥)

100 𝑛𝜋
0 −( 𝑛𝜋 )2 𝑠𝑒𝑛 (100 𝑥)

Dando como resultado se la segunda integral y evaluando:

100 𝑛𝜋 100 2 𝑛𝜋 100

= − [− 𝑥𝑐𝑜𝑠 ( 𝑥) + ( ) 𝑠𝑖𝑛 ( 𝑥)]
𝑛𝜋 100 𝑛𝜋 100 0

100 100 2
= (100)𝑐𝑜𝑠(𝑛𝜋) − ( ) 𝑠𝑖𝑛(𝑛𝜋)
𝑛𝜋 𝑛𝜋

II.

Resolvemos la segunda integral por separación de variables:

100
𝑛𝜋 100 𝑛𝜋
∫ 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 = − ( ) cos ( 𝑥)
0 100 𝑛𝜋 100
𝑛𝜋
𝑢= 𝑥
100
𝑛𝜋
𝑑𝑢 = 𝑑𝑥
100
 100
𝑑𝑥 = 𝑑𝑢
𝑛𝜋

Dando como resultado de la segunda integral y evaluando:

100
100 𝑛𝜋
II=− ( 𝑛𝜋 ) cos (100 𝑥)
0
100
=− ( 𝑛𝜋 ) cos(𝑛𝜋)

100
=− ( 𝑛𝜋 ) cos(𝑛𝜋)

cos(𝑛𝜋) =(-1)n

III.

Resolvemos la segunda integral de igual forma por bloques:

(Polinomio) Derivamos Integramos (función trigonométrica)

𝑛𝜋
𝑥 (+) 𝑠𝑒𝑛 ( 𝑥)
100

100 𝑛𝜋
1 (−) −( 𝑛𝜋 )𝑐𝑜𝑠 (100 𝑥)

100 𝑛𝜋
0 −( 𝑛𝜋 )2 𝑠𝑒𝑛 (100 𝑥)

Dando como resultado se la segunda integral y evaluando:

100 𝑛𝜋 100 2 𝑛𝜋 100

= − [− 𝑥𝑐𝑜𝑠 ( 𝑥) + ( ) 𝑠𝑖𝑛 ( 𝑥)]
𝑛𝜋 100 𝑛𝜋 100 0

100 100 2
= (100)𝑐𝑜𝑠(𝑛𝜋) −( ) 𝑠𝑖𝑛(𝑛𝜋)
𝑛𝜋 𝑛𝜋
Simplificando términos de cada integral obtenemos:

I.

2 10,000
= 𝑛𝜋 {[− ] (−1)𝑛 }
2 𝑏
100[1−𝑒 100 ] 𝑛𝜋
II.
2 100 𝑛
= 2
𝑛𝜋
𝑏
[− ( 𝑛𝜋 ) (−1) ]
100[1−𝑒 100 ]

III.
2 10,000
= 𝑛𝜋 {[ ] (−1)𝑛 }
100 [1 − 𝑒 2100𝑏 ] 𝑛𝜋

Juntando los tres términos resultantes de cada integral obtenemos:

2 100
𝑛𝜋 {[( − )] (−1)𝑛 }
100 [1 − 𝑒 2100𝑏 ] 𝑛𝜋

1 100
= 𝑛𝜋 {[( − )] (−1)𝑛 }
50 [1 − 𝑒 2100𝑏 ] 𝑛𝜋

Entonces como solución obtenemos la siguiente:

𝑛𝜋 𝑛𝜋 𝑛𝜋
1 100 𝑛𝜋
𝑢(𝑥, 𝑦) = ∑∞
𝑛=1 [ 2
𝑛𝜋
𝑏
{[− ( 𝑛𝜋 )] (−1)𝑛 }] 𝑠𝑒𝑛 (100 𝑥) [𝑒 100𝑦 − 𝑒 2100𝑦 𝑒 −100𝑦 ]
50[1−𝑒 100 ]

Por lo tanto, nuestra solución final sería la sumatoria de los dos casos dando como resultado del
problema:
1 1009900 2,000,000 𝑛𝜋 𝑛𝜋
𝑢(𝑥, 𝑦) = ∑∞
𝑛=1 [50𝑠𝑒𝑛ℎ(𝑛𝜋) [( )+( )] (−1)𝑛 ] 𝑠𝑒𝑛 (100 𝑥) 𝑠𝑒𝑛ℎ (100 𝑦) +
𝑛𝜋 𝑛3 𝜋 3
𝑛𝜋 𝑛𝜋 𝑛𝜋
1 100 𝑛𝜋
[ 2
𝑛𝜋
𝑏
{[− ( 𝑛𝜋 )] (−1)𝑛 }] 𝑠𝑒𝑛 (100 𝑥) [𝑒 100𝑦 − 𝑒 2100𝑏 𝑒 −100𝑦 ]
50[1−𝑒 100 ]

Gráficas