Derivation of Lapce
Derivation of Lapce
Derivation of Lapce
y = r sin() sin(),
z = r cos(),
(1)
(2)
2f
2f
2f
+
+
.
x2
y 2
z 2
(3)
=
=
=
f
r
f
r
f
r
r
f
f
+
+
,
x
x x
r
f
f
+
+
,
y
y y
r
f f
+
+
.
z
z
z
(4)
(5)
(6)
The next step is to convert the right-hand side of each of the above three equations so that it only has partial derivatives in terms of r, and . We can do
hDerivationOfTheLaplacianFromRectangularToSphericalCoordinatesi
created: h201303-21i by: hswapnizzlei version: h39376i Privacy setting: h1i hTopici h53A45i
This text is available under the Creative Commons Attribution/Share-Alike License 3.0.
You can reuse this document or portions thereof only if you do so under terms that are
compatible with the CC-BY-SA license.
1 Readers should note that we do not have to define the Laplacian this way. A more rigorous
approach would be to define the Laplacian in some coordinate free manner.
this by substituting the following values (which are easily derived from (2)) in
their respective places in the above three equations
r
1
= sin() cos(),
= cos() cos(),
x
x
r
r
1
= sin() sin(),
= cos() sin(),
y
y
r
r
= cos(),
= sin(),
= 0.
z
z
r
z
1 sin()
=
,
x
r sin()
1 cos()
=
,
y
r sin()
(7)
sin() cos()
f
1
f
1 sin() f
+ cos() cos()
.
r
r
r sin()
sin() cos()
(8)
f
x
h f i 1
h f i 1 sin() h f i
+ cos() cos()
(9)
.
r x
r
x
r sin() x
Now the trick is to substitute equation (8) into equation (9) in order to eliminate
any partial derivatives with respect to x. The result is the following equation
2f
x2
1
f
1 sin() f i
h
f
+ cos() cos()
sin() cos()
+
r
r
r
r sin()
h
1
f
1 sin() f i
1
f
cos() cos()
+ cos() cos()
sin() cos()
r
r
r sin()
1 sin() h
f
1
f
1 sin() f i
sin() cos()
+ cos() cos()
.
r sin()
r
r
r sin()
sin() cos()
In the hopes of simplifying the above equation, we operate the derivates on the
operands and get
2f
x2
h
2f
1
f
1
2f
sin() cos() sin() cos() 2 2 cos() cos()
+ cos() cos()
+
r
r
r
r
h
1 sin() f
1 sin() 2 f i 1
f
+
cos()
cos()
cos() cos()
+
r2 sin() r sin() r
r
r
2f
1
f
1
2f
1 sin() cos() f
sin() cos()
sin() cos()
+ cos() cos() 2 +
r
r
r sin2 ()
1 sin() 2 f i 1 sin() h
f
2f
sin() sin()
+ sin() cos()
r sin()
r sin()
r
r
1
f
1
2f
1 cos() f
1 sin() 2 f i
cos() sin()
+ cos() cos()
.
r
r
r sin() r sin() 2
After further simplifying the above equation, we arrive at the following form
2f
x2
2f
1
f
1
2f
2 cos() sin()cos2 ()
+ cos() sin()cos2 ()
+
2
r
r
r
r
f
1
2f i h 1 2
1
f
sin()
cos()
sin()
cos()
+
+ cos ()cos2 ()
r2
r
r
r
r
1
2f
1
f
1
2f
sin() cos()cos2 ()
2 sin() cos()cos2 ()
+ 2 cos2 ()cos2 () 2 +
r
r r
1 sin() cos()cos2 () f
1 cos() sin() cos() 2 f i h 1 2
f
+ sin ()
r2
r2
sin()
r
r
sin2 ()
sin2 ()cos2 ()
1
2f
1 cos()sin2 () f
1 cos() cos() sin() 2 f
sin() cos()
+ 2
2
+
r
r r
sin()
r
sin()
1 sin() cos() f
1 sin2 () 2 f i
+ 2
.
2
2
r
r sin2 () 2
sin ()
Notice that we have derived the first term of the right-hand side of equation
2
(3) (i.e. xf2 ) in terms of spherical coordinates. We now have to do a similar
2
arduous derivation for the rest of the two terms (i.e. yf2 and zf2 ). Lets do it!
After we substitute the values of (7) into equation (5) we get
f
y
sin() sin()
f
1
f
1 cos() f
+ cos() sin()
+
.
r
r
r sin()
(10)
sin() sin()
h f i 1
h f i 1 cos() h f i
+ cos() sin()
+
(11)
.
r y
r
y
r sin() y
Now we substitute equation (10) into equation (11) in order to eliminate any
partial derivatives with respect to y. The result is the following
2f
y 2
h
f
1
f
1 cos() f i
sin() sin()
+ cos() sin()
+
+
r
r
r
r sin()
h
f
1
f
1 cos() f i
1
cos() sin()
sin() sin()
+ cos() sin()
+
+
r
r
r
r sin()
1 cos() h
f
1
f
1 cos() f i
sin() sin()
+ cos() sin()
+
.
r sin()
r
r
r sin()
sin() sin()
2
y
r
r
r
r
h
1 cos() f
1 cos() 2 f i 1
f
+ cos() sin() sin() sin()
+
+
2
r sin() r sin() r
r
r
1 cos()cos() f
f
1
f
1
2f
cos() sin()
sin() sin()
+ cos() sin() 2
+
r
r
r sin2 ()
f
1 cos() 2 f i 1 cos() h
2f
+
sin() cos()
+ sin() sin()
+
r sin()
r sin()
r
r
1
f
1
2f
1 sin() f
1 cos() 2 f i
,
cos() cos()
+ cos() sin()
+
r
r
r sin() r sin() 2
and after some simplifications
h
2f
1
f
2f
1
2f
=
sin2 ()sin2 () 2 2 sin() cos()sin2 ()
+ sin() cos()sin2 ()
2
y
r
r
r
r
1
f
1
2f i h 1
f
sin()
cos()
+
sin()
cos()
+
cos() sin()sin2 ()
+
2
r
r
r
r
r
1 2
f
1
f
1
2f
cos ()sin2 ()
2 sin() cos()sin2 ()
+ 2 cos2 ()sin2 () 2
r
r
r
r
r
sin2 ()
1
2f
1 cos()cos2 () f
1 cos()cos() sin() 2 f
cos() sin()
+ 2
+ 2
r
r r
sin()
r
sin()
1 sin() cos() f
1 cos2 () 2 f i
+
.
2
r2
r2 sin2 () 2
sin ()
Now its time to derive
equation (6) we get
f
z
2f
z 2 .
cos()
f
1
f
sin() .
r
r
(12)
cos()
h f i 1
h f i
sin()
.
r z
r
z
(13)
Now we substitute equation (12) into equation (13) in order to eliminate any
partial derivatives with respect to z and we arrive at
2f
f
1
f i
h
=
cos()
cos()
sin()
z 2
r
r
r
h
1
f
1
f i
sin()
cos()
sin()
.
r
r
r
h
2f
1
f
1
2f i
cos() cos() 2 + 2 sin()
sin()
r
r
r
r
h
1
f
2f
1
f
sin() sin()
+ cos()
cos()
r
r
r
r
2 i
1
f
sin() 2 ,
r
2f
1
f
1
2 f i h 1 2 f
+
cos()
sin()
+ sin ()
cos()
sin()
r2
r2
r
r
r
r
2f
1
f
1
2f i
1
sin() cos()
+ 2 sin() cos()
+ 2 sin2 () 2 .
r
r r
cos2 ()
Now that we have all three terms of the right hand side of equation (3)(i.e.
2
and zf2 ), we add them all together (because of equation (3)) to get
2f 2f
x2 , y 2
r
r
i
h
2
f
1
1
1
f
f
sin() cos()
+ cos2 ()cos2 ()
sin() cos()
+
2
r
r
r
r
r
f
2f
1
2f
1
1
sin() cos()cos2 ()
2 sin() cos()cos2 ()
+ 2 cos2 ()cos2 () 2 +
r
r r
1 sin() cos()cos2 () f
1 cos() sin() cos() 2 f i h 1 2
f
+ sin ()
r2
r2
sin()
r
r
sin2 ()
1
2f
1 cos()sin2 () f
1 cos() cos() sin() 2 f
sin() cos()
+ 2
2
+
r
r r
sin()
r
sin()
1 sin2 () 2 f i h 2
2f
1 sin() cos() f
2
+
+
sin
()sin
()
r2
r2 sin2 () 2
r2
sin2 ()
1
f
1
2f
1
f
sin() cos()sin2 ()
+ sin() cos()sin2 ()
2 sin() cos()
+
2
r
r
r r
2f i h 1
f
1
f
1
sin() cos()
cos() sin()sin2 ()
+ cos2 ()sin2 ()
+
r
r
r
r
r
r
1
f
1
2f
1 sin() cos()cos2 () f
2
2
2
sin()
cos()sin
()
+
cos
()sin
()
+
r2
r2
2
r2
sin2 ()
1 cos() sin() cos() 2 f i h 1 2
f
1
2f
+ cos ()
+ cos() sin()
+
2
r
sin()
r
r
r
r
1 cos()cos2 () f
1 cos()cos() sin() 2 f
1 sin() cos() f
+ 2
2
+
2
r
sin()
r
sin()
r
sin2 ()
1 cos2 () 2 f i h 2 2 f
1
f
1
2f i
+ cos () 2 + 2 cos() sin()
cos() sin()
+
2
2
2
r sin ()
r
r
r
r
h1
f
1
2f
1
f
1
2f i
sin2 ()
sin() cos()
+ 2 sin() cos()
+ 2 sin2 () 2 .
r
r
r
r r
It may be hard to believe but the truth is that the above expression, after
some miraculous simplifications of course, reduces to the following succinct form
and we finally arrive at the Laplacian in spherical coordinates!
2 f
2f
1 cos() f
1 2f
1
1 2f
2 f
+ 2
.
+
+
+
r2
r2 2
r2 sin2 () 2
r r
r sin()
(14)
r2 sin2 () 2
2 Readers might be surprised how I got from expression (14) to expression (15) so fast!
Even though it is true that expression (14) is equivalent to expression (15), I wouldnt have
been able to write the expression (15) if I hadnt known about it beforehand.