SEMESTER II MATHEMATICS GENERAL LECTURE NOTE INTEGRAL CALCULUS AND DIFFERENTIAL EQUATIONS TARUN KUMAR BANDYOPADHYAY, DEPARTMENT OF MATHEMATICS In case you face difficulty in understanding the following material, you may e-mail to me at [email protected]stating your Name and Roll No. 1
106
Embed
sem ii - general - integral calculus and differential equations
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
SEMESTER II MATHEMATICS GENERAL
LECTURE NOTE
INTEGRAL CALCULUS AND DIFFERENTIAL EQUATIONS
TARUN KUMAR BANDYOPADHYAY, DEPARTMENT OF MATHEMATICS
In case you face difficulty in understanding the following material,
you may e-mail to me at [email protected] stating your Name and Roll No.
Text: An Introduction to Analysis: Integral Calculs— K. C. Maity,
R. K. Ghosh
CHAPTER 1: DEFINITE INTEGRALS—A REVISION
Let f be a real valued continuous function defined
on a closed and bounded interval [a,b]. Let us
choose a partition (collection of finite number of
points of [a,b] including a and b) P = {a = x0,x1,x2,
…,xn=b} of [a,b] (for example: {0,1/2,1} is a partition
for [0,1].How many partitions are there for a given interval—finite or infinite?).
2
Fig 1: Area under a curve is approximated by sum of areas of rectangles
Let δr = xr-xr-1, r = 1,…,n and δ=max {δr| r = 1,2,…,n}. Choose
an arbitrary point cr∈(xr-1,xr) for all r and consider sum
of areas of rectangles ∑1
n
f (cr)δ r. It can be seen that this
sum approaches more closely the actual area under the
curve if we make width of the rectangles smaller , that
is, if we increase number n of points of subdivision
(sum of areas of two rectangles on δ rgives a better approximation to the area under the curve than area of a single rectangle).
Definition 1.1 ∫a
b
f ( x )dx= limn→∞
∑1
n
f (cr)δ r, provided the limit
exists independent of choice of points of subdivision xi and that of ci , for all i. It can
be proved that for a continuous function f defined over a closed bounded interval
[a,b], ∫a
b
f ( x )dx exists in above sense.
3
Simpler equivalent expression for
calculating ∫a
b
f ( x )dx:
We can make choices of xi and ci suitably so as to
obtain equivalent simpler expression of ∫a
b
f ( x )dx.
Let us choose xi’s equi-spaced, that is , δ1 = δ2 = … = δn =
(b-a)/n. Then∫a
b
f ( x )dx= limn→∞
b−an ∑
1
n
f (cr).
Let us choose cr = a+rh, r = 1,…,n, where h = (b-a)/n.
Then∫a
b
f ( x )dx= limn→∞
b−an ∑
1
n
f (a+r b−an
) = limh→ 0
h∑1
n
f (a+rh).
As a special case,∫0
1
f ( x )dx= limn→∞
1n∑1
n
f ( rn )=limh→0h∑
1
n
f (rh)
Example 1.1 From definition, calculate ∫0
1
x2dx.
» ∫0
1
x2dx = limh→0
h∑1
n
(rh)2= limh→ 0
(nh ) (nh+h )(2nh+h)6 = 1/3, since nh = 1
holds, for every positive integer n and the
corresponding h.
Fundamental Theorem of Integral Calculus
4
Theorem 1.1 If ∫a
b
f ( x )dxexists and if there exists a
function g:[a, b]→R such that g1(x) = f(x)(suffix denotes
order of differentiation) on [a, b], then ∫a
b
f ( x )dx = g(b) –
g(a).
NOTE: g is called a primitive of f. A function f may not
possess a primitive on [a,b] but ∫a
b
f ( x )dx may exist ; in that
case,∫a
b
f ( x )dxcan not be calculated using fundamental
theorem. Primitives of f on [a, b] are given by the
indefinite integral ∫ f ( x )dx : that is the reason why we consider indefinite integrals.
Example ∫0
1
x2dx exists, since x2 is continuous on [0,1].
Also g(x)=∫ x2dx= x3
3 +c is a primitive of x2 on [0,1]. Hence
Fundamental Theorem gives ∫0
1
x2dx=¿+c)-c=13.
Note ∫0
1
x2dxis independent of c though ∫ x2dx involves c.
PROPERTIES OF DEFINITE INTEGRALS
5
We assume below that the definite integrals exist and
whenever we consider ∫a
b
f ( x )dx, a primitive g to f over
[a,b] exists, so that we can apply Fundamental
Theorem. For a<b, we define ∫b
a
f ( x )dx = - ∫a
b
f ( x )dx.
1. ∫a
b
f ( x )dx=∫a
c
f ( x )dx+∫c
b
f ( x )dx(irrespective of relative algebraic
magnitude of a,b,c)
Example ∫3
1
x2dx+∫1
4
x2dx=
−∫1
3
x2dx+∫1
4
x2dx ¿−(9−13 )+(643
−13)=377 =∫
3
4
x2dx.
2. ∫0
a
f ( x )dx=∫0
a
f (a−x ) dx
Example ∫0
π2
sin x dx=∫0
π2
cos xdx
6
3. ∫0
a
f ( x )dx=∫0
a/2
f ( x )dx+∫0
a/2
f (a−x ) dx. In particular, if f(a-x) = f(x)for all
x in [0,a], then ∫0
a
f ( x )dx=2∫0
a /2
f ( x )dx and if f(a-x) = - f(x) for all x
in [0,a], then ∫0
a
f ( x )d x= 0.
4. ∫0
na
f ( x )dx=n∫0
a
f ( x )dx, if f(a+x) = f(x), n natural.
5. ∫−a
a
f (x ) dx=∫0
a
{f ( x )+ f (−x ) }dx. If f is even,∫−a
a
f (x )dx=2∫0
a
f ( x ) dx. If f is odd,
∫−a
a
f (x)dx= 0.
Example ∫−1
1
x3cos 2x dx=0, ∫−1
1
x4 cos2 xdx= 2∫0
1
x4 cos2 x dx
PRACTICE SUMS
1. Evaluate:∫0
π2 √sin x
√sin x+√cos xdx.
2. Evaluate:∫−2
2 x2 sin xx6+12
dx
3. If f(x) = f(x+ kp) for all integer values of k, show that
∫0
np
f ( x )dx=n∫0
p
f ( x )dx.
CHAPTER 2: REDUCTION FORMULA
7
In this chapter, we study how to decrease complexity
of some integrals in a stepwise manner by the use of
recurrence relation that we derive generally using
integration by parts formula.
1.Let In = ∫sinnx dx , n natural.
In = ∫sin n−1 x sin x dx = sinn-1x(-cos x)- (n-1)∫sin n−2 x cos x¿¿¿= sinn-
1x(-cos x)+(n-1)∫sin n−2x (1−sin2 x )dx= -sinn-1x cos x+(n-1)In-2-
(n-1)In. hence In=−sinn−1 xcos xn
+ n−1n In-2.
If we denote Jn = ∫0
π2
sinnx dx , then Jn = −sinn−1 xcos xn
∨¿0π2 ¿+n−1
n Jn-2 =
n−1n Jn-2. By repeated application of the reduction
formula, it can be proved that Jn= n−1n Jn-2=…=n−1n
n−3n−2
… 12. J 0=
n−1n
n−3n−2
… 12. π2 , if n is even natural and Jn =n−1
nn−3n−2
… 23J1=
n−1n
n−3n−2
… 23 , if n is odd natural.
2.Let In = ∫ tann xdx , n natural.
Then In = ∫ tann−2 x . tan2 xdx =∫ tann−2 x .(sec2 x−1)dx = tann−1 x
n−1−I n−2.
Also, Jn = ∫0
π4
tann xdx= tann−1 x
n−1∨¿0
π4 ¿-Jn-2= 1
n−1−¿ Jn-2.
8
3.Let In = ∫ secn x dx , n natural. Then In = ∫ secn−2 x . sec2x dx = secn-2x
tanx – (n-2)∫ secn−2 x .(sec2 x−1)dx = secn-2x tan x – (n-2)(In-In-2).
Hence In = secn−2 x tan xn−1 +n−2
n−1I n−2.
4.Let Im,n = ∫sinm x cosn x dx = ∫(sinm xcos x )cosn−1x dx = sinm+1 x cosn−1 x
m+1+ n−1m+1∫cos
n−2x sinm xsin 2 x dx=
sinm+1 x cosn−1 xm+1
+ n−1m+1∫cos
n−2x sinm x (1−cos2 x )dx = s∈¿m+1 xcosn−1 xm+1
+ n−1m+1
¿Im,n-2- n−1m+1
Im, n. Transposing and simplifying, we get a reduction
formula for Im,n.
5.Let Im,n = ∫cosm xsin nx dx= -cosm x cosnxn
−mn∫cos
m−1x sin xcos nxdx = -cosm x cosnx
n−mn∫cos
m−1 x sin nxcos xdx+mn∫ cos
m−1 x sin (n−1 ) x dx [ since cos nx
sin x = sin nx cos x – sin(n-1)x] = - cosm x cosnxn
−mn Im,n+m
nIm−1 , n−1.
Transposing and simplifying, we get a reduction
formula for Im,n.
Illustrative examples
1. I n = ∫0
π2
xn sin xdx and n>1, show that In+ n(n-1)In-2 = n¿
9
»I n =(−xn cos x )∨¿0π2+n∫
0
π2
xn−1 cos xdx ¿ = n¿ = n( π2)n−1
-n(n-1)I n−2. Hence.
2. Im,n = ∫ xm(1−x)ndx = xm+1
m+1(1−x )n+ n
m+1∫ xm[1−(1−x )] (1−x )n−1dx = xm+1
m+1(1−x )n+ n
m+1(Im,n−1−Im, n¿. Hence Im,n can be obtained.
3. Im = ∫0
π2
cosm x sinmxdx. From 5 above, Im = 12m +12Im−1 =
12m
+12 [ 12(m−1)
+ 12Im−2] = 12m + 1
22(m−1)+ 122
Im−2. Repeating the use of
the reduction formula, it can be proved that Im = 12m+1 [2+ 222 +
23
3+…+
2m
m ].4. Get a reduction formula for each of the following:∫ xm¿¿¿,
∫ sin nxsin xdx.
CHAPTER 3: IMPROPER INTEGRAL
When we consider the definite integral ∫a
b
f ( x )dx in earlier
standards, we implicitly assume two conditions to
hold: (a) f is continuous on [a, b] or , to that
10
matter, at least the limit ∫a
b
f ( x )dx= limn→∞
∑1
n
f (cr)δ r exists
independent of choice of points of subdivision xi and
that of ci, for all i and (b) the interval [a,b] is bounded. We want to extend the definition of ∫a
b
f ( x )dx when either (a) or (b) or both are not
met. This extended definition of definite integral is
referred to as Improper Integrals. Improper
integrals can be of two types: (a) Type 1: interval of integration is unbounded, (b) Type 2: integrand has a finite number of infinite discontinuities in the interval of integration.
Definition of TYPE I improper integral∫a
∞
f ( x )dx,∫−∞
a
f ( x )dx
and ∫−∞
∞
f ( x )dx
Let the function f be integrable in [a ,B], for every B>a.
If limB→∞
∫a
B
f ( x )dx exist finitely, we define ∫a
∞
f ( x )dx = limB→∞
∫a
B
f ( x )dx and
we say ∫a
∞
f ( x )d x exists or converges; otherwise ∫a
∞
f ( x )dx
11
diverges. Similarly, ∫−∞
a
f ( x )dx= limB→−∞
∫B
a
f ( x )dx (provided the
limit exists) and ∫−∞
∞
f ( x )dx=∫−∞
a
f ( x )dx+∫a
∞
f ( x )dx, a is any real,
provided ∫a
∞
f ( x )dx and ∫−∞
a
f ( x )dx exist separately.
Example 3.1 ∫1
∞ dxx2 ,∫1
∞ dx√ x . The range of integration of the
integrals are unbounded. For a>1, ∫1
a dxx2 = 1-1a and ∫
1
a dx√ x =
2(√a-1). Since lima→∞ (1−1a )= 1 exists but lima→∞
(2 (√a−1 ) ) does not
exist, hence the improper integral ∫1
∞ dxx2 converges and
∫1
∞ dx√ x
diverges.(compare areas below the curves y =
1/x2 and y = 1/√ x in diagram below)
Fig.2:Comparison of areas under y=1/x2 and y=1/√ x in (0,1] and [1,∞)
12
Definition of TYPE II improper integral
Let f have an infinite discontinuity only at the point a
(that is, limx→a+¿ f ( x )=±∞¿
¿or limx→a−¿ f ( x )=±∞¿
¿ and is continuous in (a,
b].Then we define ∫ab
f ( x )dx= limc→0+¿∫
a+c
b
f ( x ) dx
¿¿,0<c<b-a, provided
the limit exists. Similarly, if f has an infinite
discontinuity only at the point b and is continuous in
[a, b), then we define ∫ab
f ( x )dx= limc→0+¿ ∫
a
b−c
f ( x )d x
¿¿, 0<c<b-a,
provided the limit exists. If f has an infinite
discontinuity at d, a<d<b, and is otherwise continuous
in [a,b], we define ∫a
b
f ( x )dx=∫a
d
f ( x )dx+∫d
b
f ( x )dx, provided both of
∫a
d
f ( x )dx and ∫d
b
f ( x )dx exist separately.
Example 3.2 ∫0
1 dx√ x ,∫0
1 dxx2 . The integrands have an infinite
discontinuity at x=0. For 0<a<1, ∫a
1 1√ x
dx=2(1−√a) and
∫a
1 dxx2
= 1a−1. Since lim
a→0+¿2 (1−√a )¿¿ = 2 exists but lim
a→ 0+¿( 1a−1)¿¿ does not
exist, so ∫0
1 dx√ x converges whereas ∫
0
1 dxx2 diverges.
13
(Compare areas between x = a,0<a<1, and x = 1, below
the curves y = 1/x2 and y = 1/√ x in diagram above)
Example 3.3 ∫0
∞ dx4+9 x2
= π12
»The integrand 1
4+9 x2 is continuous everywhere but the
interval of integration is unbounded. Let a>o be fixed.
∫0
a dx4+9 x2 =
19∫0
a dx
x2+( 23 )2 = 16 tan
−1 3a2 . Thus lim
a→∞∫0
a dx4+9 x2 = 16 .
π2 = π12.
Example 3.4 ∫0
3 dx√9−x2
=π2 .
» The integrand has an infinite discontinuity at x = 3
and is continuous on [0, 3). Let 0<a<3. Then ∫0
a dx√9−x2 =
sin-1(a/3) . So lim
a→ 3−¿∫0
a dx√9−x2
¿
¿ =
lima→ 3−¿sin−1( a3 )=π
2 ¿
¿. Hence ∫
0
3 dx√9−x2=
π2 .
Note : we can apply standard methods of integration, in particular method of substitution, only to a proper integral and not directly to an improper integral. Thus if we
14
substitute z=1/x directly in the improper integral∫−1
1 dxx2 ,
we get a value -2 of the integral whereas it can be
checked from definition that the improper integral
diverges.
PRACTICE SUM
Verify that (1 )∫−∞
∞
x e− x2dx=0 ,(2)∫0
∞ dx(x2+a2) (x2+b2 ) = π
2ab(a+b), (3)∫0
∞ xdx(x2+a2 ) (x2+b2 )=
1a2−b2
ln( ab )(a,b>0),(4)∫0
∞
e−xsin x dx= 12, (5)∫
0
π dx4+5 sinx
=23 ln 2.
TESTS FOR CONVERGENCE OF IMPROPER INTEGRALS
TYPE I INTEGRAL
Theorem 3.1 (Comparison test) Let f and g be
integrable in [a, B], for every B>a. Let g(x)>0, for all x
≥a. If limx→∞
f (x )g (x) = c≠0, then the integrals ∫
a
∞
f ( x )dx and ∫a
∞
g ( x ) dx
either both converge or both diverge. If c = 0 and ∫a
∞
g ( x )dx
converges , then ∫a
∞
f ( x )dx converges.
15
Theorem 3.2 (μTest ¿ Let f be integrable in [a, B], for
every B>a. Then ∫a
∞
f ( x )dx converges if limx→∞xμ f (x) exists with μ
>1 and ∫a
∞
f ( x )dx diverges if limx→∞xμ f (x) exists and ≠0 with μ ≤ 1.
Example 3.5 ∫0
∞ dxex+1 converges by comparison test
, since 0≤1ex for all x≥0, lim
x→∞
e x
ex+1=lim
x→∞
11+e−x =1 and ∫
0
∞ dxex converges
(need to prove!).
Example 3.6 ∫1
∞ dxx √1+x2 converges by μTest since
1x√1+ x2
is integrable∈[1 ,B ] , for B>1, limx→∞
x2{ 1x √1+ x2 } = 1, μ = 2>1. Note that
1x√1+ x2
is continuous and hence integrable in [1, B]for
B>1.
Example 3.7 ∫0
∞
e−x2dx converges by μTest sincelimx→∞¿ = 0
(verify using L’Hospital’s rule) ,μ = 2>1 and e− x2 is
continuous, and hence integrable, in [0,B] for B>0.
16
Example 3.8 ∫0
∞ x32
3 x2+5dx diverges, since lim
x→∞(x
12 ( x
32
3 x2+5 )) = 1/3, μ
= ½<1 and x32
3x2+5 is continuous, and hence integrable, in
[0,B] for B>0.
TYPE II INTEGRAL
Theorem 3.3 (Comparison test) Let f and g be
integrable in [c, b], for every c, a<c<b. Let g(x)>0 , for
all x, a<x≤b. If lim
x→a+¿f (x)g (x) ¿
¿≠0, then ∫
a
b
f ( x )dx and ∫a
b
g ( x )dx both
converge or both diverge. If lim
x→a+¿f (x)g (x) ¿
¿ = 0 and ∫
a
b
g ( x ) dx
converges, then ∫a
b
f ( x )dx converges.
Theorem 3.4 (μTest ¿ Let f be integrable in[c, b], for
every c, a<c<b. Then ∫a
b
f ( x )dx converges if limx→a+¿( x−a )μ ¿
¿f(x)
exists for 0<μ<1 and ∫a
b
f ( x )dx diverges if limx→a+¿( x−a )μ¿
¿f(x) exists
(≠0) for μ≥1.
17
Example 3.9 ∫0
1 dx(1+x )√x converges , since
limx→ 0+¿ (x−0 )
12 1(1+ x)√x
¿
¿ = 1,
for μ<1 and 1
(1+x)√ x is continuous, and hence integrable,
in [c,1]for 0<c<1.
Example 3.10 ∫12
1 dx√ x(1−x) converges, since
limx→ 1−¿(1−x )
12 1√x(1−x)
¿
¿ = 1,
for μ=12 <1 and 1
√x (1−x) is continuous, and hence integrable,
in [1/2, c]for 1/2<c<1.
THE GAMMA AND BETA FUNCTIONS
Definition (Gamma function) For n>0, Γ(n) = ∫0
∞
e−x xn−1dx.
NOTE: Gamma function is an improper integral of
type I. If 0<n<1, Γ(n) is also an improper integral of
type II. We shall assume convergence of the gamma
function in our course of study.
18
Definition (Beta function) For m, n>0,β(m,n) =
∫0
1
xm−1(1−x )n−1dx
NOTE: Beta function is an improper integral of type II
if either m or n or both lies between 0 and 1 strictly;
otherwise it is a proper integral.
Properties of Gamma and Beta functions
1. For any a>0, ∫0
∞
e−ax xn−1dx = Γ(n)/an.
» let 0<c<d. consider the proper integral I = ∫c
d
e−ax xn−1dx.
Let y = ax. Then I = ∫ac
ad
e− y yn−1
an−1 .1ady = 1an∫
ac
ad
e− y yn−1dy. Thus lim I = Γ (n)an as c→0+¿ and d→∞.
2. Γ(n+1) = nΓ(n)
» Let 0<c<d. using integration by parts on the proper
integral I = ∫c
d
e−x xndx, we get I = (−xn
e x )∨¿cd+n∫
c
d
e− x xn−1dx ¿ = (cn
ec−dn
ed ¿+
n∫c
d
e− x xn−1dx, which tends to nΓ(n) as c→ 0+ and d→∞ (by use
of L’Hospital’s rule). Hence the result.
3 Γ(1) = 1 (can be verified easily)
19
4 Γ(n+1) = n!, for a natural n (follows from property 2
and 3)
4 β(m,n) = β(n,m) (follows using a substitution y = 1-x
after passing to a proper integral)
5 β(m,n) =2∫0
π2
sin2m−1θ cos2n−1θdθ(follows using a substitution x =
sin2θ after passing to a proper integral)
6 β(12,12)=π (follows from definition)
7 β (m,n )=Γ (m )Γ (n)Γ (m+n)
8 Γ ( 12 )=√π
9 For 0<m<1, Γ(m) Γ(1-m) =π cosec(mπ ¿
Example 3.11 ∫0
∞
e−x2dx=√ π2 .
» The range of integration of the given integral is
unbounded but the integrand is continuous
everywhere. For 0<a, ∫0
a
e−x2dx = 12∫0a2
y−12 e− ydy (substituting y =
20
x2 in the proper integral ). Thus lima→∞
∫0
a
e− x2dx = 12 lima→∞∫0
a2
y−12 e− ydy=
12∫0
∞
e− y y12−1dy =12 Γ ( 12 ) = √π /2.
Example 3.12 ∫01 dx
(1−x6)16
dx.
» The integrand has an infinite discontinuity at x =
1.Let 0<c<1. Substituting x3 = sinθ in the proper
integral ∫0c dx
(1−x6)16 , ∫0
c dx
(1−x6)16 = 13 ∫
0
sin−1c3
cos53−1θ sin
13−1θdθ. Since
limc→1−¿∫
0
c dx
(1−x6)16
¿
¿
= 16 β(5/3,1/3) = 16Γ ( 53 )Γ ( 1
3)
Γ (53+ 13)
= 16 .23Γ ( 23 )Γ ¿¿ =
19Γ (13 )Γ (1−13 )= 1
27πcosec( π
3).
PRACTICE SUMS
1. ∫0
π2
√ tan x dx = π√2 .
2. ∫0
1
xm−1(1−x2)n−1dx=12β (m2 , n)
3. ∫0
1
√1−x4dx=[Γ ( 14 )]
2
6√2 π
21
CHAPTER 4: DOUBLE INTEGRAL
Fig. 3:Volume enclosed by the surface z=f(x,y) is
approximated by volume of parallelopipeds
Integrals over rectangles
Let f(x,y) be a bounded function of two independent
variables x and y defined over a closed rectangular
region R: a≤x≤b; c≤y≤d. we take partitions {a = x0,x1,
…,xr-1,xr,…,xn = b} of [a,b] and {c = y0,y1,…,ys-1,ys,…,ym =
d}. These partitions divides the rectangle R into mn
number of subrectangles Rij(1≤i≤n, 1≤j≤m). Let us
choose arbitrarily (α i , β j ¿∈R ijwhere α i∈[ x i−1 , x i] and β j∈[ y j−1 , y j],
22
1≤i≤n, 1≤j≤m. the volume of the parallelepiped with
base Rij and altitude f(α i , β j ¿ is f(α i , β j ¿(xi-xi-1)(yj-yj-1). ∑i , j
f (α i, β j ) (x i−x i−1 )( y j− y j−1), sum of the volumes of all the
parallelepipeds erected over all of the Rij’s, gives an
approximation of the volume enclosed by the curve
and the planes x = a,x = b, y = c, y = d and z = 0. The
approximation can be improved by increasing number
of subrectangles into which R is divided into. Thus the
limit limm→∞,n→∞
∑i , j
f (αi , β j ) (x i−x i−1) ( y j− y j−1), provided it exists, gives
the volume and is represented by ∬R
❑
f ( x , y )dxdy.
NOTE: Every continuous function is integrable over
any rectangle.
Theorem4.1 (equivalence of double integrals with
repeated integrals) If ∬R
❑
f ( x , y )dxdy exists over a rectangle
R: a≤x≤b; c≤y≤d and ∫a
b
f ( x , y )dx exists for each value of y in
23
[c,d], then the repeated integral ∫c
d
dy∫a
b
f (x , y )dxexists and is
equal to ∬R
❑
f ( x , y )dxdy .
Example 4.1 Evaluate ∬R
❑
sin ( x+ y )dx dy over R:0≤x≤π2 , 0≤y≤π
2.
# sin(x+y) is continuous on R, so the double integral
∬R
❑
sin ( x+ y )dx dyexists . Evaluating given double integral in
terms of repeated integrals,
∬R
❑
sin ( x+ y )dx dy = ∫0
π2
dx∫0
π2
sin ( x+ y )dy = ∫0
π2
¿¿¿= 2.
Example 4.2 Evaluate ∬R
❑
(x2+ y2)dx dy over R bounded by y =
x2, x = 2, y = 1.
Fig.4:Fixing up of range of integration of independent
variables x and y
24
¿∬R
❑
(x2+ y2 )dx dy = ∫1
2
dx∫1
x2
(x2+ y2 )dy = ∫1
2
[ x2 y+ y3
3 ]1
x2
dx = ∫1
2
( x63 +x4−x2− 13 )dx = 1006105 .
Example 4.3 Prove that ∬R
❑
xy ( x+ y )dx dy= 536 ,where R is the
region bounded by y = x and y = x2.
CHAPTER 5: EVALUATION OF AREA
Cartesian co-ordinate
It has already been seen that area of the region
bounded by the curve y = f(x), lines x = a, x = b and y =
0 is given by ∫a
b
f (x )dx, provided it exists. Similarly area of
the region bounded by the curve x = g(y), lines y = c, y
= d and x = 0 is given by ∫c
d
g ( y )dy, provided it exists. We
can define F:[a,b]→ R by F(t)=∫a
t
f ( x )dx ,a≤t≤b.
25
Observe that for the function f[-2,2]
→R , f (x )=x2−x+3 , F : [−2,2 ] →Rcanbe depictedas..\Documents\30.gif.
Fig.5: Sign convention of area calculation and area
bounded by two curves
Example 5.1 find the area of the bounded region
bounded by the curve y = x(x-1)(x-2) and the x-axis...\
Documents\31.gif
# y value is positive for 0<x<1 and for x>2 and is
negative for 1<x<2; the curve cuts the x-axis at points
whose abscissa are 0,1,2. Required area = ∫0
1
y dx+|∫12
y dx|.Example 5.2 ..\Documents\x5.mw
26
Example 5.2 Find the area of the bounded region
bounded by the curves y = x2 and x = y2.
# On solving the given equations of the curves, the
point of intersection of the two curves are (0,0) and
(4,4). Thus required area = ∫0
4
√4 xdx−∫0
4 x2
4dx.
Example 5.3 Find the area of the loop formed by the
curve y2 = x(x-2)2
# The abscissa of points of intersection of the curve
with the x-axis are given by y = 0, that is, x = 0,2,2. For
x<0, no real value of y satisfy the equation. Hence no
part of the curve exist corresponding to x<0.
Corresponding to each x-value satisfying 0<x<2, there
exist two values of y, equal in magnitude and opposite
in sign. Thus between x = 0 and x = 2, the curve is
symmetric about the x-axis and a loop is formed
thereby. For x>2, y→∞ as x→∞. the required area = 2
∫0
2
(x−2)√ x dx (by symmetry of the curve about x-axis).
27
Example 5.4 Prove that area included in a circle of
radius r unit is πr2 square unit.
# We can choose two perpendicular straight lines
passing through the centre of the circle as co-ordinate
axes. With reference to such a co-ordinate system,
equation of the circle is y = ±√r2−x2.Curve is symmetric
about the axes .Thus required area = 4∫0
r
√r2− x2dx.
Polar co-ordinate
Fig.6:Area under a polar curve and that under a
cardioide
The area of the region bounded by the curve r = f(θ) ,
the radius vector θ=α,θ=β is given by 12∫αβ
r2dθ.
28
Example 5.5 Find the area enclosed by the cardioide r
= a(1+ cosθ)
# As θvaries from 0 to π2 , r decreases continuously from
2a to a. When θfurther increases from π2 to π, r
decreases further from a to 0. Also the curve is
symmetric about the initial line (since the equation of
the curve remains unaffected on replacing θ by –θ ¿.
Hence the area enclosed by the curve = 2.a2
2 ∫0π
(1+cosθ)2dθ.
Fig.7:Area between two cardioides
Example 5.6 Find the area enclosed by the cardioide r
= a(1+ cosθ) and r = a(1- cosθ)
# The vectorial angle corresponding the points of
intersection of the curves are θ=π2 and θ=-π2 . Because of
29
the symmetry of the curves about the initial line,θ=π2 ,θ
=-π2 and θ=π , required area is 4.a22 ∫0
π2
¿¿¿.
Notation: For the rest of our discussion , suffix of dependent variable or of a function denotes order of differentiation.
CHAPTER 6: LENGTH OF ARC OF A PLANE CURVE
The length of the arc of a curve y = f(x) between two
points whose abscissae are a and b, when f1 (x) is
continuous on the interval [a,b] is given by ∫a
b
√1+{f 12(x) }dx.
We can define the arc length function s:[a,b]→ R by
30
s(t)=∫a
t
√1+{f 12(x) }dx ,a≤t≤b. ..\Documents\x6.mw .
Example 6.1 Find the length of circumference of a
circle of radius 5 units.
# With suitable choice of co-ordinates, equation of the
curve is x2+y2 = 25. Here 1+y12 = 25/(25-x2). Required
length = 4∫0
5 5√25−x2 = 10π.
Example 6.2 Find the length of the arc of the parabola
y2 = 16 x measured from the vertex to an extremity of
the latus rectum.
Example 6.3 Let f(x)=x3+1,x∈[-2,2].In adjoining figure,
f is plotted in red, integrand in blue and the arc
function s(x)=∫−2
x
√1+{f 12(t)}dt, x∈[-2,2], in green...\
Documents\32.gif
Example 6.3 Find the length of the loop of the curve
3y2 = x(x-1)2.
31
The length of an arc of the curve r = f(θ) between points
whose vectorial angles are θ=α and θ=β is given by
∫α
β
√r2+( drdθ )2
dθ.
Example 6.4 Find the length of the perimeter of the
cardioide r = a(1+cosθ ¿
# As we can see from fig.7 above, using symmetry of
the curve about the initial line (justified by the
invariance of the equation on replacing θ by –θ), the
required length is 2x[∫0π2
√a2(1−cosθ)2+a2 sin2θdθ+∫π2
π
√a2(1−cosθ )2+a2 sin2θdθ ].
32
CHAPTER 7: SURFACE AREA AND VOLUME OF
SURFACE OF REVOLUTION
The volume generated by revolving about the x-axis an
area bounded by the curve y = f(x), the x-axis and two
ordinates x = a and x = b is given by V=π∫a
b
{ f (x)}2dx.
Similarly, The volume generated by revolving about
the y-axis an area bounded by the curve x = g(y), the y-
axis and two ordinates y = c and y = d is given by
V=π∫c
d
{g( y )}2dy. The formula for surface area of the surface
generated is given by S= 2π∫a
b
y √1+ y12dx and 2π∫
c
d
x √1+x12dy.
Example 7.1 Find the volume and surface area of a
right circular cylinder of base radius r and altitude h.
Fig.8:Volume and surface area of surface of revolution
33
# A right circular cylinder of radius r and altitude h is
obtained by revolving y = r,z = 0 about the x-axis. Thus
volume = π∫o
h
r2dx = πr2h and surface area = 2π∫0
h
r √1+02dx = 2π
rh.
Example 7.2 Find the volume and surface area of a
sphere of radius r.
# A sphere of radius r is obtained by revolving y=√r2−x2
,-r≤x≤r, about x-axis. Thus V=π∫−r
r
(r2−x2)dx and S=2
π∫−r
r
√r2−x2 r√r2−x2
dx=4πr2.
Example 7.3 Surface area of the surface generated by
y=cos (x3) between x=0, x=2 for rotation about x-axis
and y-axis(frustum view)...\Documents\33.gif
Example 7.4 Surface area of the surface generated by
revolving y=cos x between x=0 and x=2 around x-axis:
..\Documents\x3.mw
34
Example 7.5 Surface area of the surface generated by
revolving y=1+ cos x between x=0 and x=2 around x-
axis:..\Documents\x4.mw
SECTION II (ORDINARY DIFFERENTIAL
EQUATIONS)
NOTE: suffix denotes order of differentiation
TEXT : (1) Differential equations and their
applications—Zafar
Ahsan
(2) Differential Equations—Richard
Bronson (Schaum)
35
CHAPTER 1: ORDER,DEGREE AND FORMATION OF
ORDINARYDIFFERENTIAL EQUATION
Definition 1.1 An ordinary differential equation(ODE)
is an equation involving derivative(s) or differentials
w.r.t. a single independent variable.
Example: y2=a2y, 3x3y3+4x2y1 = 7x2+9, x dy – ydx+2xy
dy = 0.
Definition 1.2 the order of an ODE is the order of the
heighest ordered derivative occurring in the equation.
The degree of an ODE is the largest power of the
heighest ordered derivative occurring in the equation
after the ODE has been made free from the radicals
and fractions as far as the derivatives are concerned.
Example: y = xdydx
+ dxdy : first order, second degree
(1+y12)3/2 = y1: first order, sixth degree.
36
Definition 1.3 An ODE in which the dependant
variables and all its derivatives present occur in first
degree only and no products od dependent variables
and /or derivatives occur is known as a linear ODE. An
ODE which is not linear is called nonlinear ODE. Thus
y1 = sin x+cos x is linear while y = y1+1/y1 is non-linear.
Definition 1.4 A solution of an ODE is a relation
between the dependent and independent variables,
not involving the derivatives such that this relation
and the derivatives obtained from it satisfies the given
ODE. For example, y = ce2x is a solution of the ODE y1 =
2y, since y1 = 2ce2x and y = ce2x satisfy the given ODE.
Some of the solutions of the DE y2=x2-x+1 are depicted
as follows: ..\Documents\x1.mw . Some of the
solutions y=c cos x+d sin x (c varying between -20 to
10, d between -10 to 20) of y2+y=0 are shown here: ..\
Documents\x2.mw
FORMATION OF ODE
37
GEOMETRIC PROBLEMS
Example 1.1: Let a curve under Cartesian coordinate
system satisfy the condition that the sum of x- and y-
intercepts of its tangents is always equal to a. Find the
ODE corresponding to the curve.
# Equation of tangent at any point (x,y) on the curve is
Y-y = dydx (X-x). Thus the differential equation expressing
the given condition is (x-yx1)+(y – xy1) = a.
PHYSICAL PROBLEMS
Example 1.2: Five hundred grammes of sugar in water
are being converted into dextrose at a rate which is
proportional to the amount unconverted. Form an ODE
expressing the rate of conversion after t minutes.
# Let y denote the number of grammes converted in t
minutes. Then the ODE is dydt = k(500 – y), k constant.
38
Example 1.3 A person places Rs.20,000 in a savings
account which pays 5% interest per annum,
compounded continuously. Let N(t) denotes the
balance in the account at any time t. N(t) grows by the
accumulated interest payments, which are
proportional to the amount of money in the account.
The resulting ODE is dNdt =0.05 N.
Example 1.4 Five mice in a stable population of 500
are intentionally infected with a contagious disease to
test a theory of epidemic spread that postulates the
rate of change in the infected population is
proportional to the product of the number of mice who
have the disease with the number that are disease free.
Let N(t) denote the number of mice with the disease at
time t. resulting ODE is dNdt = kN(500-N).
Example 1.5 Newton’s law of cooling states that the
time rate of change of the temperature of a body is
39
proportional to the temperature difference between
the body and the surrounding medium. Let T denote
the temperature of the body and let Tm denote the
temperature of the surrounding medium. Resulting
ODE is dTdt = -k(T-Tm).
ELIMINATION OF ARBITRARY CONSTANTS
Suppose we are given an equation(not a differential
equation) containing n parameters (arbitrary
constants).by differentiating the given equation
successively n times , we get n equations more
containing n parameters and derivatives. By
eliminating n parameters from the above (n+1)
equations and obtaining an equation which involves
derivatives upto the n th order, we get an ODE of
order n.
40
Example 1.6: Obtain an ODE by eliminating
parameters a,b from (x-a)2+(y-b)2 = c2, c constant.
(x-a)+(y-b)y1 = 0, 1+y12+(y-b)y2 = 0. Obtain values of x-
a and y-b in the given equation, obtain the ODE.
Example 1.7: Find the ODE corresponding to the
family of curves y = c(x-c)2 (2.1), where c is a
parameter.
Y1 = 2c(x-c) . so y12 = 4c2(x-c)2 (2.2). from (2.1) and
(2.2), c = y12/(4y). Substituting in (2.1), required
equation is 8y2 = 4xyy1-y12.
Example 1.8: Obtain the ODE corresponding to the
family of all circles each of which touches the axis of x
at the origin.
# The equation of an arbitrary circle of the family is of
the form (x-r)2+y2 = r2, that is , x2+y2 – 2rx = 0, r is a
parameter. Obtain the ODE by eliminating the
parameter r.
41
GENERAL, PARTICULAR AND SINGULAR
SOLUTIONS
A solution which contains a number of independent
parameters equal to the order of the ODE is called the
general solution or complete integral of the ODE. A
solution obtained from the general solution by putting
particular values to at least one of the parameters
present in the general solution is called a particular
solution. A solution which is neither the general
solution nor is a particular solution of an ODE is called
a singular solution of the ODE.
NOTE (1) In counting the parameters in the general
solution, we must check whether they are independent
and are not equivalent to a lesser number of
parameters. Thus a solution of the form c1 cosx+c2
sin(x+c3) appears to have three parameters but
42
actually they are equivalent to two; for, c1 cosx+c2
sin(x+c3) = (c1+c2 sin c3) cos x+c2 cos c3 sinx = A cos x
+B sin x.
(2) The general solution of an ODE may have more
than one form, but parameters in one form will be
related to parameters in another form. Thus y = c1
cos(x+c2)and y = c3 sin x+c4 sin x are both solutions of
y2+y = 0. Here c4 = c1 cos c2,c3 = -c1 sin c2.
Example 1.9 y = ln x is a solution of xy2+y1 = 0 on
(0,∞) but is not a solution on (-∞,∞) since ln x is not
defined on (-∞,0].
EXERCISES
1. Find the ODE of the family of circles whose centres are
on the y-axis and touch the x-axis.
2. Find the ODE of the family of parabolas whose axes are
parallel to the y-axis.
43
3. A canonical tank of height 20 cm. loses water out of an
orifice at its bottom. If the cross-sectional area of the
orifice is ¼ cm2, obtain the ODE representing the
height h of water at any time t.
CHAPTER 2: ODE OF FIRST ORDER AND FIRST
DEGREE
An ODE of the first order and first degree is of the form
y1 = f(x,y) which is sometimes conveniently written as
M(x,y) dx+N(x,y) dy = 0 (2.1), where M(x,y) and N(x,y)
are either functions of x,y possessing partial
derivatives or constants.
NOTE: First order first degree ODE whose variables
are separable or which are homogeneous or which are
44
reducible to homogeneous equations have been
studied earlier and will not be repeated here. Students
are advised to revise the same.
Exact equations
Definition 2.1 An ODE of the form (2.1) is called exact
if and only if there exists a function f(x,y) with
continuous partial derivatives such that (2.1) can be
written in the form df(x,y) = 0, that is, in the form fx
dx+fy dy = 0. In this case, f(x,y) = c will be the general
solution of given ODE.
For example, xdx+ydy = 0 is exact since the given
equation can be written in the form d(x2+y2) = 0.
Theorem2.1 (2.1) is exact iff My = Nx holds.
Example 2.1 Verify whether the differential equation
(sin x cos y+e2x) dx+(cos x sin y+tan y)dy = 0 is exact.
Also find the general solution of the equation.
45
# Comparing with the form Mdx+Ndy = 0, M = sin x cos
y+e2x, N = cos x sin y+tan y. so My=- sin x siny = Nx.
hence given equation is exact. Let f(x,y) = 0 be the
solution. Given equation will be identical with
fxdx+fydy = 0. Comparing, fx = sin x cosy+e2x(2.2), fy =
cos x sin y+tan y(2.3).From (2.2), f(x,y) = - cos x cos
y+1/2 e2x+g(y), where g(y) is the integration constant.
Substituting in (2.3), cos x sin y+ g1(y) = fy = cosx sin y+
tan y. Thus g1(y) = tan y so that g(y) = ln sec y. Hence
the general solution is f(x,y) = - cos x cos y+1/2 e2x+ln
sec y = c, c parameter.
Integrating factors
In theory, a non-exact ODE can always be made exact
by multiplying the equation by some function μ(x,y) of
x and y. such a function μ(x,y) is called an Integrating
factor (I.F.) of the ODE. Although there is always an
46
I.F. for a non-exact ODE, there is no general method of
finding the I.F. we shall now discuss methods of finding
I.F. in some particular cases .
Method 1 (by inspection) use is made of exact
differentials like
d ( yx )= xdy− ydxx2 , d(xy) = x dy+y dx ,d( ln ( yx ))= x dy− y dx
xy , d
( tan−1 yx )= x dy− y dx
x2+ y2.
Example 2.2: (1+xy)y dx+(1-xy)x dy = 0.
# d(xy)+ xy(y dx-x dy) = 0 . thus d (xy )(xy )2
− x dy− y dxxy
=¿ 0. So
d [−1xy −ln ( yx )]= 0. General solution is 1/(xy)+ln(y/x) = c, c
parameter.
Example 2.3: x dy-y dx = a(x2+y2)dy
d(tan-1(y/x)) = a dy.
Method 2 if the ODE Mdx+Ndy = 0 is homogeneous
and Mx+Ny≠0, then 1/(Mx+Ny) is an IF.
47
Example 2.4: x2y dx – (x3+y3)dy = 0
# Here Mx+Ny = -y4≠0 and the given equation is
homogeneous. Thus an IF is -1/y4.
Method 3 if M = y f(xy), N =x g(xy), then 1/(Mx-Ny) is
an IF of the ODE Mdx+Ndy = 0.
Example 2.5: (xy+2x2y2) y dx+(xy – x2y2) x dy = 0.
# Here M = y[(xy)+2(xy)2] and N = x[(xy)- (xy)2]. Thus
an IF = 1/[Mx-Ny] = 1/(3x3y3).
Method 4 e∫ f ( x ) dx is an IF of Mdx+Ndy =0,if M y−N x
N is a
function of x.
Example 2.6: (x2+y2)dx – 2xy dy = 0
# M = x2+y2, N = -2xy, M y−N x
N =-2/x. Thus an IF = e∫ f ( x )dx=e-2 ln
x = x-2.
Example 2.7 (a very useful form of ODE: linear
equation in y) y1+P(x) y = f(x), P(x), f(x) are functions
48
of x alone or constants. M = P(x)y-f(x), N = 1. ThusM y−N x
N
= P(x). An IF = e∫ P ( x )dx.
Method 5 e∫ f ( y ) dy is an IF of Mdx+Ndy =0,if N x−M y
M is a
function of y.
Example 2.8 (xy3+y) dx+2(x2y2+x+y4) dy = 0
¿N x−M y
M=1/y. hence an IF = e∫
1y dy=y.
Method 6 if the equation Mdx +Ndy = 0 is of the form
xayb(my dx+nx dy)+xcyd(py dx+qx dy) = 0, where
a,b,c,d,m,n,p are constants,then xhyk is an IF, where h,k
are constants and can be obtained by applying the
condition that after multiplication by xhyk the given
ODE becomes exact.
Example 2.9 (y2+2x2y)dx+(2x3-xy)dy = 0
# Let xhyk be an IF. Multiplying the given equation by
this factor, we have (xhyk+2+2xh+2yk+1)dx+(2xh+3-
xh+1yk+1)dy = 0 (2.4). Since (2.4) is to be exact, from the