Laplace
Laplace
Laplace
Ecuación de Laplace
Problema 7.1 Sea u(x, y) una función armónica para x, y ∈ (0, π) que es nula
en el borde de dicho cuadrado salvo u(x, 0) = 3 sin x, x ∈ (0, π).
Solución:
3 sinh(π − y)
u(x, y) = sin x. ✷
sinh π
Problema 7.2 Hallar en forma de serie la solución del problema para la ecua-
ción de Laplace ∆u(x, y) = uxx + uyy = 0, u(0, y) = 0, u(L, y) = 0, u(x, 0) =
f (x), uy (x, 0) = g(x) para x ∈ (0, L), y > 0. Aplicarlo al caso L = π, u(x, 0) =
0, uy (x, 0) = N −1 sin N x, donde N es un número natural. Tomar el lı́mite de la
solución cuando N tiende a infinito. Interpretar el resultado. Aplicarlo también
al caso L = π, u(x, 0) = 0, uy (x, 0) = x.
Solución:
163
ecuación,
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
X(0) = 0, X(L) = 0.
164
Por tanto, si f (x) tiene desarrollo en serie de senos en el intervalo [0, L],
∞ L
nπx 2 nπx
X Z
f (x) = fn sin , fn = f (x) sin dx,
n=1
L L 0 L
∞
X nπ nπx
g(x) = uy (x, 0) = Bn sin .
n=1
L L
Ası́ pues, si g(x) tiene desarrollo en serie de senos en el intervalo [0, L],
∞ L
nπx 2 nπx
X Z
g(x) = gn sin , gn = g(x) sin dx,
n=1
L L 0 L
∞ L
nπx 2 nπx
X Z
g(x) = gn sin , gn = g(x) sin dx. ✷
n=1
L L 0 L
∞
X ∞
X
f (x) = fn sin nx, g(x) = gn sin nx.
n=1 n=1
165
Si f (x) = 0, es obvio que todos los coeficientes fn son nulos. Si g(x) =
N −1 sin N x, la propia función es su desarrollo en serie, ası́ que gN = 1/N y el
resto de coeficientes son nulos. Por tanto, la solución en este caso es
sinh N y sin N x
u(x, y) = .✷
N2
Observamos que, si tomamos el lı́mite N tendiendo a infinito, el dato inicial
es trivial u(x, 0) = 0, uy (x, 0) = 0, ya que el seno está acotado por la unidad,
ası́ que g(x) decrece como 1/N . La solución en ese caso es la trivial u(x, y) ≡
0, como es fácil comprobar. Sin embargo, para y 6= 0, la solución u(x, y) =
N −2 sinh N y sin N x tiende a infinito como una exponencial eN y , en lugar de
tender a cero.
Por tanto, las condiciones iniciales conducen a malos problemas para la ecua-
ción de Laplace. Esto es razonable, ya que la ecuación de Laplace aparece en
problemas estacionarios y, por ello, no tiene mucho sentido considerar la coor-
denada y como un tiempo. ✷
En el caso f (x) = 0 todos los coeficientes fn son nulos. Y g(x) = x tiene
coeficientes
π
2 π 2(−1)n+1
2 sin nx x cos nx
Z
gn = x sin nx dx = − = ,
π 0 π n2 n 0 n
con lo cual la serie de la solución es
∞
X 2(−1)n+1
u(x, y) = sinh ny sin nx. ✷
n=1
n2
Solución:
166
Tenemos un problema de Dirichlet simplificado, en el que sólo hay condicio-
nes de contorno no triviales en uno de los lados del rectángulo,
u(x, M ) = f (x), u(x, 0) = 0, u(0, y) = 0 = u(L, y),
para x ∈ (0, L), y ∈ (0, M ).
Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,
X(0) = 0 = X(L), Y (0) = 0.
Resolvemos primero el problema para X(x),
X ′′ (x) + λX(x) = 0, X(0) = 0 = X(L).
Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores
y autofunciones,
nπ 2 nπx
λn = , Xn (x) = sin , n∈N
L L
Abordamos la otra ecuación,
nπ 2 nπy nπy
Yn′′ − Y = 0 ⇒ Yn (y) = An sinh + Bn cosh ,
L L L
para n ∈ N.
Imponemos la condición de contorno Yn (0) = 0,
0 = Yn (0) = Bn , n ∈ N,
y, por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,
∞
X nπy nπx
u(x, y) = An sinh sin .
n=0
L L
167
si la función f (x) admite desarrollo de Fourier en el intervalo (0, L) en serie de
senos”
∞
nπx 2 L nπx
X Z
f (x) = fn sin , fn = f (x) sin dx, n ∈ N.
n=1
L L 0 L
π
2
Z
fn = f (x) sin nx dx, n ∈ N.
π 0
2 π
2 h cos nx iπ 2
Z
0 n par
fn = sin nx dx = − = (1 − (−1)n ) =
π 0 π n 0 nπ 4/nπ n impar
∞
4 X sinh ny
u(x, y) = sin nx ✷.
π n impar
n sinh nπ
Solución:
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
168
uy
M
y
-uX uX
Ω
x L
-uy
X ′ (0) = 0 = X ′ (L), Y ′ (M ) = 0.
0 = Y0′ (M ) = B0 ⇒ Y0 (y) = A0 ,
nπ nπM/L nπ
0 = Yn′ (M ) = An e − Bn e−nπM/L ⇒ Bn = An e2nπM/L , n ∈ N.
L L
169
Yn (y) = An enπy/L + e−nπy/L e2nπM/L
nπ(M − y)
= An enπM/L enπ(M−y)/L + e−nπ(M−y)/L = an cosh
L
denotando a0 = A0 , an = 2An enπM/L , n ∈ N, y, por tanto, a falta de la
condición inhomogénea, la solución se puede expresar en forma de serie,
∞
X nπ(M − y) nπx
u(x, y) = an cosh cos .
n=0
L L
2 L
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, 0) = −g(x), uy (x, M ) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución
∞ nπ(M−y)
L X gn cosh L nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L
170
Problema 7.5 Hallar en forma de serie la solución del problema de Neu-
mann para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M )
con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ). Aplicar el resultado al caso en el que
L = π = M , g(x) = 2 sin2 3x. Lo mismo para g(x) = 2 sin2 3x − 1.
Solución:
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
171
y observamos que todos los coeficientes Bn tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,
∞
X nπy nπx
u(x, y) = An cosh cos .
n=0
L L
L
2
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución
∞
L X gn cosh nπ
L y nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L
172
la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g0 = 2,
g6 = −1. Por tanto, este dato de contorno no cumple la condición necesaria
g0 = 0, por lo que el problema de Neumann en este caso no tiene solución. ✷
En cambio, en el otro caso,
∞
g0 X
g(x) = 2 sin2 3x − 1 = − cos 6x = + gn cos nx,
2 n=1
la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g6 = −1.
Por tanto, la solución de problema de Neumann en este caso es
cosh 6y cos 6x
u(x, y) = K − .✷
6 sinh 6π
Problema 7.6 Sea u(x, y) una función armónica para x, y ∈ (0, π) cuya de-
rivada normal es nula en el borde de dicho cuadrado salvo uy (x, 0) = 5 cos 3x,
x ∈ (0, π).
Solución:
173
Solución:
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
0 = Y0 (0) = A0 ⇒ Y0 (y) = B0 y,
0 = Yn (0) = An
y observamos que todos los coeficientes An tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,
∞
X nπ nπx
u(x, y) = B0 y + Bn sinh y cos .
n=1
L L
174
Falta por imponer la condición inhomogénea u(x, M ) = f (x),
∞
X nπM nπx
f (x) = u(x, M ) = B0 M + Bn sinh cos ,
n=1
L L
L
2
Z
f (x) dx = f0 = 2B0 M,
L 0
∞ π
f0 X 2
Z
f (x) = + fn cos nx, fn = f (x) cos nx dx.
2 n=1
π 0
y sinh 4y cos 4x
u(x, y) = + .✷
2π 2 sinh 4π
175
En el otro caso, f (x) = cos(x/2), sı́ que nos vemos abocados a calcular el
desarrollo de Fourier,
2 π 1 π
Z
1 1
Z
fn = cos(x/2) cos nx dx = cos n + x + cos n − x dx
π 0 π 0 2 2
x=π
(−1)n (−1)n
1 sin (n + 1/2) x sin (n − 1/2) x 1
= + = −
π n + 1/2 n − 1/2 x=0 π n + 1/2 n − 1/2
1 (−1)n+1
= , n = 0, 1, . . .
π n2 − 1/4
Solución:
que en nuestro caso no es preciso realizar, ya que es patente que todos los
coeficientes bn son nulos, lo mismo que los an , salvo a0 = 1, a2 = 1/2. Por
tanto, la solución del problema es
1 r2
u(r, φ) = + cos 2φ. ✷
2 2R2
Problema 7.9 La función u(r, φ) es armónica en el interior del cı́rculo de
radio a centrado en el origen, siendo u(a, φ) = 5 + sin φ cos φ. Hallar u(r, ϕ) en
el interior del cı́rculo.
Solución:
176
propiedad. Sabemos que las funciones armónicas en el disco de radio a se pueden
expresar como
∞
a0 X r n
u(r, φ) = + (an cos nφ + bn sin nφ) ,
2 n=1
a
que en nuestro caso no es preciso realizar, ya que es patente que todos los coe-
ficientes son nulos, salvo a0 = 10, b2 = 1/2. Por tanto, la solución del problema
es
r2
u(r, φ) = 5 + 2 sin 2φ. ✷
2a
Problema 7.10 Hallar en forma de serie la solución del problema de Dirichlet
para la ecuación de Laplace ∆u(r, φ) = 0 en el semicı́rculo de radio R, y ≥ 0,
con condiciones de contorno u(R, φ) = f (φ), u(r, 0) = 0 = u(r, π), r ∈ (0, R),
φ ∈ (0, π). Aplicarlo al caso f (φ) = 6 sin 4x cos 4x.
Solución:
Φ(0) = 0, Φ(π) = 0,
177
Queda por resolver la ecuación en la coordenada radial,
que es una ecuación de Euler, que se reduce, como hemos visto, a una ecuación
lineal con coeficientes constantes con el cambio de variable independiente r =
es ⇔ s = ln r,
En nuestro caso particular f (φ) = 3 sin 8φ, con lo cual el desarrollo de Fourier
en senos sólo tiene no nulo el término f8 = 3 y, por tanto, la solución del
problema de Dirichlet queda
3r8
u(r, φ) = sin 8φ. ✷
R8
178
Problema 7.11 Sea u una función armónica en el disco de radio a centrado
en el origen y que se anula en el diámetro y = 0. Demostrar que u(x, −y) =
−u(x, y).
Solución:
serie nula de términos positivos, de donde se deduce que los coeficientes an son
todos nulos. Con la condición u(r, π) se llega a la misma conclusión.
Por tanto, la condición de anulación impone que el desarrollo sea sólo de
senos,
∞
a0 X r n
u(r, φ) = + bn sin nφ.
2 n=1
a
Problema 7.12 Sea u(x, y) una función armónica para x ∈ (0, π), y > 0. Se
tiene que u(x, 0) = 1 para x ∈ (0, π), u(0, y) = 0 = u(π, y) para y > 0. Hallar
u(x, y).
Solución:
179
y por la unicidad del desarrollo de Fourier, obtenemos un sistema infinito de
ecuaciones ordinarias lineales homogéneas,
Y ′′ − n2 Y = 0 ⇒ Yn (y) = an cosh ny + bn sinh ny, n ∈ N.
Por tanto, la función más general que es solución de la ecuación de Laplace
con dichas condiciones de contorno es
∞
X
u(x, y) = (an cosh ny + bn sinh ny) sin nx.
n=1
y los coeficientes an no son más que los del desarrollo de la función unidad en
serie de senos,
2 π 2 1 − (−1)n
2 h cos nx iπ
Z
0 n par
an = sin nx dx = − = = .
π 0 π n 0 π n 4/nπ n impar
Pero esto fija sólo los coeficientes an , no los bn . Luego el problema está mal
planteado, ya que es indeterminado, como ya sabı́amos, ya que es un proble-
ma mixto de valores iniciales para la ecuación de Laplace, que sólo permite
problemas de contorno.
El problema tendrı́a solución única si impusiéramos la condición adicional
de que u sea acotada en infinito. Dado que
an + bn ny an − bn −ny
an cosh ny + bn sinh ny = e + e ,
2 2
eliminar los términos divergentes en infinito supone eliminar las exponenciales
crecientes eny , lo que se consigue si bn = −an .
Por tanto, la función buscada tiene desarrollo en serie
∞
X 4 X e−ny
u(x, y) = an e−ny sin nx = sin nx
n=1
π n impar
n
4 X e−(2n−1)y
= sin(2n − 1)x. ✷
π n=1 2n − 1
180
por separación de variables, buscando soluciones de la forma u(r, φ) = R(r)Φ(φ),
que sustituidas en la ecuación,
que es una ecuación de Euler, que se reduce, como hemos visto, a una ecuación
lineal con coeficientes constantes con el cambio de variable independiente r =
es ⇔ s = ln r,
R̈0 = 0 ⇒ R0 (s) = A0 + B0 s = A0 + B0 ln r,
denotando por un punto la derivada con respecto a s.
Llegado este punto, observamos que aparecen funciones singulares, las poten-
cias inversas r−n y ln r. Imponiendo que las soluciones sean funciones acotadas
y regulares en el cı́rculo, es patente que todos los coeficientes Bn tienen que ser
nulos. Por tanto, las funciones radiales son Rn (r) = rn , n = 0, 1, . . .
Y podemos escribir la solución de la ecuación de Laplace en el cı́rculo como
∞
X
u(r, φ) = A0 + rn An cos nφ + Ãn sin nφ .
n=1
181
y si la función g(φ) admite desarrollo de Fourier en serie de senos y cosenos,
∞
a0 X
g(φ) = + (an cos nφ + bn sin nφ) ,
2 n=1
2π 2π
1 1
Z Z
an = g(φ) cos nφ dφ, bn = g(φ) sin nφ dφ,
π 0 π 0
2π 2π
1 1
Z Z
an = g(φ) cos nφ dφ, bn = g(φ) sin nφ dφ.
π 0 π 0
Esta solución en serie de Fourier se puede sumar, al igual que se hizo para
el problema de Dirichlet,
∞
X r0n
u(r0 , φ0 ) − K = (an cos nφ0 + bn sin nφ0 )
n=1
nRn−1
∞
1 2π r0n
Z X
= g(φ) (cos nφ0 cos φ + sin nφ0 sin φ) dφ
π 0 n=1
nRn−1
∞
1 2π r0n
Z X
= g(φ) cos n(φ − φ0 ) dφ,
π 0 n=1
nRn−1
ya que las series de Fourier se pueden integrar bajo condiciones muy generales.
182
El integrando es la integral de una serie geométrica, como se comprueba
desarrollando el coseno en exponenciales imaginarias,
∞ ∞ ∞
X r0n X r0n inγ
X rn
0
2 cos nγ = e + e−inγ
n=1
nRn−1 n=1
nR n−1
n=1
nR n−1
iγ n−1
Z r0 ∞ ∞ n−1 !
reiγ
iγ
X re −iγ
X
= e +e dr
0 n=1
R n=1
R
Z r0
Reiγ Re−iγ
= + dr
0 R − reiγ R − re−iγ
r0
−R ln R − reiγ + ln R − re−iγ 0
=
2R ln R − R ln R2 + r02 − 2Rr0 cos γ ,
=
sin 2φ
g(φ) = ,
2
tiene todos los coeficientes nulos, salvo b2 = 1/2. Por tanto, la solución del
problema es
r2
u(r, φ) = K + sin 2φ. ✷
4R
Problema 7.14 Hallar en forma de serie la solución del problema de Neu-
mann para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M )
con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ). Aplicar el resultado al caso en el que
L = π = M , g(x) = 2 sin2 3x. Lo mismo para g(x) = 2 sin2 3x − 1.
Solución:
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
183
nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones
ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.
Tenemos un problema de Neumann simplificado, en el que sólo hay condi-
ciones de contorno no triviales en uno de los lados del rectángulo,
184
L
2
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución
∞ nπy
L X gn cosh L nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L
la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g6 = −1.
Por tanto, la solución de problema de Neumann en este caso es
cosh 6y cos 6x
u(x, y) = K − .✷
6 sinh 6π
Lo comprobamos,
cosh 6y sin 6x 6 cosh 6y cos 6x
ux (x, y) = , uxx (x, y) = ,
sinh 6π sinh 6π
185
sinh 6y cos 6x 6 cosh 6y cos 6x
uy (x, y) = − , uyy (x, y) = − .
sinh 6π sinh 6π
Obviamente, ux (0, y) = 0 = ux (π, y), porque el seno se anula en cero y π.
Igualmente uy (x, 0) = 0, porque se anula el seno hiperbólico. A su vez,
Solución:
u(0, φ) = 0, u(L, φ) = f (φ), u(z, −π) = u(z, π), uφ (z, −π) = uφ (z, π),
kn = n2 , n = 0, 1, . . . ,
186
para n ∈ N.
Imponemos la condición de contorno Zn (0) = 0,
0 = Z0 (0) = A0 ⇒ Z0 (z) = B0 z,
1 π 1 π
Z Z
fn = f (φ) cos nφ dφ, f˜n = f (φ) sin nφ dφ,
π −π π −π
comparando ambos desarrollos, concluimos
f0 fn f˜n
B0 = , Bn = , B̃n = , n = 1, 2, . . .
2L sinh nL sinh nL
“El problema de la ecuación de Laplace uzz + uφφ = 0 para z ∈ (0, L),
φ ∈ (−π, π) con condiciones de contorno u(0, φ) = 0, u(L, φ) = f (φ), u(z, −π) =
u(z, π), uφ (z, −π) = uφ (z, π) tiene solución
∞
f0 sinh nz
fn cos nφ + f˜n sin nφ
X
u(z, φ) = z+ ,
2L n=1
sinh nL
187
observamos que este problema es un problema de Dirichlet para un cilindro de
radio unidad y altura L con datos de contorno sobre las circunferencias situadas
en los planos z = 0 y z = L.
Usando trigonometrı́a elemental,
1 1
f (φ) = 4 sin 2φ cos 2φ + cos2 2φ = 2 sin 4φ + + cos 4φ,
2 2
observamos que la función f (φ) tiene desarrollo de Fourier con coeficientes nulos
salvo f0 = 1, f4 = 1/2, f˜4 = 2. Por tanto, la solución de problema de Dirichlet
en este caso es
1 1 1
u(z, φ) = z+ 2 sin 4φ + cos 4φ sinh 4z. ✷
2L sinh 4L 2
Problema 7.16 Hallar en forma de serie la solución del problema para la ecua-
ción de Laplace R2 uzz + uφφ = 0, uz (0, φ) = 0, uz (h, φ) = g(φ), u(z, −π) =
u(z, π), uφ (z, −π) = uφ (z, π) para z ∈ (0, h), φ ∈ (−π, π), siendo constantes
h, R > 0. Aplicarlo al caso h = 1 = R, g(φ) = φ. Aplicarlo al caso h = 1 = R,
g(φ) = φ2 .
Solución:
188
Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,
kn = n2 , n = 0, 1, . . . ,
n2 nz nz
Zn′′ − Zn = 0 ⇒ Z0 (z) = A0 + B0 z, Zn (z) = An cosh + Bn sinh ,
R2 R R
para n = 1, 2, . . .
Necesitamos las derivadas,
nAn nz nBn nz
Z0′ (z) = B0 , Zn′ (z) = sinh + cosh ,
R R R R
para imponer la condición de contorno Zn′ (0) = 0,
nBn nz
0 = Zn′ (0) = ⇒ Zn (z) = An cosh ,
R R
y observamos que todos los coeficientes Bn tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, haciendo uso del principio
de superposición lineal, la solución se puede expresar en forma de serie,
∞
X nz
u(z, φ) = A0 + An cos nφ + Ãn sin nφ cosh ,
n=1
R
que es un desarrollo en serie de senos y cosenos para g(φ), sin término constante,
lo que impone una condición sobre el dato de contorno, como sucede en los
problemas de Neumann.
189
Por tanto, si la función g(φ) admite desarrollo de Fourier en serie de senos
y cosenos,
∞ ∞
g0 X X
g(φ) = + gn cos nφ + g̃n sin nφ,
2 n=1 n=1
1 π 1 π
Z Z
gn = f (φ) cos nφ dφ, g̃n = g(φ) sin nφ dφ,
π −π π −π
comparando ambos desarrollos, concluimos
Rgn Rg̃n
g0 = 0, An = , Ãn = , n = 1, 2, . . .
n sinh nh
R n sinh nh
R
190
Problema 7.17 Sea u(x, y) una función tal que uxx + uyy = −2 sin x sin y para
x, y ∈ (0, π) que satisface u = 0 en el borde de dicho cuadrado. Calcular u(x, y).
Solución:
191
Problema 7.18 Sea u una función armónica y no negativa en el cı́rculo de
radio R centrado en el origen, y continua en su borde. Demostrar las desigual-
dades
R − ||x|| R + ||x||
u(0) ≤ u(x) ≤ u(0).
R + ||x|| R − ||x||
Utilizar la fórmula integral de Poisson.
Solución:
Para que se anule Gx0 (x, y, 0) habrá que añadir a ψx0 una función armónica
Ux0 cuya contribución sobre el eje X cancele la de la solución fundamental. Por
tanto, es preciso que
1
Ux0 (x, y, 0) = p .
4π (x − x0 ) + (y − y0 )2 + z02
2
192
La manera de conseguirlo es situar una carga negativa en el punto simétrico
de x0 respecto al plano XY , que es x̃0 = (x0 , y0 , −z0 ). Ası́, como los puntos
del plano XY son equidistantes de x0 y de x̃0 , las contribuciones de ambos se
cancelarán, al tener signo opuesto,
1 1
Ux0 (x, y, z) = = p ,
4πkx − x̃0 k 4π (x − x0 ) + (y − y0 )2 + (z + z0 )2
2
que es una función armónica, ya que el polo (x0 , y0 , −z0 ) no pertenece al semi-
espacio z > 0.
Por tanto, la función de Green del semiespacio z > 0 es
1 1
Gx0 (x, y, z) = − .✷
4πkx − x̃0 k 4πkx − x0 k
z0 ∞
Z ∞
f (x, y)
Z
u(x0 , y0 ) = dx dy 3/2
.
2π −∞ −∞ ((x − x0 ) + (y − y0 )2 + z 2 )
2
193