ENGINEERING - MATHEMATICS - 2 VTU Syllabus PDF
ENGINEERING - MATHEMATICS - 2 VTU Syllabus PDF
ENGINEERING - MATHEMATICS - 2 VTU Syllabus PDF
com
ENGINEERING MATHEMATICS-II
SUBJECT CODE: 14MAT21
Subject Code: 14MAT21 IA Marks: 25
Hours/Week: 04 Exam. Hours: 03
Total Hours: 50 Exam Marks: 100
Objectives:
To enable students to apply knowledge of Mathematics in various engineering fields by making
m
them to learn the following
co
• Ordinary differential equations
• Partial differential equations
• Double and triple integration
s.
bu
• Curvilinear co-ordinates
la
• Laplace transform
yl
Module – 1
lls
Differential equations-1:
.a
Linear differential equations with constant coefficients: Solutions of second and higher order
w
Module – 2
Differential equations-2:
Solutions of simultaneous differential equations of first order.
Linear differential equations with variable coefficients: Solution of Cauchy’s and Legendre’s
linear differential equations.
Nonlinear differential equations - Equations solvable for p, equations solvable for y, equations
solvable for x, general and singular solutions, Clairauit’s equations and equations reducible to
Clairauit’s form.
10 hrs
Page 1
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Module – 3
Partial Differential equations:
Formulation of PDE by elimination of arbitrary constants/functions, solution of non-
homogeneous PDE by direct integration, solution of homogeneous PDE involving derivative
with respect to one independent variable only. Derivation of one dimensional heat and wave
equations and their solutions by variable separable method.
Double and triple integrals: Evaluation of double integrals. Evaluation by changing the order
of integration and changing into polar coordinates. Evaluation of triple integrals. 10hrs
Module – 4
m
Integral Calculus:
co
Application of double and triple integrals to find area and volume. Beta and Gamma functions,
definitions, Relation between beta and gamma functions and simple problems.
s.
Curvilinear coordinates
bu
Orthogonal curvilinear coordinates - Definition, unit vectors and scale factors. Expressions for
gradient, divergence and curl. Cylindrical and spherical coordinate systems. 10hrs
la
yl
Module –5
lls
Laplace Transform:
.a
f (t)
eat f (t), t n f (t) and (without proof), periodic functions, unit-step function and Impulse
t
w
function - problems
w
Page 2
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Course Outcomes:
On completion of this course, students are able to,
• Use ordinary differential equations to model engineering phenomena such as electrical circuits,
forced oscillation of mass spring and elementary heat transfer. .
• Use partial differential equations to model problems in fluid mechanics, electromagnetic theory
and heat transfer.
• Evaluate double and triple integrals to find area, volume, mass and moment of inertia of plane
and solid region.
• Use curl and divergence of a vector function in three dimensions, as well as apply the Green's
m
Theorem, Divergence Theorem and Stokes' theorem in various applications like electricity,
co
magnetism and fluid flow.
• Use Laplace transforms to determine general or complete solutions to linear ODE.
s.
bu
la
yl
lls
.a
w
w
w
Page 3
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
INDEX SHEET
m
co
Module4: Integral Calculus……………………………….86-114
s.
bu
Module5: Laplace Transforms…………………………..116-167
la
yl
lls
.a
w
w
w
Page 4
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
MODULE - 1
DIFFERENTIAL EQUATIONS –I
INTRODUCTION:
m
LINEAR DIFFERENTIAL EQUATIONS OF SECOND AND
HIGHER ORDER WITH CONSTANT COEFFICIENTS
co
s.
bu
la
yl
lls
.a
w
w
w
Page 5
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 6
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 7
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
INVERSE DIFFERENTIAL OPERATOR AND PARTICULAR INTEGRAL
yl
lls
.a
w
w
w
Page 8
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 9
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 10
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 11
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 12
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 13
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Type2:
1.
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 14
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 15
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 16
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 17
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 18
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 19
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 20
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 21
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 22
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
The trial solution to be assumed in each case depend on the form of X. Choose PI from the
following table depending on the nature of X.
m
3 K e ax sin (ax+b) c1 e ax sin (ax+b)+ c2 e ax cos (ax+b)
or
co
ax
K e cos (ax+b)
4 K x n where n=0,1,2,3….. c0 c1 x c2 x2 .... cn 1 xn 1
cn x n
s.
5 K x n e ax where n=0,1,2,3….. eax c0 c1 x c2 x 2 .... cn x n
bu
6 K x n sin (ax+b) a0 sin(ax b) b0 cos(ax b)
or a1 .x.sin(ax b) b1x cos(ax b)
la
K x n cos (ax+b)
a 2 .x 2 .sin(ax b) b 2 x 2 cos(ax b)
yl
...........
lls
a n .x n .sin(ax b) b n x n cos(ax b)
7 K x n e dx sin (ax+b) e dx a0 sin(ax b) b 0 cos(ax b)
.a
n dx
K x e cos (ax+b)
a 2 .x 2 .sin(ax b) b 2 x 2 cos(ax b)
w
...........
w
a n .x n .sin(ax b) b n x n cos(ax b)
Assume PI y p c1e 3 x substituting this in the given d.e we determine the unknown coefficient as
Page 23
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
( D 2 3D 2) y 4e 3 x
9ce 3 x 9ce 3 x 2ce 3 x 4e 3 x
2ce 3 x 4e 3 x c 2
3x
yp 2e
d2y dy
2. Solve 2
2 4 y 2 x 2 3e x by the method of undetermined coefficients.
dx dx
Sol: We have ( D 2 2 D 4) y 2x 3e x
2 12 2 2 3i
A.E is m 2 2m 4 0 m 1 3i
2 2
x
m
yc e c1 cos 3x c2 sin 3x
co
Assume PI in the form y a1 x 2 a2 x a3 a4 e x
x
Dy 2a1 x a2 a4 e
s.
D2 y 2a1 a4 e x
bu
Substituting these values in the given d.e
la
x
We get 2a1 a4 e 2(2a1 x a2 a4 e x ) 4(a1 x 2 a2 x a3 a4 e x ) 2 x 2 3e x
yl
1
x 2 : 4a1 2 a1
2
.a
1
x : 4a1 4a2 0 4 4a 2 0
w
2
w
1
2 4a 2 0 4a 2 2 a2
2
w
c : 2a1 2a2 4 a3 0
1 1
2 2 4a3 0 a3 0
2 2
x
e : a4 2a 4 4a 4 3
3a4 3 a4 1
1 2 1 x
PI : y p x x e
2 2
x 1 2 1 x
y e c1 cos 3x c2 sin 3x x x e
2 2
Page 24
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
d2y
3. Solve by using the method of undetermined coefficients 9y x 3 e 2 x sin 3x
dx2
Sol: We have ( D 2 9) y x 3 e 2 x sin 3x
A.E is m 2 9 0 m2 9 m 3
yc c1e3 x c2 e 3x
y 3 Ax 2 2 Bx C 2 Ee 2 x 3F cos3x 3G sin 3x
y 6 Ax 2 B 4 Ee 2 x 9 F sin 3x 9G cos3x
m
Substituting these values in the given d.e, we get
co
6 Ax 2 B 4 Ee 2 x 9 F sin 3x 9G cos3x 9 Ax 3 Bx 2 Cx D Ee 2 x F sin 3x G cos3x
x 3 e 2 x sin 3x
Equating the coefficient of
s.
bu
1
x3 : 9 A 1 A
la
9
x2 : 9B 0 B 0
yl
1
lls
x : 6 A 9C 0 6 9C 0
9
.a
2 2 2
9C 0 9C C
3 3 27
w
C : 2B 9D 0 D 0
w
1
e 2 x : 4E 9E 1 5E 1 E
w
5
1
sin 3 x : 9 F 9G 0 F
18
cos3 x : 9G 9G 0 G 0
1 3 2x 1 2x 1
yp x e sin 3 x
9 27 5 18
Complete solution y yc yp
1 3 2x 1 2x 1
y c1e 3 x c2 e 3x
x e sin 3x
9 27 5 18
Page 25
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 26
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 27
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 28
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 29
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 30
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
MODULE - 2
DIFFERENTIAL EQUATIONS –II
SOLUTION OF SIMULTANEOUS DIFFERENTIAL EQUATIONS:
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 31
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 32
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 33
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 34
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 35
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
PROBLEMS:
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 36
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 37
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 38
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 39
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
Supposing that the LHS of (1) is expressed as a product of n linear factors, then the
co
equivalent form of (1) is
p f1 ( x, y ) p f 2 ( x, y ) ... p f n ( x, y ) 0 ....(2)
s.
p f1 ( x, y ) 0, p f 2 ( x, y ) 0... p f n ( x, y )
bu 0
All these are differential equations of first order and first degree. They can be solved by
the known methods. If F1 ( x, y, c) 0, F2 ( x, y, c) 0,... Fn ( x, y, c) 0 respectively represents the
la
solution of these equations then the general solution is given by the product of all these solution.
yl
Note: We need to present the general solution with the same arbitrary constant in each factor.
lls
2
.a
dy dy
1. Solve : y x y x 0
dx dx
w
yp 2 ( x y ) p x 0
w
(x ( x y)2
y) 4 xy
p
2y
( y x) ( x y )
p
2y
y x x y y x x y
ie., p or p
2y 2y
ie., p 1 or p x/ y
We have,
Page 40
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
dy
1 y x c or ( y x c) 0
dx
dy x
Also, or ydy xdx 0 y dy x dx k
dx y
y2 x2
ie., k or y 2 x 2 2k or ( x 2 y 2 c) 0
2 2
Thus the general solution is given by (y-x-c) (x 2 y 2 c) 0
2. Solve : x( y ' ) 2 (2 x 3 y ) y 6y 0
m
(2 x 3 y ) (2 x 3 y ) 2 24 xy
p
co
2x
(2 x 3 y ) (2 x 3 y ) 3y
p 2 or
s.
2x x
We have
bu
dy
2 dy 2 dx c or y 2 x c or ( y 2 x c ) 0
dx
la
dy 3 y dy dx dy dx
Also or 3 3 k
yl
dx x y x y x
lls
p( p y) x( x y)
w
3) Solve
2
Sol: The given equation is, p py x( x y) 0
w
y y2 4 x( x y)
p
2
y 4 x2 4 xy y2 y (2 x y)
p
2 2
2( y x)
ie., p x or p (y x)
2
We have,
Page 41
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
dy x2
x y k
dx 2
dy
Also, y x
dx
dy
ie., y x, is a linear d .e ( similar to the previous problem)
dx
P dx
P 1, Q x; e ex
Hence ye x xe x dx c
ie., ye x ( xe x e x ) c, int egrating by parts.
Thus the general solution is given by (2 y x 2 c ) e x ( y x 1) c 0
m
co
Equations solvable for y:
s.
We say that the given differential equation is solvable for y, if it is possible to express y
in terms of x and p explicitly. The method of solving is illustrated stepwise.
bu
Y=f(x, p)
We differentiate (1) w.r.t x to obtain
la
dy dp
p F x, y ,
yl
dx dx
Here it should be noted that there is no need to have the given equation solvable for y in
lls
the explicit form(1).By recognizing that the equation is solvable for y, We can proceed to
.a
differentiate the same w.r.t. x. We notice that (2) is a differential equation of first order in p
and x. We solve the same to obtain the solution in the form. ( x, p, c) 0
w
By eliminating p from (1) and (3) we obtain the general solution of the given
w
Remark: Suppose we are unable to eliminate p from (1)and (3), we need to solve for x and y
from the same to obtain.
x F1 ( p, c), y F2 ( p, c)
Which constitutes the solution of the given equation regarding p as a parameter.
We say that the given equation is solvable for x, if it is possible to express x in terms of y
and p. The method of solving is identical with that of the earlier one and the same is as follows.
x = f(y, p )
Page 42
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
(2) Being a differential equation of first order in p and y the solution is of the form.
( y, p, c ) 0
By eliminating p from (1) and (3) we obtain the general solution of the given d.e in the form
G(x, y, c) = 0
Note: The content of the remark given in the previous article continue to hold good here also.
Solve: y 2 px tan 1 ( xp 2 )
m
1.
co
Sol : By data, y 2 px tan 1 ( xp 2 )
s.
The equation is of the form y = f (x, p), solvable for y. bu
Differentiating (1) w.r.t.x,
la
dp 1 dp
yl
p 2p 2 x 2 4
x.2 p p2
dx 1 x p dx
lls
dp 1 dp
ie., p 2x 2 4
2 xp p2
.a
dx 1 x p dx
w
p2 dp p
w
ie., p 2x 1
1 x2 p4 dx 1 x 2 p 4
w
1 x2. p4 p dp p 1 x 2 p 4
ie., p 2x
1 x2 p4 dx 1 x 2 p 4
ie., log x 2log p k
consider y 2 px tan 1`( xp 2 )
and xp 2 c
Using (2) in (1) we have,
y = 2 c / x . x tan 1(c)
Thus y 2 cx tan 1 c, is the general solution.
Page 43
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2. Obtain the general solution and the singular solution of the equation y px p2 x4
Consider , y px p2 x4
x p c or x 2 p c or p c / x2
m
Using (2) in (1) we have, y (c / x 2 ) x (c 2 / x 4 ) x 4
co
Thus xy c c 2 x is the general solution.
s.
bu
Now, to obtain the singular solution, we differentiate this relation partially w.r.t c,
treating c as a parameter.
la
yl
1 1
xy x
w
2x 4 x2
w
Page 44
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
ie.,sin p x c or x sin p c
Thus we can say that y p sin p cos p and x sin p - c constitutes
the general solution of the given d.e
Note :sin p x c p sin 1 ( x c).
We can as well substitute for p in (1) and present the solution in the form,
y ( x c)sin 1( x c) cos sin 1 ( x c)
m
p
Differentiating w.r.t y we get
co
2 1 y dp dp
2
p y
p p p dy dy
s.
1 y dp
p 1 0
bu
p p dy
1 dp
la
Ignoring p which does not contain , this gives
p dy
yl
y dp dy dp
lls
1 0 or 0
p dy y p
.a
Integrating we get
yp c........(2)
w
y2 2cx c 2
w
5) Solve p 2 2 py cot x y2 .
Page 45
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
y
cot x cos ecx
p
y
cot x cos ecx
dy / dx
dy sin x dy sin x
dx and
y cos x 1 y cos x 1
Integrating these two equations we get
y (cos x 1) c1 and y (cos x 1) c2
general solution is
y (cos x 1) c y (cos x 1) c 0
2
6) Solve: p 4 x5 p 12 x 4 y 0 , obtain the singular solution also. .
m
Sol: The given equation is solvable for y only.
co
p 2 4 x5 p 12 x 4 y 0 ...........(1)_
s.
2 5
p
4x p
y f ( x, p )
12 x 4
bu
Differentiating (1) w.r.t.x,
la
dp dp
2p 4 x5 20 x 4 p 12 x 4 p 48 x 3 y 0
dx dx
yl
dp p2 4 x5 p
(2 p 4 x 5 ) 8 x 3 ( xp ) 0
lls
dx 2x4
dp 2p
.a
( p 2 x5 ) ( p 2 x5 )
dx x
w
dp 2 p
0
dx x
w
2 2
p c x equation (1) becomes
4 2 3
c 4c x 12 y
Setting c 2 k the general solution becomes
k 2 4kx3 12 y
Differentiating w.r.t k partially we get
2k 4 x 3 0
Using k 2 x3 in general solution we get
x6 3 y 0 as the singular solution
Page 46
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
12 y
dy p y
co
dp 2 p 3 8 y 2 p3 4y2
dy p y
s.
2 dp 3 p3 4 y 2
( p 4y2 )
bu
p dy y
2 dp 1
la
p dy y
yl
cy cy 4 xy cy 8 y 2 0
3
w
2
Dividing throughoutby y y y we have,
w
c c 4x c 8 y 0
w
c (c 4 x ) 8 y
Thus the generalsolutionisc(c 4 x) 2 64 y
Clairaut’s Equation
The equation of the form y px f ( p) is known as Clairaut’s equation.
This being in the form y = F (x , p), that is solvable for y, we differentiate (1) w.r.t.x
dy dp dp
p p x f ( p)
dx dx dx
Page 47
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
dp
This implies that 0 and hence p=c
dx
Using p c in (1) we obtain the genertal solution of clairaut's equation in the form
y cx f (c)
a
1. Solve: y px
p
Sol: The given equation is Clairaut’s equation of the form y px f ( p) , whose general solution
is y cx f (c)
a
Thus the general solution is y cx
m
c
co
Singular solution
Differentiating partially w.r.t c the above equation we have,
s.
a
0 x
c2
bu
a
c
la
x
yl
Hence y cx (a / c) becomes,
lls
y a / x. x a x / a
Thus y 2 4ax is the singular solution.
.a
2. Modify the following equation into Clairaut’s form. Hence obtain the associated general
w
2 py kp a 0
and singular solutions xp
w
w
Sol : xp 2 py kp a 0, by data
ie., xp 2 kp a py
p ( xp k ) a
ie., y
p
a
ie., y px k
p
Here (1) is in the Clairaut’s form y=px+f(p) whose general solution is y = cx + f (c)
Page 48
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
a
Thus the general solution is y cx k
c
Now differentiating partially w.r.t c we have,
a
0 x
c2
c a/x
Hence the general solution becomes,
y - k = 2 ax
Thus the singular solution is (y - k)2 4ax.
Remark: We can also obtain the solution in the method: solvable for y.
m
co
3. Solve the equation (px – y) (py + x) = 2p by reducing into Clairaut’s form, taking the
substitutions X = x2 , Y = y2
s.
dX
Sol : X x2 2x
dx
bu
dY
Y y2 2y
la
dy
yl
dy dy dY dX dY
Now, p and let p
dx dY dX dx dx
lls
1
ie., p .P.2 x
.a
2y
w
X
ie., p P
w
Y
w
Consider ( px y )( py x) 2p
X X X
ie., P X Y P Y X 2 P
Y Y Y
ie., ( PX Y ) ( P 1) 2P
2P
ie., Y PX is in the Clairaut's form and hence the associated genertal solution is
P 1
2c
Y cX
c 1
2c
Thus the required general solution of the given equation is y 2 cx 2
c 1
Page 49
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
X X X
P X Y P Y X 2 P
co
Y Y Y
( PX Y )( P 1) 2 P
s.
2P
Y PX
bu
P 1
Is in the Clairaut’s form and hence the associated general solution is
la
2c
Y cX
c 1
yl
2c
Thus the required general solution of the given equation is y 2 cx 2
lls
c 1
5) Obtain the general solution and singular solution of the Clairaut’s equation xp 3 yp 2 1 0
.a
y 2
y px i s in the Clairaut ' s form y px f ( p )
p p2
w
1
Thus general solution is y cx
c2
Differnetiating partially wr
. .t. c we get
1/3
2 2
0 x 3 c
c x
Thus general solution becomes
1/3 2/3
2 x
y x 22/3 y 3 x 2/3
x 2
or 4 y 3 27 x 2
Page 50
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 51
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
MODULE – 3
PARTIAL DIFFERENTIAL EQUATIONS
Introduction:
Many problems in vibration of strings, heat conduction, electrostatics involve two or more variables.
Analysis of these problems leads to partial derivatives and equations involving them. In this unit we first
discuss the formation of PDE analogous to that of formation of ODE. Later we discuss some methods of
solving PDE.
Definitions:
m
An equation involving one or more derivatives of a function of two or more variables is called a partial
co
differential equation.
s.
The order of a PDE is the order of the highest derivative and the degree of the PDE is the degree of highest
order derivative after clearing the equation of fractional powers.
bu
A PDE is said to be linear if it is of first degree in the dependent variable and its partial derivative.
la
In each term of the PDE contains either the dependent variable or one of its partial derivatives, the PDE is
yl
P p + Q q =R
Page 52
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2 z 2y 1 1 z q
2
or 2
y b b y x y
m
z z
(x a ) z 0, ( y b) z 0
co
x y
s.
Substituting for (x- a) and (y- b) from these in (i), we get bu
2 2
2 z z
z 1 k 2 This is the required partial differential equation.
la
x y
yl
z
w
a cy..(ii)
x
w
z
w
b cx..(iii)
y
2
z
c Using this in (ii) and (iii)
x y
2
z z
a y
x x y
Page 53
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2
z z
b x
y x y
2 2 2
z z z z z
z x y y x xy
x x y y x y x y
2
z z z
z x y xy
x y x y
x2 y2 z2
m
(5) 2 1
a b2 c2
co
Sol: Differentiating partially w.r.t. x,
2x 2z z x z z
s.
bu
0, or
a2 c2 x a2 c2 x
la
Differentiating this partially w.r.t. x, we get
yl
2 2 2
1 1 z z c2 z z 2
lls
z 2 or 2 z
a2 c2 x x a x x2
.a
2 2 2 2
z z z z z z z z
z 2 z
w
y y y y x x x x2
Note:
1) z = f(x2 + y2)
Page 54
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
z z
p f ' (x2 y 2 ).2 x, q f ' ( x2 y 2 ).2 y
x y
m
co
2 2
z z
2 2
0
t x
(3) x + y + z = f(x2 + y2 + z2)
s.
bu
Sol:Differentiating partially w.r.t. x and y
la
z z
yl
1 f ' (x2 y2 z 2 ) 2x 2z
x x
lls
.a
z z
1 f ' (x2 y2 z 2 ) 2 y 2z
y y
w
1 ( z / x) 1 ( z / y)
2 f ' ( x2 y2 z2 )
w
x z ( z / x) y z ( z / y)
w
z z
( y z) ( z x) x y is the required pde
x y
(4) z = f ( xy / z ).
z xy y xy z
f'
x z z z2 x
Page 55
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
z xy x xy z
f'
y z z z2 x
xy z/ x z/ y
f'
z ( y / z ){1 ( x / z )( z / x} ( x / z ){1 ( y / z )( z / y}
z z
x y
x y
or xp = yq is the required pde.
m
z 1
Sol : 2 y 2 f '(1/ x log y)
co
y y
s.
z 1
2 f ' (1 / x log y )
x2
x
bu
la
z z
2 f ' (1 / x log y) x2 y 2y
x y
yl
lls
2 z z
Hence x y 2y 2
x y
.a
w
z z
Sol : ( y) y '( x); x '( y) ( x)
w
x y
Page 56
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
7) Form the partial differential equation by eliminating the arbitrary functions from
z = f(y-2x) + g(2y-x) (Dec 2011)
Sol: By data, z = f(y-2x) + g(2y-x)
z
p 2 f ( y 2 x) g (2 y x)
x
z
q f ( y 2 x) 2 g (2 y x)
y
2
z
r 4 f ( y 2 x) g (2 y x)...............(1)
x2
m
2
z
s 2 f ( y 2 x) 2 g (2 y x).........(2)
co
x y
s.
2
z
t f ( y 2 x) 4 g (2 y x)................(3)
bu
y2
la
(1) 2 (2) 2r s 6 f ( y 2 x)..............(4)
yl
lls
2r s
2 or 2r 5s 2t 0
w
2s t
w
2 2 2
z z z
w
Thus 2 2
5 2 0 is the required PDE
x x y y2
dx dy dz
Sol :
yz zx xy
From the first two and the last two terms, we get, respectively
Page 57
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
dx dy dy dz
or xdx ydy 0 and or ydy zdz 0.
y x z y
Integrating we get x2 - y2 = a, y2 – z2 = b.
Φ(x2-y2, y2 –z2) = 0
dx dy dz
Sol :
y2 xy x( z 2 y )
m
From the first two ratios we get
co
x2 + y2 = a from the last ratios two we get
dz z
s.
bu
2
dy y
la
dz z
2 ordinary linear differential equation hence
dy y
.a
yz – y2= b
w
w
dx dy dz
Sol : 2
zx zy y x2
dx dy
, or xdy ydx 0 or d(xy) 0,
x y
on integration, yields xy = a
Page 58
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Φ(xy,x2+y2+z2)=0
y z z x x y
(4) Solve: p q
yz zx xy
yz zx xy
Sol : dx dy dz
y z z x x y
x dx + y dy + z dz = 0 …(i)
m
Integrating (i) we get
co
x2 + y2 + z2 = a
yz dx + zx dy + xy dz = 0 …(ii)
s.
bu
Dividing (ii) throughout by xyz and then integrating,
la
yl
we get xyz = b
lls
Φ( x2 + y2 + z2, xyz ) = 0
.a
dx dy dz
w
Sol : ..(i)
x 2z 4 zx y 2x2 y
w
xy – z2 = b ….(iii)
5) Solve zxy sin x sin y for which z y 2sin y when x 0 and z 0
Page 59
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Sol: Here we first find z by integration and apply the given conditions to determine the arbitrary
functions occurring as constants of integration.
z
The given PDF can be written as sin x sin y
x y
Integrating w.r.t x treating y as constant,
z
sin y sin x dx f ( y) sin y cos x f ( y)
y
Integrating w.r.t y treating x as constant
z cos x sin y dy f ( y ) dy g ( x)
z cos x ( cos y ) F ( y ) g ( x),
m
where F ( y ) f ( y ) dy.
co
Thus z cos x cos y F ( y ) g ( x )
s.
z
bu
Also by data, 2sin y when x 0. U sin g this in (1)
y
la
2sin y ( sin y ).1 f ( y ) (cos 0 1)
yl
2
w
Page 60
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
u u
3 ( xy ) 2 ( xy ) 0
x x
dX dY 3 dX 2 dY
3Y 2X 0
dx dy X dx Y dy
3 dX 3dX
Let K kdx
X dx X
Kx
3log X kx c1 log X c1
3
m
kx
c1
co
3
X e
2 dY dY Kdy
s.
Let k
Y dy Y 2 bu
ky
Kdy 2
c2
log Y c2 Y e
2
la
x y
K c1 c2
lls
3 2
U e
x
Also u ( x1o) 4e
.a
2x kx
w
k
x 6 x 3
i.e., 4e Ae 4e Ae
w
x y
3
3 2
U 4e is required solution.
u u
2) Solve by the method of variables 4 3u, giventhat u (0, y) 2e5 y
dx y
u u
Solution: Given 4 3u
x y
Assume solution of (1) as
u XY where X X ( x); Y Y ( y)
Page 61
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
4 ( XY ) ( XY ) 3 XY
x y
dX dY 4 dX 1 dY
4Y X 3 XY 3
dx dy X dx Y dy
4 dX 1 dY
Let k, 3 k
X dx Y dy
Separating var iables and int egrating we get
kx
log X c1 , log Y 3 k y c2
4
kx
c1
X e 4
and Y e3 k y c2
m
kx kx
c1 c2 3 k y 3 k y
Hence u XY e e 4
Ae 4
where A e c1 c2
co
put x 0 and u 2e5 y
s.
The general solution becomes bu
2e5 y Ae 3 k y A 2 and k 2
Particular solution is
la
x
5y
yl
2
u 2e
lls
variables.
w
We need to obtain the solution of the ODEs by taking the constant k equal to
w
Various possible solutions of the one dimensional heat equation ut =c2uxx by the method of
separation of variables.
2
u 2 u
Consider c
t x2
Page 62
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2
XT 2 XT dT 2d2X
c or X c
t x2 dt dx 2
1 dT 1 d2X
Dividing by c2XT we have
c 2T dt X dx 2
1 d2X 1 dT
=k and =k
X dx 2 c 2T dt
dT
d2X c 2 kT 0
kX 0 and dt
dx 2
m
co
D2 k X 0 and D c 2 k T 0
s.
d2 d
Where D2 = 2
in the first equation and D = in the second equation
dx dt
bu
Case (i) : let k=0
la
AEs are m=0 amd m2=0 amd m=0,0 are the roots
yl
lls
T = c1e0t c2 x c3 e0 x
.a
c1 and X c2 x c3
w
U= XT= c1 c2 x c3
w
Page 63
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2 2
Or u(x,t) = c'1ec p t
(A’ e px +B e px
) where c1’c2’=A’ and c1’c3’=B’
m
c 2 p 2t
T c''1e and X c''2 cos px c''3 sin px
co
Hence the solution of the PDE is given by
s.
c 2 p 2t
u XT c''1e .(c''2 cos px c''3 sin px) bu
c2 p2t
u( x, t ) e ( A'' cos px B'' sin px)
la
yl
2
u u
lls
1. Solve the Heat equation c2 given that u(0,t)=0,u(l,0)=0 and u(x,0)= 100x/l
t x2
.a
l l
2 100 x n x 200 n x
Soln: bn sin dx = x sin dx
w
2
l 0 l l l 0 l
w
l
n x n x
w
x. cos sin
200 l l
bn 2
1 2
l n /l n /l
0
n n 1
200 1 200 1 200 1
bn . l cos n .
l2 n n n
n 1 n2 2 2
c t
200 1 n x
Thus u ( x, t ) e 2
sin
n 1 n l l
Page 64
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2
u u
2. Obtain the solution of the heat equation c2 given that u(0,t)=0,u(l,t)and
t x2
2Tx l
in 0 x
l 2
u(x,0) =f(x)where f ( x)
2T l
l x in x l
l 2
l
2 n x
Soln: bn f ( x) sin dx
l 0 l
l
2 l
2 2Tx n x 2Tx n x
bn sin dx (l x) sin dx
l l l l l
m
0 l
2
co
l
2 l
4T n x n x
s.
x sin dx (l x) sin dx
l 0
l l l bu
2
8T n
la
bn 2 2
sin
n 2
yl
n2 2 2
c t
8T 1 n n x
.a
Thus u ( x, t ) 2 2
sin e 2
sin
n 1 n 2 l l
w
2
u u
w
u(x,0) =3sin x
p 2t
Soln: u( x, t ) e ( A cos px B sin px)............................(1)
p 2t
0= e (A) thus A=0
p 2t
0= e (Bsinp)
Page 65
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
n2 2c2t
u ( x, t ) e ( B sin n x)
n2 2c2t
In general u ( x, t ) bn e sin n x
n 1
m
We substitute these values in the expanded form and then get
co
2 t
u ( x, t ) 3e (sin x)
s.
Various possible solutions of the one dimensional wave equation utt =c2uxx by the method of
bu
separation of variables.
2 2
la
u u
Consider 2
c2
t x2
yl
2 2
w
XT 2 XT d 2T d2X
c or X c2
t2 x2 dt 2 dx 2
w
w
1 d 2T 1 d2X
Dividing by c2XT we have
c 2T dt 2 X dx 2
1 d2X 1 d 2T
=k and =k
X dx 2 c2T dt 2
d 2T
d2X c 2 kT 0
kX 0 and dt 2
dx 2
Page 66
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
D2 k X 0 and D 2 c 2 k T 0
d2 d2
Where D2 = in the first equation and D 2
= in the second equation
dx 2 dt 2
AEs are m=0 amd m2=0 amd m=0,0 are the roots
T = c1e0t c1 and X c2 x c3 e0 x c2 x c3
m
U= XT= c1 c2 x c3
co
Or u(x,t) =Ax+B where c1c2=A and c1c3=B
2
p 2t
T c'1ec andX c'2e px c'3e px
w
2
p 2t
u XT c'1ec .( c'2e px c ' 3e px
)
w
2 2
Or u(x,t) = c'1ec p t
(A’ e px +B e px
) where c1’c2’=A’ and c1’c3’=B’
c 2 p 2t
T c''1e and X c''2 cos px c''3 sin px
Page 67
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
c 2 p 2t
u XT c''1e .(c''2 cos px c''3 sin px)
c2 p2t
u( x, t ) e ( A'' cos px B'' sin px)
1. Solve the wave equation utt=c2uxx subject to the conditions u(t,0)=0 ,u(l,t)=0,
u
x, 0 0 and u(x,0) =u0sin3( x/l)
t
n x n ct
Soln: u x, t bn sin cos
m
n 1 l l
co
n x
s.
u x,0 bn sin
n 1 l
bu
x n x
u0 sin 3 bn sin
la
l n 1 l
yl
3 3 x 1 3 x n x
u0 sin sin bn sin
lls
4 l 4 l n 1 l
.a
3u0 x u0 3 x x 2 x 3 x
sin sin b1 sin b2 sin b3 sin
4 l 4 l l l l
w
3u0 u0
b1 , b2 0 , b3 , b4 0 b5 0 ,
w
4 4
3u0 x ct u0 3 x 3 ct
u( x, t ) sin cos sin cos
4 l l 4 l l
2. Solve the wave equation utt=c2utt subject to the conditions u(t,0)=0 ,u(l,t)=0,
u
x, 0 0 when t=0and u(x,0) =f(x)
t
Page 68
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
n x n ct
Soln: u x, t bn sin cos
n 1 l l
n x
Consider u(x,0) = bn sin
n 1 l
n x
F(x) = bn sin
n 1 l
The series in RHS is regarded as the sine half range Fourier series of f(x) in (0,l) and hence
l
2 n x
bn f ( x) sin dx
m
l 0 l
co
Thus we have the required solution in the form
s.
n x n ct
u x, t bn sin cos
n 1 l l
bu
DOUBLE INTEGRAL
la
yl
lls
.a
w
w
w
Page 69
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
PROBLEMS:
la
yl
lls
.a
w
w
w
Page 70
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 71
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
Page 72
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 73
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 74
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 75
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 76
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 77
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Triple Integrals:
The treatment of Triple integrals also known as volume integrals in R 3 is a simple and straight
extension of the ideas in respect of double integrals.
Let f(x,y,z) be continuous and single valued function defined over a region V of space. Let V be
divided into sub regions v1 , v2 ....... vn in to n parts. Let ( xk , yk , zk ) be any arbitrary point
Page 78
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
n
within or on the boundary of the sub region vk . From the sum s f ( xk , yk , zk ) vk
k 1
…………(1)
Sub region approaches zero the sum (1) has a limit then the limit is denoted by f ( x, y, z )dv
V
For the purpose of evolution the above triple integral over the region V can be expressed as an
iterated integral or repeated integral in the form
m
b h( x) ( x, y )
f ( x, y, z )dxdydz f ( x, y, z )dz dy dx
co
V a g ( x) ( x, y )
s.
Where f(x,y,z) is continuous in the region V bounded by the surfaces z= z bu ( x, y), z ( x, y) ,
y g ( x), y h( x), x a, x b . the above integral indicates the three successive integration to be
performed in the following order, first w.r.t z, keeping x and y as constant then w.r.t y keeping x
la
as constant and finally w.r.t.x.
yl
Note:
lls
If the limits are not constants the integration should be in the order in which dx, dy, dz is
w
Evaluation of the integral may be performed in any order if all the limits are constants.
If f(x,y,z) = 1 then the triple integral gives the volume of the region.
w
1 2 2
1. Evaluate xyz 2 dxdydz
0 0 1
1 2 2 1 2 2
2 x 2 yz 2
Sol : xyz dxdydz dydz
0 0 1 0 0
2 1
1 2 2
2 y2 z2
2 yz dydz
0 0
2 1
Page 79
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 2
2 y2 z2 y2 z2
dz
0
2 4 0
1 2
2 2 y2 z2
y z dz
0
4 0
1 2
3y2 z2
dz
0
4 0
1
a a a
m
2. Evaluate ( x2 y2 z 2 )dxdydz
co
0 0 0
a a a a a a
2 2 2 x3 2 2
s.
Sol : (x y z )dxdydz y x z x dydz
0 0 0 0 0
bu 3 0
a a
a3
y 2a z 2 a dydz
la
0 0
3
a a
yl
a3
[ y 2a z 2 a dy ]dz
lls
0 0
3
a a
a3 y y 3a
.a
2
z ay dz
3 3
w
0 0
a 4 4
w
a a
a 2 z 2 dz
3 3
w
0
a
a4 z a4 z a2 z3
3 3 3 0
5 5 5
a a a
3 3 3
a5
Page 80
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
a a2 x2 a2 x2 y 2
3. Evaluate
x y z dx dy dz
0 0 0
a a2 x2 a2 x2 y 2
Sol : I x y z dz dy dx
0 0 0
a2 x2 y 2
a a2 x2 2
xyz
dy dx
0 0
2 0
a a2 x2
xy 2 2 2
m
(a x y ) dy dx
2
co
0 0
a a2 x2
1
(xya 2 x 3 y x y3 )dy dx
s.
20 0
bu
a 2 2 3 2 2 a2 x2
1 xy a x y xy
dx
la
20 2 2 2 0
yl
a
1
(a 4 x x 5 2 a 2 x 3 )
lls
80
.a
a
1 4 x2 x6 2a 2 x 4 a6
a
w
8 2 6 4 0
48
w
w
4. Evaluate xyz dx dy dz over the region R enclosed by the coordinate planes and the
R
plane x + y + z=1
Sol: In the given region, z varies from 0 to 1 – x – y
For z-=0, y varies from 0 to 1 – x. For y=0,x varies from 0 to 1.
1 1 x1 x y
x y z dx dy dz x y z dx dy dz
R x 0y 0 z 0
1 1 x
1
x y (1 x y) 2 dy dx
0 0
2
Page 81
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 1 x
1
x (1 x) 2 y 2(1 x) y 2 y 3 dy dx
20 0
1
1 1 2 1
x (1 x) 2 (1 x) 2 (1 x)(1 x)3 (1 x) 4 dx
20 2 3 4
1 1
1 1 (1 x)6
x(1 x) 4 dx
24 0 24 30 0
1
720
Change of variable in triple integrals
m
Computational work can often be reduced while evaluating triple integrals by changing
co
the variables x, y, z to some new variables u, v, w, which related to x,y,z and which are
s.
such that the bu
x x x
u v w
la
( x, y , z ) y y y
Jacobian J 0
yl
(u , v, w) u v w
z z z
lls
u v w
.a
f ( x, y, z )dxdydz
R
w
(u, v, w) Jdudvdw........(1)
w
R is the region in which (x,y,z) vary and R is the corresponding region in which
(u,v,w)vary and (u, v, w) f x(u, v, w), y (u, v, w), z (u, v, w)
Once the triple integral wrt (x,y,z) is changed to triple integral wrt (u,v,w) by using the
formula(1), the later integral may be evaluated by expressing it in terms of repeated
integrals with appropriate limit of integration
Page 82
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
Hence dxdydz has to be changed to R dR d dz
co
Thus we have
s.
f ( x, y, z )dxdydz
R
bu
( R, , z ) RdRd dz
R
la
R is the region in which ( R, , z) vary, as (x,y,z) vary in R
yl
( R, , z) f ( R cos , R sin , z)
lls
Suppose (x,y,z) are related to three variables (r, , ) through the relations
w
x r sin cos , y r sin sin , z r cos . Then (r, , ) are called spherical polar
w
coordinates.
w
PROBLEMS:
Page 83
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 1
R3dR 2
sin cos zdz
0 0 0
1 1 cos 2 2
R 3dR
4 0 2 0
1 1
R 3dR
4 0
1
16
m
co
Sol: In the region, r varies from 0 to a, varies from 0 to and varies from 0 to.
2
s.
The relations between Cartesian and spherical polar coordinates are
x r sin cos , y r sin sin , z r cos .....(1)
bu
Also dxdydz r 2 sin drd d
la
We have x 2 y2 z2 a 2 .....(2)
yl
lls
a
xyzdxdydz 2 2
r sin cos r sin sin r cos r 2 sin drd d
0 0 r 0
R
.a
2
0
w
a6
cos cos 0
w
96
a6
48
Page 84
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 85
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
MODULE-4
INTEGRAL CALCULUS
Application of double integrals:
Introduction: we now consider the use of double integrals for computing areas of plane and
curved surfaces and volumes, which occur quite in science and engineering.
m
f ( x, y )dA f ( x, y )dxdy f ( x, y )dydx f ( x, y )dxdy
co
A R a y1 ( x ) c x1 (y)
b y2 ( x ) d x2 (y)
dA dxdy dydx dxdy..........(1)
s.
A R a y1 ( x ) c x1 (y) bu
The integral dA represents the total area of the plane region R over which the iterated integral
A
la
are taken . Thus (1) may be used to compute the area A. nNote that dx dy is the plane area
yl
Also dxdy rdrd , rdrd is the plane area element in polar form.
R R
.a
Let the curves AB and CD be y1 f1 ( x)andy2 f 2 ( x) . Let the ordinates AC and BD be x=a and
w
x=b. So the area enclosed by the two curves and x=a and x=b is ABCD. Let p(x,y) and be
Q(x x, y y) two neighbouring points, then the area of the small rectangle PQ= x y
y2 y2
Area of the vertical strip = lim x y x dy
y 0 y1 y1
Since x the width of the strip is constant throughout, if we add all the strips from x=a to x=b
we get
h y2 b y2
The area ABCD =
limy 0 a
x dy dx dy
y1 a y1
Page 86
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
b y2
Area= dxdy
a y1
x2 y2
1. Find the area of the ellipse 1 by double integration .
a2 b2
b
Soln: For the vertical strip PQ, y varies from y =0 to y a2 x 2 when the strip is slided
a
from CB to A, x varies from x=0 to x=a
b 2 2
a x
a a
m
Therefore Area of the ellipse=4 Area of CAB 4 dydx
co
x 0 y 0
b 2 2
a x
s.
a a a b 2 2
a x
a
4 dy dx 4 y 0 dx
0 0 0
bu
a
la
a
b b x a2 x2 a2 x
4 a2 2
x dx 4 sin 1
a a 2 2 2
yl
0 0
2 2
lls
b a b a
4 sin 1 1 4 . . ab
a 2 a 2 2
.a
Solving (1) and (2) we get the point of intersections (0,0) and (4a,4a) . The shaded portion
in the figure is the required area divide the arc into horizontal strips of width y
y2
x varies from p, to Q 4ay and then y varies from O, y=0 to A, y=4a .
4a
4a 4 ay 4a 4 ay
dy dx dy x
y2
0 y2 0
4a
4a
Page 87
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
4a
3
4a
y2 y 2
1 y3
4ay dy 4a .
0
4a 3 4a 3
2 0
3
4 a 1 3
4a 2 4a
3 12a
32 2 16 2 16 2
a a a
3 3 3
The double integral can made use in evaluating the surface area of a surface.
m
Consider a surface S in space .let the equation of the surface S be z=f(x,y) . it can be that surface
co
area of this surface is
1
s.
2 2 2
z z
Given by s 1 dxdy
x y
bu
A
la
Where A the region representing the projection of S on the xy-plane.
yl
1
2 2 2
z z
w
s 1 dydz
A
z y
w
and
1
w
2 2 2
z z
s 1 dzdx
A
z x
Soln: the required surface arc is twice the surface are of the upper part of the given sphere,
whose equation is
1
2 2 2 2
z a x y z 0
1
z
this, gives, a2 x2 y2 2 2x
x
Page 88
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
x
1
a2 x2 y2 2
z y
similarly , 1
x
a2 x2 y2 2
2 2
z z a2
1
x x a2 x2 y2
hence, the, required , surface, area , is
1 1
2 2 2
z z a2 2
s 2 1 dxdy 2 dxdy
A
x y A
a2 x2 y2
m
Where A the projection of the sphere on the xy-plane . we note that this projication is the area
co
bounded by circle x2+y2=a2.hence in A ,Ѳ varies from 0to2
s.
And r varies from 0to a, where (r, Ѳ) are the polar coordinates. put x=cos θ ,y=sin θ dxdy=rdrd θ bu
2 a 2 a
a r
s 2 rdrd 2 d rdr
a2 r2 a2 r2
la
0r 0 0 0
2 a 2
2 2 2
2a 2 4 a2
yl
2a d a r 0 2a d a 0
0 0
lls
2) Find the surface area of the portion of the cylinder x2+z2=a2 which lies inside the
.a
cylinder x2+y2+=a2.
w
Soln: Let s1 be the cylinder x2+z2=a2 and s2 be the cylinder x2+z2=a2 for the cylinder
w
z x z
s1 , 0
x z y
w
2 2
z z x2 z2 x2 a2
so that ,1 1 0
x y z2 z2 a2 x2
The required surface area is twice the surface area of the upper part of the cylinder S 1 which lies
inside the cylinder x2+y2=a2. Hence the required surface area is
1
2 2 2
z z a
s 2 1 dA 2 dA,
A
x y a a2 x2
Page 89
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Where A is the projection of the cylinder S1 on the x y plane that llies with in the cylinder
S2:x2+y2=a2. In Ax varies from –a toa and for each x,y varies from a2 x 2 to a 2 x2
a a2 x2
a
s 2 dydx
x ay a 2
x 2 a2 x2
a
1 a2 x2
2a y dx
a2 x2
a a2 x2
a
1
2a 2 a2 x 2 dx
2 2
a a x
a
a
4a dx 4a x a 4a a a 8a 2
m
a
co
Volume underneath a surface:
Let Z=(x,y)be the equation of the surface S. let P be a point on the surface S.let A denote the
s.
orthogonal projection of S on the xy- plane . divide it into area elements by drawing thre lines bu
parallel to the axes of x and y on the elements x y as base ,erect a cylinder having generators
parallel to QZ and meeting the surface S in an element of area s .the volume underneath the
la
surface bounded by S, its projection A on xy plane and the cylinder with generator through the
boundary curve of A on the xy plane and parallel to OZ is given by,
yl
lls
v f x, y dxdy Zdxdy
A A
.a
x2 y2 z2
w
Sol: Let S denote the surface of the ellipsoid above the xy-plane .the equation of this surface
w
x2 y2 z2
1z 0
a2 b2 c2
is 1
x2 y2 2
or, z c1 f x, y
a2 b2
The volume of the region bounded by this surface and the xy-plane gives the volume v1of the
upper half of the full ellipsoid .this volume is given by v1 x, y dxdy
A
Page 90
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
x2 y2
Note that A is the area bounded by the ellipse 1
a2 b2
1
2 2
x y 2 2
v1 1 dxdy c ab
A a2 b2 3
2
ac
3
2 4
The volume of the full ellipsoid is 2v1.thus the required volume is v 2. abc abc
3 3
m
Let y=f(x) be a simples closed plane curves enclosing an area A. suppose this curve is revolved
co
about the x-axis. Then it can be proved that the volume of the solid generated is given by the
formula .
s.
v 2 ydA 2 ydxdy
A
bu
r 2 sin d dr
la
In polar form this formula becomes v
A
yl
1) Find the volume generated by the revolution of the cardioids r =a (1+cosθ) about the
lls
intial line.
.a
Sol: The given cardioids is symmetrical about the initial line θ=0.therfore the volume generated
by revolving the upper part of the curve about the initial line is same as the volume
w
generated by revolving the whole the curve .for the upper part of the curve θ varies form 0
w
a (1 cos )
v 2 r 2 sin drd
0 r 0
a 1 cos
r3
2 sin d
0
3 0
3
2 a3
1 cos sin d
3 0
4
2 a3 1 cos 8 3
a
3 4 0
3
Page 91
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
b h x x, y
f x, y, z dv f x, y, z dxdydz x, y, z dz dy dx
v R a g x x, y
dv dxdydz dzdydx………………………..(1)
v R a g x x, y
The integral dv represents the volume V of the region R. thus expression (1)may be used to
m
v
co
compute V.
If(x,y,z) are changed to (u,v,w)we obtained the following expression for the volume,
jdudvdw…………………………(2)
s.
dv dxdydz
bu
v R R*
la
Taking (u,v,w)= (R,φ,z) in (2)
yl
v R
v R
w
coordinates.
w
PROBLEMS:
a a2 x2 a2 x2
v dzdydx
x ay a 2 2
x z a 2
x 2
Page 92
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
a2 x2
a a2 x2
z dydx
a a 2
x 2
a 2
x 2
2 2
a a x
2 a2 x 2 dydx
a a2 x2
a2 x2
a
a2 x 2 dy dx
2
a2 x2
a
a a2 x2
2 2
2 a x y dx
m
a a2 x2
a
co
2 a2 x2 2 a2 x 2 dx
a
s.
a a
2 2 2 x3
4 a x dx 4 a x
bu
a
3 a
3 3
a a
la
4 a3 a3
3 3
yl
2a 3 16a 3
lls
3
4 2a
3 4
.a
2) Find the volume bounded by the cylinder X2+Y2=4 and the planes y+z=3 and z=0
w
Soln: Here z varies from 0 to 3-y, y varies from () to () and x varies from -2 to 2
w
Required volume
w
2 4 x2 3 y
V dzdydx
x 2y 4 x z 0
2
3 y
2 4 x2 3 y 2 4 x2
dx dy dz dx dy z
2 4 x2 0 2 4 x2 0
2 4 x2 2 4 x2
y2
dx 3 y dy dx 3
2 4 x 2 2
2 4 x2
Page 93
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2
4 x2 4 x2
3 4 x2 3 4 x2 dx
2
2 2
2 2
x 4 x
6 4 x 2 dx 6 4 x2 sin 1
2
2 2 2 2
1 2 1 2
6 2sin 2sin 12 12
2 2 2 2
Curvilinear coordinates:
Introduction: the cartesian co-ordinate system is not always convenient to solve all sorts of
problems. Many a time we come across a problem having certain symmetries which decide the
choice of a co ordinates systems .our experience with the cylindrical and spherical polar co-
m
ordinates systems places us in a good position to analyse general co-ordinates systems or
curvilinear coordinates. Any suitable set of three curved surface can be used as reference surface
co
and their intersection as the reference axes. Such a system is called curvilinear system.
s.
Definition: bu
The position of a point P(x,y,z)in Cartesian co-ordinates system is determined by intersection of
three mutually perpendicular planes x=k1, y=k2, and z=k 3 where ki (i=1,2,3)
la
Are constants in curvilinear system, the axes will in general be curved. Let us the denote the
yl
It should be noted that axis is the intersection of two surfaces u1= constant and u2=constant and
.a
so on.
w
Cartesian coordinates (x,y,z) are related to (u1,u2,u3) by the relations which can be expressed as
w
Equation (1) gives the transformation equation from 1 coordinates system to another.
The inverse transformation equation can be written as u1= u1 (x,y,z), u2= u2 (x,y,z), u3= u3
(x,y,z)……(2).
Page 94
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
r
co
So that where h1
u1
r
s.
h1eˆ1
u1
bu
Similarly if e2 ande3 and are unit tangent vector to the u and u curves at p respectively.
Than
la
r r
yl
u2 u2
e2
lls
r h2
u2
.a
r
So that h2 eˆ2
w
u2
w
r r
w
The quantities h1, h2 and h3 are called scale factors. The unit vectors are in the directions of
increasing u1, u2, and u3 respectively.
Page 95
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 x ˆ y ˆ r ˆ
i.e;; e1 i j k
h1 u1 u1 u1
1 x ˆ y ˆ r
e2 i j k
h2 u2 u2 u2
1 x ˆ y ˆ r ˆ
e3 i j k
h3 u3 u3 u3
r xiˆ yjˆ zk
1 x ˆ 1 x ˆ 1 x
e1.iˆ ; e2 .i ; e3 .i
h1 u1 h2 u2 h3 u3
m
1 y ˆ 1 y ˆ 1 y
e1. ˆj ; e2 . j ; e3 . j
h1 u1 h2 u2 h3 u3
co
1 z ˆ 1 z ˆ 1 z
e1.kˆ ; e2 .k ; e3 .k
s.
h1 u1 h2 u2 h3 u3
bu
Elementary arc length:
la
Let r r u1 , u2 , u3
yl
r r r
lls
i.e; dr eˆ1h1du1 eˆ2 h2 du2 eˆ3h3du3
w
r
On the curve u1 cure u2 and u3 are constants du2 du3 0 ds h1 du1 du1
u1
Page 96
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Volume of the parallelepiped =dv=h1 h2 h3 du1 du2 du3 is called the volume element.
( x, y , z )
v= du1 du2 du3
(u1 , u 2 , u 3 )
xyz
= j du1 du2 du3
m
u1u 2 u 3
co
Jacobian is positive since each h1, h2, h3 of are positive.
s.
2
bu
Expression for , divF , curlFand in orthogonal curvilinear coordinates:
la
Suppose the transformations from Cartesian coordinates x,y,z to curvilinear coordinates u1 , u2 , u3
yl
be x=f( u1 , u2 , u3 ), y=g( u1 , u2 , u3 ), z=h( u1 , u2 , u3 ) where f,g,h are single valued function with
lls
continuous first partial derivatives in some given region. The condition for the function f,g,h to
be independent is if the jacobian
.a
x x x
w
u1 u2 u3
w
( x, y , z ) y y y
= 0
(u1 , u2 , u3 ) u1 u2 u3
w
z z z
u1 u2 u3
When this condition is satisfied, u1 , u2 , u3 can be solved as single valued functions odf x, y and z
with continuous partial derivatives of the first order.
Let p be a point with position vector op xiˆ yˆj zk in the Cartesian form. The change of
r r r
coordinates to u1 , u 2 , u 3 makes r a function of u1 , u 2 , u 3 . The vectors , , are along
u1 u2 u3
Page 97
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
tangent to coordinate curves u1 c 2 , u 2 c 2 , u 3 c3 . Let ê1 , ê2 , ê3 denote unit vector along
r r r
these tangents. Then = eˆ1h1 , eˆ2 h2 eˆ3 h3
u1 u2 u3
r r r
Where h1= , h2= , h3=
u1 u2 u3
If ê1 , ê2 , ê3 are such that ê1 . ê2 =0, ê2 . ê3 =0, ê3 . ê1 =0
Then the curvilinear coordinates will be orthogonal and ê1 = ê2 x ê3 , ê2 = ê3 x ê1 ê3 = ê1 x ê2
r r r
Now r = r ( u1 , u 2 , u 3 ) dr du1 + du2 du3
u1 u2 u3
m
co
Gradient in orthogonal curvilinear coordinates:
s.
Let Φ(x,y,z) be a scalar point function in orthogonal curvilinear coordinates.
bu
letgrad eˆ1 eˆ2 eˆ3 where 1 2, 3 are functions of u1 , u 2 , u3
la
dr eˆ1 h1 du1
eˆ2 h2 du2 eˆ3 h3 du3
lls
alsod grad dr ˆ
1 e1
ˆ
2 e2 eˆ
3 3 eˆ1 h1 du1 eˆ2 h2 du2 eˆ3 h3 du3
.a
i.e, d h du1
1 1 h du2
2 2 3 3 h du3 .......... .......( 2)
w
1 1 1
w
Page 98
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
f f ( f1eˆ1 ) ( f 2 eˆ2 ) ( f 3 eˆ3 )
consider, , ( f1eˆ1 ) ( f 1eˆ2 eˆ3 ) ( f1 h2 h3 u 2 u 3 )(u sin g ....( 4))
co
( f1eˆ1 ) ( f 1 h2 h3 ) ( u 2 u3 ) f1 h2 h3 ( u2 u3 )
s.
(u sin g ) ( A) A A) bu
also.. u2 0 . u 3 sin ce..curl..grad 0
eˆ eˆ
.( f1eˆ1 ) ( f 1 h2 h3 ).( u 2 u3 ) ( f1 h2 h3 ). 2 3 from(5)
la
h2 h3
eˆ1
yl
( f 1 h2 h3 ).
h2 h3
lls
1
( f1eˆ1 ) ( f 1 h2 h3 )
.a
h1 h2 h3 u1
1
( f1eˆ2 )
w
similarly ( f 2 h3 h1 )
h1 h2 h3 u 2
w
1
( f1eˆ3 ) ( f 3 h1 h2 )
w
h1 h2 h3 u 3
1
f ( f1 h2 h3 ) ( f 2 h3 h1 ) ( f 3 h1 h2 )
h1 h2 h3 u1 u2 u3
Expression for curlF in orthogonal curvilinear coordinates
Let F (u1 , u 2 , u 3 ) be a vector point function such that f f1eˆ1 f 2 eˆ2 f 3 eˆ3
curlF = curl( f1eˆ1 ) + curl( f 2 eˆ2 ) + curl ( f 3 eˆ3 )
Page 99
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
= grad f1 h1 u1
1 1 1 eˆ
[ ( f1 h1 )eˆ1 ( f1 h1 )eˆ2 ( f1 h1 )eˆ3 ] 1 u sin g (3)and (5)
h1 u1 h2 u 2 h3 u 3 h1
1
{ ( f1 h1 )}eˆ2 h2 { ( f1 h1 )}eˆ3 h3
h1 h2 h3 u3 u2
similarly
1
curl ( f 2 eˆ2 ) { ( f 2 h2 )}eˆ3 h3 { ( f 2 h2 )}eˆ1 h1
h1 h2 h3 u1 u3
1
curl ( f 3 eˆ3 ) { ( f 3 h3 )}eˆ1 h1 { ( f 3 h3 )}eˆ2 h2
h1 h2 h3 u2 u1
m
( f 3 h3 ) ( f 2 h2 ) eˆ1 h ( f1 h1 ) ( f 3 h3 ) eˆ2 h2
co
1 u2 u3 u3 u1
curlf
h1 h2 h3
( f1 h1 ) eˆ3 h3
s.
( f 2 h2 )
u1 u2 bu
eˆ1 h1 eˆ2 h2 eˆ3 h3
la
1
Thus curlf is the expression for curlf in orthogonal curvilinear
h1 h2 h3 u1 u2 u3
yl
f1 h1 f 2 h2 f 3 h3
lls
coordinates.
.a
2
Expression for in orthogonal curvilinear coordinates
w
We know
w
1 1 1
eˆ1 eˆ2 eˆ3
h1 u1 h2 u 2 h3 u 3
1 1 1
2
eˆ1 eˆ2 eˆ3
h1 u1 h2 u 2 h3 u 3
2 1 h2 h3 h1 h3 h2 h1
h1 h2 h3 u1 h1 u1 u2 h2 u 2 u3 h3 u 3
2
This is the expression for in orthogonal curvilinear coordinates.
Page 100
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Definitions
m
co
s.
bu
la
Page 101
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 102
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 103
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 104
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 105
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 106
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 107
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1.
Sol:
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 108
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2.
m
co
s.
bu
la
yl
lls
w
dA1 dxdy, dA2 dydz , dA3 dzdx the elementary volume element is given by dv dxdydz
Page 109
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
r r r
The scalar factors are given by h1 1, h2 1, h3 1,
z
The elementary arc length is given by (ds)2 h12 (du1 )2 h 22 (du2 )2 h32 (du3 )2
i.e; (ds ) 2 (d ) 2 2
(d ) 2 ( dz ) 2
m
unit vectors at P in the direction of the tangents to , and z curves respectively, then we have
co
r r r
h1eˆ , h2eˆ , h3eˆz
z
s.
For cylindrical coordinate system h1 1, h2
bu
, h3 1
r 1 r r
la
eˆ , eˆ , eˆz eˆ cos iˆ sin ˆj; eˆ sin iˆ cos ˆj; eˆz kˆ
z
yl
.a
Hence the unit vectors e , e , ez are mutually perpendicular, which shows that the cylindrical
polar coordinate system is orthogonal curvilinear coordinate system.
w
w
In this case u1 r , u2 , u3 . Also x r sin cos , y r sin cos , z r cos . In this system
unit vectors e1 , e2 , e3 are denoted by e , e , ez respectively. These unit vectots are extended
respectively in the directions of r increasing, increasing and increasing.
Let r be the position vector of the point P. Then
r (r sin cos )iˆ (r sin sin ) ˆj (r cos )kˆ
r ˆ r
ˆ ˆ
sin cos i sin sin j cos k ; r cos cos iˆ r cos sin ˆj r sin kˆ
r
Page 110
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
r
r sin sin iˆ r sin cos ˆj
r r r
The scalar factors are h1 1, h2 r , h3 r sin
r z
The elementary arc length is given by (ds)2 h12 (du1 )2 h 22 (du2 )2 h32 (du3 )2
Show that the spherical coordinate system is orthogonal curvilinear coordinate system and
m
also prove that (er , e , e ) form a right handed basis.
co
Proof: We have for spherical Polar coordinate system
s.
r (r sin cos )iˆ (r sin sin ) ˆj (r cos )kˆ bu
Let er , e , e be the base vectors at P in the directions of the tangents to r, , curves respectively
la
then we have
yl
r r r
h1eˆ1 ; 1, h2eˆ2 r , h3eˆ3
lls
r
r r r
i.e h1eˆr ; 1, h2eˆ r , h3eˆ
.a
r
w
w
We know that for spherical polar coordinate the scalar factors h1 1, h2 2, h3 r sin
w
r
eˆr sin cos iˆ sin sin ˆj cos kˆ reˆ
r
r cos cos iˆ r cos sin ˆj r sin kˆ
r
r sin eˆ r sin sin iˆ r sin cos ˆj
Page 111
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
This shows that eˆr , eˆ and eˆ are mutually perpendicular. Hence spherical polar coordinates are
also orthogonal curvilinear coordinates.
iˆ ˆj kˆ
Further eˆr eˆ sin cos sin sin cos sin iˆ cos ˆj eˆ
cos cos cos sin sin
Similarly we can show that eˆ eˆ eˆr and eˆ eˆr eˆ which shows that (er , e , e ) form a right
handed basis.
m
We can use from matrix algebra, if Y=AX then X=A-1Y provided A is non singular.
co
1) Cylindrical polar coordinates (eρ, eφ, ez)
s.
bu
We have for cylindrical coordinate system
la
e cos sin 0 i
.a
e 0 0 1 k
w
i cos sin 0 e
w
i
sin cos cos cos sin er
we get j = sin sin cos sin cos e ………..(b)
k cos sin 0 e
m
i sin cos cos cos sin
j = sin sin cos sin cos
co
k cos sin 0
s.
This gives the transformation of (i,j,k) in terms of the base vectors (er eθ,,eφ).
bu
3) Relation between cylindrical and spherical coordinates
la
Now from (a) and (b)
yl
e sin cos 0 er
e = 0 0 1 e
e cos sin 0 e
Page 113
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
sin 0 cos e
= cos 0 sin e
0 1 0 e
These two results give us the relation between cylindrical and spherical coordinates bases and
vice versa.
PROBLEMS:
m
1. Express vector f=2yi-zj+3xk in cylindrical coordinates and find fρ, fφ fz.
co
Sol:The relation between the Cartesian and cylindrical coordinates given by
X=ρcosφ,y=ρsinφ,z=z
s.
i=cosφeρ-sinφeρ ; j=sineρ+ cosφeφ, k=ez.
bu
We have f=2yi-zj+3xk
la
Therefore
w
2) Express the vector f=zi-2xj+yk in terms of spherical polar coordinates and find fr, fθ, fφ,
Page 114
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 115
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
MODULE – 5
LAPLACE TRANSFORM
INTRODUCTION
Laplace transform is an integral transform employed in solving physical problems.
Many physical problems when analysed assumes the form of a differential equation
subjected to a set of initial conditions or boundary conditions.
By initial conditions we mean that the conditions on the dependent variable are specified
m
at a single value of the independent variable.
co
If the conditions of the dependent variable are specified at two different values of the
s.
independent variable, the conditions are called boundary conditions.
bu
The problem with initial conditions is referred to as the Initial value problem.
la
The problem with boundary conditions is referred to as the Boundary value problem.
yl
d2y dy
Example 1: The problem of solving the equation y x with conditions y(0) = y
lls
dx 2 dx
.a
d2y dy
Example 2: The problem of solving the equation 3 2 2 y cos x with y(1)=1,
w
dx dx
w
Laplace transform is essentially employed to solve initial value problems. This technique
is of great utility in applications dealing with mechanical systems and electric circuits.
Besides the technique may also be employed to find certain integral values also. The
transform is named after the French Mathematician P.S. de’ Laplace (1749 – 1827).
Page 116
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
LAPLACE TRANSFORMS
Definition:
m
Let f(t) be a real-valued function defined for all t 0 and s be a parameter, real or
co
st
complex. Suppose the integral e f (t )dt exists (converges). Then this integral is called the
0
s.
Laplace transform of f(t) and is denoted by L[f(t)]. bu
st
Thus, L[f(t)] = e f (t )dt (1)
la
0
yl
We note that the value of the integral on the right hand side of (1) depends on s. Hence
lls
Consider relation (2). Here f(t) is called the Inverse Laplace transform of F(s) and is
denoted by L-1 [F(s)].
w
Thus, (3)
f1(t), 0<t<a
f3(t), t>b
Note that f(t) is piecewise continuous. The Laplace transform of f(t) is defined as
Page 117
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
st
L[f(t)] = e f (t )
0
a b
st st st
= e f1 (t )dt e f 2 (t )dt e f 3 (t )dt
0 a b
NOTE: In a practical situation, the variable t represents the time and s represents frequency.
Hence the Laplace transform converts the time domain into the frequency domain.
Basic properties
The following are some basic properties of Laplace transforms:
m
1. Linearity property: For any two functions f(t) and (t) (whose Laplace transforms exist)
co
and any two constants a and b, we have
s.
L [a f(t) + b (t)] = a L[f(t)] + b L[ (t)]
bu
Proof :- By definition, we have
la
yl
st st st
L [af (t) + b (t)] = e af (t ) b (t ) dt = a e f (t )dt b e (t )dt
0 0 0
lls
= a L[f(t)] + b L[ (t)]
.a
1 s
2. Change of scale property: If L L[f(t)] = F(s), then L[f(at)] = F , where a is a
a a
positive constant.
st
L[f(at)] = e f (at )dt (1)
0
Page 118
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
s
1 a
x
L f(at) = e f ( x)dx
a0
1 s
F
a a
m
Proof :- By definition, we have
co
st
L [eatf (t)] = e e at f (t ) dt
s.
0
bu
(s a)
la
= e f (t )dt
yl
0
lls
= F(s-a)
This is the desired property. Here we note that the Laplace transform of eat f(t) can be written
.a
L[(eat)] = e st e at dt e ( s a )t
dt
0 0
e ( s a )t 1
= , s>a
( s a) 0 s a
Thus,
Page 119
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
L[(eat)] =
s a
1
L(1) = , s>0
s
1 1
L1 e at L 1 e at
s a s
e at at
1
m
e st
2. L(cosh at) = L = e e at e at
dt
2 20
co
1
s.
( s a )t ( s a )t
= e e dt
20
bu
Let s > |a|. Then,
la
yl
1 e ( s a )t e ( s a )t
L(cosh at )
lls
s
2 ( s a) ( s a) 0 = 2
s a2
.a
w
s
w
and so
s
L1 2
cosh at
s a2
e at e at
a
3. L (sinh at) = L , s > |a|
2 s2 a2
Thus,
Page 120
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
a
L (sinh at) = , s > |a|
s2 a2
and so,
1 sinh at
L1 2 2
s a a
st
4. L (sin at) = e sin at dt
0
Here we suppose that s > 0 and then integrate by using the formula
m
e ax
co
e ax sin bxdx a sin bx b cos bx
a2 b2
s.
Thus,
bu
a
L (sinh at) = , s>0
la
2
s a2
yl
and so
lls
1 sinh at
L1
.a
2 2
s a a
w
w
w
st
5. L (cos at) = e cos atdt
0
e ax
e ax cosbxdx a cos bx b sin bx
a2 b2
s
Thus, L (cos at) = , s>0
s2 a2
Page 121
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
s
and so L1 cos at
s2 a2
L(tn) = e st t n dt
0
m
n x x dx 1
L(t ) e n 1
e x x n dx
s s s
co
0 0
s.
The integral e x x n dx is called gamma function of (n+1) denoted by (n 1) . Thus
bu
0
(n 1)
la
L(t n )
sn 1
yl
lls
n!
L(t n )
sn 1
w
w
and so
w
1 1 tn tn
L or as the case may be
sn 1
(n 1) n!
Application of shifting property:-
The shifting property is
Page 122
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
s s a
1. L(e at cosh bt ) L(coshbt ) s s a 2 2 =
s b s s a ( s a) 2 b 2
Thus,
s a
L(eatcoshbt) =
( s a) 2 b 2
and
s a
L1 e at cosh bt
( s a) 2 b 2
m
a
co
at
2. L(e sinh bt )
( s a) 2 b2
s.
and
bu
1
L1 e at sinh bt
la
( s a) 2 b 2
yl
lls
s a
.a
and
w
s a
L1 e at cos bt
( s a) 2 b 2
b
4. L(e at sin bt)
( s a) 2 b2
and
Page 123
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 1 e at sin bt
L
( s a) 2 b2 b
(n 1) n!
5. L(e at t n ) or as the case may be
( s a) n 1 ( s a) n 1
Hence
1 1 e at t n n!
L or as the case may be
( s a) n 1
(n 1) ( s a) n 1
m
Examples :-
co
1. Find L[f(t)] given f(t) = t, 0 < t < 3
s.
4, t>3 bu
Here
3
la
st
L[f(t)]= e f (t )dt e st tdt 4e st dt
yl
0 0 3
lls
1 3s 1 3s
L[f(t)] = e (1 e )
w
s s2
w
0, t>
Here
st st st
L[f(t)] = e f (t )dt e f (t )dt = e sin 2tdt
0 0
Page 124
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
st
e 2
= 2
s sin 2t 2 cos 2t = 1 e s
s 4 0 s 2
4
1
(i) Here L(sin3t sin4t) = L [ (cost cos 7t )]
2
m
1
= L(cost ) L(cos7t ) , by using linearity property
co
2
s.
1 s s 24s
= 2 2 2
2 s 1 s 49 (s 1)(s 2 49)
bu
(ii) Here
la
1 1 1 s
yl
L(cos24t) = L (1 cos8t ) 2
2 2 s s 64
lls
.a
(iii) We have
w
1
w
1
sin3 2t 3sin 2t sin 6t
4
so that
1 6 6 48
L(sin3 2t ) 2 2 2
4 s 4 s 36 (s 4)(s 2 36)
Page 125
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
Here cos2t cos3t = [cos 5t cost ]
2
so that
1
cost cos2t cos3t = [cos 5t cos t cos 2 t ]
2
1
= [cos 6t cos 4t 1 cos 2 t ]
4
1 s s 1 s
Thus L(cost cos2t cos3t) =
m
2 2 2
4 s 36 s 16 s s 4
co
L(cosh22t)
s.
5. Find bu
We have
1 cosh 2
la
cosh2
2
yl
1 cosh 4t
.a
cosh2 2t
2
w
Thus,
w
1 1 s
w
L(cosh2 2t ) 2
2 s s 16
1
6. Evaluate (i) L( t ) (ii) L (iii) L(t-3/2)
t
(n 1)
We have L(tn) =
sn 1
1
(i) For n= , we get
2
Page 126
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
( 1)
L(t1/2) = 2
s3 / 2
1 1 1
Since (n 1) n (n) , we have 1
2 2 2 2
Thus, L( t ) 3
2
2s
1
(ii) For n = - , we get
2
m
1
1 2
co
2
L(t ) 1
s 2 s
s.
bu
3
(iii) For n = - , we get
2
la
1
yl
3 2 2
lls
2
L(t ) 1 1
2 s
2 2
s s
.a
We have,
w
n!
L (tn) =
w
sn 1
2! 2
L (t2) =
s3 s3
3! 6
L (t3) =
s4 s4
Page 127
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Given =
s 3 15
=2 , by using shifting property
(s 3) 25 ( s 3)2 25
2
2s 9
= 2
, on simplification
s 6s 34
m
eat e at
co
s.
1 a a
=
2 ( s a) 2 a 2
( s a) 2 a2
bu
la
a ( s 2 2a 2 )
, on simplification
yl
[( s a ) 2 a 2 ][( s a ) 2 a2 ]
lls
Given
w
w
et e t
3 sin 2t sin 6t
L
w
2 4
1
= 3 L et sin 2t L(et sin 6t ) 3L(e t sin 2t ) L(e t sin 6t )
8
1 6 6 6 6
=
8 ( s 1)2 4 ( s 1) 2 36 ( s 1) 2 4 ( s 1)2 36
3 1 1 1 1
=
4 ( s 1)2 4 ( s 1)2 36 ( s 1)2 24 ( s 1)2 36
Page 128
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
5
4t 2
11. Find L(e t )
We have
(n 1)
L(tn) = Put n= -5/2. Hence
sn 1
( 3 / 2) 4
L(t-5/2) = Change s to s+4.
s 3/ 2 3s 3/ 2
4
Therefore, L (e 4 t t 5/ 2
) 3/ 2
3( s 4)
m
Transform of tn f(t)
co
Here we suppose that n is a positive integer. By definition, we have
s.
e st
f (t )dt
bu
F(s) =
0
la
Differentiating ‘n’ times on both sides w.r.t. s, we get
yl
dn n
lls
st
F ( s) e f (t )dt
ds n s n
0
.a
dn
F ( s) ( t)n e st
f (t )dt
ds n
w
dn
( 1) n F ( s) (t n f (t )e st dt L[t n f (t )] , by definition
ds n 0
Thus,
n dn
L [t f(t)]= ( 1) n F ( s)
ds n
Also,
1 dn
L F ( s) ( 1) n t n f (t )
ds n
In particular, we have
d
L[t f(t)] = F (s) , for n=1
ds
2 d2
L [t f(t)]= 2 F (s) , for n=2, etc.
ds
m
d
Also, L1 F ( s) tf (t ) and
ds
co
d2
s.
L1 F ( s) t 2 f (t )
ds 2
bu
la
f(t)
Transform of
yl
t
lls
st
e f (t )dt
.a
We have, F(s) =
0
w
Therefore,
w
w
F ( s)ds e st f (t )dt ds
s s 0
= f (t ) e st ds dt
0 s
st
e
= f (t ) dt
0 t s
Page 130
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
st f (t ) f (t )
= e dt L
0
t t
f (t )
Thus, L F ( s)ds
t s
f (t )
This is the transform of
t
1 f (t )
Also, L F ( s)ds
s
t
m
Examples :
co
1. Find L [te-t sin4t]
s.
4
t
bu
We have, L[e sin 4t ]
(s 1) 2 16
la
So that,
yl
d 1
L [te-t sin4t] = 4
lls
ds s 2 2s 17
.a
8( s 1)
=
w
( s 2s 17) 2
2
w
3
We have L (sin3t) = 2
s 9
So that,
2 d2 3
L (t sin3t) = 2 2
ds s 9
d s
= 6
ds ( s 9) 2
2
18( s 2 3)
=
( s 2 9) 3
Page 131
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
e t sin t
3. Find L
t
We have
1
L(e t sin t )
( s 1) 2 1
e t sin t ds
Hence L = 2
tan 1 ( s 1) s
t 0 ( s 1) 1
m
sin t sin at
co
4. Find L . Using this, evaluate L
t t
s.
1
We have L (sint) =
s 2
1
bu
la
sin t ds
So that L [f (t)] = L = 2
tan 1 s S
t s 1
yl
s
lls
= tan 1 s cot 1 s F ( s)
2
.a
Consider
w
w
sin at sin at
L =aL aLf (at)
t at
w
1 s
=a F , in view of the change of scale property
a a
1 s
= cot
a
Page 132
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
cos at cosbt
5. Find L
t
s s
We have L [cosat – cosbt] = 2 2 2
s a s b2
cos at cosbt s s
So that L = 2 2 2
ds
t s s a s b2
1 s2 a2
= log 2
2 s b2 s
m
1 s2 a2 s2 a2
= Lt log 2 log
co
2 s s b2 s2 b2
s.
1 s2 b2
= 0 log 2
2 s a2
bu
la
1 s2 b2
= log 2
2 s a2
yl
lls
3
6. Prove that e 3t t sintdt
50
.a
0
w
We have
w
d d 1
e st t sin tdt L(t sin t )
w
= L(sint ) = 2
0 ds ds s 1
2s
=
( s 1) 22
3
e 3t t sintdt
0
50
Page 133
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Consider
st
L f (t ) = e f (t )dt
0
st st
= e f (t ) 0 ( s)e f (t )dt , by using integration by parts
0
st
= Lt (e f (t ) f (0) sLf (t )
t
= 0 - f (0) + s L[f(t)]
m
Thus
co
L f (t ) = s L[f(t)] – f(0)
s.
Similarly, bu
L f (t ) = s2 L[f(t)] – s f(0) - f (0)
la
yl
In general, we have
lls
t
w
Transform of f(t)dt
0
w
w
st
Now, L (t) = e (t )dt
0
st st
e e
= (t ) (t ) dt
s 0 0
s
Page 134
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
= (0 0) f (t )e st dt
s0
t
1
Thus, L f (t )dt L[ f (t )]
0
s
t
1 1
Also, L L[ f (t )] f (t )dt
s 0
Examples:
m
1. By using the Laplace transform of sinat, find the Laplace transforms of cosat.
co
a
Let f(t) = sin at, then Lf(t) = 2
s a2
s.
We note that
bu
f (t ) a cos at
la
1 1
.a
1 sa
= 0
w
a s a2
2
w
Thus
s
L(cosat) = 2
s a2
t 1 1 1
2. Given L 2 3/ 2
, show that L
s t s
t 1
Let f(t) = 2 , given L[f(t)] = 3/ 2
s
Page 135
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2 1 1
We note that, f (t )
2 t t
1
Lf (t ) L
t
Hence
1
L Lf (t ) sLf (t ) f (0)
t
m
1
= s 3/ 2
0
s
co
1 1
s.
Thus L
t s
bu
This is the result as required.
la
yl
t
cos at cos bt
lls
3. Find L dt
0
t
.a
cos at cos bt 1 s2 b2
w
t
w
1
Using the result L f (t )dt Lf (t )
0
s
t
cos at cos bt 2 2
We get, L dt = 1 log s 2 b 2
0
t 2s s a
t
t
4. Find L te sin 4tdt
0
Page 136
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
8( s 1)
Here L te t sin 4t
( s 2 2s 17) 2
t
8( s 1)
Thus L te t sin 4tdt =
0 s( s 2s 17) 2
2
m
Proof :By definition, we have
co
st su
e f (t )dt e f (u )du
s.
L f (t) =
0 0
bu
T 2T ( n 1)T
su su su
= e f (u )du e f (u )du ....... e f (u )du ....
la
0 T nT
yl
( n 1)T
lls
su
= e f (u )du
n 0 nT
.a
T
s ( t nT )
L f(t) = e f (t nT )dt
w
n 0t 0
Here
Hence
T
sT n st
Lf (t ) (e ) e f (t )dt
n 0 0
Page 137
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
T
1 st
= ST
e f (t )dt , identifying the above series as a geometric series.
1 e 0
T
1 st
Thus L[ f(t)] = sT
e f (t )dt
1 e 0
Examples:-
1. For the periodic function f(t) of period 4, defined by f(t) = 3t, 0 < t < 2
6, 2 < t < 4
m
find L [f(t)]
co
Here, period of f(t) = T = 4
s.
We have, bu
T
1 st
L f(t) = sT
e f (t )dt
la
1 e 0
yl
4
1 st
lls
= 4s
e f (t )dt
1 e 0
.a
2 4
1
w
st
= 4s
3te dt 6e st dt
1 e
w
0 2
w
2 4
st 2 st st
1 e e e
= 4s
3 t 1. dt 6
1 e s 0 0
s s 2
2s 4s
1 31 e 2se
= 4s
1 e s2
Thus,
3(1 e 2 s 2se 4s
)
L[f(t)] =
s 2 (1 e 4 s )
Page 138
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
2
3. A periodic function of period is defined by
2
0, t
E
where E and are positive constants. Show that L f(t) =
(s 2 w 2 )(1 e s/w
)
2
Sol: Here T = . Therefore
m
co
2 /
1 st
L f(t) = s(2 / )
e f (t )dt
1 e
s.
0
bu
/
1 st
= s(2 / )
Ee sin tdt
1 e
la
0
yl
st /
E e
lls
= s(2 / )
s sin t cos t
1 e s2 2
0
.a
E (e s /
1)
w
= s(2 / )
1 e s2 2
w
w
E (1 e s / )
= s /
(1 e )(1 e s / )(s 2 2
)
E
= s /
(1 e )(s 2 2
)
f (t)= E, 0 t a
Page 139
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
-E, a < t 2a
E as
show that L [f (t)] = tanh
s 2
2a
1 st
Sol: Here T = 2a. Therefore L [f (t)] = 2 as
e f (t )dt
1 e 0
a 2a
1 st
= 2 as
Ee dt Ee st dt
1 e 0 a
m
sa 2 as as
= 2 as
1 e (e e )
s(1 e )
co
E as 2
1 e
s.
2 as
s(1 e ) bu
E (1 e as ) 2
la
s (1 e as )(1 e as
)
yl
E e as / 2 as / 2
lls
e
=
s e as / 2 e as / 2
.a
w
E as
tanh
s 2
w
w
Step Function:
In many Engineering applications, we deal with an important discontinuous function H
(t-a) defined as follows:
0, t a
H (t-a) = 1, t > a
Page 140
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
This function is known as the unit step function or the Heaviside function. The function is
named after the British electrical engineer Oliver Heaviside.The function is also denoted by
u (t-a). The graph of the function is shown below:
H (t-1)
Note that the value of the function suddenly jumps from value zero to the value 1 as t a
from the left and retains the value 1 for all t>a. Hence the function H (t-a) is called the unit step
function.
0, t 0
m
H(t) = 1, t>0
co
s.
Transform of step function bu
By definition, we have L [H(t-a)] = e st H (t a)dt
0
la
a
yl
= e st 0dt e st (1)dt
lls
0 a
.a
as
e
=
w
s
w
1
In particular, we have L H(t) =
w
as
e 1
Also, L1 H (t a) and L 1 H (t )
s s
Proof: - We have
Page 141
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
st
= e f (t a)dt
a
s(a u)
L[f(t-a) H(t-a)] = e f (u )du
0
= e-as L [f(t)]
m
co
This is the desired shift theorem.
s.
Examples:
bu
1. Find L [et-2 + sin(t-2)] H(t-2)
la
Sol: Let
lls
1 1
so that L f(t) = 2
s 1 s 1
w
w
Thus,
1 1
L[e (t 2)
sin(t 2)]H (t 2) e 2s
s 1 s2 1
so that
Hence
6 8 8
L[ f (t )]
s3 s2 s
Thus
= e-s L [f(t)]
m
s 6 8 8
=e
co
s3 s2 s
s.
bu
3. Find Le-t H (t-2)
e2
Thus, L [f(t)] =
lls
s 1
2 ( s 1)
2s e
L[ f (t 2) H (t 2)] e Lf (t )
w
s 1
w
Thus
2 ( s 1)
t e
L e H (t 2)
s 1
f1 (t), t a
Sol: Consider
Page 143
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
0, t a
= f2 (t), t > a
5. Express the following functions in terms of unit step function and hence find their
Laplace transforms.
m
1. f(t) = t2 , 1 < t 2
co
4t , t>2
s.
bu
Sol: Here, f(t) = t2 + (4t-t2) H(t-2)
2
la
Hence, L f(t) = L(4t t 2 ) H (t 2) (i)
s3
yl
Let (t-2) = 4t – t2
lls
2 4
w
Now, L[ (t )]
s3 s
w
2
L f(t) = L (t 2) H (t 2)
s3
2 2s
= e L (t )
s3
2 2s 4 2
= e
s3 s s3
Page 144
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
s
Hence, L[ f(t)] = 2
L(sin t cos t ) H (t ) (ii)
s 1
1 s
so that L[ (t)] = 2 2
m
s 1 s 1
co
s
Expression (ii) reads as L [f(t)] = 2
L (t ) H (t )
s 1
s
s.
bu
s
= 2
e L (t )
s 1
la
Page 145
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Solution:
Page 146
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
3.
m
co
s.
bu
4.
la
yl
lls
.a
w
w
w
Page 147
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
F(s) L 1 F (s)
m
f (t )
co
1 1
,s 0
s
s. e at
1
,s a
bu
s a
la
s Cos at
2 2
,s 0
s a
yl
lls
1 Sin at
2 2
,s 0
s a a
.a
1 Sin h at
w
2 2
,s a
s a a
w
s
w
2
,s a
s a2
Cos h at
1 tn
n 1
,s 0
s n!
n = 0, 1, 2, 3, . . .
1 tn
,s 0
sn 1
n 1
n > -1
Page 148
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Example
1. Find the inverse Laplace transforms of the following:
1 s b 2s 5 4s 9
(i) (ii) 2 2
(iii) 2
2s 5 s a 4s 25 9 s 2
Here
5t
1 1 1 1 1 1 2
(i) L L e
2s 5 2 s 5 2
2
s b s 1 b
(ii) L 1 L1 bL1 cos at sin at
s a2
2 2 2 2 2
m
s a s a a
co
s.
5 s 9
2s 5 4s 8 2 1 s 2 2
(iii) L 1 L 4L 1
bu
4s 25 9 s 2
2
4 25 s 2
9
s2
4
la
yl
1 5t 5t 3
cos sin 4 cos h3t sin h3t
lls
2 2 2 2
.a
w
Examples
3s 1
1. Evaluate : L 1 4
s 1
3 s 1 -1 1 1 1
Given L-1 4
3L1 3
2L1 4
s 1 s 1 s 1
Page 149
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
t 1 1
3e L1 2e t
L1
s3 s4
1 1 tn
L and taking n 2 and 3, we get
sn 1
n!
3e t t 2 e tt 3
Given
2 3
s 2
2. Evaluate : L-1 2
s - 2s 5
m
co
s 2 s 1 3
Given L-1 2
L1 2
s 1 4 s 1 4
s.
s 1 1
L1 3L 1
2
bu 2
s 1 4 s 1 4
la
s 1
e t L-1 3 et L 1
s2 s2
yl
4 4
lls
3 t
et cos 2t e sin 2t
2
.a
2s 1
w
3.Evaluate :L 1 2
s 3s 1
w
w
3 3
-1
s 1 1
s 1
Given 2L 2
2 L 2
L1
3 2 5 3 2 5 3 2 5
s s s
2 4 2 4 2 4
3t 3t
s 1
2 e 2
L1 e 2
L1
s 2 5 s 2 5
4 4
3t
2
5 2 5
2e cos h t sin h t
2 5 2
Page 150
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
-1 2s 2 5s 4
4. Evaluate : L
s 3 s 2 2s
we have
2 s 2 5s 4 2 s 2 5s 4
s 3 s 2 2s s s2 s 2
2 s 2 5s 4
ss 2 s 1
A B C
s s 2 s 1
m
co
For s = 0, we get A = 2, for s = 1, we get C = 1 and for s = -2, we get B = -1. Using these values
in (1), we get
s.
2s 2 5s 4 2 1 1
bu
s 3 s 2 2s s s 2 s 1
la
Hence
yl
1 2 s 2 5s 4 2t
L 2 e et
lls
s 2 s 2 25
.a
w
4s 5
w
5. Evaluate : L 1 2
s 1 s 2
w
Let us take
4s 5 A B C
2 2
s 1 s 2 s 1 s 1 s 2
Page 151
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
4s 5 1 3 3
2 2
(1), we get s 1 s 2 s 1 s 1 s 2
1 4s 5 1 1 1
Hence L 2
e tL 1 3e t L 1 3e 2t
L1
s 1 s 2 s2 s s
t t 2t
te 3e 3e
1 s3
6. Evaluate :L
s4 a4
m
s3 A B Cs D
co
4 4
(1)
Let s a s a s a s2 a2
s.
Hence s3 = A(s + a) (s2 + a2) + B (s-a)(s2+a2)+(Cs + D) (s2 – a2) bu
For s = a, we get A = ¼; for s = -a, we get B = ¼; comparing the constant terms, we get
la
D = a(A-B) = 0; comparing the coefficients of s3, we get
yl
s3 1 1 1 1 s
.a
4 4
s a 4 s a s a 2 s a2
2
w
s3 1 at 1
w
1 at
L 4 4
e e cos at
s a 4 2
1
cos hat cos at
2
s
7. Evaluate :L 1 4
s s2 1
s s 1 2s
4
Consider
s s2 1 s 2
s 1 s 2
s 1 2 s s 1 s2 s 1
2
Page 152
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1 s2 s 1 s2 s 1
2 s2 s 1 s2 s 1
1 1 1
2 2
2 s s 1 s s 1
1 1 1
2 1 2 3 1 2 3
s 2
s 2
4 4
m
Therefore
co
s.
1 1
1 s 1 2t 1 1 t 1
L e L e 2
L1
4
bu
s s2 1 2 3 3
s2 s2
4 4
la
3 3
yl
1 sin t 1 sin t
1 2t 2 2
t
2
e e
lls
2 3 3
2 2
.a
w
2 3 t
sin t sin h
w
3 2 2
w
Page 153
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Examples
ss
1 e
(1) Evaluate: L 4
s 2
Here
1
a 5, F ( s ) 4
s 2
1 1 e 2t t 3
Therefore f (t ) L 1F ( s ) L1 4
e 2t L 1
s 2 s4 6
Thus
e 5s
m
L1 4
f (t a ) H (t a)
s 2
co
3
e2 t 5
t 5
H t 5
s.
6
bu
e s se 2 s
la
1
(2) Evaluate: L
s2 1 s2 4
yl
lls
Given f1 t H t f2 t 2 H t 2 (1)
.a
1
Here f1 (t ) L1 2
sin t
s 1
w
s
f 2 (t ) L1 cos 2t
w
2
s 4
w
Page 154
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
d
L1 F s tf (t )
Hence ds
Examples:
s a
(1) Evaluate: L 1 log
s b
s a
Let F ( s ) log log s a log s b
m
s b
co
d 1 1
s.
Then F s
ds s a s b
bu
d
la
So that L 1 F s e at
e bt
ds
yl
lls
bt at
or t f t e e
.a
w
bt at
e e
Thus f t
w
b
w
a
(2) Evaluate L 1 tan 1
1 a
Let F ( s) tan
s
d a
Then F s 2
ds s a2
Page 155
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
d
or L1 F s sin at so that
ds
or t f t sin at
sin at
f t
a
F s
Inverse transform of
m
s
co
1
(1) Evaluate: L 1
s.
s s2 a2 bu
1
Let us denote F s so that
s2 a2
la
sin at
f (t ) L 1F s
yl
a
lls
t
1 1 F s
1 sin at
Then L L dt
.a
2 2
ss a s 0
a
w
1 cos at
w
a2
w
Convolution Theorem:
Page 156
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
la
yl
lls
.a
w
w
w
Page 157
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
m
co
s.
bu
Using Convolution
la
yl
lls
inverse laplace
w
w
transforms
Page 158
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
(2) Evaluate : L 1 2 2
s s a
1
Solution : we have L 1 2
e at t
s a
t
1 1
Hence L 2
e at t dt
s s a 0
1 at
1 e 1 at , on integration by parts.
a2
m
co
Using this, we get
s.
t
1 1
L-1 at
1 e 1 at dt
bu
2
2
s s a a2 0
la
1
yl
at at
at 1 e 2 e 1
a3
lls
.a
t
1
L F ( s) G ( s) f (t ) g (t ) f t u g u du
0
This expressionis called the convolution theorem for inverse Laplace transform
Page 159
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
Examples
Employ convolution theorem to evaluate the following:
1
(1) L 1
s a s b
1 1
Sol:Let us denote F(s) , G ( s)
s a s b
m
Therefore, by convolution theorem,
co
t t
s.
1 a t u
-1
L e e bu
du e at
ea b u
du
s a s b 0
bu 0
at ea bt 1
e
la
a b
yl
bt at
e e
lls
a b
.a
s
( 2) L 1 2
w
s2 a2
w
1 s
Sol: Let us denote F(s) , G ( s) Then
w
2 2 2
s a s a2
sin at
f(t) , g(t) cos at
a
t
-1 s 1
L 2
sin a t u cos au du
s2 a2 0
a
Page 160
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
t
1 sin at sin at 2au
du, by using compound angle formula
a0 2
t
1 cos at 2au t sin at
u sin at
2a 2a 0
2a
s
(3) L 1
s 1 s2 1
Sol: Here
m
1 s
F(s) , G( s)
co
2
s 1 s 1
Therefore
s.
f(t) e t , g(t) sin t
bu
la
By convolution theorem, we have
yl
t
-1 1 t u te u
lls
et
w
t 1 t
e sin t cos t 1 e sin t cos t
2 2
w
w
problems. The solution of such a problem is obtained by using the Laplace Transform of the
The following are the expressions for the derivatives derived earlier.
Page 161
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
t
1. Solve by using Laplace transform method y y t e , y(o) 2
Sol: Taking the Laplace transform of the given equation, we get
1
sL y t yo Lyt 2
s 1
m
1
s 1L yt 2 2
co
s 1
so that
s.
2s 2 4s 3
Lyt 3
bu
s 1
Taking the inverse Laplace transform, we get
la
yl
1 2s 2 4s 3
Y t L
lls
3
s 1
.a
2
1 2s 1 1 4s 1 1 3
w
L 3
s 1
w
w
2 1
L1 3
s 1 s 1
1 t
e t2 4
2
Page 162
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
s 2 Ly (t ) sy(0) y (0) 2 s Ly (t ) y (0) 3 L y (t ) 2
s 1
1
L y (t ) s 2 2s 3
s2 1
or
1
L y (t )
s 1 s 3 s2 1
or
m
1
y (t ) L1
co
s 1 s 3 s2 1
s.
A B
Cs D
L1
bu
s 1 s 3 s2 1
la
s 1
yl
1 1 1 1 10 5
L1
lls
8 s 1 40 s 3 s2 1
.a
w
1 t 1 3t 1
e e cos t 2 sin t
8 40 10
t
f(t) l f u sin t u du
0
Page 163
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
t
1
L f(t) L f u sin t u du
s 0
1 1 L f (t )
L f(t) Lf (t ) L sin t
s s s2 1
Thus
s2 1 s2 1 t2
L f (t ) or f (t ) L1 1
m
s3 s3 2
co
This is the solution of the given integral equation.
s.
bu
d2x dx
(4) A particle is moving along a path satisfying, the equation 2
6 25 x 0 where
dt dt
la
x denotes the displacement of the particle at time t. If the initial position of the particle is at x 20
yl
and the initial speed is 10, find the displacement of the particle at any time t using Laplace transforms.
lls
L x (t) s 2 6s 25 20 s 130 0 or
20 s 130
L x (t)
s 2 6 s 25
so that
20 s 130
x(t) L1 2
s 3 16
Page 164
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
20 s 3 70
L1 2
s 3 16
s 3 1
20 L 1 2
70 L 1 2
s 3 16 s 3 16
3t
3t e sin 4t
20 e cos 4t 35
2
m
(5) A voltage Ee -at is applied at t
co
0 to a circuit of inductance L and resistance R. Show that the
Rt
E
s.
at L
current at any time t is e e
R - aL
bu
la
Sol: The circuit is an LR circuit. The differential equation with respect to the circuit is
yl
di
L Ri E (t )
lls
dt
.a
Here L denotes the inductance, i denotes current at any time t and E(t) denotes the E.M.F.
Thus, we have
di at
L Ri Ee or
dt
at
Li ' (t ) R i (t ) Ee
at
L LT i' (t) R LT i' (t) E LT e or
Page 165
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
1
L s LT i(t ) i(0) R LT i(t ) E
s a
E
Since i(o) o, we get LT i(t ) sL R or
s a
E
LT i (t )
s a sL R
E
Taking inverse transform L, we get i (t ) LT1
( s a )(sL R)
m
E 1 1
LT1 L LT1
co
R aL s a sL R
Thus
s.
Rt
bu
E at L
i (t ) e e
R aL
la
dx
(6) Solve the simultaneous equations for x and y in terms of t given 4y 0,
dt
.a
dy
9x 0 with x(o) 2, y(o) 1.
dt
w
w
s Lx(t) - x(o) 4 Ly (t ) 0
9 L x(t ) s Ly (t ) y (o ) 0
s L x(t ) 4 L y (t ) 2
9 L x(t ) 5 L y (t ) 1
s 18
L y(t)
s 2 36
Page 166
www.allsyllabus.com
vtu.allsyllabus.com www.allsyllabus.com
so that
s 18
y(t) L1 2 2
s 36 s 36
cos 6t 3 sin 6t (1)
dy
Using this in 9x 0, we get
dt
1
x(t) 6 sin 6t 18 cos 6t
9
or
m
2
co
x(t) 3 cos 6t sin 6t
3 (2)
s.
(1) and (2) together represents the solution of the given equation.
bu
la
yl
lls
.a
w
w
w
Page 167
www.allsyllabus.com