Mathematics I PDF
Mathematics I PDF
Mathematics-I
(R18A0021) Mathematics -I
UNIT I: Matrices
TEXT BOOKS:
i) Higher Engineering Mathematics by B V Ramana ., Tata McGraw Hill.
ii) Higher Engineering Mathematics by B.S. Grewal, Khanna Publishers.
iii) Advanced Engineering Mathematics by Kreyszig, John Wiley & Sons.
REFERENCE BOOKS:
i)Advanced Engineering Mathematics by R.K Jain & S R K Iyenger, Narosa Publishers.
ii)Advanced Engineering Mathematics by Michael Green Berg, Pearson Publishers .
iii)Engineering Mathematics by N.P Bali and Manish Goyal.
Course Outcomes: After learning the concepts of this paper the student will be able to
1.Analyze the solution of the system of linear equations and to find the Eigen values and
Eigen vectors of a matrix.
2.Find the extreme values of functions of two variables with / without constraints.
3.Solve first and higher order differential equations.
4.Solve first order linear and non-linear partial differential equations.
5.Solve differential equations with initial conditions using Laplac
UNIT-I
MATRICES
Matrix : A system of mn numbers real (or) complex arranged in the form of an ordered set
of ‘m’ rows, each row consisting of an ordered set of ‘n’ numbers between [ ] (or) ( ) (or) || ||
is called a matrix of order m xn.
a11 a12 .........a1n
a
21 a12 .........a 2 n
Eg: ......................... = [aij ]mxn where 1≤i≤m, 1≤j≤n.
.......................
a m1 a m 2 .......a mn
mxn
Order of the Matrix: The number of rows and columns represents the order of the matrix. It
is denoted by mxn, where m is number of rows and n is number of columns.
Types of Matrices:
Row Matrix: A Matrix having only one row is called a “Row Matrix”.
Eg: 1 2 31x3
Column Matrix: A Matrix having only one column is called a “Column Matrix” .
1
Eg: 1
2 3 x1
Null Matrix: A= aij mxn such that aij=0 i and j. Then A is called a “Zero Matrix”. It is
denoted by Omxn.
0 0 0
Eg: O2x3=
0 0 0
Rectangular Matrix: If A= aij mxn , and m n then the matrix A is called a “Rectangular
Matrix”
1 1 2
Eg :
3 4
is a 2x3 matrix
2
Square Matrix: If A= aij mxn and m = n then A is called a “Square Matrix”.
1 1
Eg : is a 2x2 matrix
2 2
Lower Triangular Matrix: A square Matrix AnXn aij nxn is said to be lower Triangular
of aij = 0 if i <j i.e. if all the elements above the principle diagonal are zeros.
4 0 0
Eg: 5 2 0 is a Lower triangular matrix
7 3 6
aij = 0 if i.>j. i.e. all the elements below the principle diagonal are zeros.
1 3 8
Eg: 0 4 5 is an Upper triangular matrix
0 0 2
Triangle Matrix: A square matrix which is either lower triangular or upper triangular is
called a triangle matrix.
Principal Diagonal of a Matrix: In a square matrix, the set of all aij, for which i = j are
called principal diagonal elements. The line joining the principal diagonal elements is called
principal diagonal.
Note: Principal diagonal exists only in a square matrix.
Diagonal elements in a matrix: A= [aij]nxn, the elements aij of A for which i = j.
i.e. a11, a22….ann are called the diagonal elements of A
1 2 3
Eg: A= 4 5 6 diagonal elements are 1, 5, 9
7 8 9
Diagonal Matrix: A Square Matrix is said to be diagonal matrix, if aij = 0 for i j i.e. all
3 0 0
Eg : A = diag (3,1,-2)= 0 1 0
0 0 2
Scalar Matrix: A diagonal matrix whose leading diagonal elements are equal is called a
2 0 0
“Scalar Matrix”. Eg : A= 0 2 0
0 0 2
Unit/Identity Matrix: If A = aij such that aij=1 for i = j, and aij=0 for i j then A is
nxn
Trace of Matrix: The sum of all the diagonal elements of a square matrix A is called Trace
of a matrix A, and is denoted by Trace A or tr A.
a h g
Eg : A = h b f then tr A = a+b+c
g f c
Singular & Non Singular Matrices: A square matrix A is said to be “Singular” if the
determinant of │A│= 0, Otherwise A is said to be “Non-singular”.
Note: 1. Only non-singular matrices possess inverse.
2. The product of non-singular matrices is also non-singular.
Inverse of a Matrix: Let A be a non-singular matrix of order n if there exist a matrix B
such that AB=BA=I then B is called the inverse of A and is denoted by A-1.
If inverse of a matrix exist, it is said to be invertible.
Note: 1. The necessary and sufficient condition for a square matrix to posses inverse is that
| |≠ .
2 .Every Invertible matrix has unique inverse.
3. If A, B are two invertible square matrices then AB is also invertible and
AB
1
B 1 A1
AdjA
4. A1 where detA 0 ,
det A
Theorem: The inverse of a Matrix if exists is Unique.
Note: 1. (A-1)-1 = A 2. I-1 = I
Theorem: If A, B are invertible matrices of the same order, then
(i). (AB)-1 = B-1A-1
(ii). (A1)-1 = (A-1)1
Sub Matrix: - A matrix obtained by deleting some of the rows or columns or both from the
given matrix is called a sub matrix of the given matrix.
1 5 6 7
Eg: Let A = 8 9 10 5 . Then 8 9 10 is a sub matrix of A obtained by deleting first
3 4 5
2 x3
3 4 5 1
2 1 1
3 1 2
Eg: A = be a 4x3 matrix
1 2 3
5 6 7
det C = 2(7-12)-1(21-10)+(18-5) = -9
Properties of trace of a matrix: Let A and B be two square matrices and be any scalar
1) tr ( A) = (tr A) ; 2) tr(A+B) = trA + trB ; 3) tr (AB) = tr(BA)
Idempotent Matrix: A square matrix A Such that A2=A then A is called “Idempotent
Matrix”.
Eg: A = [ ]
Eg: A = [ ]
Eg: A = [ ]
Transpose of a Matrix: The matrix obtained by interchanging rows and columns of the
given matrix A is called as transpose of the given matrix A. It is denoted by AT or A1
Eg: A = [ ] Then AT = [ ]
Properties of transpose of a matrix: If A and B are two matrices and AT , B T are their
transposes then
1) AT A ; 2) A B AT BT ; 3) KA KAT ; 4) AB BT AT
T T T T
a h g
Eg: h b f is a symmetric matrix
g f c
Note: All the principle diagonal elements of a skew symmetric matrix are always zero.
Since aij = -aij aij = 0
Theorem: Every square matrix can be expressed uniquely as the sum of symmetric and
skew symmetric matrices.
A A = A AT A AT =
1 1
Proof: Let A be a square matrix, A =
2 2
1
2
A AT A AT = P + Q, where P = A AT ; Q = A AT
1
2
1
2
1
2
Thus every square matrix can be expressed as a sum of two matrices.
T
Consider PT 1 A AT 1 A AT 1 AT AT = 1 A AT =P, since PT P ,
T T
2 2 2 2
P is symmetric
T
Consider QT A AT A AT AT AT = - 1 A AT = - Q
1 1 T 1 T
2 2 2 2
Since QT Q , Q is Skew-symmetric.
To prove the representation is unique: Let A= R+S 1 be the representation, where R is
Consider AT R S RT S T R S 2
T
1 2 A AT 2S S
1
2
A AT Q
Therefore every square matrix can be expressed as a sum of a symmetric and a skew
symmetric matrix
Ex. Express the given matrix A as a sum of a symmetric and skew symmetric matrices
2 4 9
where A= 14 7 13
9 5 11
2 14 3
Solution: A 4 7 5
T
9 3 11
4 10 12 2 5 6
A A 10 14 18 P A A 5 7 9 ; P is symmetric
T 1 T
2
12 18 22 6 9 11
0 18 6 0 9 3
A A 18 0 8 Q A A 9 0 4 ; Q is skew symmetric
T 1
2
6 8 0 3 4 0
2 5 6 0 9 3
Now A =P+Q = 5 7 9 + 9 0 4
6 9 11 3 4 0
Solved Problems :
cos sin
1. Show that A =
cos
is orthogonal.
sin
cos sin cos sin
Sol: Given A = AT =
cos cos
then
sin sin
cos sin cos sin
Consider A.AT =
sin cos sin cos
1 2 2 1 2 2
Sol: Given A = 2 1 2
1 Then AT = 1 2 1 2
3 3
2 2 1 2 2 1
1 2 2 1 2 2 9 0 0 1 0 0
Consider A .AT = 1 2 1 2 2 1 2 = 1 0 9 0 = 0 1 0 = I
9 9
2 2 1 2 2 1 0 0 9 0 0 1
⇒ A.AT = I
Similarly AT .A = I
Hence A is orthogonal matrix
0 2b c
3. Determine the values of a, b, c when a b c is orthogonal.
a b c
4b 2 c 2 2b 2 c 2 2b 2 c 2
2 1 0 0
2b c a2 b2 c2 a2 b2 c2 =I= 0 1 0
2
2b 2 c 2 a 2 b 2 c 2 a 2 b 2 c 2
0 0 1
2 3 1 2 4 3
Sol:- Given A= 4 3 1 3 3 1
3 1 9 1 1 9
2 3 1 2 4 3 14 0 0
4 3 1 3 3 1 = 0 26 0 I3
3 1 9 1 1 9 0 0 91
I3
Matrix is not orthogonal.
Complex matrix: A matrix whose elements are complex numbers is called a complex
matrix.
Conjugate of a complex matrix: A matrix obtained from A on replacing its elements by the
corresponding conjugate complex numbers is called conjugate of a complex matrix. It is
denoted by A
2 3i 5 2 3i 5
Eg: If A= then A =
6 7i 5 i 6 7i 5 i
Note: If A and B be the conjugate matrices of A and B respectively, then
(i) A = A
(ii) A B A B
3i 2 i
3i 2 i
Eg: Let A= then A =
3i
2 i
2 i
i
and A
T
2 i i
2 i i
A =-A A is skew-Hermitian matrix.
T
Note: 1. In Skew-Hermitian matrix the principal diagonal elements are either Zero or Purely
Imaginary.
2. The Skew- Hermitian matrix over the field of Real numbers is nothing but real
Skew - Symmetric matrix.
1 4 2 4i
6 2 4i 4
Eg: B
Theorem1: Every square matrix can be uniquely expressed as a sum of Hermitian and
skew – Hermitian Matrices.
1 1 1
A(2 A) ( A A) ( A A A A )
Proof: - Let A be a square matrix write 2 2 2
1 1
A ( A A ) ( A A )i.eA P Q
2 2
Let P
1
2
A A ; Q A A
1
2
1 1
Consider P ( A A ) ( A A ) ( A A ) P
2 2
Consider A ( R S ) R S R S . Ie A R A (2)
1 2 A A
1
2
A A P
2 R ie R
1 2 A A 2S ie S A A Q
1
2
Thus every square matrix can be uniquely expressed as a sum of Hermitian & skew
Hermitian matrices.
Solved Problems :
3 7 4i 2 5i
1) If A= 7 4i 2 3 i then show that A is Hermitian and iA is skew-
2 5i 3 i 4
Hermitian.
3 7 4i 2 5i
Sol:
Given A= 7 4i 2 3 i then
2 5i 3 i 4
3 7 4i 2 5i 3 7 4i 2 5i
A 7 4i 2
3 i And A 7 4i
T
2 3 i
2 5i 3 i 4 2 5i 3 i 4
T
A A Hence A is Hermitian matrix.
Let B= iA
3i 4 7i 5 2i
i.e B= 4 7i 2i 1 3i then
5 2i 1 3i 4i
3i 4 7i 5 2i
B 4 7i 2i 1 3i
5 2i 1 3i 4i
3i 4 7i 5 2i 3i 4 7i 5 2i
B 1 3i (1) 4 7i 2i 1 3i B
T
4 7i 2i
5 2i 1 3i 4i 5 2i 1 3i 4i
T
B =-B
Now AB BA AB BA
T T
T
AB BA
BA B A A B
T T T T T T
AB
BA AB (By (1))
AB BA
a ic b id
Sol: Given A=
b id a ic
a ic b id
Then A
b id a ic
a ic b id
T
Hence A A
a ic
b id
a ic b id a ic b id
AA
b id a ic b id a ic
a 2 b2 c 2 d 2 0
=
0 a 2 b2 c 2 d 2
0 1 2i
, show that I AI A is a unitary matrix.
1
4)Given that A=
1 2i 0
1 0 0 1 2i
Sol: we have I A
0 1 1 2i 0
1 1 2i
1
And
1 2i
1 0 0 1 2i
I A
0 1 1 2i 0
1 1 2i
=
1 2i 1
1 1 1 2i
( I A) 1
1 4i 1 1 2i
2
1
1 1 1 2i
6 1 2i 1
Let B I AI A
1
1 1 1 2i 1 1 2i 1 1 (1 2i)( 1 2i) 1 2i 1 2i
B
6 1 2i 1 1 2i 1 6 1 2i 1 2i (1 2i)(1 2i) 1
1 4 2 4i
B
6 2 4i 4
1 4 2 4i 1 4 2 4i
T
Now B
4
and B
6 2 4i 4 6 2 4i
1 4 2 4i 4 2 4i
T
2 4i 4
36 2 4i 4
B B
1 36 0 1 0
I
36 0 36 0 1
B
T
B 1
Sol: Let A be a unitary matrix. Then AA I
i.e AA 1
I 1
A A1 I
1
A1 A1 I
1
Thus A is unitary.
Rank of a Matrix:
Let A be mxn matrix. If A is a null matrix, we define its rank to be ‘0’. If A is a non-zero
matrix, we say that ‘r’ is the rank of A if
i. Every (r+1)th order minor of A is ‘0’ (zero) &
ii. At least one rth order minor of A which is not zero.
It is denoted by (A) and read as rank of A.
Note: 1. Rank of a matrix is unique.
2. Every matrix will have a rank.
3. If A is a matrix of order mxn, then Rank of A ≤ min (m,n)
4. If ρ (A) = r then every minor of A of order r+1, or minor is zero.
5. Rank of the Identity matrix In is n.
6. If A is a matrix of order n and A is non-singular then ρ (A) = n
7. If A is a singular matrix of order n then (A) < n
Important Note:
1. The rank of a matrix is ≤ r if all minors of (r+1)th order are zero.
2. The rank of a matrix is ≥ r, if there is at least one minor of order ‘r’ which is not equal to
zero.
1 2 3
1. Find the rank of the given matrix 3 4 4
7 10 12
1 2 3
Sol: Given matrix A = 3 4 4
7 10 12
Note: 1. The corresponding column transformations will be denoted by writing ‘c’. i.e
ci ↔cj, ci k cj cj cj + kci
2. The elementary operations on a matrix do not change its rank.
Equivalance of Matrices: If B is obtained from A after a finite number of elementary
transformations on A, then B is said to be equivalent to A.It is denoted as B~A.
Note : 1. If A and B are two equivalent matrices, then rank A = rank B.
2. If A and B have the same size and the same rank, then the two matrices are equivalent.
Elementary Matrix or E-Matrix: A matrix is obtained from a unit matrix by a single
elementary transformation is called elementary matrix or E-matrix.
Notations: We use the following notations to denote the E-Matrices.
1) Eij Matrix obtained by interchange of ith and jth rows (columns).
2) Ei k Matrix obtained by multiplying ith row (column) by a non- zero number k.
3) Eij k Matrix obtained by adding k times of jth row (column) to ith row (column).
1 0 0 0
Eg: 1. 0 1 0 0
is a row echelon form.
0 0 1 1
0 0 0 0
1 3 1
2. 0 1 1 is a row echelon form.
0 0 0
Solved Problems :
2 3 7
1. Find the rank of the matrix A = 3 2 4 by reducing it to Echelon form.
1 3 1
2 3 7
Sol: Given A = 3 2 4 Applying row transformations on A.
1 3 1
1 3 1
R1 ↔ R3 A ~ 3 2 4
2 3 7
1 3 1
R2 → R2 –3R1; R3→ R3 -2R1 ~ 0 7 7
0 9 9
1 3 1
R2 → R2/7,R3→ R3/9 ~ 0 1 1
0 1 1
1 3 1
R3 → R3 –R2 ~ 0 1 1
0 0 0
4 4 3 1
1 1 1 0
A=
k 2 2 2
9 9 k 3
Applying R2 → 4R2-R1, R3 →4R3 – kR1, R4 → 4R4 – 9R1
4 4 3 1
0 0 1 1
We get A ~
0 8 4k 8 3k 8 k
0 0 4k 27 3
Since Rank A = 3 det A =0
0 1 1
4 8 4k 8 3k 8 k 0
0 4k 27 3
2 1 3 5
4 2 1 3
Sol: Given A
8 4 7 13
8 4 3 1
2 1 3 5
0 0 5 7
By applying R2 R2 2R1 ; R3 R3 4R1 ; R4 R4 4 R1 ~
0 0 5 7
0 0 15 21
2 1 3 5
R R R 0 7
R1 1 , R2 2 , R3 3 ~
0 5
1 1 3 0 0 5 7
0 0 5 7
2 1 3 5
0 7
R3 R3 R2 , R4 R4 R2 ~
0 5
0 0 0 0
0 0 0 0
1 2 0 1
Sol: By applying R2 R2 2R1; R3 R3 3R1; R4 R4 R1 A ~ 0 3 1 2
0 3 1 2
0 3 1 2
1 2 0 1
R3 R3 R2 A ~ 0 3 1 2
0 0 0 0
0 3 1 2
1 2 0 1
R3 R4 A ~ 0 3 1 2
0 3 1 2
0 0 0 0
1 2 0 1
R3 R3 R2 A ~ 0 3 1 2
0 0 0 0
0 0 0 0
Solved Problems :
1 2 3 4
1. By reducing the matrix 2 1 4 3 into normal form, find its rank.
3 0 5 10
1 2 3 4
Sol: Given A = 2 1 4 3
3 0 5 10
1 2 3 4
R2 → R2 – 2R1; R3 → R3 – 3R1
A ~ 0 3 2 5
0 6 4 22
1 2 3 4
R3 → R3/-2 A ~ 0 3 2 5
0 3 2 11
1 2 3 4
R3 → R3+R2
A ~ 0 3 2 5
0 0 0 6
0 0 0 0
c2→ c2 - 2c1, c3→c3-3c1, c4→c4-4c1
A ~ 0 3 2 5
0 0 0 6
1 0 0 0
c3 → 3 c3 -2c2, c4→3c4-5c2 A ~ 0 3 0 0
0 0 0 18
1 0 0 0
c2→ c2 /-3, c4→c4/18 A ~ 0 1 0 0
0 0 0 1
1 0 0 0
c4 ↔ c3 A~ 0 1 0 0
0 0 1 0
normal form.
1 1 1 1
Sol: Given A = 1 2 3 4
2 3 5 5
3 4 5 8
1 1 1 1
By applying R2 R2 R1 , R3 R3 2R1 , R4 R4 3R1 0 1 2 5
~
0 1 3 7
0 7 8 5
1 1 1 1
R3 R3 R2 , R4 R4 7 R2 0 1 2 5
~
0 0 1 2
0 0 6 30
1 1 1 1
R4 R4 6R3 0 1 2 5
~
0 0 1 2
0 0 0 18
1 1 1 1
R 0
R4 4 1 2 5
18 ~
0 0 1 2
0 0 0 1
1 0 0 0
0 2 5
Apply C2 C2 C1 , C3 C3 C1 , C4 C4 C1 1
~
0 0 1 2
0 0 0 1
1 0 0 0
0 0 0
C3 C3 2C2 ; C4 C4 5C2 1
~
0 0 1 2
0 0 0 1
1 0 0 0
0 0
C4 C4 2C3 ~
1 0
0 0 1 0
0 0 0 1
3). Define the rank of the matrix and find the rank of the following matrix
2 1 3 5
4 2 1 3
8 4 7 13
8 4 3 1
2 1 3 5
2 1 3
Sol: Let A= 4
8 4 7 13
8 4 3 1
R2 R2 2 R1 2 1 3 5
0 0 5 7
R3 R3 4 R1 A~
0 0 5 7
R4 R4 4 R1
0 0 15 21
2 1 3 5
R3 R3 R2 0 0 5 7
R4 R4 3R2 A~
0 0 0 0
0 0 0 0
It is in echelon form. So, rank of matrix = no. of non zero rows in echelon form.
Rank ( A) 2
2 1 3 4
0 3 4 1
4). Reduce the matrix A to normal form and hence find its rank A
2 3 7 5
2 5 11 6
2 1 3 4
1
Sol: Given A 0 3 4
2 3 7 5
2 5 11 6
1 1 3 4
1 0 3 4 1
C1 C1 A~
2 1 3 7 5
1 5 11 6
1 1 3 4
R3 R3 R2 0 3 4 1
R4 R4 R1 A~
0 2 4 1
0 4 8 2
1 1 3 4
0 0
R2 R2 R3 A~
1 0
0 2 4 1
0 4 8 2
1 0 0 0
R3 R3 2 R2 0 1 0 0
R4 R4 4 R2 A~
0 0 4 1
0 0 8 2
1 0 0 0
0 0
R4 R4 2R3 A~
1 0
0 0 4 1
0 0 0 0
1 0 0 0
0 0
C4 4C4 C3 A~
1 0
0 0 4 1
0 0 0 0
1 0 0 0
1 0 1 0 0 I 0
C3 C3 A~ ⇒ A~ 3
3 0 0 1 0 0 0
0 0 0 0
This is in normal form. Thus Rank of matrix = Order of identify matrix. Rank ( A) 3
0 1 2 2
5). Reduce the matrix A = 4 0 2 6 into canonical form and then find its rank.
2 1 3 1
1 0 2 2
Sol: Apply C1 C2 A ~ 0 4 2 6
1 2 3 1
1 0 2 2
R3 R3 R1 A ~ 0 4 2 6
0 2 1 3
1 0 0 0
C3 C3 2C1; C4 C4 2C1 A ~ 0 4 2 6
0 2 1 3
R2 1 0 0 0
R2 A ~ 0 2 1 3
2
0 0 0 0
1 0 0 0
C2 C3 A ~ 0 1 2 3
0 0 0 0
1 0 0 0
C3 C3 2C2 ; C4 C4 3C2 A ~ 0 1 0 0
0 0 0 0
I 0
Which is in the normal form 2 (A) = 2
0
,
0
Note: .If A is an mxn matrix of rank r, there exists non-singular matrices P and Q such that
I r 0
PAQ =
0 0
Suppose we want to find P and Q we have procedure.
Let order of matrix ‘A’ is ‘3 i.e. A = I3 A I3
1 0 0 1 0 0
A = 0 1
0 A 0 1 0
0 0 1 0 0 1
Now we go on applying elementary row operations and column operations on the matrix A
I 0
(L.H.S) until it is reduced to the normal form r
0 0
Every row operations will also be applied to the pre-factor of on R.H.S
Every column operation will also be applied to the post –factor of on R.H.S.
Solved Problems :
1 0 2
1. Find the non-singular matrices P and Q is of the normal form where A = 2 3 4
3 3 6
Sol: Write A = I3 A I3
1 0 2 1 0 0 1 0 0
~ 2 3 4 = 0 1
0 A 0 1 0
3 3 6 0 0 1 0 0 1
1 0 2 1 0 0 1 0 0
R2 → R2-2R1, R3→R3-3R1 ~ 0 3 0 = 2 1 0 A 0 1
0
0 3 0 3 0 1 0 0 1
1 0 2 1 0 0 1 0 0
R3 → R3-R2, R2 →1/3 R2 ~ 0 1 0 2 / 3 1 / 3 0 A 0 1
0
0 0 0 1 1 1 0 0 1
1 0 0 1 0 0 1 0 2
c3 → c3-2c1 ~ 0 1 0 = 2 / 3 1 / 3 0 A 0 1
0
0 0 0 1 1 1 0 0 1
1 0 0 1 0 2
~[
]= PAQ where P = 2 / 3 1 / 3 0 Q = 0 1 0
1 1 1 0 0 1
2. Find the non-singular matrices P and Q such that the normal form of A is P A Q.
1 3 6 1
Where A = 1 4 5 1 . Hence find its rank.
1 5 4 3
1 3 6 1 1 0 0 1 0 0 0
0 1 0 0
~ 1 4 5 1 =
0 A 1 0
0 0 1 0
1 5 4 3 0 0 1
0 0 0 1
1 3 6 1 1 0 0 0
1 0 0 0
R2 →R2- R1, R3→R3-R1 ~ 0 1 1 2 = 1 1 0 A 1 0 0
0 0 1 0
0 2 2 4 1 0 1
0 0 0 1
1 3 6 1 1 0 0 0
1 0 0 0 0
R3 →R3-2R2 ~ 0 1 1 2 = 1 1 0 A 1 0
0 0 1 0
0 0 0 0 1 2 1
0 0 0 1
I 0 1 3 9 7
1 0 0 0
2 = P A Q Where P = 1 1 0 Q= 1 1 2
0 0 0 0 1 0
1 2 1
0 0 0 1
I 0
Here A ~ 2 , ∴ Hence ρ(A) =2
0 0
System of linear equations: In this chapter we shall apply the theory of matrices to study the
existence and nature of solutions for a system of m linear equations in ‘n’ unknowns.
The system of m linear equations in ‘n’ unknowns x1, x2, x3 xn given by
Hence eqn (2) Reduces to AX = 0 which are known as homogeneous linear equations
Non-Homogeneous Linear equations:
It at least one of b1, b2----bm is non zero. Then B 0, the system Reduces to AX = B is known
as Non-Homogeneous Linear equations.
Solutions: A set of numbers x1 , x2 xn which satisfy all the equations in the system is known
as solution of the system.
Consistent: If the system possesses a solution then the system of equations is said to be
consistent.
Inconsistent: If the system has no solution then the system of equations is said to be
Inconsistent.
Augmented Matrix: A matrix which is obtained by attaching the elements of B as the last
column in the coefficient matrix A is called Augmented Matrix. It is denoted by [A|B]
a11 a12 a13 a1n : b1
A B C a21 a22 a23 a2 n : b2
am1 am 2 am3 amn : b3
b). If ( A / B) (A) = r < n (no. of unknowns) then the system of equations AX = B will
have an infinite no. of solutions. In this case (n-r) variables can be assigned arbitrary values.
Let us first consider n equations in n unknowns ie. m=n then the system will be of the form
a11x1+a12x2+a13x3+…..+a1nxn =b1
a21x1+a22x2+a23x3+…..+a2nxn =b2
…………………………………..
an1x1+an2x2+an3x3+…..+annxn =bn
The above system can be written as AX = B --------- (1)
Where A is an n n matrix.
Solving AX = B using Echelon form:
Consider the system of m equations in n unknowns given by
a11x1+a12x2+a13x3+…..+a1nxn =b1
a21x1+a22x2+a23x3+…..+a2nxn =b2
…………………………………..
am1x1+am2x2+am3x3+…..+amnxn =bm
We know this system can be we write as AX = B
a11 a12 a1n b1
a a22 a2 n b2
The augmented matrix of the above system is [A / B] = 21
am1 am 2 amn bm
1 1 1 4
Consider the Augment matrix is [A /B] [A/B] = 2 5 2 3
1 7 7 5
1 1 1 4
Applying R2 →R2-2R1 and R3 → R3-R1, we get [A/B] ~ 0 3 4 5
0 6 8 1
1 1 1 4
Applying R3→ R3-2R2, we get [A/B] ~ 0 3 4 5
0 0 0 11
1 3 8 10
Augmented matrix [A/B] is [A/B] = 3 1 4 0
2 5 6 3
1 3 8 10
R2 → R2-3R1 R3 → R3 -2R1 ~ 0 10 20 3 0
0 11 22 23
1 3 8 10
R2 →1/10 R2 ~ 0 1 2 3
0 11 22 23
1 3 8 10
R3 →R3-11R2 ~ 0 1 2 3
0 0 0 10
1 1 2 4
The Augmented matrix [A/B] = 2 1 3 9
3 1 1 2
1 1 2 4
Applying R2→ R2-2R1 and R3→R3-3R1, we get [A/B] ~ 0 3 1 1
0 4 7 10
1 1 2 4
Applying R3 →3R3-4R2, we get [A/B] ~ 0 3 1 1
0 0 17 34
-17z = -34 z = 2
-3y-z =1 -3y =z+1 -3y =3 y= -1
and x+y+2z =4 x=4-y-2z =4+1-4=1
x=1, y=-1, z=2 is the solution.
4). Show that the equations x+y+z=6, x+2y+3z=14, x+4y+7z=30 are consistent and solve
them.
1 1 1 x 6
Sol: We write the given equations in the form AX=B i.e. 1 2 3 y 14
1 4 7 z 30
1 1 1 6
The Augmented matrix [A/B] = 1 2 3 14
1 4 7 30
1 1 1 6
Applying R2→ R2-R1 and R3→R3-R1, we get [A/B] ~ 0 1 2 8
0 3 6 24
1 1 1 6
Applying R3 →R3-3R2, we get [A/B] ~ 0 1 2 8
0 0 0 0
1 2 1 3
3 1 2 1
The Augmented matrix is AB
2 2 3 2
1 1 1 1
R2 R2 3R1 1 2 1 3
0 7 5 8
R3 R3 2 R1 ~
0 6 5 0
R4 R4 R1
0 3 2 4
1 2 1 3
0 1 0 4
R2 R2 R3 ~
0 6 5 4
0 3 2 4
1 2 1 3
R3 R3 3R1 0 1 0 4
R4 R4 3R2 ~
0 0 5 20
0 0 2 8
1 2 1 3
R 0 1 0 4
R3 3 ~
5 0 0 1 4
0 0 2 8
1 2 1 3
R4 R4 2R3 0 1 0 4
~
0 0 1 4
0 0 0 0
∴ρ (A) = 3 = ρ (A/B
∴ρ (A) =ρ (A/B) = No. of unknowns = 3
The given system has unique solution.
The systems of equations equivalent to given system are
x 2y z 3 y 4; z 4
x84 3 y 4; z 4
x 3 4 1
x 1, y 4, z 4.
6). Solve x y z 3; 3x 5 y 2 z 8; 5 x 3 y 4 z 14
1 1 1 x 3
Sol: - 3 5 2 y 8
5 3 4 z 14
1 1 1 3
Augmented Matrix is AB 3 5 2 8
5 3 4 14
1 1 1 3
R2 R2 3R1
~ 0 8 1 1
R3 R3 5R1
0 8 1 1
1 1 1 3
R3 R3 R2
~ 0 8 1 1
0 0 0 0
Let z k y
1 k
and x 3
1 k k 24 1 k 8k 23 7k
8 8 8 8
x 238 87 k 238
X y 18 8k X 81 where k is any real number.
z 0 k 1
7). Find whether the following system of equations is consistent. If so solve them.
x 2 y 2 z 2, 3x 2 y z 5, 2 x 5 y 3z 4, x 4 y 6 z 0.
1 2 2 2
3 2 1 x
y 5
Sol: In Matrix form it is 2 5 3 4
z
1 4 6 0
AX B
1 2 2 2
R2 R2 3R1 0 8 5 1
~
R3 R3 2 R1 0 9 1 8
R4 R4 R1 0 2 4 2
1 2 2 2
0 1 4 7
~
0 9 1 8
R2 R2 R3 0 2 4 2
1 2 2 2
0 1 4 7
~
R3 R3 9 R2 0 0 37 55
R4 R4 2 R2 0 0 12 16
1 2 2 2
1 4 7
R4
1
R4 ~ 0
4 0 0 37 55
0 0 3 4
1 2 2 2
1 4 7
R4 37 R4 3R3 ~ 0 is in echelon form
0 0 37 55
0 0 0 17
8). Discuss for what values of , the simultaneous equations x+y+z = 6, x+2y+3z =10,
x+2y+ z = have
(i). No solution
(ii). A unique solution
(iii). An infinite number of solutions.
1 1 1 x 6
Sol: The matrix form of given system of Equations is A X = 1 2 3 y 10 = B
1 2 z
1 1 1 6
The augmented matrix is [A/B] = 1 2 3 10
1 2
1 1 1 6
R2 → R2 – R1, R3 → R3-R1 [A/B] ~ 0 1 2 4
0 1 1 6
1 1 1 6
R3 →R3 – R2 ~ 0 1 2 4
0 0 3 10
Case (i): let ≠ 3 the rank of A = 3 and rank [A/B] = 3
Here the no. of unknowns is ‘3’ ∴ρ (A) =ρ (A/B)= No. of unknowns
The system has unique solution if ≠3 and for any value of ‘ ’.
1 1 1 3
Augmented matrix AB 1 2 2 6
1 a 3 b
1 1 1 3
R2 R2 R1
~ 0 1 1 3
R3 R3 R1
0 a 1 2 b 2
1 1 1 3
R3 R3 R2
~ 0 1 1 3
0 a 3 0 b 9
1 1 1 3
For a 3 & b 9 ~ 0 1 1 3
0 0 0 0
1 1 1 3
For a 3 & b any value ~ 0 1 1 3
0 a 3 0 b 9
1 1 1 3
For a 3&b 9 ~ 0 1 1 3
0 0 0 b 9
1 1 1 x 1
1 2 4 y
Sol: The above system in matrix notation is
1 4 10 z 2
A X B
1 1 1 1
Augmented Matrix is AB 1 2 4
1 4 10 2
R2 R2 R1 1 1 1 1
~ 0 1 3 1
R3 R3 R1
0 3 9 2 1
1 1 1 1
~ 0 1 3 1
R3 R3 3R2 0 0 0 2 3 2
1 1 1 1
The equivalent matrix is 0 1 3 0
0 0 0 0
x 1 2k 1 2
X y 0 3k X 0 k 3 where k is any arbitrary constant.
z 0 k 0 1
1 1 1 1
The equivalent matrix is ~ 0 1 3 1
0 0 0 0
The system of equations equivalent to the given system is x + y + z = 1; y + 3z = 1
Let z k y 1 3k and x (1 3K ) k 1 x 2k
x 0 2k 0 2
X y 1 3k X 1 k 3
z 0 k 0 1
where k is any arbitrary constant.
11). Show that the equations 3x 4 y 5 z a; 4 x 5 y 6 z b;5 x 6 y 7 z c don’t have
3 4 5 a
Augment Matrix is AB 4 5 6 b
5 6 7 c
3 4 5 a
R2 3R2 4 R1
~ 0 1 2 3b 4a
R3 3R3 5R1
0 2 4 3c 5a
3 4 5 a
R3 R3 2R2
~ 0 1 2 3b 4a
0 0 0 3a 6b 3c
3 4 5 1
The equivalent matrix is 0 1 2 1
0 0 0 0
The system has infinite number of solutions. The system of equations equivalent to the
given system 3x 4 y 5 z 1; y 2 z 1 y 2 z 1
Let z k y 1 2k and 3x 4 8k 5k 1 x 1 k
x 1 k 1 1
X y 1 2k 1 k 2
z 0 k 0 1
1 3 1
Sol: The Given Vector X 1 2 X 2 2 X 3 6
3 1 5
The Vectors X1, X2, X3 from a square matrix.
1 3 1 1 3 1
Let A 2 2 6 Then A 2 2 6
3 1 5 3 1 5
= 1(10+6)-2(15-1) + 3(-18+2)
= 16+32-48 = 0
The given vectors are linearly dependent |A| = 0
2). Show that the Vector X1=(2,2,1), X2=(1,4,-1) and X3=(4,6,-3) are linearly dependent.
Sol: Given Vectors X1=(2,-2,1) X2=(1,4,-1) and X3=(4,6,-3) The Vectors X1, X2, X3 form a
square matrix.
2 1 4 2 1 4
Let A 2 4 6 Then A 2 4 6
1 1 3 1 1 3
0
a11 a12 a13 ....a1n x1
x 0
i.e. a 21 a 22 a 23 ....a 2 n 2 Here A is called Co –efficient matrix.
a ...
m1 a m 2 a m3 ....a mn ..
x n .0
Note: 1. Here x1= x2 = --------- xn = 0 is called trivial solution or zero solution of AX = 0
2. A zero solution always linearly dependent.
Solved Problems :
1). Solve the system of equations x+3y-2z = 0, 2x-y+4z = 0, x-11y+14z = 0
1 3 2 x 0
Sol: We write the given system is AX = 0 i.e. 2 1 4 y 0
1 11 4 z 0
1 3 2
R2 →R2 -2R1; R3→R3-R1 A ~ 0 7 8
0 14 16
1 3 2
R3 →R3 -2R2 A ~ 0 7 8
0 0 0
2 1 3 0
The Augumented matrix A O 3 2 1 0
1 4 5 0
1 4 5 0
R1 R3 ~ 3 2 1 0
2 1 3 0
1 4 5 0
R2 R2 3R1
~ 0 14 14 0
R3 R3 2 R1
0 7 7 0
1 4 5 0
R3 2 R3 R2 ~ 0 14 14 0 it is echelon form.
0 0 0 0
x k 1
Thus, the solution set is X y k k 1 K .
z k 1
3). Show that the only real number for which the system x+2y+3z = x, 3x+y+2z= y,
2x+3y+z = z, has non-zero solution is 6 and solve them.
1 2 3 x 0
Sol: Above system can we expressed as AX = 0 i.e. 3 1 2 y 0
2 3 1 z 0
Given system of equations possess a non –zero solution i.e. ρ (A) < no. of unknowns.
For this we must have det A = 0
1 2 3
3 1 2 0
2 3 1
1 1 1
(6 ) 3 1 2 0
2 3 1
→ − , → − − | − − − |=
− −
(6- )[(-2- )(-1- )+1] =0
(6- ) ( 2+3 +3) = 0
= 6 only real values.
When = 6, the given system becomes
5 2 3 x 0
3 5 2 y 0
2 3 5 z 0
5 2 3 x 0
R2 → 5R2+3R1, R3→5R3+2R1 ~ 0 19 19 y 0
0 19 19 z 0
5 2 3 x 0
R3 →R3+R2 ~ 0 19 19 y 0
0 0 0 z 0
into a scalar multiple of itself i.e. AX = Y = X, Then the unknown scalar is known as an
“Eigen value” of the Matrix A and the corresponding non-zero vector X is known as “Eigen
Vector” of A. Corresponding to Eigen value . Thus the Eigen values (or) characteristic
values (or) proper values (or) latent roots are scalars which satisfy the equation.
AX = X for X 0, AX IX 0 ( A I ) X 0
Which represents a system of ‘n’ homogeneous equations in ‘n’ variables x1, x2, ----, xn this
system of equations has non-trivial solutions If the coefficient matrix (A- I) is singular i.e.
a11 a12 a1n
a21 a22 a2 n
A I 0 0
an1 an 2 ann
characteristic vector (or) Eigen vector of A if there exists a scalar such that AX = X.
Method of finding the Eigen vectors of a matrix.
Let A=[aij] be a nxn matrix. Let X be an eigen vector of A corresponding to the eigen value .
Then by definition AX = X.
⇒ AX = IX
⇒ AX – IX = 0
⇒ (A- I)X = 0 ------- (1)
This is a homogeneous system of n equations in n unknowns.
Will have a non-zero solution X if and only |A- I| = 0
A- I is called characteristic matrix of A
|A- I| is a polynomial in of degree n and is called the characteristic polynomial of A
|A- I|=0 is called the characteristic equation
Solving characteristic equation of A, we get the roots , 𝜆 , 𝜆 , 𝜆 , … … . 𝜆 , These are
called the characteristic roots or eigen values of the matrix.
Corresponding to each one of these n eigen values, we can find the characteristic
vectors.
Procedure to find Eigen values and Eigen vectors
……
…..
Let A = [… … . ……. … … .] 𝑖 𝑖
…….
ℎ 𝑖 𝑖 𝑖 𝑖 − λI
Solved Problems
−
1. Find the eigen values and the corresponding eigen vectors of [ ]
−
Sol: = [ ]
−λ −
Characteristic matrix = [ − 𝜆 ] = [ ]
−λ
Characteristic equation is | − I|=0
−λ −
⟹ | |=
−λ
−λ −λ + =
⟹ +λ − λ+ =
⟹ λ − λ+ =
⟹ λ− λ− =
⟹ λ = , are eigen values of A
−λ −
Consider the system [ ] =
−λ
Eigen vector corresponding to 𝛌 =
−
Put λ = in the above system, we get =
−
⟹ − = −−− − = −−−
ℎ =
Let x1 =
x 1
Eigen vector is 1
x2 1
[ ]𝑖 𝑖 𝑖 , 𝑖 𝑖 𝜆=
2 2
Eigen vector
1
[ ]𝑖 𝑖 𝑖 𝑖 𝑖 𝜆=
2. Find the eigen values and the corresponding eigen vectors of matrix [ ]
Sol: Let A = [ ]
𝛌=
[ ] [ ]=[ ]
+ =
=
+ =
=− , =
=𝛼
⇒ = −𝛼 = , =𝛼
−𝛼 −
[ ]= [ ] = 𝛼[ ]
𝛼
−
[ ]𝑖 𝑖
𝛌=
[ ] [ ]=[ ]
= = 𝑖 𝑖. =𝛼
[ ] = [𝛼 ] = 𝛼 [ ]
𝑖 𝑖 [ ]
𝛌=
−
[ − ] [ ]=[ ]
−
− + =
− =
− =
ℎ 𝑖 = , = =∝
=∝ , = , =∝
x1 1
x 0 0
2
x3 1
𝑖 𝑖 [ ]
3 6 2
3. Find the Eigen values and Eigen vectors of the matrix is 6 7 4
2 4 3
8 6 2
Sol: Let A 6 7 4
2 4 3
8 6 2
Consider characteristic equation is A I 0 i.e. 6 7 4 0
2 4 3
Case (i): If =
−
[− − ] =
−
−
[− − ][ ] = [ ]
−
8 x1 6 x2 2 x3 0 (1)
6 x1 7 x2 4 x3 0 (2)
2 x1 4 x2 3 x3 0 (3)
x1 k 1
Eigen Vector is x2 2k k 2
3
x 2 k
2
Case (ii): If =
8 6 2 0
6 7 4 x 0
2 4 3 0
5 6 2 x1 0
6 4 4 x2 0
2 4 0 x3 0
5 x1 6 x2 2 x3 0 (1)
6 x1 4 x2 4 x3 0 (2)
2 x1 4 x2 0 0 (3)
Consider (2) & (3)
x1 x2 x3
6 4 4
2 4 0
x1 x2 x3
k
0 16 08 24 8
x1 x2 x
3 k
16 8 16
x1 x2 x3
k
2 1 2
x1
k , x2 k , x3 2k
2
x1 2k , x2 k , x3 2k
x1 2k 2
Eigen Vector is x2 k k 1
x3
2k
2
Case (iii): If =
7 6 2
6 8 4 x 0
2 4 12
7 6 2 x1 0
6 8 4 x2 0
2 4 12 x3 0
7 x1 6 x2 2 x3 0 (1)
6 x1 8 x2 4 x3 0 (2)
2 x1 4 x2 12 x3 0 (3)
x1 2k 2
Eigen Vector is x2 2k 2 k
x3 k 1
4. Find the Eigen values and the corresponding Eigen vectors of the matrix.
2 2 3
6
2 1
1 2 0
2 2 3
A 6
Sol: Let
2 1
1 2 0
2 2 3
The characteristic equation of A is A I 0 i.e. 2 1 6 0
1 2
2 2 2 6 3
1 12 2 2 1
0
3 2 21 45 0
3 3 5 0
3, 3, 5
The Eigen values are -3,-3, and 5
Case (i): If =−
2 3 2 3 x1 0
We get 1 3 6 x2 0
2
1 2 0 3 x3 0
1 2 3 0
The augment matrix of the system is 2 4 6 0
1 2 3 0
1 2 3 0
Performing R2 2 R1 , R3 R1 , we get 0 0 0 0
0 0 0 0
x1 2 3
Hence x2 k1 1 k2 0
x3 0 1
2 3
So 1 and 0 are the Eigen vectors corresponding to 3
0 1
Case (ii): If =
7 2 3 x1 0
We get 2 4 6
x2 0
1 2 5 x 0
3
7 x1 2 x2 3x3 0 (1)
2 x1 4 x2 6 x3 0 (2)
x1 2 x2 5 x3 0 (3)
x1 x2 x3
k3
20 12 10 6 4 4
x x2 x3
1 k3
8 16 8
x x2 x3
1 k3
1 2 1
x1 1
Eigen vector is x2 2 k3
x3 1
5. Find the Eigen values and Eigen vectors of the matrix A and it’s inverse where
1 3 4
A 0 2 5
0 0 3
1 3 4
Sol: Given A 0 2 5
0 0 3
1 3 4
0 2 5 0
0 0 3
1 2 3 0
1, 2,3 i.e. EigenValues are 1, 2,3
Note: In upper le (or) Lower lar of a square matrix the Eigen values of a diagonal matrix
are just the diagonal elements of the matrix.
Case (i): If =
∴ A I x 0
1 3 4 x1 0
0 2 5 x2 0
0 0 3 x3 0
0 3 4 x1 0
0 1 5 x2 0
0 0 2 x3 0
x1 k1 1
X x2 0 0 k1
x3 0 0
Case (ii): If =
1 3 4 x1 0
0 2 5 x 2 0
0 0 3 x3 0
1 3 4 x1 0
0 0 5 x2 0
0 0 1 x3 0
1 3 4 x1 0
0 0 5 x2 0
0 0 1 x3 0
x1 3x2 4 x3 0; 5 x2 0; x3 0
x1 3k 4(0) 0 x1 3k 0 x1 3k
x1 3k 3
X x2 k k 1
x3 0 0
Case (iii): If =
1 3 4 x1 0
0 2 5 x2 0
0 0 3 x3 0
2 3 4 x1 0
0 1 5 x2 0
0 0 0 x3 0
2 x1 3 x2 4 x3 0; x2 5 x3 0; x3 0
Let x3 k
x2 5 x3 0 x2 5k
and 2 x1 3 x2 4k 0 2 x1 15k 4k 0
19
2 x1 19 x 0 x1 k
2
x1 19 2 k 19
2
X x2 5k k 5
x3 k 1
1 1 1 1 1
Note: Eigen Values of A-1 are , , i.e1, , and the Eigen vectors of A-1 are same as
1 2 3 2 3
Eigen vectors of the matrix A
6. Determine the Eigen values and Eigen vectors of
1 8 4
B 2 A2 A 3I where A
2 2 2
1 8 4
Sol : Given that B 2 A A 3 A
2 2 2
8 4 8 4 56 40
we have A2 A. A 2 2 20 4
2 2
B 2 A2 1 A 3I
2
56 40 1 8 4 1 0
2 3
20 4 2 2 2 0 1
112 80 4 2 3 0
40 8 1 1 0 3
111 78
39 6
Characteristic equation of B is B I 0
111 78
0
39 6
i.e. 2 105 2376 0
33 72 0
33 or 72
Eigen Values of B are 33 and 72.
Case (i): If = 33
111 78
0
39 6
78 78 x1 0
39 39 x2 0
39 x1 78 x2 0 x1 2 x2
x1 x2
k ( say )
2 1
x 2
1 k
x2 1
Case (ii): If = 72
111 78
0
39 6
111 72 78
X 0
39 6 72
39 78
X 0
39 78
39 78 x1 0
39 78 x2 0
39 x1 78 x2 0 x1 2 x2
x1 x2
k ( say )
2 1
x 2
1 k
x2 1
a11 a22 ann a12 (a polynomial of degree n – 2)+ a13 (a polynomial of degree n
-2) + … = 0
(1) n n (a11 a22 .... ann )n1 a polynomial of deg ree (n 2) 0
(1) n n (1) n1 (Trace A)n1 a polynomial of deg ree (n 2) in 0
λ ,λ …..λ ℎ ℎ𝑖 𝑖
− 𝑛+
ℎ = − 𝑛
= s
ℎ |A − λI| = − λ + ⋯ . +a
𝜆= ℎ | |=
− −
− λ +a − λ +a − λ + ….+ a =
−
ℎ = =
−
= | | = det
ℎ
Theorem 2: If is an eigen value of A corresponding to the eigen vector X, then is eigen
value An corresponding to the eigen vector X.
Proof: Since is an eigen value of A corresponding to the eigen value X, we have
AX= ----------(1)
Pre multiply (1) by A, A(AX) = A( X)
(AA)X = (AX)
A2X= ( X)
A2X= 2
X
2
is eigen value of A2 with X itself as the corresponding eigen vector. Thus the
theorm is true for n=2
Let we assume it is true for n = k
i.e. AKX = K
X------------(2)
Premultiplying (2) by A, we get
A(AkX) = A( KX)
(AAK)X= K
(AX)= K
( X)
AK+1X= K+1
X
K+1 K+1
is eigen value of with X itself as the corresponding eigen vector.
n
Thus, by Mathematical induction. is an eigen value of An.
Theorem 3: A Square matrix A and its transpose AT have the same eigen values.
Proof: We have (A- I)T = AT- IT
= AT - I
|(A- I)T|=|AT- I| (or)
|A- I|=|AT- I| A A
T
|A- I|=0 if and only if |AT- I|=0
Hence the theorem.
Theorem 4: If A and B are n-rowed square matrices and If A is invertible show that A-1B and
B A-1 have same eigen values.
Proof: Given A is invertble i.e, A-1 exist
We know that if A and P are the square matrices of order n such that P is non-singular
then A and P-1 AP have the same eigen values.
Taking A=B A-1 and P=A, we have
B A-1 and A-1 (B A-1 ) A have the same eigen values
ie.,B A-1 and (A-1 B)( A-1 A) have the same eigen values
ie.,B A-1 and (A-1 B)I have the same eigen values
ie.,B A-1 and A-1 B have the same eigen values
Theorem 5: If 1 , 2 ,...........n are the eigen values of a matrix A then k 1, k 2, ….. k n
are the eigen value of the matrix KA, where K is a non-zero scalar.
Proof: Let A be a square matrix of order n. Then |KA- KI| = |K (A- I)| = Kn |A- I|
Since K≠0, therefore |KA- KI| = 0 if and only if A I 0
= ( +K) X
+k is an eigen value of the matrix A+KI.
X = A1 X A1 X 1 X ( 0)
Hence 1 is an eigen value of A1
| |
Theorem 10: If is an eigen value of a non-singular matrix A, then is an Eigen value of
(adjA) A A I
A A
X (adj A) X or (adj A) X X
Since X is a non – zero vector, therefore the relation (1)
A
it is clear that is an eigen value of the matrix Adj A
Theorem 11: If is an eigen value of an orthogonal matrix A, then is also an Eigen value A
Proof: We know that if is an eigen value of a matrix A, then is an eigen value of A–1
But the matrices A and A1 have the same eigen values, since the determinants |A- I|
and | A1 - I| are same.
Theorem 12: If is eigen value of A then prove that the eigen value of B = a0A2+a1A+a2I is
2
a0 +a1 +a2
Proof: If X be the eigen vector corresponding to the eigen value , then AX = X --- (1)
Premultiplying by A on both sides
We have B = a0A2+a1A+a2I
BX = a0A2+a1A+a2I) X
= a0A2 X+a1AX+a2 X
2 2
= a0 X+a1 X+a2X = (a0 +a1 +a2 ) X
2
(a0 +a1 +a2 ) is an eigen value of B and the corresponding eigen vector of B is X.
Theorem 13: Suppose that A and P be square matrices of order n such that P is non singular.
Then A and P-1AP have the same eigen values.
Proof: Consider the characteristic equation of P-1AP
1
It is | ( P-1AP)- I| = | P-1AP- P-1IP| ( I P P)
= | P-1 (A- I) P| = | P-1 | |A- I| |P|
= |A- I| since |P-1 | |P| = 1
Thus the characteristic polynomials of P-1AP and A are same. Hence the eigen values of
P-1AP and A are same.
Corollary 1: If A and B are square matrices such that A is non-singular, then A-1B and BA-1
have the same eigen values.
Corollary 2: If A and B are non-singular matrices of the same order, then AB and BA have
the same eigen values.
Theorem 14: The eigen values of a triangular matrix are just the diagonal elements of the
matrix.
i.e.,
Proof: Let be an eigen value of a real symmetric matrix A and Let X be the corresponding
eigen vector then AX=
Take the conjugate
Taking the transpose
Since
Post multiplying by X, we get ------- (2)
T T
Premultiplying (1) with , we get X AX X X ------ (3)
– (3) gives X X 0 but
T
(1) 0
X 2T AT X 2T T
1 X 1T X 2T
T
-------- (3) s
Premultiplying (2) by
Hence from (3) and (4) we get
( 1 2 )
Note: If is an eigen value of A and f(A) is any polynomial in A, then the eigen value of
f(A) is f( ) .
Theorem17: The Eigen values of a Hermitian matrix are real.
Proof: Let A be Hermitian matrix. If X be the Eigen vector corresponding to the eigen value
of A, then AX = X -------------------- (1)
We get X AX X X
i.e X A X
X X ABC C B A and KA K A
X
X , A A -------------------- (3)
(or) X A X X X
From (2) and (3), we have
X X X X
i.e X
X 0 0
X X 0
Hence is real.
Note: The Eigen values of a real symmetric are all real
Corollary: The Eigen values of a skew-Hermitian matrix are either purely imaginary (or)
Zero
Proof: Let A be the skew-Hermitian matrix
If X be the Eigen vector corresponding to the Eigen value of A, then
AX X (or ) iA X i X
AX X X A X 2
T T
T
A I 3
T
Since A is unitary, we have A
AX X
T T T
(1) and (2) given X A X
T T
i.e X X X X From (3)
T
X X 1 0
T
Since X X 0 ,we must have 1 0
1
Since =
We must have =1
Note 1: From the above theorem, we have “The characteristic root of an orthogonal matrix is
of unit modulus”.
2. The only real eigen values of unitary matrix and orthogonal matrix can be 1
Theorem 19: Prove that transpose of a unitary matrix is unitary.
Proof: Let A be a unitary matrix, then A. A A . A I
AA A A I
T T T
AA A A I
T T T
A
T
AT AT A I
T
A
T
AT AT AT I
Solved Problems:
1 2 3
1. For the matrix A 0 3 2 find the Eigen values of 3 A 5 A 6 A 2 I
3 2
0 0 2
1 2 3
Sol: The Characteristic equation of A is A I 0 i.e. 0 3 2 0
0 0 2
⇒ 1 3 2 0
Let f A 3 A 5 A 6 A 2I
3 2
2. Find the eigen values and eigen vectors of the matrix A and its inverse, where
A=
Sol: Given A =
and x1
1
X 0 0 is the solution where is arbritary constant
0 0
1
X 0 is the eigen vector corresponding to 1
0
Case (i): If =
5 x3 0 x3 0
x1 3x2 0 x1 3x2
Let x2 k
x1 3k
3k 3
X k k 1
0 0
3
X 1 is the eigen vector corresponding to 2
0
Case (iii): If =
2 3 4 x1 0
For 3, becomes 0 1 5 x2 0
0 0 0 x3 0
19
X 10 is the eigen vector corresponding to 3
2
A = A
T
3i 2 i
AT 0
2i i
2 2i 8 0
4i, 2i are the Eigen values of A
1 3
4.Find the eigen values of A 2
i
2
3 1
i
2 2
1 3 1 3
i i
Now A 2 2 and T
A 2 2
3 1 i 3 1 i
2 2 2 2
T 1 0
We can see that A . A I
0 1
Thus A is a unitary matrix
1 3
i
2 2 0
3 1
i
2 2
3 1 3 1
Which gives i and i and
2 2 2 2
Hence above values are Eigen values of A.
Cayley-Hamilton Theorem: Every Square Matrix satisfies its own characteristic equation
To find Inverse of matrix: If A is non-singular Matrix, then A-1 exists, Pre multiplying (1)
above by A-1 we have a0 An 1 a1 An 2 an 1I an A1 0 ,
1
A1 a0 An 1 a1 An 2 an 1 I
an
To find the powers of A: - Let K be a +ve integer such that K n
a0 AK a1 AK 1 an Ak n 0
Pre multiplying (1) by Ak-n we get ,
1
AK a1 Ak 1 a2 AK 1 an AK n
a0
Solved Problems :
1. S.T the matrix A = satisfies its characteristic equation and hence find A-1
C2 → C2+C3
⇒ − + − =
By Cayley – Hamilton theorem, we have A3-A2+A-I = 0
1 2 2 1 0 0 1 2 2
A 1 2 3 A 1 1 2 A 2 2 1
2 3
1 2 2 1 0 0 1 2 2 1 0 0
A A A I 2 2 1 1 1 2 1 2 3 0 1 0 =
3 2
2. Using Cayley - Hamilton Theorem find the inverse and A4 of the matrix
−
A =[− − ]
−
-
Sol: Let A = [- - ]
-
7 2 2
The characteristic equation is given by |A- I|=0 i.e., 6 1 2 0
6 2 1
⇒ − + − =
By Cayley – Hamilton theorem we have A3-5A2+7A-3I=0…..(1)
Multiply with A-1 we get
−
= [ − + ]
25 8 8 79 26 26
A 24 7 8 A 78 25 26
2 3
24 8 7 78 26 25
3 2 2
A 6 5 2
11
3
6 2 5
2 1 2
3. If A 5 3 3 Verify Cayley-Hamilton theorem hence find A-1
1 0 2
2 1 2
Sol: - Given that A 5 3 3
1 0 2
2 1 2 2 1 2 7 5 3
A 5 3 3 5 3 3 22 14 13
2
1 0 2 1 0 2 0 1 2
7 5 3 2 1 2 36 22 23
A A . A 22 14 13 5 3 3 101 64 60
3 2
0 1 2 1 0 2 7 3 7
Now A 3 A 7 A I 0
3 2
36 22 23 21 15 9 14 7 14 1 0 0 0 0 0
101 64 60 66 42 39 35 21 21 0 1 0 0 0 0 0
7 3 7 0 3 6 7 0 14 0 0 1 0 0 0
⇒ A3 3 A2 7 A I 0
Multiply A-1, we get
A1 A3 3 A2 7 A I 0
A2 3 A 7 I A1 0
A1 A2 3 A 7 I
7 5 3 6 3 6 7 0 0
A 22 14 13 15 9 9 0 7 0
1
0 1 2 3 0 6 0 0 7
6 2 3
A 7 2 4
1
3 1 1
Check A.A-1= I
2 1 2 6 2 3 1 0 0
A.A = 5 3 3 7 2 4 0
-1
1 0 I
1 0 2 3 1 1 0 0 1
1 2
4. Using Cayley – Hamilton theorem, find A8, if A
2 1
1 2
Sol: Given A
2 1
Characteristic equation of A is A I 0
1 1 4 0
2 5 0 (1)
Substitute A in place of
A2 5 I 0 A 2 5 I
find A8
A8 5 A6 5( A2 )( A2 )( A2 )
5 5 I 5 I 5 I
625I
A8 625I
Diagonalization of a Matrix:
Let A be a square Matrix. If there exists a non-singular Matrix P and a diagonal Matrix D
such that P-1AP=D, then the Matrix A is said to be diagonalizable and D is said to be
“Diagonal” form (or) canonical diagonal form of the Matrix A
Modal Matrix:The modal matrix which diagonalizes A is called the modal Matrix of A and
is obtained by grouping the Eigen vectors of A into a Square Matrix.
Spectral Matrix: The resulting diagonal Matrix D is known as Spectral Matrix.
In this spectral Matrix D whose principal diagonal elements are the Eigen values of the
Matrix.
Calculation of powers of a matrix:
We can obtain the power of a matrx by using diagonalization
Let A be the square matrix then a non-singular matrix P can be found such that D = P-1AP
D2= (P–1AP) (P–1AP)
= P–1A (PP–1) AP
= P–1A2P (since PP–1=I)
Simlarly D3 = P–1A3P
In general Dn = P–1AnP……..(1)
To obtain An, Premultiply (1) by P and post multiply by P–1
n 1
Then PDnP–1 = P(P–1AnP)P–1 = (PP–1)An (PP–1) = An A PD P
n
1n 0 0 0
n
Hence A = P 0
n
2 0 0 1
P
0 0 0 nn
Diagonalization of a matrix:
Theorem: If a square matrix A of order n has n linearly independent eigen vectors
(X1,X2…Xn) corresponding to the n eigen values 1, 2…. n respectively then a matrix P can
be found such that P-1AP is a diagonal matrix.
Note: 1. If X1,X2…Xn are not linearly independent this result is not true.
2. Suppose A is a real symmetric matrix with n pair wise distinct eigen values 1 , 2 n then
the corresponding eigen vectors X1,X2…Xn are pairwise orthogonal.
Hence if P = (e1,e2…en)
Where e1 = (X1 / ||X1||), e2 = (X2 / ||X2||)….en = (Xn)/ ||Xn|| then P will be an orthogonal
matrix.
i.e, PTP=PPT=I
Hence P–1 = PT
−
∴ =
Solved Problems :
matrix.
When =5
Similarly, for the given eigen value = -3 we can have two linearly independent eigen vectors
X2 =
P=[ ]
−
−
= [ ]
=[ − ] = diag [5,-3,-3].
−
−
∴ = diag [5,-3,-3].
−
2. Find a matrix P which transform the matrix A = [ ] to diagonal form. Hence
calculate A4.
=> (1- [ − − − ]− −[ − − ]=
=> 9 − − − =
=> = , = , =
Thus the eigen values of A are 1, 2, 3.
If x1, x2, x3 be the components of an eigen vector corresponding to the eigen value , we have
[A- I] X =
Case (i): If =
P=
1
0 1 1 1 1 2 1
2
0
Now P 1 AP 1 1 0 1 2 1 1 1 1
1
1 1 2 2 3 0 2 2
2
A4 = PD4P-1
1
1 2 1 1 0 0 0 1
2
1 1 1 0 16 0 1 1 0
0 2 2 0 0 81 2 2 1
1 1 3
The characteristic equation of A is A I 0 i.e. 1 5 1 0
3 1 1
1 1 3 x1 0
⇒ 1 5 1 x2 0
3 1 1 x3 0
3 1 3 x1 0
⇒ 1 7 1 x 0
2
3 1 3 x3 0
x1 20k 1
X x2 0k 20k 0
x3 20k 1
Case (ii): If =
[A- I]X =
2 1 3 x1 0
⇒ 1 2 1 x2 0
3 1 2 x3 0
2 x1 x2 3 x3 0 (1)
x1 2 x2 x3 0 (2)
3 x1 x2 2 x3 0 (3)
Case (iii): If =6
[A- I]X =
1 1 3 x1 0
1 5 1 x2 0
3 1 1 x3 0
5 1 3 x1 0
1 11 1 x2 0
3 1 5 x3 0
5 x1 x2 3x3 0 (1)
x1 x2 x3 0 (2)
3x1 x2 5 x3 0 (3)
Consider (2) & (3)
x1 x2 x3
1 1 1
3 1 5
x1 x2 x
3 k
5 1 5 3 1 3
x x x
1 2 3 k
4 8 4
x x x
1 2 3 k
1 2 1
x1 1
X x2 2 k
x3 1
1 1 1
0 1 2
1 1 1
3 2 1
0 2 2
3 2 1
3 2 1
Adj ( ) 0 2 2
3 2 1
Adj
1
3 0 3 12 0 1
2
1
Cofactor of 2 2 1
2 1 1
6 3 3 3
1 2
1
6
1 1
3
1
6
D 1 A
12 0 1
2 1 1 3 1 1 1
D 13 13 1
3
1 5 1 0 1 2
16 1
3
1
1 1 1 1
6 3 1
12 0 2 1 0 3 1 1 3
1
1 2 3
13 13 1
3 1 0 1 1 5 1 1 10 1
16 1
3
1
6 3 0 1 3 1 1 3 2 1
12 0 1
2 2 3 6 2 0 0
13 13 1
3
0 3 12 0 3 0
16 1
3
1
6
2 3 6
0 0 6
1 2 3 4 1 0 4 1 0 4 2 0
1 6 2 3 2 4 4 8 4 x 0
1 2 5 2 4 4 0
2 5 2 3 5 2 2 4 4 0
3 6 2 11 6 0
1, 2,3
=>[ ][ ] = [ ]
−
x1+x2=0, x1+x2=0, -4x1+4x2+2x3=0
Let x3 = k, x2 + k=0, x2 = -k
k
4 x1 4( k ) 2 K 0 4 x1 2k 0 4 x1 2k x1
2
x1 2k 12 1
x k 1 k 2 k
2 2
x3 k 1 2
1
X 1 2
2
Case (ii): If =
=> [ − I]X =
−
=> [ − ] =
− −
−
=> [ ][ ] = [ ]
−
x1 x2 x3 0 (1)
x3 0 (2)
4 x1 4 x2 x3 0 (3)
x1 k 1 1
x2 k 1 k X 2 1
x3 0 0
0
Case (iii): If =
=> [ − I]X =
−
=> [ − ] =
− −
−
=> [ − ][ ] = [ ]
−
2 x1 x2 x3 0
x2 x3 0
4 x1 4 x2 0
Let x1 k and x3 k
2 x1 x2 x3 0 2k x2 k 0 x2 k
x1 k 1
X 3 x2 k 1 k
x3 k 1
1 1 1
p X 1 X 2 X 3 2 1 1
2 0 1
1 1 0
1
4 3 1
2 2 1
1 0 0 1 0 0
D 1
AP 0 2 0 0 2 0
0 0 3 0 0 3
18 0 0 1 0 0
D8 0 28 0 0 256 0
0 0 38 0 0 6561
( a ). A8 PD8 P 1
1 1 1 1 0 0 1 1 0
2 1
1 0 256 0 4 3 1
2 0 1 0 0 6561 2 2 1
12099 12355 6305
12100 12356 6305
13120 13120 6561
14 0 0 1 0 0
(b). D 4 0 24 0 0 16 0
0 0 34 0 0 81
1 1 1 1 0 0 1 1 0
A =PD P 2 1 1 0 16 0 4 3 1
4 4 -1
UNIT-II
1. Evaluate lim
→ + +
→
Sol. lim + +
= lim {lim [ + +
]}
→ → →
→
= lim
→ +
=
=
(or)
lim = lim {lim [ ]}
→ + + → → + +
→
= lim
→ +
=
=
−
2. If f(x ,y) = +
show that lim {lim , } ≠ lim {lim , }
→ → → →
−
Sol. lim {lim , } = lim {lim }
→ → → → +
= lim
→
=
−
lim {lim , } = lim {lim }
→ → → → +
−
= lim
→
= -1
Hence the result follows.
= lim +
→
= +
Which is different for the different m selected.
∴ lim
→
, does not exist.
→
Consider
lim , = lim +
=lim = = ,
→ → →
lim , = lim =lim = = ,
→ → + →
∴ f(x,y) is continuous for given values of x and y but it is not continuous at (0,0)
Partial Differentiation:
+∆ , − ,
Let z = f x, y be a function of two variables x and y. Then lim ∆
, if it exists ,
→
i.e, for the partial derivative of = , w.r.t. ‘x’ , ‘y’ is kept constant.
Similarly, the partial derivative of = , w.r.t. ‘y’ , ‘x’ is kept constant and is defined
, +∆ − ,
as lim and is denoted by or or .
→ ∆
In general the first order partial derivatives and are also functions of x and y and they
can be differentiated repeatedly to get higher order partial derivatives,
∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂
So ∂
=∂ , =∂ , =∂ ,∂ =∂ ,
∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂
∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂
∂ ∂
= ∂
, ∂ ∂
= ∂
, ∂ ∂
= ∂ ∂
, and so on.
Then the derivative of z w.r.t. ‘t’ i.e, is called the total differential coefficient or total
derivative of z.
∂ ∂ ∂ ∂ ∂
∴ ∂
= ∂ ∂
+ ∂ ∂
∂ ∂
Note: In the differential form , this result can be written as = ∂
. dx + ∂
. dy
x + y + z − xy − yz − zx
=
x + y + z x + y + z − xy − yz − zx
∂U ∂U ∂U
⟹ + + = ……..(1)
∂ ∂ ∂ + +
∂ ∂ ∂ ∂ ∂ ∂ ∂U ∂U ∂U
Now ∂
+∂ +∂ U = +∂ +∂ +∂ +∂
∂ ∂
∂ ∂ ∂
= +∂ +∂ [from ]
∂ + +
∂ ∂ ∂
= + +
∂ + + ∂ + + ∂ + +
=- + +
- + +
− + +
=- + +
𝛛 −
2. If = e show that at x = y = z , 𝛛 𝛛
=− 𝐥
When x = y = z , we have
= -1 and = -1
− ∂ +l ∂
= - (1+ log y)[− + log z ]= …….(3)
∂ +l ∂
=- +l
=- l +l
(since log e =1)
−
= - = − x log ex
l
=− , = , =
∂ ∂ ∂
and =1, ∂ = -1 , ∂ = 0
∂
= + +
= + − + =− + …(1)
= + +
= + − + − = − …(2)
and
∂u ∂u ∂r ∂u ∂s ∂u ∂t
= + +
∂z ∂r ∂z ∂s ∂z ∂t ∂z
∂ ∂ ∂ ∂ ∂
= ∂
− + + =−∂ +∂ …(3)
∂ ∂
Jacobian :
Let u = u (x , y) , v = v(x , y) are two functions of the independent variables x , y.
The jacobian of ( u , v ) w .r .t (x , y ) or the jacobian transformation is given by the
∂ ∂
∂ ∂
determinant |∂ ∂
| (or)
∂ ∂
J( ) = =
Properties of Jacobians
(u, v) ( x, y )
1. If J and J 1 then JJ 1 1
( x, y ) (u, v)
, , ,
2. If u,v are functions of r, s and r, s are functions of x ,y , then = .
, , ,
Solved Problems:
( x, y , z )
1. If x + y2 = u , y + z2 = v , z + x2 = w find
(u , v, w)
Sol : Given x + y2 = u , y + z2 = v , z + x2 = w
We have = =
= =
=6
c1 c1 c2
c2 c2 c3
0 0 1
6 2x 2 y 2 y 2z z yx
x yz y 2 xz
2
y xz z xy z 2 xy
2 2
Sol: x + y + z = u
y + z = uv
z = uvw
y = uv – uvw = uv (1 – w)
x = u – uv = u (1 – v)
R2 R2 R3
= uv [ u –uv +uv]
= u2v
= =
= (2r)(2r)
We have
ux = yz(-1/x2) = , uy = , uz =
= , xz(-1/y2) = ,
= , = , = xy (-1/z2) =
= . .
= , =
= = =
= , = sin-1( )
1 1 y
x y 2
y x x x2 y2
2
1 2
x
= (1/x) =
x y
r , (x y2 ) (x y2 )
2 2
y 1
( x, y )
x x2 y2 x2 y2
=[ - ]
= = =
. =1
Functional Dependence:
Two functions u and v are functionally dependent if their Jacobian i.e,
J( ) = = =0
If the Jacobian of u, v is not equal to zero then those functions u, v are functionally
independent.
Solved Problems :
+ − − ,
1. If = −
, = 𝐚 + 𝐚 . Find ,
. Hence prove that u and v are
∂ ∂ + +
∂ , ∂ ∂ − −
∴ ∂ ,
= |∂ ∂
| =| |= −
- −
=0
∂ ∂ + +
∂ , u u e sin y e cos y
Sol: i) Jacobian = ∂ = |v v | =| |
, e cos y − e sin y
= e ( -sin y − cos y = − ≠
∴ u,v are functionally independent .
+
ii) u= ,v=
−
−
, ∂ ,
∴ J( ) = =| − | = - =0
, ∂ , − −
− −
+
∴v= −
is the functional relation between u and v .
+ + +
=2| | ( Applying → + )
= 2(x + y + z) | |
∂ ∂ ∂
∂ ∂ ∂
−
, , |∂ ∂ ∂ |
Now J( = ∂ =| − − | = 2(-1) | − | = (-2)(0) =0
, | ∂ ∂ |
∂ ∂ ∂ −
∂ ∂ ∂
, , | | − − − − − −
Now = =| |
, , | |
− − −
=6| − − − − − − |
− − −
=6 | − − − − | C1 → − and C2
− + + − + + −
→ −
, , − −
∴ , ,
= 12 | − + + − + + |
= 12 (x-y)(y-z) | |
+ + + +
= 12 (x-y)(y-z) (0) [C1 and C2 are identical]
=0
Hence the functional relationship exists between u, v, w.
−
6. Prove that u = ,v= are functionally dependent and find the relation between
+ +
them.
−
Sol. We are given u = ,v=
+ +
( + ). −( − ). + − +
∴ = = =
+ + +
( + ). − −( − ). − + − + −
= = =
+ + +
( + ). − . −
=2 y [ ]= and
+ +
( + ). − . −
=2 x [ ]=
+ +
−
, + +
Thus =| |=| − −
|
,
+ +
− −
= +
+ +
( − )− −
= +
=0
∴ u,v are functionally dependent.
− +
u2 + v2 = + = =1
+ + +
2 2
Hence u + v = 1 is the functional relation between u and v.
= 0 i.e the pairs (a1, b1), (a2, b2) ………….. are called
Stationary.
Working procedure:
1. Find and Equate each to zero. Solve these equations for x & y we get the pair
3 i) If l n –m2 > 0 and l < 0 at (a1,b1) then f(x ,y) is maximum at (a1,b1) and maximum
value is f(a1,b1)
ii) If l n –m2 > 0 and l > 0 at (a1,b1) then f(x ,y) is minimum at (a1,b1) and minimum
value is f(a1,b1) .
iii) If l n –m2 < 0 and at (a1, b1) then f(x, y) is neither maximum nor minimum at (a1, b1).
In this case (a1, b1) is saddle point.
iii) If l n –m2 = 0 and at (a1, b1) , no conclusion can be drawn about maximum or
minimum and needs further investigation. Similarly we do this for other stationary
points.
Solved Problems:
1. Locate the stationary points & examine their nature of the following functions.
u =x4 + y4 -2x2 +4xy -2y2, (x > 0, y > 0)
Sol: Given u(x ,y) = x4 + y4 -2x2 +4xy -2y2
u
For maxima & minima = 0, u = 0
x y
(1) x3 – 2x x = 0, 2, 2
Hence (3) y = 0, - 2, 2
2u 2u 2u
l= = 12x 2
– 4, m = = ( ) =4&n= = 12y2 – 4
x 2 xy y 2
ln – m2 = (12x2 – 4 )( 12y2 – 4 ) -16
At ( , ), ln – m2 = (24 – 4)(24 -4) -16 = (20) (20) – 16 > 0 and l=20>0
The function has minimum value at ( , )
At (0,0) , ln – m2 = (0– 4)(0 -4) -16 = 0
(0,0) is not a extreme value.
2. Investigate the maxima & minima, if any, of the function f(x) = x3y2 (1-x-y).
2x = 1 => x = ½
4 ( ½) + 3y – 3 = 0 => 3y = 3 -2 , y = (1/3)
2 f
l= = 6xy2-12x2y2 -6xy3
x 2
2 f
2 (1/2,1/3) = 6(1/2)(1/3)2 -12 (1/2)2(1/3)2 -6(1/2)(1/3)3 = 1/3 – 1/3 -1/9 = -1/9
x
2 f f
m= = = 6x2y -8 x3y – 9x2y2
xy x y
2 f
(1/2 ,1/3) = 6(1/2)2(1/3) -8 (1/2)3(1/3) -9(1/2)2(1/3)3 = =
xy
2 f
n= = 2x3 -2x4 -6x3y
y 2
2 f
2 (1/2,1/3) = 2(1/2)3 -2(1/2)4 -6(1/2)3(1/3) = - - = -
y
1
ln- m2 =(-1/9)(-1/8) –(-1/12)2 = - = = > 0 and l = 0
9
The function has a maximum value at (1/2 , 1/3)
1 1 1 1 1 1 1 1 1 1
Maximum value is f , 1
2 3 8 9 2 3 72 2 3 432
3. Find three positive numbers whose sum is 100 and whose product is maximum.
Sol: Let x ,y ,z be three +ve numbers.
Then x + y + z = 100
z = 100 – x – y
Let f (x,y) = xyz =xy(100 – x – y) =100xy –x2y-xy2
100 -2x –y = 0
200 -2x -4y =0
-----------------------------
-100 + 3y = 0 => 3y =100 => y =100/3
100 – x –(200/3) = 0 => x = 100/3
2 f
= =- 2y
x 2
2 f
2 (100/3 , 100/3 ) = - 200/3
x
2 f f
= = = 100 -2x -2y
xy x y
2 f
(100/3 , 100/3 ) = 100 –(200/3) –(200/3) = -(100/3)
xy
2 f
= = -2x
y 2
2 f
2 (100/3 , 100/3 ) = - 200/3
y
ln -m2 = (-200/3) (-200/3) - (-100/3)2 = (100)2 /3
The function has a maximum value at (100/3 , 100/3)
100 100 100
i.e. at x = 100/3, y = 100/3 z = 100
3 3 3
The required numbers are x = 100/3, y = 100/3, z = 100/3
4. Find the maxima & minima of the function f(x) = 2(x2 –y2) –x4 +y4
Sol: Given f(x) = 2(x2 –y2) –x4 +y4 = 2x2 –2y2 –x4 +y4
2 f
l = 2 = 4-12x2
x
2 f
m = = f
=0
x y x y
2 f
n = 2 = -4 +12y2
y
we have ln – m2 = (4-12x2)( -4 +12y2 ) – 0
= -16 +48x2 +48y2 -144x2y2
= 48x2 +48y2 -144x2y2 -16
i) At ( 0 , ± 1 )
ln – m2 = 0 + 48 - 0 -16 =32 > 0
l = 4-0 = 4 > 0
F
= 0 => + = 0 --------------- (4)
z
Step 3: Solving the equations (1) (2) (3) & (4) we get the stationary point (x, y, z).
Step 4 : Substitute the value of x , y , z so obtained in equation (1) we get the extremum.
Solved Problems:
1. Find the minimum value of x2 +y2 +z2, given x + y + z =3a
Sol: u = x2 +y2 +z2
= x + y + z - 3a = 0
Using Lagrange’s function
F(x , y , z) = u(x , y , z) + ( x , y , z)
For maxima or minima
F
= + = 2x + = 0 ------------ (1)
x
F
= + = 2y + = 0 ------------ (2)
y
F
= + = 2z + = 0 ------------ (3)
z
(1) , (2) & (3)
= -2x = -2y = -2z
= x + x + x - 3a = 0
=a
= y =z = a
Minimum value of u = a2 + a2 + a2 =3 a2
2. Find the minimum value of + + , given that xyz = 𝐚
Sol: Let u = x + y + z …… (1)
And ∅ = xyz - =0 …… (2)
Consider the lagrangean function F(x,y,z) = u( x,y,z ) + λ ϕ x, y, z
i.e , F(x,y,z) = x + y + z + λ xyz − a ) ……(3)
∂F ∂ ∂∅
Now ∂ = 0 ⇒ ∂ + λ ∂ = 2x + λyz = 0 ……(4)
∂F ∂ ∂∅
∂
= 0 ⇒ ∂ + λ ∂ = 2y + λxz = 0 ……(5)
∂F ∂ ∂∅
= 0 ⇒ ∂ + λ ∂ = 2z + λyx = 0 ……(6)
∂
λ
From (4) , (5) and (6) ,we have = = = − ……(7)
4. Show that the rectangular solid of maximum volume that can be inscribed in a
sphere is a cube.
Sol: Let 2x ,2y, 2z are the length , breadth and height of rectangular solid
Then its volume V = 8 xyz ……..(1)
Let the sphere have a radius of ‘r’ so that x + y + z = r ……(2)
Consider the lagrangean function F(x,y,z) = u( x,y,z ) + λ ϕ x, y, z
i.e, F(x,y,z) = V + 𝜆 x + y + z − r
= 8xyz + 𝜆 x + y + z − r …….(3)
∂F ∂F ∂F
For maxima or minima ∂ = 0 , ∂ = 0 , ∂
=0
∂F
= 0 ⇒ 8yz +2 λx = 0 ……(4)
∂
∂F
∂
= 0 ⇒ 8zx +2 λy = 0 ……(5)
∂F
∂
= 0 ⇒ 8xz +2 λz = 0 ……(6)
Thus for a maximum value x = y = z which shows that the rectangular solid is a cube.
y−b f a, b + (y − b f a. b ] + [ x − a f a, b + x−a y−
!
f a, b = 0 ⇒ f ,− =0
f a, b = 0 ⇒ f ,− =0
f a, b = 2 ⇒ f ,− =2
f a, b = 0 ⇒ f ,− =0
f a, b = 0 ⇒ f ,− =0
All other partial derivatives of higher order will vanish
By substituting above values in taylor ′ s series , we get
x y+ y− = +[ x− − + y+ ]+ [ x− − + x− y+
+ y+ + [ x− +3(x − y+ + x− y+ +
y+ ]
= 10 - 4(x - 1) + 4( y + 2) - 2(x − + 2 (x - 1) (y + 2) + (x − (y + 2)
UNIT – III
E. g: 1. 4z
x y
z z
2. x y 2z
x y
Order of a Differential equation: It is the order of the highest derivative occurring in
the Differential equation. Differential equation is said to be of order ‘n’ if the
derivative is the highest derivative in that equation.
dy
E. g : (1). (x2+1) . + 2xy = 4x2
dx
Order of this Differential equation is 1.
d2 y dy
(2) x 2 2x 1 x 1 y e x
dx dx
1 x 2 dx
dy
2xy. 1 y 0 . Degree = 2
dy
dx
2
32
d 2 y dy
2
2) a. 1 on solving . we get
dx 2 dx
3
d 2 y dy
2 2
a . 2 1 . Degree = 2
2
dx dx
dy
The general form of first order, first degree differential equation is f x, y or
dx
′
, , = [i.e + = Where are functions of and ]. There
is no general method to solve any first order differential equation The equation which
belong to one of the following types can be easily solved.
In general the first order first degree differential equation can be classified as:
(1) Exact equations
(2) Non exact equations (reducible to exact equations).
Exact Differential Equations
Def: Let M(x, y) dx +N(x, y) dy =0 be a first order and first degree Differential Equation
where M & N are real valued functions of x, y . Then the equation M dx + N dy =0 is said to
be an exact Differential equation if a function f .
f f
d f x, y dx dy
x y
Eg : d ( x 2 y) 2 xydx x 2 dy
Condition for Exactness: If , & , are two real functions which have
continuous partial derivatives then the necessary and sufficient condition for the Differential
equation + = is to be exact if
M N
y x
Hence solution of the exact equation , + , = is
Mdx Ndy c .
(y is taken as constant) (terms free from x are taken).
Solved Problems :
dy y cos x sin y y
1. Solve 0.
dx sin x x cos y x
Sol : Given equation can be written as
( y cos x sin y y )dx (sin x x cos y x)dy 0 …(1)
It is of the form + = .
Here
M y cos x sin y y
N sin x x cos y x
M
cos x cos y 1
y
N
cos x cos y 1
x
M N
Clearly
y x
⟹Equation is exact.
The general solution is given by Mdx Ndy c
y sin x (sin y y ) x c .
x
x
x
2. Solve 1 e y dx e y 1 dy 0
y
x
x x
y
x N 1 x 1
x x x
M
ey 2 & e y 1 e y
y y x y y y
x N x
x x
M
ey 2 & ey 2
y y x y
M N
equation is exact
y x
General solution is
Mdx Ndy c.
ey
x c
1
y
x
x ye c y
Sol : Here M x y 2 ;
N x y 1
M N
1; 1
y x
M N
Clearly
y x
Thus the equation is exact.
General solution is Mdx Ndy c.
x2 y2
xy 2 x yc
2 2
x 2 y 2 2 xy 4 x 2 y c1
M N
e y cos x; e y cos x
y y
M N
e y c os x
y x
Equation is exact.
Gen Sol. is Mdx Ndy c.
e 1 cos x dx 0 dy c
y
e y sin x c
1
5. Solve y 1 cos y dx x log x x sin y dy 0 .
x
1
Sol : Here M y 1 cos y, N x log x x sin y
x
M 1 N 1
1 sin y 1 sin y
y x x x
M N
so the equation is exact
y x
x sin y y cos x c
ycosx - xsiny = c
M N
We have sin cos ; cos sin
r
M N
Clearly
r
The given equation is exact.
x2 y 2
4. d = x dx + y dy
2
x xdy ydx
5. d log =
y xy
x ydx xdy
6. d log =
y xy
x ydx xdy
7. d tan 1 =
y x2 y 2
x xdy ydx
8. d tan 1 =
y x2 y 2
xdy ydx
9. d log xy =
xy
2 xdx ydy
10. d log x 2 y 2 =
x2 y 2
ex ye x dx e x dy
11. d =
y y2
Solved Problems
y( xy e x )dx e x dy
1) Solve 0.
y2
( xy 2 ye x )dx e x dy
Sol : It can be written as 0.
y2
xy 2 e x ydx e x dy
2 dx 0
y y2
ye x dx e x dy
xdx 0
y2
ex
xdx d 0
y
x2 ex
On integrating, we get c
2 y
This is the required solution.
xdy ydx
3) Solve xdx ydy 0
x2 y 2
xdy ydx
Sol : Given equation is xdx ydy 0
x2 y 2
x2 y 2 1 y
d d tan 0 on Integrating we get
2 x
x2 y 2 y
tan 1 c .
2 x
4) Solve y x3 .e xy y dx x y x 3 .e xy dy 0 .
Dividing by x3
y xdy ydx
e xy ydx xdy . 0
x x2
y y
d e xy .d 0
x x
on Integrating
2
y
e xy 1 c is required G.S.
2 x
1
On integrating we get log y log x log c
xy
1
log x log y log c .
xy
ydx xdy
a dx
x 2 y2
x
d tan 1 a dx
y
x
On Integrating tan 1 ax c where c is an arbitrary constant.
y
1
Mx Ny 0 then is an integrating factor of M dx+ N dy =0.
Mx Ny
Solved Problems :
1 . Solve x 2 ydx x 3 y3 dy 0
Where M x 2 y & N x 3 y3
M N
Consider x2 & 3x 2
y x
M N
y x
equation is not exact .
But given equation (1) is homogeneous differential equation then
Mx Ny x x 2 y y x 3 y3 y 4 0 .
1 1
I.F = 4
Mx Ny y
1
Multiplying equation (1) by
y4
x2 y x 3 y3
dx dy 0 ---------------------- (2)
y4 y4
x2 x 3 y3
dx dy 0
y3 y4
x 2 x 3 y3
For M1 & N1
y3 y4
M1 3x 2 N1 3x 2
4 & 4
y y x y
M1 N1
equation (2) is an exact D.E.
y x
x3
log y c
3y3
2. Solve y y 2 2x 2 dx x 2y 2 x 2 dy 0
M N
Consider 3y2 2x 2 & 2y 2 3x 2
y x
M N
equation is not exact .
y x
Since equation(1) is Homogeneous differential equation then
3xy y 2 x 2 0 .
1
I.F.
3xy y 2 x 2
1
Multiplying equation (1) by we get
3xy y 2 x 2
y y2 2x 2 x 2y 2 x 2
dx dy 0
3xy y 2 x 2 3xy y 2 x 2
y(y2 - 2x 2 ) x(2y2 - x 2 )
M1 = ; N1 =
3xy(y2 - x 2 ) 3xy(y2 - x 2 )
y 2
x2 x2
dx
y2 x y2 x 2
dy 0
3x y 2 x 2 3y y 2 x 2
dx xdx ydy dy
2 2 0.
x y x 2
y x 2
y
dx dy 2ydy 2xdx
0
x y 2 y x 2 y2 x 2
2 2
On integrating we get
Problems :
1 . Solve xysin xy cos xy ydx xysin xy cos xy xdy 0 .
Sol : Given equation xysin xy cos xy ydx xysin xy cos xy xdy 0 -------(1).
N xysin xy cos xy x
M N
y x
xysin xy y dx xysin xy cos xy x dy 0
2xy cos xy 2xy cos xy
1 1
y tan xy dx y tan xy dy 0
x y
M1dx N1dy 0
M1 N1
= tanxy + xysec2 xy =
y x
Now the equation is exact.
General solution is M dx N dy c. .
1 1
y.log sec xy
log x log y log c
y
x
log sec xy log log c.
y
x
.sec xy c.
y
2. Solve 1 xy ydx 1 xy xdy 0
Sol : Here M 1 xy y : N 1 xy x
M N
= 1+ 2xy; = 1- 2xy
y x
Hence, the equation is not exact
Also Mx - Ny= 2 x 2 y 2 0
1 1
I.F 2 2 0
Mx Ny 2x y
Multiply the given equation by I.F, we get
1 1 1
2 dx dy 0
2x y 2x 2y
M1 -1 N
= 2 2 = 1
y 2x y x
⟹Equation is exact.
On integrating, we get
1 1 1
2x y 2x dx 2y dy c
2
1 1 1
log x log y c .
2xy 2 2
1 x
log c1 where = 2c.
xy y
Method 4: If there exists a continuous single variable function such that
M N
y x
f x , then I.F. of Mdx Ndy 0 is e
f x dx
N
Solved Problems :
1 . Solve 3xy 2ay 2 dx x 2 2axy dy 0
M N
y x 1
f x .
N x
1
e x x is an Integrating factor of (1)
dx
3xy 2ay xdx x 2 2
2axy
xdy 0
1 1
3x y 2ay x dx x
2 2 3
2ax 2 y dy 0
M1 3x 2 y 2ay 2 x, N1 x3 2ax 2 y
M1 N1
3x 2 4axy, 3x 2 4axy
y x
M1 N1
Equation is exact
y x
x 3 y ax 2 y2 c ..
2. Solve ydx xdy 1 x 2 dx x 2 sin ydy 0
Sol: Given equation is y 1 x 2 dx x 2 sin y x dy 0.
M y 1 x 2 & N x 2 sin y x
M N
1, 2x sin y 1
y x
M N
the equation is not exact.
y x
M N
y x 1 2x sin y 1 2x sin y 2 2 x sin y 1 2
So consider 2
N x 2 sin y x x sin y x x x sin y 1 x
1
2 dx
I.F e f x dx 1
e x e2log x 2
x
y 1 x2 x 2 sin y x
Equation (1) x I.F gives dx dy 0
x2 x2
It is the form of M1dx+ N1dy =0.
M1 1 N1
= 2=
y x x
Equation is exact
y 1
Gen. sol. is thus 2 2 1dx sin ydy 0
x x
y 1
x cos y c .
x x
x 2 y 1 x cos y cx.
⇒ ∫ − − − + ∫ =
(y constant)
⇒ -x + + = c
This is the general solution o (2) and hence of (1)
M N
y x
Method -5: For the equation Mdx + Ndy = 0 if g ( y ) (is a function of y alone)
M
Solved Problems :
1 . Solve (3x2y4+2xy)dx +(2x3y3-x2) dy = 0
Sol : Given equation (3x2y4+2xy)dx +(2x3y3-x2) dy =0 -----------------(1).
Equation is of the form Mdx + Ndy =0.
where M =3x2y4+2xy & N = 2x3y3-x2
M N
12x 2 y3 2x ; 6x 2 y3 2x
y x
M N
equation(1) is not exact.
y x
M N
y x 2
So consider g ( y) g y
M y
1
2 dy
g ( y )dy 2 log y
I .F e
y 1
e e .
y 2
3x 2 y4 2xy 2x 3 y3 x 2
Equation (1) x I.F dx dy 0
y2 y2
2x x2
3x 2 y2 dx 2x 3 y 2 dy 0
y y
M1 x N1
= 6x 2 y - 2 2 =
y y x
Equation is exact
3x 3 y2 2x 2
c.
3 2y
x2
x 3 y2 c.
y
2 . Solve (xy3+y) dx + 2(x2y2+x+y4) dy =0
Sol : Here M xy3 y ; N 2 x 2 y 2 x y 4
M N
3xy 2 1; 4xy 2 2
y x
We see equation is not exact.
M N
y x
M N
y x
Also g y
M
xy 2 1 1
g(y)
y xy 1
2
y
1
dy
Thus I .F e
g ( y ) dy
e y y.
M1 N
= 4xy3 + 2y = 1 where M1 = xy 4 + y 2 ; N1 = 2x 2 y3 + 2xy + 2y 5
y x
xy y 2 dx 2y5 dy c
4
Gen Sol:
x 2 y4 2y6
y2 x c.
2 6
3. Solve y 4 2y dx xy3 2y 4 4x dy 0
y
dy dx
2) Replace by in (2)
dx dy
Then the Differential Equation of family of O.T is
dx
F(x, y, ) =0 ---------------- (3).
dy
3) Solve equation (3) to get the equation of family of O.T’s of equation(1).
Solved Problems
1 . Find the O.T’s of family of semi-cubical parabolas ay2=x3 where a is a parameter.
Sol : The given family of semi-cubical parabola is ay2=x3 ---------------(1)
dy
Differentiating with respect to ‘x’ ⟹ a 2y = 3 x 2 -------------(2)
dx
x3 dy
Eliminating ‘a’ from (1) and (2) ⟹ y 2 2y dx 3x
2
2 x 3 dy
⟹ 3x 2
y dx
dy dx 2x 3
dx
Replace by ⟹ 3x 2
dx dy y dy
-2 dx
x =y
3 dy
-2
xdx - ydy = c
3
-x 2 y 2
- =c
3 2
x 2 y2
+ =1
3c 2c
2. Find the O.T of the family of circles x2+y2+2gx+c =0, Where g is the parameter
Sol : x2+y2+2gx+c =0. --------- (1)
represents a system of co- axial circles with g as parameter.
dy
Differentiating with respect to ‘x‘ ⟹ 2x+ 2y + 2g =0--------------(2)
dx
Substituting equation from (2) in (1)
dy
⟹x2+y2 -(2x+ 2y )x +c =0.
dx
dy
⟹y2-x2-2xy + c=0
dx
dy dx
Replace by
dx dy
dx
⟹y2-x2-2xy + c=0
dy
dx
⟹y2-x2+2xy + c=0
dy
This can be written as
dx 1 2 -(c + y 2 )
2x x =
dy y y
This is a Bernoulli’s equation in x
dx du
So put x2=u 2x =
dy dy
du 1 -(c + y2 )
- u=
dy y y
which is a linear equation in ‘u’
1
y dy 1
I.F = e e log y
y
General solution is u(I.F) = Q(y).I.Fdy + k
1 -(c + y2 ) 1
x2 = dy k
y y y
-2
= -c - y + k
y
x2 c
= -y+k
y y
2 2 2
3. Find orthogonal trajectories of the family of curves x + y = a 3 3 3
2 2 2
Sol : The equation of the given family of curves is x 3 + y 3 = a 3 …(1)
-1
-1 -1 3
2 2 dy dy y
Differentiating (1) w.r.t x, we get x + y 3
=03
= -1 …(2)
3 3 dx dx
x3
This is the differential equation of the given family of curves.
-1 -1
3 3
dy dx dx y dx y
Changing to - in (2), we get - = - -1 = -1 …(3)
dx dy dy 3
dy
x x3
1 1
Separating the variables in (3), we get x dx = y dy
3 3
1 1
Integrating, we get x 3dx = y 3dy + c
4 4
3 4 3 4
x 3 = y 3 + c or x 3 - y 3 = c1
4 4
This is the orthogonal equation.
4.Find the O.T’s of the family of parabolas through origin and foci on y –axis.
Sol : The equation of the family of parabolas through the origin and foci on y-axis is
x2=4ay where ‘a’ is parameter.
dy
Differentiating with respective ‘x’ 2x =4a
dx
dy x
=
dx 2a
dx x
O.T ⟹ - =
dy 2a
dx dy
- =
x 2a
y
On integrating, we get -logx = +c
2a
y + 2alogx = c1 is the equation of required O.T.
5.Find the O.T of the one parameter family of curves e x + e-y =c.
Sol : Given equation is e x + e-y =c.
dy
Differentiating with respective ‘x’ e x + e-y (- )=0
dx
dx
Its O.T e x + e-y ( ) = 0 ⟹ e-x dx + e ydy = 0
dy
dr
1.) Differentiating (1) with respect to obtain F [r, θ, ] =0 by eliminating the
dθ
parameter c.
dr dθ
2.) Replace by r2 then the Differential Equation of family of O.T is
dθ dr
dθ
F [ r, θ,- r2 ] =0
dr
3.) Solve the above equation to get the equation of O.T of (1)
Self-orthogonal
A family of curves is said to be self orthogonal when the differential equation of the
family of O.T is same as that of the original family.
Solved Problems
1. Find the orthogonal trajectories of the family of curves r n cosθ = a n .
2. Find the Orthogonal trajectories of the family of cardioids r = a(1- cosθ) where
‘a’ is the parameter.
Sol : Given r = a(1- cosθ) …(1)
Equations (3) and (4) are same, hence the given system is self orthogonal.
x2 y2
4. Show that the system of confocal conics 2 + = 1 where is a
a + λ b2 + λ
parameter,
is self orthogonal.
x2 y2
Sol : Given equation of family of confocal conics is 2 + = 1 …(1)
a + λ b2 + λ
Differentiating (1) w.r.t x, we get
2x 2y dy
+ 2 = 0. For convenience, we write =
a + λ b + λ dx
2
x y
+ 2 p = 0 x(b 2 + λ) + py(a 2 + λ) = 0
a +λ b +λ2
-(b 2 x + a 2 yp)
λ=
x + yp
-(a 2 - b2 )yp
a + λ =
2 (a 2 - b 2 )x and b2 + λ = …(2)
x + yp x + yp
Eliminating from (1) and (2), we get
x(x + yp) y(x + yp) x + yp y
- 2 2 = 1 2 2 x - = 1
a - b p
2 2
a -b (a - b )p
…(3)
y
(x + yp) x - = a 2 - b 2
p
This is the differential equation of the family of curves (1). We get the differential
equation of the family of orthogonal trajectories by replacing
dy dx 1 1
p= with - =- =- .
dx dy dy p
dx
Hence the differential equation of orthogonal trajectories is
y
x - (x + py) = a - b …(4)
2 2
p
which is same as (3). Thus we see that the differential equation of the family of
orthogonal trajectories is same as that of the original family. Hence the given family of
curves is orthogonal to itself. Hence it is a self orthogonal family of curves.
NEWTON’S LAW OF COOLING
STATEMENT: The rate of change of the temperature of a body is proportional to the
difference of the temperature of the body and that of the surrounding medium.
Let be the temperature of the body at time ‘t’ and θ0 be the temperature of its
surrounding medium(usually air). By the Newton’s law of cooling , we have
dθ dθ
α(θ - θ 0 ) = -k(θ - θ 0 ) k is +ve constant
dt dt
dθ
= -k dt
(θ - θ 0 )
c = log(θ1 - θ0 )
θ - θ0
log = -kt.
θ1 - θ 0
θ - θ0 -kt
=e
θ1 - θ 0
θ = θ 0 + (θ1 - θ 0 ) = e -kt
which gives the temperature of the body at time‘t’.
Solved Problems
1. A pot of boiling water 1000 C is removed from the fire and allowed to cool at
30 0 C room temperature. 2 mins later, the temperature of the water in the pot is
90 0 C .What will be the temperature of water after 5 mins?
k = 0.0771
When t = 5,θ - 30 = 70e-5k
θ = 77.460
0
2. A body is originally at 800C and cools down to 60 C in 20 min . If the
temperature of the air is 400 C find the temperature of body after 40 min.
Sol : By Newton’s law of cooling we have
dθ
= -k (θ - θ 0 ) where θ0 is the temperature of the air.
dt
dθ
= -k dt log(θ - θ 0 ) = -kt + logc
(θ - θ 0 )
Here θ 0 = 400 C
log(θ - 40) = -kt + log c
θ - 40
log = -kt
c
θ - 40
=e
-kt
c
θ = 40 + ce-kt --------------(1)
When t=0, θ = 800 C
80 = 40 +c c = 40 ------------(2).
When t=20, θ = 600 C 60 = 40+ ce-20k ------------(3).
Solving (2) & (3) ce-20k =20
40e 2 k =20
k = -20log2
1
- log2 40
When t = 40 C then equation (1) is θ = 40 + 40e
0 20
= 40 +40 e -2log2
1
= 40 + 40X
4
θ = 500 C
dθ
It is given that = -kt
dt
θ = ce -kt ------------(1).
Initially when t=0 θ = 750 - 250 = 500
c= 500
Hence c = 50 θ = 50e-kt -----------------(2)
When t= 10 min θ = 650 - 250 = 400
40= 50 e-10k
4
e -10k = ---------------------(3).
5
The value of θ when t=20 θ = ce -kt
θ = 50e -k
θ = 50(e-10k )2
2
4
50
5
3
log 5
t1 = 10 = 22.9min
log 4
5
4. A body kept in air with temperature 250C cools from 1400C to 800C in 20 min.
Find when the body cools down in 350C.
dθ dθ
Sol : By Newton’s law of cooling = -k(θ - θ 0 ) = -kdt
dt θ - θ0
dθ
Sol : = -k(θ - θ 0 )
dt
log(θ - 20) = -kt + logc
4096
log5 - log80
when θ = 250 C t = 10 = 74.86min
(log11- log18)
6. If the air is maintained at 150 C and the temperature of the body drops from
700 C to 400 C in 10 minutes. What will be its temperature after 30 minutes?
dθ
Sol : If θ be the temperature of the body at time t, then = -k(θ -15) , where k is
dt
constant
dθ
Integrating, θ -15 = -k dt + logc
i.e log(θ -15) = -kt + logc i.e , θ -15 = ce-kt …(1)
When t 0, 70 C and when t 10, 40 C
0 0
70 15 ce 0 c 55 40 15 ce 10 k
and
25 5
= e-10k or e -10k = …(2)
55 11
Then (1) becomes θ -15 = 55e-kt
When t=30 min, θ = 15 + 55e-30k
3
5
θ = 15 + 55(e ) = 15 + 55 using (2)
-10k 3
11
625 2441
= 15 + = = 20.160 C .
121 121
7. In a pot of boiling water 1000 C is removed from the fire and allowed to cool at
30 0 C room temperature. Two minutes later, the temperature of the water in the pot
is 90 0 C . What will be the temperature of the water after 5 minutes?
dθ
= -k(θ - θ 0 ) ;k>0
dt
θ 0 = 240 C log(θ - 24) = -kt + log c--------------(1).
When t=0; θ = 92 c =68
68
When t =1; θ = 800 C e-k =
56
56
k = log .
68
65x41
When θ = 650 C , t = = 0.576min
682
LAW OF NATURAL GROWTH OR DECAY
Statement : Let x(t) or x be the amount of a substance at time ‘ t’ and let the substance be
getting converted chemically . A law of chemical conversion states that the rate of change
of amount x(t) of a chemically changed substance is proportional to the amount of the
substance available at that time
dx
αx
dt
Note: a) In case of Natural growth we take
dx
= kx (k > 0)
dt
dx
b) In case of Natural decay, we take = -kx (k>0)
dt
where k is a constant of proportionality
RATE OF DECAY OF RADIO ACTIVE MATERIALS
Statement : The disintegration at any instant is proportional to the amount of material
present in it.
du
If u is the amount of the material at any time ‘t’ , then = - ku , where k is any
dt
constant (k >0). i.e Law of Natural Decay is applied.
Solved Problems
dN
Sol : The differential equation to be solved is = kN
dt
dN
= kdt
N
dN
N
= kdt
logN = kt + log c
N = c e -kt ------------(1).
When t= 0sec, N =100 100 = c c =100
When t =3600sec, N =332 332 =100 e3600k
332
e3600k =
100
3
Now when t = hours = 5400 sec then N=100 e5400k
2
3
N =100 e 3600k
2
3
332 2
N=100 = 605.
100
N = 605.
2. A bacterial culture, growing exponentially, increases from 100 to 400 gms in 10
hrs. How much was present after 3 hrs, from the initial instant?
1 3
-3 - log2
5 5
= 100e = 100e(log2)
3 1
= 100X(2) 5 = 100X8 5 = 100X1.414
= 141.4gms
3.If a radioactive Carbon-14 has a half life of 5750 years, what will remain of one
gram after 3000years?
Sol : Let mass of radioactive Carbon-14 at any time be denoted by .
dx
Then it is known that = -kt where k is a constant
dt
x = Ae -kt where A is also a constant.
It is known that at t=0, we have 1gm of Carbon-14
1 = Ae0 A = 1
x = e-kt
However when t=5750 years, we have 1/2gm of Carbon-14.
1 1
= e-k(5750) k = log2
2 5750
Suppose t=3000years, we have to find x.
3000
- log2
x = e-kt = e-3000k = e 5750
3000
- gms
x = (2) 5750
4. If 30% of a radioactive substance disappears in 10 days, how long will it take for
90% of it to disappear?
Required time at is
10m1 1
= ce-kt = m1e-kt = e-kt
100 10
1
t = log(10)
k
10log(10)
64.5 days.
log10 log 7
dy d2 y dn y
Dy= ; D2y= 2 ;…………………… Dn y= n
dx dx dx
1
2. Operator Q = Qdx i e D-1Q is called the integral of Q.
D
irrational say
and rest
are real and
distinct.
Solved Problems
d3y dy
1. Solve 3
-3 + 2y = 0
dx dx
Sol : Given equation is of the form f(D).y = 0
Where f(D) = (D3 - 3D + 2)y = 0
Now consider the auxiliary equation f (m) = 0
f(m) = (m3 - 3m + 2)y = 0 (m-1) (m-1) (m+2) = 0
m =1, 1, -2
Since m1 and m 2 are equal and m 3 is -2
2. Solve (D4-2D3-3D2+4D+4)y = 0
Sol : Given f(D) = (D4 -2 D3 - 3 D2 + 4D +4) y = 0 …(1)
Auxiliary equation is f(m)=0
m 4 - 2m3 - 3m 2 + 4m + 4 = 0 …(2)
By inspection m+1 is its factor.
(m +1)(m3 - 3m2 + 4) = 0 …(3)
(m +1) 2 (m - 2) 2 = 0
m = -1, -1, 2, 2
Hence general solution of (1) is
y = (c1 + c 2 x)e-x + (c3 + c 4 x)e 2x
(m2 + 4)2 = 0
(m+2i)2 (m+2i)2 = 0
m= 2i ,2i , -2i , -2i
Here roots are complex and repeated
Hence general solution is
It is given by y = yc + yp
i.e. y = C.F+P.I
Where the P.I consists of no arbitrary constants and P.I of f (D) y = Q(x)
1
Is evaluated as P.I = Qx
f(D)
Depending on the type of function of Q(x) , P.I is evaluated .
1 2
1. Find (x )
D
1 2 x3
Sol : (x ) = x 2dx =
D 3
1
2. Find Particular value of (x)
D +1
1
Sol : (x) = e-x xe x dx (By definition)
D +1
= e-x (xe x - e x )
= x-1
1
P.I of f(D)y=Q(x), when is expressed as partial fractions.
f(D)
Q. Solve (D2 + a 2 )y = secax
Let f(D) = D2 + a 2
1 cosax - isinax
secax = eiax secaxdx = eiax dx
D - ai cosax
i
= eiax (1- itanax)dx = eiax x + logcosax …(4)
a
1 i
Similarly we get secax = e-iax x - logcosax …(5)
D + ai a
Case 2: If f(a)=0 then the above method fails. Then if f D = D - a D (i.e ‘a’ is
k
Type 3.P.I for f(D)y=Q(x) where Q(x)=xk where k is a positive integer .f(D) can be
expressed as f(D) = [1 (D)]
1 1
Express [1 (D)]1
f(D) [1 (D)]
1
Hence P.I = Q(x)
[1 (D)]
= [1 (D)]1 x k
Type 4.P.I of f(D)y=Q(x) when Q(x)=eax V where ‘a’ is a constant and V is function of
x. where V =sinax or cosax or xk
1
Then P.I = Q(x)
f(D)
1 ax
= e V
f(D)
1 1
= eax V & V is evaluated depending on V.
f(D + a) f(D + a)
Type 5. P.I of f(D)y=Q(x) when Q(x)=xV where V is a function of x.
1
Then P.I = Q(x)
f(D)
1
= V
f(D)
1 1
= x - f'(D) V
f(D) f(D)
Type 6. P.I. of f(D)y = Q(x) where Q(x)= xm v where v is a function of x.
1 1 m
When P.I. = ×Q(x) = x v, where v = cosax or sinax
f(D) f(D)
1 m 1 m iax
i. P.I. = x sinax = I.P.of x e
f(D) f(D)
1 m 1 m iax
ii. P.I. = x cosax = R.P.of x e
f(D) f(D)
Formulae
1. = (1 – D)-1 = 1 + D + D2 + D3 + ------------------
2. = (1 + D)-1 = 1 - D + D2 - D3 + ------------------
Solved Problems
1. Solve (4D2 - 4D +1)y = 100
1 -1
Sol : A.E is 4m 2 - 4m +1 = 0 (2m -1) 2 = 0 m = ,
2 2
x
C.F = (c1 + c2 x)e 2
Sol : A.E is m 2 + m +1 = 0
-1± 1- 4 -1± 3i
m= =
2 2
-x
x 3 x 3
yc = e 2 c1cos + c 2sin …(1)
2 2
To find P.I :
sin2x sin2x
yp = 2 =
D + D +1 -4 + D +1
sin2x (D + 3)sin2x (D + 3)sin2x
= = =
D-3 D2 - 9 -4 - 9
Dsin2x + 3sin2x 2cos2x + 3sin2x
= =
-13 -13
-x
x 3 x 3 1
y = y c + y p = e 2 c1cos + c 2sin - (2cos2x + 3sin2x)
2 2 13
x +3 1 x3 + 3
yp = 3 =
2
(D - 3D -10D + 24) 24 D3 - 3D 2 -10D
1+
24
-1
1 1+ D3 - 3D 2 -10D
24 (x + 3)
=
24
1 D3 - 3D 2 -10D
= 1- (x + 3)
24 24
1 10 24x + 82
= x +3+ =
24 24 576
General solution is y = yc + y p
24x + 82
y = e 2x (c1 + c 2 x) + c3e -x +
576
8. Solve the differential equation (D2 - 4D + 4)y = e2x + x2 + sin3x .
2D sinx
= x - 2 2
D + 4 D + 4
xsinx 2(Dsinx)
= -
3 3(D 2 + 4)
xsinx 2cosx
= -
3 9
Hence the general solution is
1 2
y = y c + y p = c1cos2x + c 2sin2x + xsinx - cosx
3 3
m=-2 or m=-3
The roots are real and distinct
C.F=yc= c1e-2x + c 2e-3x
1
Particular Integral=yp= Q(x)
f(D)
1 1
= 2
ex = ex
D + 5D + 6 (D + 2)(D + 3)
Put D = 1 in f(D)
1
P.I= ex
(3)(4)
1 x
Particular Integral = yp = e
12
General solution is y = yc + y p
ex
y = c1e-2x + c 2e-3x +
12
11. Solve y''- 4y' + 3y = 4e3x , y(0) = -1, y'(0) = 3
d 2 y dy
i.e - 4 + 3y = 4e3x it can be expressed as
dx 2 dx
D2 y - 4Dy + 3y = 4e3x
1
P.I= y p = Q(x)
f(D)
1
= 2
4e3x
D - 4D + 3
1
= 4e3x
(D -1)(D - 3)
Put D=3
given y'(0) = 0
A.E is m3 + 2m 2 - m - 2 = 0
(m2 -1)(m + 2) = 0
m2 = 1 or m = -2
m = 1, -1, -2
C.F= c1e x + c2 e-x + c3e-2x
1 -1
P.I = (1- 4x 3 ) =
3 2
(D + 2D - D - 2) (D + 2D 2 - D)
3
2 1- (1- 4x 3 )
2
-1
-1 (D3 + 2D 2 - D)
= 1- (1- 4x 3 )
2 2
-1 5 3 5 2 1
= 1- D + D - D (1- 4x 3 )
2 8 4 2
-1 5 5 1
= (1- 4x 3 ) - (-24) + (-24x) - (-12x 2 )
2 8 4 2
-1
= -4x 3 + 6x 2 - 30x +16
2
= [2x 3 - 3x 2 +15x -8]
The general solution is
y= C.F + P.I
y = c1e x + c 2e-x + c3e-2x +[2x 3 - 3x 2 +15x - 8]
ex cos2x
P.I=
(D3 - 7D2 +14D -8)
1
= ex cos2x
(D +1) - 7(D +1) 2 +14(D +1) -8
3
1 ax 1
P.I = e v = eax v
f(D) f D + a
1
= ex cos2x
(D - 4D2 + 3D)
3
1
= ex 3 cos2x
(D - 4D2 + 3D)
1
= ex cos2x (Replacing D2 with -22)
(-4D + 3D +16)
1
= ex cos2x
(16 - D)
16 + D
= ex cos2x
(16 - D)(16 + D)
16 + D
= ex cos2x
256 - D 2
16 + D
= ex cos2x
256 - (-4) 2
ex
= (16cos2x - 2sin2x)
260
2ex
= 8cos2x - sin2x
260
ex
= 8cos2x - sin2x
130
General solution is y = yc + yp
ex
y = c1e x + c 2e2x + c3e 4x + 8cos2x - sin2x
130
Variation of Parameters :
This method is applied when P,Q in above equation are either functions of ‘x’ or real
constants but R is a function of ‘x’.
Working Rule :
1. Find C.F. Let C.F= y c = c1u(x) + c2 u(x)
vRdx uRdx
2. Take P.I= y p =Au+Bv where A= - and B =
uv' - vu' uv' - vu'
3. Write the G.S. of the given equation y = yc + y p
d2y
1. Apply the method of variation of parameters to solve + y = cosecx
dx 2
Sol : Given equation in the operator form is (D2 +1)y = cosecx …(1)
A.E is (m2 +1) = 0
m = ±i
The roots are complex conjugate numbers.
C.F is yc = c1cosx + c2sinx
Let yp= A cosx+ B sinx be P.I. of (1)
dv du
u -v = cos 2 x + sin 2 x = 1
dx dx
A and B are given by
vRdx sinxcosecx
A= - = - dx = - dx = -x
uv' - vu' 1
uRdx
B= 1
uv - vu1
= cosx.cosecxdx = cotxdx = log(sinx)
y p = -xcosx + sinx.log(sinx)
Sol : A.E is m 2 - 2m + 2 = 0
2 ± 4 - 8 2 ± i2
m = = = 1± i
2 2
UNIT-IV
Ex: + = +
Now eliminate the constants a and b from (1), (2) and (3). We get a partial differential
equation of the first order of the form. x, y, z, p, q 0
Note : 1. If the number of arbitrary constants is equal to the number of variables, a partial
differential equation of first order can be obtained.
2.If the number of arbitrary constants is greater than the number of variables, a partial
differential equation of order higher than one can be obtained.
Solved Problems
Putting the values of a and b from equation (2) and (3) in (1), we get
= + +
This is the required partial differential equation
2. Form the partial differential equation by eliminating the arbitrary constants a and b
a
from (a) z ax by a 2 b 2 (b) z ax by b
b
Sol. (a) we have z ax by a 2 b 2 ....... (1)
Differentiating (1) partially w.r.t. ‘x’ and ‘y’, we get
z z
a p a ……… (2) and b q b ………..(3)
x y
Putting the values of a and b from equation (2) and (3) in (1), we get
z px qy p 2 q 2
x a y b z2 r 2
2 2
(OR )
Find the differential equation of all spheres of fixed radius having their centre on the
xy-plane.
Sol. The equation of sphere of radius r having their centers on xy-plane is
x a y b z 2 r 2 ……….(1)
2 2
zp zq z 2 r 2
2 2
P= = 1. + −−−−−−−−−−
q= = 1. + −−−−−−−−−
from (2 = +
= +
Substituting in (1) we get
= .
This is the required partial differential equations
5. Form the partial differential by eliminating the arbitrary constants from
log az 1 x ay b
q
q z 1or pq qz p or p q 1 q 2
p
is the required partial differential equation.
1 1
6.Form the differential equation by eliminating a and b from 2z x a 2 y a 2 b
1 1
Sol: We have 2z x a 2 y a 2 b …………..(1)
1
yq ………… (3)
16q 2
Adding (2) and (3), we get
1 1 1
x y 2 2
16 p q
or 16 x y p2q2 p2 q2
3z px qy
8.Form the partial differential equation by eliminating the arbitrary constants a and b
from = + +
Sol:-The given equation = + + −−−−−−−−−−
Differentiating (1) w.r.t.,
P= = + −−−−−−−−−−
∴ + =
q= = 2y + −−−−−−−−−
∴ + =
− =
9.Form the partial differential equation by eliminating the arbitrary constants from
− + − = 𝜶
Sol:Given
− + − = 𝜶……..(1)
Differentiating (1) w.r.t.,
− = 𝜶
Differentiating (1) w.r.t.,
− = 𝜶
Substituting (2),(3) in (1),we get
𝜶 + 𝜶 = 𝜶
∴The required Partial differential equation is
+ = 𝐚 𝜶
u u v v
i.e., p p 0 ……(2)
u x z v x z
4 4 v v
and q q 0 ……(3)
4 y z v y z
Eliminating and from (2) and (3)
u v
u u v v u u v v
p q q p
y z y z y z x z
u v u v u v u v u v u v
i.e. p q
y z z y z x x z x y y x
u, v u, v u, v
p q
y, z z, x x, y
1 1 p2 (4)
. zr 0
a2 c2 c2
Multiplying (4) by ‘x’ and then subtracting (2) from it, we get
xp 2 z
.r 2 2 . p 0 or 2 xzr xp 2 zp 0
xz 1
2
c c c c
pz xp 2 xzr
is the required partial differential equation.
2. Form a partial differential equation by eliminating the arbitrary
function𝝋 + , − =
Sol:-Given 𝜑 + , − =
This can be written as − = f( + --------------(1)
Now we have to eliminate f from (1)
Differentiating (1) w.r.t., x
′
− = +
′
− = + --------(2)
Differentiating (2) w.r.t., y
′
− = + --------(3)
Dividing (2) by (3)
py–qx= −
is the required partial differential equation.
3. Form a partial differential equation by eliminating the arbitrary function
from z f x 2 y 2
Similarly we get
=− (u) (4)
3 4 , gives =
−
∴ + =
2 3 , gives
yz xyp 2 x 2 zp
xz xyq 2 y 2 zq
yz xyp y zq xz xyq x zp
y 2 z z 2 yq xy 2 p xyzpq x 2 z x 2 zp x 2 yq xyzpq
x y2 z 2 p y z 2 x2 q x2 y 2 z
Y xp + z + q = x yq + z ( +
− + − = −
− + − = −
z x y . p z
y y. p z x. p f 1 (2)
x y x y
2
z x y q z
And x z yq xq f 1 (3)
x y x y
2
x y p y z x y p z
x y q x z x y q z
from z f x e y .g x
2 z y z
e .g x [using (3)]
y 2
y
2 z z
2 0
y y
t q 0
z u
And f 1 u . f 1 u .2 y
y y
q f 1 u 2 y (4)
p f u .2 x x
1
3 4 , gives 1
q f u 2 y y
py qx 0
+ = + + + ….
implies
+ +
=
+ +
Let x, y, z, a, b 0 (2)
1
Differentiating (5) partially w.r.t. ‘a’, we get . f a 0 (6)
a f
Eliminate ‘a’ between (5) and (6), when it exists is called the general integral or general
solution of (1).
(i) Method of grouping : In some problems, it is possible that two of the equations dx = dy or
P Q
From this we may express y, as a function of x. Using this in dy dz and integrating we get
Q R
, = . ℎ 𝑖 = , = 𝑖 ℎ 𝑖
2. Method of multipliers: This is base on the following elementary result.
Consider dx dy dz
P Q R
Where + + = Then + + =
Integrating this we get , , = .
Similarly we get another solution , , = 𝑖ndependent of the earlier one.
We have the complete solution of (1) constituted = =
sin x sin x
or log log c1 or c1 (2)
sin y sin y
dy dt
Taking the last two members, we have
tan y tan z
Integrating, 𝑖 = 𝑖 +
sin y sin y
or log log c2 or c2 (3)
sin z sin z
From (2) and (3). The general solution of (1) is
c1 , c2 0
sin x sin y
i.e. , 0
sin y sin z
is the required Complete Solution.
2. Find the general solution of y 2 zp x 2 zq y 2 x
x3 x3 y 3
Integrating, y 3 c1 or c1 (2)
3 3 3 3
x2 z 2 x2 z 2
Integrating c2 or c 2 (3)
2 2 2 2
From (2) and (3) The general solution of (1) is
c1 , c2 0 i.e.
x3 y 3 x 2 z 2
, 0
3 3 2 2
is the required Complete Solution.
3. Solve p x q y z
Sol. The given equation can be written as
x p yq z (1)
Comparing with Pp+Qq=R, we have
P x,Q y, R z
dx dy dz
The auxiliary equations are
x y z
dx dy
From the first two members, we have
x y
Integrating, 2 x 2 y c1 or 2 x 2 y c1 or x y a (2)
dy dz
From the last two members, we have
y z
Integrating, 2 y 2 z c2 or 2 y 2 z c2
or y z b (3)
From (2) and (3). The general solution of (1) is
a, b 0 i.e.,
x y, y z 0
is the required Complete Solution.
4. Solve x y z p y z x q z x y
P x y z , Q y z x , R z x y
dx dy dz
The auxiliary equations are
x y z y z x z x y
dx dy dz 0
Integrating, x y z a (2)
1 1 1
Again using l , m , n as multipliers, we get
x y z
1 1 1
dx dy dz
x y z
Each fraction =k(say)
0
1 1 1
dx dy dz 0
x y z
Integrating, logx+logy+logz=logb. or xyz=b……. (3)
From (2) and (3). The general solution of (1) is
a, b 0 i.e.,
x y z, xyz 0
is the required Complete Solution.
5. Solve x2 y z p y 2 z x q z 2 x y
dx dy dz
The auxiliary equations are 2 2
x y z y z x z x y
2
1 1 1
Using l 2
, m 2 , n 2 as multipliers, we get
x y z
1 1 1
2
dx 2 dy 2 dz
x y z
Each fraction =k(say)
0
1 1 1
2
dx 2 dy 2 dz 0
x y z
1 1 1 1 1 1
Integrating , a or c1 (2)
x y z x y z
1 1 1
Again using l , m , n as multipliers, we get
x y z
1 1 1
dx dy dz
x y z
Each fraction =k(say)
0
1 1 1
dx dy dz 0
x y z
Integrating log x log y log z log c2
or xyz c2 (3)
1 1 1
, xyz 0
x y z
is the required Complete Solution.
6. Solve mz ny p nx lz q ly mx
x2 y 2 z 2
Integrating, a or x 2 y 2 z 2 c1 (2)
2 2 2
Again using , , as multipliers, we get
ldx mdy nz
Each fraction =k(say)
0
ldx+mdy+ndz=0
Integrating, lx my nz c2 (3)
Sol. Here P x, Q y, R y 2 x 2
dx dy dz
The auxiliary eqn’s are 2
x y y x2
dx dy
From the first two members,
x y
Integrating, log x log y log c1 or xy c1 (1)
Using l=x, m=y, n=1 as multipliers, we get
xdx ydy dz
Each fraction
0
xdx ydy dz 0
1 2 1 2
Integrating, x y z c or x 2 y 2 2 z c2 (2)
2 2
From (1) and (2), The general solution is
xy, x 2 y 2 2 z 0
1 1 1
Using l , m , n as multipliers, we get
x y z
1 1 1
dx dy dz
Each fraction = x y z
0
1 1 1
dx dy dz 0
x y z
Integrating, log x log y log z log a
or xyz a (2)
Again using l=x, m=y, n=-1 as multipliers, we get
xdx ydy dz
Each fraction =k(say)
0
xdx ydy dz 0
x2 y 2
Integrating, z c or x 2 y 2 2 z b (3)
2 2
Given that z=1, using this (2) and (3), we get
xy=a and x 2 y 2 2 b
2a b 2 0
Hence the required surface is
x 2 y 2 2 z 2 xyz 2 0
is the required Complete Solution.
9.Solve + =
Sol:Given
+ = is a Lagrange’s linear equation
The Auxillary equations are
= =
∫ =∫
∫ =∫
= …..(2)
∴ , =
10. Solve − − + − − = −
= =
− − − − −
− − =
Integrating,we get
− − = ………(1)
Taking , − , as multipliers,we get
−
=
− − − + + −
−
=
− − −
𝐥 − =𝐥
−
= ………………… )
−
∴ 𝝋 − − , =
11. Solve − − + = +
= =
− − + +
+ + =
+ + = ………………….
- + + =
Integrating, we get
= …….(2)
From (1),(2),
∴ 𝝋 , + + =
12.Solve − = −
Sol:Comparing with Pp+Qq=R, we have
=
−
Integrating,we get
+ = …..(2)
=
− −
− = −
− =
− = ………..(3)
From (2) and (3). The general solution of (1) is
i.e., ∅ − , + =0
12. Solve + + + = +
+ + − = ………..(2)
From (2) and (1). The general solution of given pde is
−
i.e., ∅ , + + − =0
−
13. Solve − = −
=
−
Integrating,we get
+ = …..(1)
Integrating,we get
+
= ….
14. Solve − + − = −
( )
=( )
=
− − −
− −
∴( − − − )
=
(− + + − )
− −
− + +
= − + +
𝑖 , , , , multipliers,we get
+ + + +
+ + −
= + + − − −
+ + + +
+ + + + − − −
=
+ + − − −
+ + + + = + +
+ + + + = + +
Integrating,we get
+ + + +
= +
∴ + + = + + +
+ + = …..(2)
−
∴ 𝑖 𝑖 𝑖 𝜑 + + , =
−
= = = =
− − − − [ + ] [ + ]
Solved Problems
1.Solve + =
Sol: Given = + −
The auxillary equations are
= = = =
− − − − [ + ] [ + ]
Here = ; = ; = ; = ; =
Substituting them in above ,we get
= = = =
− − +
By considering first and last groups,
=
− −
But we have
= +
+
= + +
= + +
2.Solve =
Sol: Given = −
= = = =
− − − − [ + ] [ + ]
Here =− ; =− ; =− ; =− ; =
Substituting them in above ,we get
= = = =
− + − +
Considering
=− +
and =− +
Taking
, 𝑖 𝑖 𝑖 , , ℎ 𝑖 ℎ
+ +
− + +
= − + +
Solving ,
Integrating,we get
log = log +log
=
Substituting in given pde,we get
=
∴ = ; = √
√
But we have
= +
= √ +
√
∫ = ∫ √ +∫
√
Integrating,we get,
∴ Required solution is
= √ √
Now let us start solving some standard forms of first order partial differential equations by
using Charpit’s method
STANDARD FORM I:
Equation of the form f(p,q)=0
i.e., equations containing p and q only.
Given partial differential equation is , = ……
The auxillary equations are
= = = =
− − − − [ + ] [ + ]
Here = ; = ; =
Substituting above and considering last group, we get
= = = =
− − − −
∴ = ,𝑖 𝑖 =
Put = , , 𝝋 .
But we have
= +
= +
= + 𝝋 ……(2)
: , 𝑖 = +
. = +
STEP3:Integrating it ,we get required complete solution of (1) .
Solved Problems
1.Solve = , where k is a constant.
Sol. Given that = ….
Since (1) is of the form , =
But we have
= +
= +
Integrating,we get ,
k
z ax yc
a
which contains two arbitrary constants a and c.
2.Solve + =
Sol : Given that + = … … … … . . 1)
But we have
= +
= + [ ±√ − ]
Integrating,we get , = ∫ + [ ±√ − ]∫
a
z ax n n2 4 y c
2
This is the complete integral of (1), which contains. two arbitrary constants a and c.
3.Find the complete integral of + =
Sol. Given that + = ……….
Since (1) is of the form , =
Put = 𝑖 , =√ −
But we have
= + …….(2)
Put the values of p,q in (2), we get
z ax m2 a 2 y c
Which is the complete integral of (1)
STANDARD FORM II :
Equation of the form , , = (i.e., not containing x and y)
PROCEDURE
Given partial differential equation is , , = ……
STEP1:Put = 𝑖 , ℎ 𝑖 , .then
: , 𝑖 = +
. = +
STEP3:Integrating it ,we get required complete solution of (1) .
Solved Problems :
Solve the following partial differential equations
1. = + 2. + = 3. = +
Put = 𝑖 , ℎ =√
+
∴ = √
+
Putting the values 𝑖 = + , we get
= + ,
√ √ +
∫ = ∫ +
√ √ +
∴ √ = +
√ +
+
∴ =
+
= +
Integrating,we get
∫ = ∫ +
+
∴ a tan− = + +
which is the required complete solution of (1)
3. Given that = + …..(1)
Since (1) is of the form , , =
+
Put = 𝑖 , ℎ =
+
∴ =
∴ +
= + +
= +
dz F1 x, k dx F2 y, k dy
Let p 2 x q 2 y k (constant)
p 2 x k and q 2 y k
p 2 k x and q 2 y k
p k x and q y k
z z
Since dz dx dy pdx qdy
x y
dz k x dx y k dy
Integrating on both sides
1 1
z k x 2 dx y k 2 dy c
3 3
2 2
z k x2 y k 2 c
3 3
xp x 2 k and yq y 2 k
k x2 k y2
p and q
x y
z z
We have dz dx dy pdx qdy
x y
k k
dz x dx y dy
x y
Integrating on both sides
k k
z x dx y dy c
x y
x2 y2
k log x k log y c
2 2
z k log xy
1 2
2
x y2 c
2 2
p q
3. Solve x y 1
2 2
Sol: Separating , the given equation can be written as.
2 2
p q
x 1 y
2 2
2 2
p q
Let x 1 y k 2 (arbitrary constant)
2 2
2 2
p q
x k 2 and 1 y k 2
2 2
2
p q q
x k and y 1 k 2 or y 1 k 2
2 2 2
p 2 k x and q 2 1 k 2 y
z z
We have dz dx dy pdx qdy
x y
dz 2(k x)dx 2 1 k 2 y dy
Integrating on both sides
z 2 (k x)dx 2 1 k 2 y dy c
z 2(kx
x2
2
y2
) 2 1 k 2 y c
2
z 2kx x 2 2
1 k 2 y y2 c
But we have
z z
dz dx dy pdx qdy
x y
Integrating,we get
= + + + +c
5. Solve − = −
Sol: Let − = - =
Then = + =√ +
But
z z
dz dx dy pdx qdy
x y Integrating,we get
= + + + +C
But
z z
dz dx dy pdx qdy
x y Integrating,we get
−
=− + [ √ + + )+k +C
√
Solved Problems:
1. Solve = + +
Sol : The given PDE is form IV
Therefore complete solution is given by
= + +
2. Find the solution of p q z px qy 1
3. Solve pqz p 2 qx p 2 q 2 py q 2
p2 q2
z p x q y
q p
p3 q3
z px qy (1)
q p
4. Solve = + + +
Sol. We have = + + + …………….
Since (1) is of the form = + + , .
Hence the complete solution of (1) is given by
= + + + …….
For singular solution, differentiating (2) partially w.r.t. a and b, we get
= , = ,
Implies that
= + …. = + + … … … … 4)
Eliminating a, b between (2), (3) and (4), we get
= − − − − +
∴ =
is the singular solution
EQUATIONS REDUCIBLE TO STANDARD FORMS:
𝑖 𝑖 = ℎ =
now the given equation reduces to the form , =
EQUATIONS OF THE FORM , , = :
This can be reduced to an equation of the form , , = by the substitutions given for the equation
, , = as above.
Solved Problems:
x2 y 2
1. Solve the partial differential equation z
p q
Sol. Given equation can be written as
x 2 p 1 y 2 q 1 z or x 2 p y 2 q z (1)
1 1
y 2 q 3Q
Now equation (1), becomes.
3P 3Q z (2)
1 1
+
= +
Integrating,we get
∫ = ∫ +∫
+
= + +
+
a 1 3
3z 2 2 x ay c1
3
a
+
,taking =
Which is the required solution of (1)
p q
2. Solve the partial differential equation z
x2 y 2
Sol. The given equation can be written as
px 2 qy 2 z (1)
Since (1) is of the form f x m p, y n q, z 0 With =− , =−
Put X x1m x3 , and Y y1n y 3
z z X z
Now p . P.3x 2 where P
x X x X
2
x p 3P
z z Y z
and q . Q3 y 2 where Q
y Y y Y
y 2 q 3Q
Equation (1) becomes, 3P 3Q z (2)
Since (2) is of the form f P, Q, z 0
Put = 𝑖 , ℎ =
+
∴ =
+
= +
+
Integrating,we get
∫ = +
∫ +∫
log = + +
+
log z
1
3 1 a
x3 ay 3 c
∴ = [−a ± √a + ]
px qy z 2z
2 2 2 2
(1)
∴ = [ ±√ a + ]
Putting the values of Pand Q in = + ,we get
= [ ±√ a + ] +
Integrating,we get
∫ = [ ±√ a + ] ∫ +∫
log = [ ±√ a + ] + +
∴ log = [ ±√ a + ] + +
= , P=a
√ + √ +
But we have
= +
Substuting P,q ,we get
= +
√ +
√ + log = + +
be the complete integral of (1)
6.Solve z p 2 x q 2 y
Sol. Given that z p 2 x q 2 y
The given equation can be written as
2 2
1 1
p x q y
2 2
z or px 2 qy 2 z (1)
=√ , P=a√ ,
+ +
But we have
= +
Substuting P,Q ,we get
= √ +
+
= +
√ √ +
√ + √ = + +
√ + √ =( √ + √ ) +
Which is the complete integral of (1)
7.Solve + =
Sol: Given + = ……..
+ =
Since (1) is of the form f x p, y q, z 0 with
m n
= =
Put X log x and Y log y
z z X 1 z
Now p . P. where P
x X x x X
xp P
z z Y 1 z
and q . Q. where Q
y Y y y Y
qy Q
Equation (1), becomes
+ = ……(2)
= 𝑖 , = ; =
√ + √ +
But we have
= +
Substuting P,Q ,we get
= +
√ +
= +
√ +
√ + log = + +
√ + log = log + log +
is the Complete solution of (1)
8. Solve + =
+ = …… 2
= 𝑖 , =√ −
But we have
= +
Substuting P,Q ,we get
= +√ −
z =( +√ − )+
z=( log + √ − log ) +
is the Complete solution of (1)
Solved Problems :
1. Solve z 2 p 2 q 2 x 2 y 2
dZ 2 x 2 k 2 dx 2 y 2 k 2 dy
Integrating on both sides
Z 2 x 2 k 2 dx 2 y 2 k 2 dy
x 2 k2 x x k2 y
2 x k 2 sinh 1 2 y 2 k 2 cosh 1 c
2 2 k 2 2 k
x y
x x 2 k 2 k 2 sinh 1 x y 2 k 2 k 2 cosh 1 c
k k
x y
or z 2 x x 2 k 2 y y 2 k 2 k 2 sinh 1 cosh 1 c
k k
x x2 k 2
or z x x k y y k k log
2 2 2 2 2 2
c
y y2 k 2
Put Z z n1 z 2
Z z P
Now 2 z. P 2 zp or pz where = ; =
x x 2
Z z Q
and 2 z. Q 2 zq or qz
y y 2
2 2
P Q
Then equation (1) becomes, sin 2 x 1 cos 2 y
2 2
2 2
P Q
i.e. sin 2 x 1 cos 2 y (2)
4 4
This is of the form f1 x, p f 2 y, q
2 2
P Q
Let sin 2 x 1 cos 2 y k 2 (constant)
4 4
2 2
P Q
sin 2 x k 2 and 1 cos 2 y k 2
4 4
P 2 sin 2 x 4k 2 and Q 2 cos 2 y 4 1 k 2
2k 2 1 k 2
P and Q
sin x cos y
Z Z
We have dZ dx dy [By total differential]
x y
dZ Pdx Qdy
2k 2 1 k 2
dZ dx dy
sin x cos y
Integrating on both sides
= ∫ csc + √ − ∫ sec
𝑖 𝑖 ℎ , , , = =
+
Put = =
Differentiating partially w.r.t ‘x’,we get = 𝑖 𝑖 ℎ =
But = = implies = = ; 𝑖 𝑖 =
Substitute in (1),we get
+ + + =1
Separating , the given equation can be written as.
+ = − + =
+ = − + =
+ =
Implies that
= √ − )−
= −
Z Z
dZ dx dy
We have x y
dz 2(k x)dx 2 1 k 2 y dy
Integrating on both sides
z 2 (k x)dx 2 1 k 2 y dy c
x2
y2
z 2(kx ) 2 1 k y c
2
2
2
z 2kx x 2 2
1 k 2 y y2 c
𝑖 𝑖 ℎ , , , = =
+
Put = =
Differentiating partially w.r.t ‘x’,we get =
𝑖 𝑖 ℎ =
− = −
Separating , the given equation can be written as.
− =− + =
Solving, we get
= √ + and = √ +
Z Z
We have dZ dx dy [By total differential]
x y
dZ Pdx Qdy
= [√ + +√ + ]
Integrating on both sides
= [∫ √ + +∫√ +
= + + + +
−
1. Solve by the method of separation of variables = + where U(x,0)=6
Solution:- Given equation is = + ----------(1)
Let U(x,t)= X(x) T(t) =XT ---------------(2)
be a solution of (1)
Differentiating (2) partially w.r.t x ant t
′ ′
= , =
Put these values in equation (1), we have
′ ′
= + Dividing by XT
′ ′
= + --------(3)
Since L.H.S is a function of ‘x’ and the R.H.S is a function of ‘t’ where x and t are
independent variables, the two sides of (3) can be equal to each other for all values of ‘x’ and
‘t’ if and only if both sides are equal to a constant.
′ ′
Therefore = + =k--------(4) where k is a constant
′ ′
Now from (4) = k----(5) and + =k-------(6)
X kX 0
′
Now consider (5) = k X C1ekx
k 1
k 1 t
Now consider (6)
′
+ =k T T 0 T C 2 e 2
(8)
2
Substituting the values of X and T in (2) we get
k 1
t
U ( x, t ) X C1e C2e kx 2
k 1
kx 2
t
U ( x, t ) X Ae e ( where A=C1 C2 )
kx
Put t=0 in the above equation ,we have U(x,0) = A e -----(9)
but given that U(x,0)=6 − --------(10)
from (9) and (10) we have A e kx =6 −
A=6 and k=-3 the solution of the given equation becomes
U ( x, t ) X 6e 3 x e 2 t 6e (3 x 2t )
′ ′ ′′
= , = =
Put these values in equation (1), we have
′′ ′
= +
′′ ′′
Dividing by XY on both sides we have = +2
′′ ′′
− = -----------------(3)
Since L.H.S is a function of ‘x’ and the R.H.S is a function of ‘y’ where x and y are
independent variables, the two sides of (3) can be equal to each other for all values of ‘x’ and
‘y’ if and only if both sides are equal to a constant.
′′ ′′
− = = --------(4)
Now from (4)
From (5) ′′ − = ′′
− + =
Which is second order differential equation
Auxiliary equation is − + = → m=±√ +
Solution of the given equation (5) is = √ +
′
′
Now consider equation (6) = → =
Integrating on both sides we get logy = ky + log
Y
log ky Y C3eky (8)
C3
Substituting the values of X and Y in (2) we have
U C1e (2 k ) x C2 e (2 k ) x C3e ky
U Ae (2 k ) x Be (2 k ) x e ky
UNIT V
LAPLACE TRANSFORMS
INTRODUCTION
Laplace Transformations were introduced by Pierre Simmon Marquis De Laplace
(1749-1827), a French Mathematician known as a Newton of French. Laplace
Transformations is a powerful technique, it replaces operations of calculus by operations of
algebra. An Ordinary (or) Partial Differential Equation together with Initial conditions is
reduced to a problem of solving an Algebraic Equation by this method.
USES
Particular Solution is obtained without first determining the general solution.
Non-Homogeneous Equations are solved without obtaining the complementary
integral.
L.T is applicable not only to continuous functions but also to piecewise continuous
functions, complicated periodic functions, step functions and impulse functions.
APPLICATIONS:
L.T is very useful in obtaining solution of linear differential equations, both ordinary
and partial, solution of system of simultaneous differential equations, solution of
integral equations, solution of linear difference equations and in the evaluation of
definite integrals.
DEFINITION:
Let f (t) be a function of‘t’ defined for all positive values of t. Then Laplace
transforms of f (t) is denoted by L {f (t)} is defined by
L f t e st f t dt f s (1)
0
provided that the integral exists. Here the parameter‘s’ is a real (or) complex number.
The relation (1) can also be written as f t L1 f s
In such a case the function f(t) is called the inverse Laplace transform of f s .The
symbol ‘L’ which transform f(t) into f s is called the Laplace transform operator. The
finite quantity.
Ex: (i). The function t2 is of exponential order
3
(ii). The function et is not of exponential order (which is not finite quantity)
Piece – wise Continuous function: A function f (t) is said to be piece-wise continuous over
the closed interval [a,b] if it is defined on that interval and is such that the interval can be
divided into a finite number of sub intervals, in each of which f (t) is continuous and has both
right and left hand limits at every end point of the subinterval.
Sufficient conditions for the existence of the Laplace transform of a function:
The function f (t) must satisfy the following conditions for the existence of the L.T.
(i).The function f(t) must be piece-wise continuous (or sectionally continuous) in any limited
interval 0 a t b .
(ii).The function f (t) is of exponential order.
Laplace Transforms of standard functions:
1
1. Prove that L 1
s
Proof: By definition
e st e e0
L 1 e .1dt
st
s s 0 s if s 0
1
0 s 0
L 1 1
s e
0
Prove that L t n
n!
3. where n is a +ve integer
s n 1
e st n 1 e st
Proof: By definition L t
n.t . s dt
st
n
e .t dt t n .
n
0
s 0 0
n st n 1
s 0
00 e t dt
L t n 1
n
s
n 1
Similarly L t n 1 L t n 2
s
n2
L t n 2 L t n 3
s
By repeatedly applying this, we get
n n 1 n 2
L t n . ..... . L t n n
2 1
.
s s s s s
n! n! 1 n!
n
L 1 n . n 1
s s s s
Note: L t n can also be expressed in terms of Gamma function.
n ! n 1
i.e., L t n n 1 n!
s n 1 s n 1
Def: If n>0 then Gamma function is defined by n e x x n1dx
0
We have L t n e st .t n dt
0
1 When t 0, x 0
s n 1 0
e x .x n dx
When t , x
L t n
1
. n 1
s n 1
If ' n 'is a +ve integer then n 1 n !
L t n n !
s n 1
Note: The following are some important properties of the Gamma function.
1. n 1 n. n if n 0
2. n 1 n! if n is a +ve integer
3. 1 1, 1
2
2 1 1 1!
3 2 2 2!
4 3 3 3!
and so on.
In general n 1 n! provided ‘n’ is a +ve integer.
Prove that L e at
1
4.
sa
Proof: By definition,
0 0
e s a t
s a 0
e e0 1
if s a
sa sa sa
Similarly L e at
1
if s a
sa
a
5. Prove that L sinh at
s a2
2
1 1 1 1 s a s a 2a a
2 s a s a 2 s 2 a 2 2 s 2 a 2 s 2 a 2
s
6. Prove that L cosh at
s a2
2
eat e at
Proof: L cosh at L
2
L e at L e at
1 1 1 1
2 2 s a s a
1 s a s a 2s s
2 s 2 a 2 2 s 2 a 2 s 2 a 2
a
7. Prove that L sin at
s a2
2
Proof: By definition, L sin at e st sin atdt
0
eax
2
e sin bxdx
ax
a sin bx b cos bx
a b
2
a
s a2
2
s
8. Prove that L cos at
s a2
2
s ia
i.e., L cosat i sin at
s2 a2
Equating the real and imaginary parts on both sides, we have
s a
L cos at and L sin at 2
s a
2 2
s a2
Solved Problems :
1. Find the Laplace transforms of (t 2 1)2
1 1 1 1
a s a s s(s a)
3. Find the Laplace transform of Sin2tcost
1 1
Sol: W.K.T sin 2t cos t [2 sin 2t cos t ] [sin 3t sin t ]
2 2
1 1
L{sin 2t cos t} L [sin 3t sin t ] L sin 3t L sin t
2 2
1 3 1 2( s 2 3)
2 s 2 9 s 2 1 ( s 2 1)( s 2 9)
4. Find the Laplace transform of Cosh22t
1
Sol: w.k.t cosh 2 2t 1 cosh 4t
2
1 1 s s2 8
2 s s 2 16 s( s 2 16)
5. Find the Laplace transform of Cos33t
Sol: Since cos9t=cos3(3t)
1
cos9t=4cos3 3t-3cos3t (or) cos3 3t= cos 9t 3cos 3t
4
1 3
L{cos3 3t} = L{cos 9t} + L{cos 3t}
4 4
∴ 1 s 3 s
= . 2 . 2
4 s 81 4 s 9
s 1 3 s s 2 63
= 2
4 s 81 s 2 9 s 2 9 s 2 81
1 2 s 2 2s 4
2
s s 4 s s2 4
1
= [1 cos 2t cos 4t cos 6t ]
4
1
L{cos t cos 2t cos 3t} L{1 cos 2t cos 4t cos 6t}
4
1
= [ L{1} L{cos 2t} L{cos 4t} L{cos 6t}]
4
= 1 1 2 s 2 s 2 s
4 s s 4 s 16 s 36
9. Find L(√ )
1
1
t L t 2
1 where n is not an integer
2
1
Sol: L
1
s2
1 1
2 2
3
3
n 1 n. n
2 2
s 2s
10. Find 𝑳 { + 𝜶 }, where a constant is
Sol: { 𝑖 +𝛼 }= { 𝑖 𝛼+ 𝑖 𝛼}
= 𝛼 { 𝑖 }+ 𝑖 𝛼 { }
= 𝛼 +
+ 𝑖 𝛼 +
constant
Proof: (i) By definition
L cf t e cf t dt c e st f t dt cL f t
st
0 0
(ii) By definition
L f t g t e st f t g t dt
0
e st
f t dt e st g t dt L f t L g t
0 0
Similarly the inverse transforms of the sum of two or more functions of ‘s’ is the sum of the
inverse transforms of the separate functions.
Thus, L1 f s g s L1 f s L1 g s f t g t
Corollary: L c1 f t c2 g t c1L f t c2 L g t , where c1, c2 are constants
a. e st f t dt b e st g t dt c e st h t dt
0 0 0
du
Put at u dt
a
when t then u and t = 0 then u = 0
s
su
du 1 a .u 1
L{ f (at )} e a
f (u ) . e f u du . f s
a a 0 a a
0
Solved Problems :
1. Find 𝑳{ }
Sol: {sinh } = =
−
1 1 3
2
1 s 9
2
3 s
3
2. Find 𝑳{ }
Sol: {cos } = =
+
f t dt
s a t
e
0
e ut f t dt where u s a
0
= f (u ) f s a
- at
Note: Using the above property, we have L{e f (t )} = f (s+ a)
Applications of this property, we obtain the following results
n! n!
1. L{e at t n } n 1
L(t n ) n 1
( s a) s
b b
2. L{e at sin bt} 2
L(sinbt) 2
( s a) b
2
s b 2
sa s
3. L{e at cos bt} 2
L(cosbt) 2
( s a) b
2
s b 2
b b
4. L{e at sinh bt} 2
L(sinhbt) 2
( s a) b
2
s b 2
sa s
5. L{e at cosh bt} 2
L(coshbt) 2
( s a) b
2
s b 2
Solved Problems :
1. Find the Laplace Transforms of t 3e 3t
−
+ +
{ cos }= =
+ + + +
3. Find L.T of
+
Sol: - L[ ]=L[ ]
= { [ ]+ [ ]}
= −
+ { [ ]} → −
−
= +
− − +
−
= −
+ − +
0
e
s (u a )
F u H u du e s (u a ) F u H u du
a 0
− + ∞ − +
= ∫− . +∫ .
− { }= −
=
Solved Problems
1.
cos t
Find the L.T. of g (t) when g t
3 if t
3
0 if t
3
∴ { }= { }=
+
=
( − 𝜋⁄ ) = ( − 𝜋⁄ ), 𝑖 > 𝜋⁄
={
, 𝑖 < 𝜋⁄
Now applying second shifting theorem, then we get
−𝜋
−𝜋
.
{ }= =
+ +
−
2. Find the L.T. of (ii) − − ) (ii) −
Sol: (i). Comparing the given function with f(t-a) u(t-a), we have a=2 and f(t)=t3
!
∴ { }= { }= = =
{ − }= { − − − − − −
(ii). − . − }= { − }
−
= then =
+
Multiplication by‘t’:
−
Theorem: If 𝑳{ }= 𝑳{ }=
d
ds
d
f s e st f t dt
ds 0
∞ −
= −∫ { } =− { }
−
Thus { }=
𝑛
∴ { }= − 𝑛
=
𝛼
{∫ ,𝛼 }=∫ 𝛼
,𝛼
{ }=− [ ]
+
− + − . −
= +
= +
2. Find t2sin at
Sol: Since { 𝑖 }=
+
{ . 𝑖 }= −
+
d 2as 2a 3s a
2 2
=
ds s 2 a 2 2 s 2
a
2 3
3. Find L.T of tet sin 3t
Sol: Since { 𝑖 }=
+
−
∴ { 𝑖 }= [ +
]= +
Now using the shifting property, we get
− + +
{ 𝑖 }= =
+ + + +
4. Find 𝑳{ }
Sol: Since { 𝑖 }=
+
L e 2t sin 3t
3 3
s 2 9 s 4s 13
2 2
− −
{ 𝑖 }= − [ ]= − [ ]
− + − +
3 2s 4 6 s 2
=
s 4s 13 s 4s 13
2 2 2 2
1 d 1 2 d 1
2
2 1 1
s ds s 1 ds 2 s 2
1 2 d 1
s s 1 ds s 2 2
2
1 2 2
s s 1 s 2 3
2
{ − }= − { − }
Sol:
− ! −
=− =
+ +
!
= +
7. Find 𝑳{ }
+ −
Sol. {cosh sin }= { . sin }
−
= [ { sin }+ { sin }]
1 a a
2 s a 2 a 2 s a 2 a 2
f t t 1 , t 1
2
− ∞ −
∞ − ∞
=∫ − =[ − −
] −∫ − −
2 st
e t 1 dt
s 1
0
2
e st
e
st
t 1 1 dt
s s s
1
2 e
st
2 1 2 st
0 e st dt 2 3 e 1
s s 1
s s 1 s
2
3
0 e s 3 e s
2
s s
9. Find the L.T of f (t) defined as f (t ) 3 , t >2
0, 0<t<2
∞ −
Sol: { }=∫
e st f t dt e st f t dt
2
0 2
− ∞ −
=∫ . +∫
3 st 3
e 2 s 0 e2s
0 e st 3dt
2 s
3
e 2 s
s
10. Find 𝑳{ + }
Sol: { + } = {cos cos − sin sin }
= cos . {cos } − sin {sin }
= cos . − sin .
+ +
−
{ . + }= [cos . − sin . ]
+ +
2 s 2 a 2 .0 a.2s
s a 2 .1 s.2s
= cos b. sin b
s a s a
2 2 2 2 2 2
1 s 2 a 2 2 cos b 2as sin b
s 2
a 2 2
= +
Division by‘t’:
∞
Theorem: If 𝑳{ }= 𝑳{ }=∫
Proof: We have f s e st f t dt
0
f t e st ds dt (
t is independent of‘s’)
0 s
e st
f t dt
0
t s
∞ −
=∫ { }
Solved Problems:
1. Find 𝑳 { }
Sol: Since { 𝑖 } = =
+
Division by t, we have
i a ∞ ∞
{ }=∫ =∫ +
1
a. Tan 1 s Tan 1 Tan 1 s
a as a
Tan1 s cot 1 s
2 a a
− 𝐚
3. Evaluate 𝑳 { }
1
1 s2 s2 2 s2 a2
l og 2 2
log 2 2
log
2 s a s a s2
− s2 1
Note: 𝑳 { } = log (Putting a=1 in the above problem)
s
− − −
4. Find 𝑳 { }
− − − ∞
Sol: { }=∫ −
+ +
s a
log s a log s b s log
s b s
a
1 s sa
l t log log
s
1 b sb
s
sb
log1 log( s a) log( s b) log
sa
1 cos t
5. Find L
t
2
1 cos t 1 1 cos t
Sol: L L . ..... 1
t t
2
t
1 cos t 1 s
s 2 ds log s log s 1
1 2
Now L
t s s 1 2 s
1 s2 1 s2 1 s2 1
log 2 log 2 log 2
2 s 1 s 2 s 1 2 s
1 cos t s 1 2
1
L
t
2
s 2
log 2 ds
s
1 s 2 1 s
2
2
log 2 .s 2 3 .sds
s s 1 s
2 s s
1 1 s2 1 ds
lt s.log 1 2 s log 2 2 2
s
2 s s s s 1
1 1 1 1 s 2 1
lt s 2 4 6 .... s log 2 2Tan 1s
2 s s 2s 3s s s
1 1 x 2 x3 x 4
0 s log 1 2 2 Tan1s
log 1 x x .....
2 s 2 2 3 4
1 1
cot 1 s s log 1 2
2 s
− − −
6.Find L.T of
− −
Sol: W.K.T L[ ]= ,L[ ]=
+ +
∞
L[ ]=∫ ̅
− − ∞
−
∴ L[ ]=∫ −
+ +
= [log + − log + ]∞
+ ∞
= log +
+
∞
= log
+
+
=log (1)-log ( )
+
+ +
=0- log ( ) = log ( )
+ +
e st f t s e st f t dt
0 0
st
lt e
= t f (t ) f (0) + . { }
∴ { }= − + { }
= −
The Laplace Transform of the second derivative f11(t) is similarly obtained.
∴ { }= . { }−
s. s f s f 0 f 1 0
s 2 f s sf 0 f 1 0
∴ { }= . { }−
= [ { }− − ]−
= { }− − −
Proceeding similarly, we have
{ }= { }− − − −
− ……
−
Note 1: { }= 𝑖 = = , = … =
Note 2: Now | | . .
− −
We have | |= | | .
− −
=M. → → ∞ if s>a
−
∴ →∞
= >
Solved Problems:
Using the theorem on transforms of derivatives, find the Laplace Transform of the
following functions.
(i). eat (ii). cosat (iii). t sin at
(i). Let f t eat Then f 1 t a.eat and f 0 1
{ }= . { }−
𝑖. . , { }= . { }−
𝑖. . , { }− . { }=−
𝑖. . , − { }=−
∴ { }=
−
(ii). = ℎ =− 𝑖 =−
∴ { }= { }− . −
Now f 0 cos0 1and f 1 0 a sin 0 0
ℎ {− cos }= {cos }− . −
⟹− {cos }− {cos }=−
⟹− + {cos } = − ⇒ {cos }=
+
Also f = and =
{ }= { }− −
𝑖. . , { cos − sin }= {tsin }− −
𝑖. . , {cos }− { sin }− {tsin }=
−
𝑖. . , − + {t sin }= ⇒ {t sin }=
+ +
Proof: Let g t f x dx
t
Then = [∫ ]= =
But { }= { }− = { }− [ 𝑖 = ]
∴ { }= { }
⟹ { }= { } ⇒ { }= { }
=∫
∴ {∫ }=
Solved Problems:
t
1. Find the L.T of sin atdt
0
Sol: L{sin }= =
+
{∫ }=
∴ {∫ sin }= +
sin t
t
2. Find the L.T of
0 t
dt
i
Sol: {sin } = →
= 𝑖
+
i ∞ ∞
∴ { }=∫ {sin } =∫ +
s 2 s
i − −
𝑖. . , { }= ( ⁄)
i − −
∴ {∫ }= ( ⁄)
−
3. Find L.T of ∫
− 𝑖
Sol: L[ ∫ ]
We know that
L {sint} = +
= ̅
𝑖 ∞ ∞
L{ }=∫ ̅ =∫ +
− ∞
=(
− − 𝜋 − −
= ∞− = − =
𝑖 −
∴ L{ }=
𝑖 −
Hence L {∫ }=
−
𝑖 −
∴ L[ ∫ ]= +
+
L f t
1
e st f t dt
a
1 e sa 0
Solved Problems:
L f t
1
e st f t dt
2b
Sol:
1 e 2bs 0
1 b st
e f t dt e st f t dt
2b
2 bs 0
1 e b
1 b st
e dt e st dt
2b
2 bs 0
1 e b
1 e st b e st 2b
1 e2bs s 0 s b
2 bs
e sb 1 e 2bs e sb
1
s 1 e
L f t
1
1 2e sb e 2bs
s 1 e
2 bs
2. Find the L.T of the function f t sin t if 0 t
2 2
0if t where f t has period
2
Sol: Since f (t) is a periodic function with period
L f t
1
e st f t dt
a
1 e sa 0
2
L f t
1
s 2
e st f t dt
1 e 0
1 st 2
0 e .0dt
st
2 s
e sin t dt
1 e
1 e st s sin t cos t
2 s
1 e
s2 2 0
eat
a sin bt b cos bt
b
a
eat sin bt
a 2 b2
1 e
1
2 s
1 s
s 2 2 e .
Laplace Transform of Some special functions:
1. The Unit step function or Heaviside’s Unit functions:
<
It is defined as − ={
>
Laplace Transform of unit step function:
−
To prove that { − }=
<
Proof: Unit step function is defined as − ={
>
∞ −
Then { − }=∫ −
e st u t a dt e st u t a dt
a
0 a
a
e st .0dt e st .1dt
0 a
e st 1 as e as
st
e dt s . e e s
a
s a
−
∴ { − }=
⁄ ∈
Proof: By the definition ∈ ={ ∈
>∈
∞ −
And Hence { ∈ }=∫ ∈
∈ − ∞ −
=∫ ∈ + ∫∈ ∈
∈ − ∞ −
=∫ + ∫∈ .
∈
− ∈
− ∈ − − ∈
= ∈ [ − ] = − ⁄∈ [ − ]=
∈
− − ∈
∴ { ∈ }=
∈
− − ∈
Now {𝛿 } =∈→ { ∈ } = ∈→
∈
∴ {𝛿 }= 𝑖 L-Hospital rule.
Properties of Dirac Delta Function:
∞
1. ∫ 𝛿 =
∞
2. ∫ 𝛿 = where G(t) is some continuous function.
∞
3. ∫ 𝛿 − = where G(t) is some continuous function.
4. G (t ) 1 t a G1 (a)
0
Solved Problems:
}= −
1. Prove that 𝑳{𝜹 −
Sol: By Translation theorem
}= − {𝛿 }
L{𝛿 −
− {𝛿 }= ]
= [sin
∞
2. Evaluate ∫ 𝜹( − 𝝅⁄ )
Sol: By using property (3) then we get
∞
∫ 𝛿 − =
Here = 𝜋⁄ , = cos
G a G 3 cos 2 3 12
∞
∴ ∫ cos 𝛿( − 𝜋⁄ ) = cos 𝜋⁄ = −𝜋⁄
3. Evaluate e 4t 1 t 2 dt
0
t a G t dt G (a)
1 1
G t e4t and a 2
G1 t 4.e4t
G1 a G1 2 4.e8
e 4t 1 t 2 dt G1 a 4.e 8
0
{ }= ⁄ −
1. { ⁄ }=
−
2.
{ }= { ⁄ − }=
−
−
3.
{ − }= { ⁄ + }=
+
4. ! −
{ }= +
𝑖 + 𝑖 { + }=
!
5. − ! −
{ − }= −
{ ⁄ }= , = , , …
− !
6. {sin }= −
+ { }= . sin
+
7. {cos }= −
{ } = cos
+ +
8. {sinh }= −
− { }= sinh
−
9. {cosh }= −
{ } = cosh
− −
{ }= −
10. sin { }= . sin
− + − +
11. − −
{ cos }= −
{ }= cos
− + − +
12. −
{ sinh }= { }= . sinh
− − − −
13. − −
{ cosh }= −
{ }= cosh
− − − −
14. − − −
{ sin }= { }= . sin
+ + + +
15. −
+ −
+ −
{ cos }= { }= cos
+ + + +
{ }= − −
16. − { − }= { }
{ − }= − − − −
17. + { + }= { }
Solved Problems :
s 2 3s 4
1. Find the Inverse Laplace Transform of
s3
− − + −
Sol: { }= { ⁄ − . ⁄ + ⁄ }
L1 1 s 3L 1 s L 4 s
1
2
1
3
t2
1 3t 4. 1 3t 2t 2
2!
+
2. Find the Inverse Laplace Transform of − +
− + − + − − +
Sol: { − +
}= { }= { }
− + − +
− − −
= { }+ . { }
− + − +
4
e 2t cos 3t e 2t sin 3t
3
2s 5
3. Find the Inverse Laplace Transform of
s2 4
− − −
Sol: { −
}= { − }
− −
− −
= { }− { }
− −
1
2.cosh 2t 5. sinh 2t
2
+
4. Find 𝑳− { }
+
− + + −
Sol: { }= { + }
+ +
− − −
= { }+ { }= +
+
−
5. Find 𝑳− { +
}
−
Sol: 𝑳− { }= −
{ }− −
{ }
+ + +
− −
= { }− { }
+( ⁄ ) +( ⁄ )
3 5 8 2 5
.cos t . sin t
4 2 4 5 2
3 5 4 5
cos t sin t
4 2 5 2
s
6. Find the Inverse Laplace Transform of
s a
2
− − + − − − −
Sol: { +
}= { }= { }
+
− −
= { − }
− − −
= [ { }− . { }]
eat 1 at
+
7. Find 𝑳− { − −
}
3s 7 A B
Sol: Let
s 2s 3 s 1 s 3
2
− + + = +
= , = ⇒ =
= − ,− = ⇒ =−
3s 7 1 4
s 2s 3 s 1 s 3
2
−
+ −
− − −
{ }= { + }=− { }+ { }
− − + − + −
e t 4.e3t
8. Find 𝑳− { (
}
+ + )
+
Sol: + +
= +
+ +
+ +
+ + + + + + + =
Equating Co-efficient of s3, A+C=0……..(1)
1 1
s 2 2
12 s 2 1
2 2
s 1 s 1 s
− − −
{ + +
}= [ { +
}− { +
}]
− −
= [sin − { }]
1
sin t te t
2
9. Find 𝑳− { +
}
= + + − +
s As B Cs D
Let 2 2
s 4a
4 4
s 2as 2a 2
s 2as 2a 2
+ − + + + + + =
1 1
Solving we get A 0, C 0, B ,D
4a 4a
−
− −
{ }= { }+ { }
+ + + − +
1 1 1 1 1
a.L 2
.. L1 2
( s a ) a 4a ( s a) a
2 2
4
1 1 at 1 1
. .e sin at . e at sin at
4a a 4a a
− + ( − )
10. Find i. 𝑳− { } ii. 𝑳− { }
Sol:
− − + − −
i. { }= { − + }= { − + }
− − −
= { }− { }+ { }
= − + = − +
!
− ( − ) − ( − ) − − +
ii. { }= { }= { }
− − − −
= { − + }+ { { }− { }+ { }}
3 t 2 4t 4 3 t4 1 4
1 4 1 2t 2
t 6t 2 6
2 2! 4! 2 6 4
11. Find 𝑳− [ −
]
Sol:
− − − −
[ ]= [ ]= [ ]= [ + ]
− − − + − +
1 at at
e e cosh at
2
4
12. Find L1
( s 1)(s 2)
4 1 1 1
Sol: L1 4 L1 4 L1 4[et e2t ]
( s 1)(s 2) ( s 1)(s 2) s 1 s 2
1
13. Find L1
( s 1) ( s 4)
2 2
1 A B Cs D
Sol: 2
( s 1) ( s 4) s 1 ( s 1)
2 2 2
s 4
2 1 2 3
A , B ,C ,D
25 5 25 25
1 2 1 1 1 1 1 2 1 s 3 1 1
L1 L L 2
L 2 L 2
( s 1) ( s 4) 25 s 1 5 ( s 1) 25 s 4 25 s 4
2 2
2 t 1 1 1 t 1 1 2 3 1
e L e L 2 cos 2t . sin 2t
25 s 5 s 25 25 2
2 t 1 t 2 3
e e .t cos 2t sin 2t
25 5 25 50
s2 s 2
14. Find L1
s( s 3)(s 2)
s2 s 2 A B C
Sol:
s s 3 s 2 s s 3 s 2
A 1 , B 4 ,C 2
3 15 5
s2 s 2
1 1 1 4 2
L L
s( s 3)(s 2) 3s 15( s 3) 5( s 2)
1 4 1 2
L1 L1 L 5( s 2)
3s 15( s 3)
1 4 2
e 3t e 2t
3 15 5
s 2 2s 4
15. Find L1 2
( s 9)(s 5)
s 2 2s 4 A Bs C
Sol: = 2
( s 9)(s 5) s 5 s 9
2
A 31 ,B 3 , C 83
34 34 34
s 2 2s 4 1 s 2 s 4
2
L1 2 L 2
( s 9)(s 5) ( s 9)(s 5)
31 1 3 1 83
L1 L L
34(s 5) 34( s 9) 34( s 9)
2 2
31 5t 1 83
e 3cos 3 t sin 3t
34 34 3
Proof: We have seen that L eat f (t ) f ( s a ) L1 f ( s a) eat f (t ) eat L1 f (s )
Solved Problems :
1
1. Find L1
1
L f ( s 2)
( s 2) 16
2
1 2 t 1 1
Sol: L1 e L 2
( s 2) 16 s 16
2
1 e2t sin 4t
e2t . sin 4t
4 4
3s 2
2. Find L1 2
s 4 s 20
3s 2 1 3s 2 1 3( s 2) 4
Sol: L1 2 L L 2
s 4s 20 ( s 2) 16 ( s 2) 4
2 2
s2 1
3L1 2
4 L1 2
( s 2) 4 ( s 2) 4
2 2
s 1
3e 2t L1 2 2
4e 2t L1 2 2
s 4 s 4
1
3e 2t cos 4t 4e 2t sin 4t
4
s3
3. Find L1 2
s 10s 29
s3 1 s3 1 s 5 8
Sol: L1 2 L 2
L 2
s 10s 29 ( s 5) 2 ( s 5) 2
2 2
s 8 5t 1
e5t L1 2 2
e cos 2t 8. sin 2t
s 2 2
then L G t e as . f (s)
L1 e as f ( s) G t
Solved Problems :
1 e s 1 e
3 s
1. Evaluate (i) L1 2 (ii) L 2
s 1 ( s 4)
1 e s 1 1 1 e
s
Sol: (i) L1 2 = L 2 + L 2
s 1 s 1 s 1
1
Since L1 sin t f (t ) , say
s 1
2
e s sin(t ) , if t
1
∴ By second Shifting theorem, we have L 2
s 1 0 , if t
e s
or L1 2 =sin(t-π)H(t-π)= -sint. H(t-π)
s 1
1 e s
Hence L1 2 =sint-sint. H (t-π) =sint [1- H (t-π)]
s 1
Where H (t-π) is the Heaviside unit step function
1 4t 1 1
(ii) Since L1 2
e L 2
( s 4) s
= e4t .t f (t ) , say
If L f t f s , Then L1 f as 1 t
f ,a 0
a a
1 t
Then f as L f , a 0
a a
L1 f as 1 t
f ,a 0
a a
Solved Problems :
s 1 8s
1. If L1 2 2
= t sin t , find L1 2
( s 1) 2 (4s 1)
2
s 1
Sol: We have L1 2
= t sin t ,
( s 1) 2
2
Writing as for s,
as 1 1 t t t t
L1 2 2 2
= . . sin = 2 .sin , by change of scale property.
(a s 1) 2 a a a 2a a
Putting a=2, we get
2s t t 8s 1 t
L1 2 2
= sin or L1 2
= sin
(4s 1) 8 (4s 1) 2
2
2 2
n
Theorem: L1 f ( s) f (t ) , then L1 f ( s ) (1) n t n f (t ) where f ( s)
n dn
f (s)
ds n
n
Proof: We have seen that L t n f (t ) (1) n
d
f (s)
ds n
n
L1 f ( s ) (1) n t n f (t )
Solved Problems :
s 1
1. Find L1 log
s 1
s 1
Sol: Let L1 log f (t )
s 1
s 1
L f (t ) log
s 1
d s 1
L tf (t ) log
ds s 1
1 1
L tf (t )
s 1 s 1
1 1
tf (t ) L1
s 1 s 1
1 1 1
tf (t ) 1. L1 L
s 1 s 1
−
= +
2sinh t
t f t 2sinh t f t
t
s 1 2sinh t
L1 log
s 1 t
t
1 s 1 e
Note: L1 log =
s t
L f (t ) cot 1 (s)
d 1 1
L tf (t ) [cot 1 ( s ) 2
1 s 1 s
2
ds
1
tf (t ) L1 2 sin t
s 1
sin t
f t
t
s t
f (t )
Proof: we have seen that L f ( s)ds
t s
f (t )
L f ( s )ds
1
s t
Solved Problems :
s 1
1. Find L1 2
( s 2s 2)
2
s 1
Sol: Let f ( s )
( s 2s 2) 2
2
s 1
Then L 1
f ( s) = L 2
1
ds
s ( s 2 s 2)
2
s 1
= L1 2
[( s 1) 1]
2
s
= et L1 2 2
, by First Shifting Theorem
( s 1)
t t s t
et sin t e t sin t L1 2 2 2
sin at
2 2 ( s a ) 2a
Theorem: L1 f ( s) f (t ) , and f (0),then L1 s f ( s) f 1 (t )
Proof: we have seen that L f 1 (t ) s f ( s ) f (0)
L f 1 (t ) s f ( s ) [ f (0) 0] or
L1 s f ( s ) f 1 (t )
Solved Problems :
s s
1. Find (i) L1 2
(ii) L1 2
( s 2) ( s 3)
1
Sol: Let f ( s) Then
( s 2) 2
1 2t 1 1 2t
L1 f ( s ) L1 2
( s 2)
= e L 2 = e .t f (t ) ,
s
Clearly f (0) =0
s 1 s 1
Thus L1
( s 2) 2
= L s.
( s 2) 2
= L s. f ( s) = f 1 (t )
d
= (te 2t ) = t(2e2t ) e2t .1= e2t (1 2 t)
dt
s 3t
Note: in the above problem put 2=3, then L1 2
= e (1 3 t)
( s 3)
Division by S:
f s
t
Theorem: If L1 f ( s) f (t ) , Then L1 f u du
s 0
t f s
Proof: We have seen that L f u du
0 s
f s t
L1 f u du
s 0
f s
t t
Note: If L 1
f ( s) f (t ) , then L 2 f u du.du
1
s
00
Solved Problems :
1
1. Find the inverse Laplace Transform of
s (s a 2 )
2 2
1 1
Sol: Since L1 2 2
= sinat , we have
(s a ) a
1 t1
L1 2 2
= sin atdt
s( s a ) 0 a
t
1 cos at 1 1
= = a 2 (cosat 1) = a 2 (1 cosat)
a a 0
1 t 1
1
Then L 2 2 2
= 2 (1 cos at )dtdt
s (s a ) 0 a
t
1 sin at 1 sin at
= t = t a
a2 a 0 a 2
1 1 sin at
L1 2 2 2
= 2 t
s (s a ) a a
Convolution Definition:
If f (t) and g (t) are two functions defined for t 0 then the convolution of f (t) and g (t) is
defined as f t * g t f u g t u du
t
Properties:
The convolution operation * has the following properties
3. Distributive f * g h t f * g t f * h t for t 0
L f t * g t L f t .L g t f s .g s
i.e., The L.T of convolution of f(t) and g(t) is equal to the product of the L.T of f(t) and g(t)
Proof: WKT L t e st
0 f u g t u dudt
0
t
f u g t u du dt
t
e
st
0 0
The double integral is considered within the region enclosed by the line
u=0 and u=t
On changing the order of integration, we get
L t e st f u g t u dt du
0 u
e su f u
0 e
u
s t u
g t u dt du
f u e
e su sv
g v dv du put t u v
0 0
e su f u g s du g s e su f u du g s . f s
0 0
L f t * g t L f t .L g t f s .g s
Solved Problems :
s
1. Using the convolution theorem find L1 2 2 2
(s a )
s 1 s 1
Sol: L1 2 2 2 L 2
= . 2 2
(s a ) s a s a
2
s 1
Let f s and g s 2
s a 2
2
s a2
s
So that L1 f ( s ) L1 2
2
s a
cos at f (t ) say
s 1
L1 g ( s ) L1 2 2
s a a
sin at g (t ) say
s t 1
L1 2 2 2
cos au. .sin a(t u)du
(s a ) 0 a
t
1
sin(au at au) sin(au at au) du
2a 0
t
1
sin at sin(2au at ) du
2a 0
t
1 1
sin at.u .cos(2au at )
2a 2a 0
1 1 1
t sin at cos 2 at at cos at
2a 2a 2a
1 1 1
t sin at cos at cos at
2a 2a 2a
t
sin at
2a
s2 1
2. Use convolution theorem to evaluate L 2 2
( s a )(s b )
2 2
s2 s s
Sol: L1 2 2 2 2 L1 2 2 . 2 2
( s a )( s b ) s a s b
s s
Let f s and g s 2
s a22
s b2
s
So that L1 f ( s ) L1 2 2
s a
cos at f (t ) say
s
L1 g ( s) L1 2 2 cos bt g (t ) say
(s b )
By convolution theorem, we have
t
s s
1
L 2 . 2 2
s a s b 0
2
cos au.cosb(t u )du
t
1
cos(au bu bt ) cos(au bu bt ) du
2 0
t
1 sin(au bu bt ) sin(au bu bt )
2 a b ab 0
1
3. Use convolution theorem to evaluate L1 2 2
s( s 4)
1 1 1 s
Sol: L1 2 2
L 2. 2 2
s( s 4) s ( s 4)
1 s
Let f s and g s
s 4
2 2
s 2
1
So that L1 g ( s ) L1 2 t g (t ) say
s
t.sin 2t ts in 2t
L1 f ( s) L1 2
s
2
( s 4)
f (t ) say L1 2
s
2 2
4 ( s a ) 2a
1 s tu
1
L 2 . 2 2
sin 2u(t u )du
s ( s 4) 0 4
t t
1 2
t
4 u sin 2udu
0
4 0
u sin 2udu
t
u 1
t cos 2u sin 2u
4 2 4 0
t
1 u 2 u 1
cos 2u sin 2u cos 2u
4 2 2 4 0
1
1 t sin 2t cos 2t
16
1
4. Find L1
( s 2)( s 1)
2
1 1 1
Sol: L1 L1 . 2
( s 2)( s 1) s 2 s 1
2
1 1
Let f s and g s 2
s2 s 1
1
So that L1 f ( s ) L1 e f (t ) say
s 2
2t
1
L1 g ( s ) L1 2 sin t g (t ) say
s 1
t
1 1
L1 . 2
s 2 s 1
f (u).g(t u) du
0
(By Convolution theorem)
t t
e sin(t u) du (or) sinu.e2(t u ) du
2u
0 0
sin ue
2u
e 2t
du
0
t
e 2u
e 2 2sin u cos u
2t
2 1 0
1
e 2t e 2t 2sin t cos t 1 1
5 5
5
e 2sin t cos t
1 2t
1
5. Find L1
( s 1)( s 2)
1 1 1 1
Sol: L1 L .
( s 1)( s 2) s 1 s 2
1 1
Let f s and g s
s 1 s2
1 t
So that L1 f ( s ) L1 e f (t ) say
s 1
1
L1 g ( s ) L1 e g (t ) say
s 2
2t
1 t u 2(t u )
e e
1
L du
( s 1)(s 2) 0
t
t t
e 3u 1 2t
e e 2 t 3u
du e 2t
e
3u
du e 2t
e e
t
0 0 3 0 3
1
6. Find L1 2 2 2
s (s a )
1 1 1
Sol: L1 2 2 2
L1 2 . 2 2
s (s a ) s s a
1 1
Let f s and g s 2
s 2
s a2
1
So that L1 f ( s ) L1 2 t f (t ) say
s
1 1
L1 g ( s ) L1 2 2
s a a
sinh at g (t ) say
1 t 1
L1 2 2 2 u. sinh a(t u )du
s (s a ) 0 a
t
1
u sinh(at au)du
a0
sin at au
t
1 u
cosh at au
a a a2 0
1 t 1
cosh(at at ) 0 2 [0 sinh at ]
=a a a
1 t 1
sinh at
a a a 2
1
at sinh at
a3
3. Using Convolution theorem, evaluate 𝑳− { }
+ +
Sol: −
{ +
. +
}= −
{ +
. +
}= −
{ ̅ . ̅ }
̅ = = { }⇒ = −
{ }= −
----------------- (1)
+ +
−
̅ = = { }⇒ = { }= ----------------- (2)
+ +
Where ∗ =∫ −
− − −
∴ { +
. +
}=∫
−
= ∫
−
= . [ − 𝑖 ]
+
−
= [ − − 𝑖 ]
−
[ − ]−
= − 𝑖
Where f 1 0 s f 0 – f 0
f 2 (0) s 2 f s – s f 0 f 1 0 , and so on
4. Step4: take the inverse Laplace transform on both sides this gives f as a function of t
which gives the solution of the given DE
Solved Problems :
1. Solve y111 2 y11 y1 2 y 0 using Laplace Transformation given that
L y111 (t ) 2L y11 (t ) L y1 2L y 0
s 3L y (t ) s 2 y (0) sy1 (0) y11 (0) 2 s 2 L y (t ) sy (0) y1 (0)
sL y(t ) y(0) 2L y(t ) 0
s 3 2s 2 s 2 L y (t ) s 2 y (0) sy1 (0) y11 (0) 2sy (0) 2 y1 (0) y (0)
0 0 6 2.0 2.0 0
s 3 2s 2 s 2 L y (t ) 6
6 6
L y (t )
s 2s s 2 ( s 1)( s 1)( s 2)
3 2
A B C
s 1 s 1 s 2
A( s 1)( s 2) B( s 1)( s 2) C ( s 1)( s 1) 6
A(s 2 3s 2) B(s 2 s 2) C (s 2 1) 6
A BC 0
2A 2B C 6
______________
3A B 6
3A B 0
_____________
6A 6 A 1
3 A B 0 B 3 A B 3
A B C 0 C A B 1 3 2
1 3 2
L y (t )
s 1 s 1 s 2
1 1 1 1 1
y (t ) L1 3.L 2.L
t
= e 3e 2.e
t 2 t
s 1 s 1 s 2
Which is the required solution
y 0 1and y1 0 1
L{ y11 (t ) 3L y1 (t ) 2L y (t ) 4L t L e3t
4 1
s 2 L{ y (t ) sy (0) y1 (0) 3 sL{ y (t ) y (0) 2L{ y (t )
s 2
s 3
4 1
( s 2 3s 2)L{ y (t ) s4
s 2
s 3
4s 12 s 4 s 2 3s3 4s3 12s 2
( s 2 3s 2)L{ y(t )
s 2 (s 3)
s 4 7s3 13s 2 4s 12
L{ y(t )
s 2 (s 3)(s 2 3s 2)
s 4 7s3 13s 2 4s 12
L{ y(t )
s 2 (s 3)(s 1)( s 2)
s 4 7s3 13s 2 4s 12 As B C D E
s ( s 3)(s 1)(s 2)
2
s 2
s 3 s 1 s 2
( As B)(s 1)(s 2)(s 3) C (s 2 )(s 1)(s 2) D(s 2 )(s 2)(s 3) E (s 2 )(s 1)(s 3)
s 2 (s 3)(s 1)(s 2)
s 4 7 s3 13s 2 4s 12 ( As B)( s3 6s 2 11s 6)
C ( s 2 )( s 2 3s 2) D( s 2 )( s 2 5s 6) E.s 2 ( s 2 4s 3)
6 A B 3C 5D 4 E 7......................(2)
1
put s 1, 2 D 1 D
2
put s 2, 4 E 8 E 2
1
put s 3,18C 9 C
2
1 1
from eq.(1) A 1 2 A3
2 2
3 5
from eq.(2) B= -7+18+ 8 3 1 2
2 2
3 2 1 1 2
y (t ) L1 2
s s 2( s 3) 2( s 1) s 2
1 1
y t 3 2t e3t et 2.e 2t
2 2
d2y dy dy
3. Using Laplace Transform Solve 2
2 3 y sin t , given that y 0 when t=0
dt dt dt
d2y dy
Sol: Given equation is 2
2 3 y sin t.
dt dt
L y11 t 2L y1 t 3L y t L sin t
s 2 L y t sy 0 y1 0 2 sL y t y 0 3.L y t
1
s 1
2
s 2 2s 3 L y t
1
s 1
2
L y t 2
1
s 1 s 2 2s 3
1
y t L1
s 1 s 3 s 2 1
Now consider
1 A B Cs D
2
s 1 s 3 s 1 s 1 s 3 s 1
2
A s 3 s 2 1 B s 1 s 2 1 Cs D s 1 s 3 1
1 1 1 1
s
y t L
1 8 40 2 5
10
s 1 s 3 s 1
1 1 1 1 1 1 1 s 1 1 1
L1 L L 2 L 2
8 s 1 40 s 3 10 s 1 5 s 1
1 1 1 1
y t et e 3t cos t sin t
8 40 10 5
dx
4. Solve x sin t, x 0 2
dt
dx
Sol: Given equation is x sin t
dt
L x1 t L x t L sin t
s.L x t x 0 L x t
s 2
2
s.L x t 2 L x t
s 2
2
s 1 L x t 2
s 2
2
2
x t L1
s 1 s s 1
2 2
1 1
2L1 L
s 1 s 1 s
2 2
(By using partial fractions)
s
2
2et L1 1 12 2 12 2
2 2
s 1 s s
1
2e t et cos t . sin t
1
2
1 2
1 2
5. Solve D 2 n 2 x a sin nt given that x=Dx=0, when t=0
x11 t n2 x t a sin nt
s 2 n 2 L x t a cos
n s
a sin . 2
s n
2 2
s n2
L x t a cos
n s
a sin
s
2 2
s n2
2
2
n 2
1 1 a sin 1
d 1
.sin nx. sin n t x dx
t
na cos L 2
s n
0 n 2
n 2 ds
a cos a sin 1
cos nt 2nx cos nt dx
t
t sin nt
2n 0 2 n
a cos t
2n 0 cos n t 2 x cos nt dx
a
sin t sin nt
2n
a cos 1 at sin
t
2n 2n .sin n t 2 x x cos nt 2n sin nt
0
a cos sin nt at
cos cos nt sin sin nt
2n 2 2n
a cos sin nt at
cos nt
2n 2 2n
−
6. Solve − + = using L.T given that y (0) = y1 (0) = 1.
−
Sol: Given equation is − + =
−
Applying L.T on both sides we get − + =
− − +
y= [ + + +
]+ [ + +
]
Let us consider
− −
[ + + +
]= [ + +
]
=
+ + + + +
= +
+ +
+ +
−
= +
+ +
+ +
−
−
= [ +
+ +
+ +
]
−
−
= [ +
+
+
+
+
]
− − −
=− [ +
]+ [ +
]+ [ +
]
− − − −
[ ]=− + + − −−→
+ + +
− + − + −
[ + +
]= [ + −
]+ [ + −
]
− − − − −
= [ ]+ + [ ]
− −
− + −
[ ] = + − − −→
+ +
𝝅
7. Solve + = using L.T. given x (0) =1, x ( ) = -1.
Sol: Given + =
L[ ]+ [ ]= [ ]
⇒ [ ]− − + [ }=
+
⇒( + [ ]− − = +
⇒( + [ ]= + +
+
[ ]= + +
+ + + +
− − −
X= [ ]+ [ ]+ [ ]
+ + + +
−
= [ − ]+ + 𝑖
+ +
− −
= [ ]− [ ]+ + 𝑖
+ +
= − + + 𝑖 ------------------→
𝜋
Given x ( ) = -1.
𝜋 𝜋 𝜋 𝜋 𝜋
∴ − = − + + + 𝑖
⇒ -1= - − + −
=− +
∴ x= + + 𝑖 From (1)
⇒ − + − [ ]− − + + + − = − [ ]
⇒ − L[y]- + =
−
=
−
⇒ − L[y] = + −
−
[ ]= + −
− − −
− − −
= [ ]+ [ ]− [ ]
− − −
− − −
= [ ]+ [ ]− [ ]
− − −
− − −
= [ ]+ −
− [ ]
−
= !
− !
+ [ −
]
−
Consider [ ]
−
− −
W.K.T [ ]= [ ]= =
− !
−
[ ]= = + = + + +
−
= + +
∴ = !
− !
− + +