A S A S A S A S A S A S A S A S S S F
A S A S A S A S A S A S A S A S S S F
A S A S A S A S A S A S A S A S S S F
1 s −a+ s +a 1æ 1
s 1 ö
fˆ (s ) = = çç + =
÷÷ Þ
( )(
s −a s +a 2 s −a s +a 2è s +a )
s −aø ( )( )
1 æç 1 1 ö
fˆ (s ) = + ÷ (1)
2 s çè 1 + a s 1 − a s ÷ø
1 a2 a3
a a4 a5 æ1ö
= 1+ + + + 2 + 2 + Oç 3 ÷ (2)
1− a s s s s s s s s ès ø
1 a a2 a3 a4 a5 æ1ö
= 1− + − + 2 − 2 + Oç 3 ÷ (3)
1+ a s s s s s s s s ès ø
ˆf (s ) = 1 é2 + 2 a + 2 a + Oæç 1 ö÷ ù = 1 + a + a + Oæç 1
2 4 2 4
ö
ê 3 ú ç 3 ÷÷ (4)
2 së s s2 è s øû s s s s2 s ès sø
f (t ) =
1
πt
+ 2a 2 t 4a 4
π
+
3
t
t
π
+ O t2 t ( ) (5)
In order to compare this with the analytical solution, we need the Taylor series of:
( )
a t a t
ò [1 − τ ( )]
2 2 2 2 2 2
ò e dτ =
−τ
ae a t erf a t = ae a t ae a t 2
+ O τ 4 dτ =
π 0 π 0
=
2a
[1 + a t + O(t )]éêa
2 a 3t t 3
ù
+ O t2 t ú = t− ( ) (6)
π ë 3 û
=
2a é
ê
π ë
2
a t + a 3t t + O t 2 t ú
3
ù
( )
û
f (t ) =
1
πt
2
+ ae a t erf a t = ( ) 1
πt
+ 2a 2
t 4a 4
π
+
3
t
t
π
+O t2 t ( ) (7)
2a
fˆ (s ) = 1 + −1 (8)
s
1 1 1
1+ x = 1+
2
x − x2 + x3 + O x4
8 16
( ) (9)
Thus we have:
2a 1 2a 1 4a 2 1 8a 3 æ 1ö a a
2
a3 æ 1ö
fˆ (s ) = 1 + −1 = − 2
+ 3
+ O ç 4 ÷ = − 2
+ 3
+ Oç 4 ÷ (10)
s 2 s 8 s 16 s è s ø s 2s 2s ès ø
a2 a3
f (t ) = a − t + t 2 + O(t 3 ) (11)
2 4