Updated Engineering Maths Year 3.first Semester
Updated Engineering Maths Year 3.first Semester
MATHEMATICS
FENG 310 LECTURER: Sallu Jaward
Phone:+23275727287
E-mail: [email protected]
[email protected]
Course Material:
ü Course material, including lecture notes, PowerPoint slides, sample exams,
Home work solutions, hints, etc … will be available to the students
throughout the semester.
ü Any student registered for the class should have access to all these materials.
2
q ADVANCED ENGINEERING MATHEMATICS by K.A. Stroud and Dexter J. Booth,
4th edition, McGraw-Hill Education, New York, NY. ISBN: 978-1-259-69653-4.
q ADVANCED ENGINEERING MATHEMATICS by Dennis G. Zill, 6 th edition, |
Burlington, MA : Jones & Bartlett Learning, [2017]. ISBN: 1284105903 (casebound).
q CALCULUS, EARLY TRANSCENDENTAL FUNCTIONS by Ron Larson, Robert
Hostetler and Bruce H. Edwards 4th edition, Houghton Mifflin , ISBN 13: 978-0-618-
73069-8 ISBN 10: 0-618-73069-9
q MULTIVARIABLE CALCULUS by Ron Larson 11th edition
Homework
ü A list of weekly problem sets along with due dates can be found throughout the
lecture slides. Although you may discuss homework problems with me, or your
classmates, you are expected to work out the problems independently
ü There is no specific format for the Home Work assignments. Please write legibly.
ü You are expected to turn in your assignments on the due date, typically by 12
noon of that day unless noted otherwise on the Assignment.
ü Late Home Work submissions will receive a zero grade unless permission is
obtained from the instructor before the due date.
4
Academic Honesty
ü You may discuss homework problems with others, but you must
not retain written notes from your conversations with other
students, or share data electronically (e.g. files, emails, etc …) to be
used in completing your homework.
ü Your written work must be completed without reference to such
notes, with the exception of class notes.
5
Classroom etiquette
üPlease try to arrive on time for each class because late arrivals are
distracting.
üIf you arrive late please enter quietly and take a seat near the back of
the room.
üPlease do not converse with your classmates during class.
üAll electronic devices (cell phones, iPods, iPads, laptops, …) must be
turned off
6
Exams
ü In addition to a 3-hour final exam, there are three tests during the
semester.
ü In case you miss a test, your score on that test will be zero without an
option for a make-up.
ü All tests and exams will be closed book.
ü The regular exams will be given at the end of the semester.
7
Grading
ü Each of the three tests will count as 15 points with the average of the scores
counting toward the final grade.
ü The final exam will count as 70 points, and the homework will count collectively
as 10 points, for a total of 95 points.
ü Regular attendance will count for 5 points.
ü Your final grade for the course will be based primarily on your total point score,
but other appropriate factors, such as participation and personal initiative, will
also be taken into account.
8
Derivatives of Functions
Derivatives of Functions
Integration Formulas
Integration Formulas
OUTLINE OF TOPICS TO BE COVERED
1. Introduction to Differential Equations 5. Applications of Differential Equations
ü Definition and classification (order, degree, linearity) ü Mechanical vibrations (damped and undamped)
ü General and particular solutions ü Electrical circuits (RLC circuits)
2. First-Order Differential Equations ü Population models
ü Separable equations 6. Laplace Transforms
ü Homogeneous equations ü Definition and basic properties
ü Exact equations and integrating factors ü Inverse Laplace transform
ü Linear differential equations ü Laplace transform of common functions
3. Higher-Order Differential Equations ü Step functions and Dirac delta function
ü Second and higher-order linear differential equation 7. Applications of Laplace Transforms
ü Solutions using characteristic equations ü Solving differential equations using Laplace transforms
ü Method of undetermined coefficients ü Initial value problems
ü Variation of parameters ü Systems of linear differential equations
4. Systems of Differential Equations 8. Introduction to Partial Differential Equations
ü Introduction to systems of first-order differential ü Brief Introduction to Partial Differential Equations
equations ü Method of Separation of Variables
LECTURE 1
Introduction to
Differential Equations
15
Definition and classification
Definition of Differential Equation: An equation containing the derivatives of one or more dependent
variables, with respect to one or more independent variables, is said to be a differential equation
(DE). To discuss these equations, we will categorize each differential equation by its type, order, and
linearity.
Classification by Type
Partial Differential Equations (PDEs): An equation involving only partial derivatives of one or
more functions with respect to two or more independent variables. Example:
The following example demonstrates both types of differential equations (ODEs and PDEs).
16
Examples: Types of Differential Equations
�� −� �2 � �� �� ��
+ 6� = � , + − 12� = 0, + = 3x+2y
�� ��2 �� �� ��
�2 � �2 � �2� �2� �� �� ��
+ = 0, = − , =− ,
��2 ��2 ��2 ��2 �� �� ��
are examples of partial differential equations. Notice that in the third equation, there are two
dependent variables and two independent variables.
This indicates that u and v must each be functions of two or more independent variables.
17
Classification by Order
The order of a differential equation (ODE or PDE) is the order of the highest derivative in the equation.
EXAMPLE: Order of a Differential Equation
�2� �� 3 �4 � �2 �
The differential equations 2 + 5 − 4� = �� , 2 4 + 2 = 0 are examples of a second-order ordinary differential equation
�� �� �� ��
and a fourth-order partial differential equation, respectively.
A first-order ordinary differential equation is sometimes written in the differential form M(�, �)�� + �(�, �)�� = 0
EXAMPLE 3 Differential Form of a First-Order ODE: Assuming y is the dependent variable in a first-order Ordinary Differential
Equation (ODE), we recall from calculus that the differential dy is defined as: �� = �’ ��
a) By dividing both sides by the differential dx, an alternative form of the equation (� − �)�� + 4��� = 0 is given by: � − � +
�� ��
4 = 0 �� 4 + � = �
�� ��
��
b) By multiplying the differential equation 6�� + �2 + � 2 = 0
��
An n-th order ordinary differential equation with one dependent variable can be generally expressed as:
� �, �, �,’ . . . . . , �� = 0
���
= �(�, �, �,’ . . . . . , �� )
�� �
Classification by Linearity
��� � � − 1� ��
�� (�) � + ��−1 (�) � − 1 + , . . . , �1 (�) + � (�)� = �(�)
0
�� �� ��
The dependent variable y and all its derivatives y, y, … , y(n) are of the first degree; that
is, the power of each term involving y is 1.
The coefficients a0, a1, … , an of y, y , … , y(n) depend at most on the independent
variable x
EXAMPLE: Linear and Nonlinear Differential Equations
5. + ���� − 5� = 0 �� 3
��2
10. �� = �3�
��
Solution of an ODE
verify that the indicated function is an explicit solution of the given differential equation.
Practice problems set 1
Verify that the indicated function is an explicit
solution of the given differential equation. Find values of m so that the function � = ��� is a solution
of the given differential equation.
−�
1) 2�′ + � = 0, �=� 2
1) �′ + 2� = 0, 2. �′′ − 5�′ + 6� = 0
2) �′′ − � = 0 � = �� 3. 2�′′ + 9�′ − 5� = 0 4. 3�′ = 4�
3) �′′ − 6�′ + 13� = 0; � = �3� ���2� find values of m so that the function y= �� is a solution
of the given differential equation.
4) �′′ + � = tan � ;
� = −(cos � )ln (sec � + tan � )
1. �2 �′′ − 7��′ + 15� = 0 2. x�′′ + 2�′ =
0
5) (1 − 2� − �2 )�′′ + 2(1 + �)�′ − 2� = 0;
�=�+1
LECTURE 2
First-Order Differential Equations
vSeparable equations
vHomogeneous equations
v Exact equations and integrating factors
vLinear differential equations
Separable equations
A first-order differential equation is called separable if all terms involving � and y can be
grouped with dx and dy , respectively. In other words, a differential equation of the form
��
= g(x). h(y)
��
is said to have separable variables.
�� ��
1. =P−� 2
6. = 4(�2 + 1), x( � 4)=1
�� ��
� −� ) �� 2 �� ��+3�−�−3
2. (� + �
��
=� , 7. = .
�� ��−2�+4�−8
2 �� �� ��+2�−�−2
3. � = � − ��, �(−1) =− 1. 8. = .
�� �� ��−3�+�−3
�� �+1 2 �� �2 −1
4. yln � = 9. = , y(2)=2
�� � �� �2 −1
2 )1 2 �� 2 )1 2
� )2 −� � )3 −� 10. � (1 + � = y (1 + � ,
5. (1 + � � �� + (1 + � � �� = 0 ��
Homogeneous equations
A differential equation is said to be a homogeneous DE if every term has
the same degree of
Homogeneous DE can be written as: M(x,y)dx + N(x,y)dy=0
If M is simpler use: x=vy, dx=vdy+ydv
If N is simpler use: y= vx, dy=vdx+xdv
��
Homogeneous DE = �(�, �)
��
�� � �
=� , let v= , y=vx
�� � �
Classwork set
Solve the following DE
Section A Section B
1. (�� + �� + �� )�� − �� �� = � ��
1. (� + �) �� = � − �
�
� ��
2. � = �� + �� � �
��
�� 2. �� = �� + ���
��
��
3. ��� �� = �� + ��� ��
3. ��� �� = �� + ���
��
4. � �� = � + � �� ��
4. �� �� = �� + � ��� + ��
��
5. (� − �) �� = � + � ��
5. � �� + � = �� + ��
Exact Equations
Exact equations are a special type of first-order differential equation where the left-hand side of the
equation is the exact differential of some function f(x,y). This allows us to integrate directly to find a
solution
Recognizing an Exact Equation
A first-order DE of the form:
�(�, �)�� + �(�, �)�� = �
�� ��
is exact if there exists a function f(x,y) such that: = �(�, �) and = �(�, �)
�� ��
Criterion for Exactness
For M(x,y) and N(x,y) with continuous first partial derivatives, the equation is exact if:
�� ��
=
�� ��
This is a necessary and sufficient condition for the equation to be exact. If this condition holds, the
equation represents the total differential of a function f(x,y)
Exact Equations: Method of Solution
1. Verify Exactness
�� ��
Check whether the condition = is satisfied.
�� ��
2. Determine f(x,y)
Ø Integrate M(x,y) with respect to x, treating y as a constant, to find f(x,y) (up to a function of y, g(y)):
�(�, �) = �(�, �)�� + �(�)
Ø Differentiate f(x,y) with respect to y and equate it to N(x,y) to find g′ (y)
��
= �′ (�) + ( ����� ���� �)
��
Ø Integrate g′ (y) to find g(y)
3. General Solution
Combine results to express the solution implicitly as
�(�, �) = �
where C is an arbitrary constant.
Practice problems set 3
Determine whether the given differential equation is exact. If it is exact, solve it.
1) (2� − 1)�� + (3� + 7)��� = 0
2) (2� − �)�� − (� + 6�)�� = 0
3) (���� − �����)�� + (���� + ����� − �)�� = 0
4) (2��2 − 3)�� + (2�2 � + 4)�� = 0
1 �� �
5) 2� − + ���3� + 2 − 4�2 + 3����3� = 0
� �� �
�
6) (1 + ln � + � )�� − (1 − ln � )�� = 0
3 �� 3
7) 1 − � + � ��
+�= −1
�
1 1 � � 1
8) �
+
�2
−
�2 +�2
�� + (�� +
�2 +�2
)�� = 0
� ��
a) Solve + ���� − 2�� = �(� + ����), y(0)=1
1+�2 ��
Integrating Factors for Nonexact Equations
Integrating factors are used to transform a nonexact differential equation into an exact one, enabling us
to solve it using the techniques for exact equations.
the equation is nonexact if the criterion for exactness does not hold:
�� ��
≠
�� ��
To make it exact, we multiply through by an integrating factor μ(x,y), such that:
2. Simplifying Assumptions
If �(�, �) is assumed to depend only on � or �, the equation simplifies:
• � = �(�):
�� �� ��−�� ��
The equation becomes: =�
�� �
• � = �(�):
�� �� ��−�� ��
The equation becomes: = �
�� �
3. Conditions for Use:
v if (�� ��−�� ��) � depends only on � ,use �(�)
v If (�� ��−�� ��) � depends only on �, use �(�)
Practice problems set 4
Verify that the given differential equation is not exact, apply the integrating factor, and solve the resulting
exact equation.
General Solution
1. Identify P(x) and Q(x)
2. Find the integrating factor �(�):
�(�) = � �(�)�� .
3. Multiply through by �(�) to make the left-hand side exact:
��
�(�) + �(�)P(�)� = �(�)Q(�)
��
4. Recognize the left-hand side as a derivative
�
[�(�)�] = �(�)Q(�)
��
5. Integrate both sides:
�(�)� = �(�)� ��
Practice problems set 5
LECTURE 3
Higher-Order Differential Equations
1, �′′ − � = 0 �1 = ��
2. �2 �′′ − 3��′ + 4� = 0 �1 = �2
4. �′′ − � = 0 �1 = ��
5. �′′ + 2�′ + � = 0 �1 = ��−�
6. �′′ + 16� = 0 �1 = cos 4�
7. �′′ + 9� = 0 �1 = sin 3�
8. �′′ − � = 0 �1 = cosh �
9. (1 − 2� − �2 )�′′ + 2(1 + �)�′ − 2� = 0 �1 =�+1
Superposition Principle for Homogeneous
Equations
The Superposition Principle of homogeneous equations states that if �� (�), �� (�), . . . , �� (�) are solutions to a
homogeneous linear differential equation
�(�) = �� �(�) + ��−1 �(�−1) + . . . . + �1 �’ + �0 � = 0
is also a solution of the homogeneous equation �(�) = �. The solutions to the homogeneous equation can be
combined linearly to form the general solution to the homogeneous equation.
Superposition Principle for Nonhomogeneous
Equations
The Superposition Principle for nonhomogeneous equations states that if ��� (�), ��� (�), . . . , ��� (�) are
particular solutions of k different nonhomogeneous linear differential equations, corresponding to different
forcing functions �� (�), �� (�), . . . , �� (�)), then the linear combination of these particular solutions:
where �� (�)is the general solution of the corresponding homogeneous equation �(�) = �.
Homogeneous Linear Equations with Constant
Coefficients
For linear higher-order differential equations with Auxiliary Equation
By assuming a solution of the form , we substitute into
constant coefficients, all solutions can be expressed
the differential equation, which leads to
as exponential functions or combinations of ��2 ��� + ����� + ���� = 0
or equivalently
exponential functions.
��� (��2 + �� + �) = 0
General Form: Since ��� ≠ �, the term in parentheses must vanish
��2 + �� + �) = 0
The homogeneous linear equation with constant
This quadratic equation is called the auxiliary equation.
coefficients is written as: Solving it provides the roots, �� and �� , which
determine the form of the general solution.
(�) (�−1) ’
�� � + ��−1 � + . . . . + �1 � + �0 � = 0
where �� are constants, �� ≠ �
For a second-order equation
��′′ + ��′ + �� = 0
Three Cases of Roots
� = �1 ��1� + �2 ��2�
2. Repeated Real Roots �� − ��� = � :
If the roots �� = �� = � are real and equal, the general solution is
� = (�1 + �2 �)���
3. Complex Conjugate Roots �� − ��� < � :
If the roots are complex �� = � + �� and �� = � − ��, the general solution is expressed in terms of real
functions:
� = ��� (�1 cos �� + �2 sin ��)
Here, � represents the exponential decay/growth, while β governs the oscillatory
behavior.
Classwork and Practice problems set 7
find the general solution or the particular solution of the given second-order differential equation
�4 � �2 � +
vi. + 2 2 �=0
6. �′′ + 4�′ + 7� = 0 ��4 ��
�4 � �2 � +
7. �′′ − 10�′ + 25� = 0 vii. 16 + 24 2 9� = 0
��4 ��
8. �′′ + 8�′ + 25� = 0 a) Two roots of a cubic auxiliary equation with real coefficients are �1 =−
1
and �2 = 3 + �. What is the corresponding homogeneous linear
9. �′′′ + 3�′′ − 4� = 0 2
differential equation?
10. 4�′′ − 5�′ + 6� = 0
LECTURE 4
Method of undetermined coefficients
The method of undetermined coefficients is a practical technique for solving nonhomogeneous linear differential
equations of the form
If �(�) = ��� guess �� = ���� Ø Substitute �� and its derivatives into the given equation.
1. Consider the differential equation ��′′ + ��′ + �� = ��� , where a, b, c, and k are constants. The auxiliary
equation of the associated homogeneous equation is ��2 + �� + � = 0
a) If k is not a root of the auxiliary equation, show that we can find a particular solution of the form �� = ���� ,
c) If k is a root of the auxiliary equation of multiplicity two, show that we can find a particular solution of the
form �� = ��2 ��� , where A = 1 (2�).
LECTURE 5
Variation of Parameters
The general solution is
�(�) = �� (�) + �� (�) Particular Solution ��
of a nonhomogeneous linear differential equation. Assume in the form:
where �� (�) is the complementary solution and �� (�) is the �� (�) = �1 (�)�1 (�) + �2 (�)�2 (�)
particular solution
where �� (�) �nd �� (�) are functions to be determined
General form Using the product rule to differentiate �� twice
�2 (�)�′′ + �1 (�)�′ + �0 (�)� = �(�)
Key Results:
or in standard form −�(�)�2 �(�)�1
�′1 = �′� =
�′′ + �(�)�′ + �(�)� = �(�) �(�1 , �2 ) �(�1 , �2 )
Assumes �(�), �(�) ��� �(�) are continuous on an interval �
where �(�� , �� ) is the Wronskian
Solve the associated homogeneous equation
�1 �2
′′ ′ W= �′
� + �(�)� + �(�)� = 0 �′ 2
1
The complementary solution is:
�� = �1 �1 (�) + �2 �2 (�) Integrate to find �� (�) �nd �� (�)then substitute back to find
where �1 (�) ��� �2 (�)are linearly independent solutions. ��
The general solution is
�(�) = �� (�) + �� (�)
Classwork and Practice problems set 9
solve each differential equation by variation of parameters.
�
� ��
�−1
� � − 1� ��
�� � + � �−1 � + , . . . , � 1 � + � � = �(�),
0
�� � �� � − 1 ��
where �� , ��−� , �� , . . . , �� are constants
2 � � − 2� �� − 4� = 0
2 1. (� + 3)2 �′′ − 8(� + 3)�′ +
1. � 2
�� �� 14� = 0
2 � � + 8� �� + � = 0
2
2. 4� 2
�� ��
2. (� + 2)2 �′′ + (� + 2)�′ + � =
0
3. 4�2 �′′ + 17� = 0, �(1) =− 1, �′ (1) =
1 3. (� − 4)2 �′′ − 5(� − 4)�′ +
− 9� = 0
2
CONT. Practice problems set 10
SECTION C SECTION D
use the substitution � = �� to transform find a homogeneous Cauchy–Euler differential
the given Cauchy–Euler equation to a equation whose general solution is given.
differential equation with constant
coefficients. Solve the original equation 1) � = �1 �4 + �2 �−2
by solving the new equation
2) � = �1 + �2 �5
1. �2 �′′ + 9��′ − 20� = 0
3) � = �1 �−3 + �2 �−3 ln �
2. �2 �′′ + 10��′ + 8� = �2
4) � = �1 cos( ln �) + �2 sin (ln � )
3. �2 �′′ − 4��′ + 6� = ln �2
1 1 1 1
5) � = �1 � 2 cos(
2
ln �) + �2 � 2 sin (
2
ln � )
LECTURE 6
Differential Equation: Power Series Method
LECTURE 7
Applications of Differential Equations: Linear Models
This section introduces linear models for solving first-order differential equations, focusing on real-world
applications like population growth, radioactive decay, Newton’s Law of Cooling, and mixing problems
��
= �(� − �� )
��
• �� : Ambient temperature
Class Work and Practice problems set 3
1. When a cake is removed from an oven, its temperature is measured at 3000F. Three minutes later its
temperature is 2000F. How long will it take for the cake to cool off to a room tem perature of 700F. ?
2. A thermometer is removed from a room where the temperature is 700F and is taken outside, where the air
temperature is 10F. After one-half minute the thermometer reads 500F. What is the reading of the
thermometer at t 1 min? How long will it take for the thermometer to reach 150F?
3. A thermometer is taken from an inside room to the outside, where the air temperature is 50F.. After 1
minute the thermometer reads 550F, and after 5 minutes it reads 300F. What is the initial temperature of
the inside room?
4. A small metal bar, whose initial temperature was 200C, is dropped into a large container of boiling water.
How long will it take the bar to reach 900C if it is known that its temperature increased 2 in 1 second?
How long will it take the bar to reach 980C?
5. Two large containers A and B of the same size are filled with different fluids. The fluids in containers A
and B are maintained at 00C and 1000C, respectively. A small metal bar, whose ini tial temperature is
1000C, is lowered into container A. After 1 minute the temperature of the bar is 900C. After 2 minutes the
bar is removed and instantly transferred to the other container. After 1 minute in container B the
temperature of the bar rises 10. How long, measured from the start of the entire process, will it take the
bar to reach 99.90C?
Mixture of two Solutions
Net rate of change of salt
��
= ����� ���� − ������ ����
��
1. A tank contains 200 liters of fluid in which 30 grams of salt is dissolved. Brine containing 1 gram of salt per
liter is then pumped into the tank at a rate of 4 L /min; the well-mixed solution is pumped out at the same
rate. Find the number A(t) of grams of salt in the tank at time t.
2. Solve Problem 1 assuming that pure water is pumped into the tank.
3. A large tank is filled to capacity with 500 gallons of pure water. Brine containing 2 pounds of salt per gallon
is pumped into the tank at a rate of 5 gal/min. The well-mixed solution is pumped out at the same rate. Find
the number A(t) of pounds of salt in the tank at time t
Series Circuits:Resistor-Inductor (R-L)
The sum of voltage drops in a series circuit equals the impressed voltage
��
� + �� = �(�)
��
LR-series circuit
1
�� + � = �(�)
�
�� 1
� + � = �(�)
�� �
��
�=
R: Resistance (Ohms), C: Capacitance (Farads)
�� RC-series circuit
Class Work and Practice problems set 3
A 100-volt electromotive force is applied to an RC-series circuit in which the resistance is 200 ohms and the
capacitance is 10–4 farad. Find the charge q(t) on the capacitor if q(0) =0. Find the current i(t).
A 200-volt electromotive force is applied to an RC-series cir cuit in which the resistance is 1000 ohms and
the capacitance is 5 × 10−6 farad. Find the charge q(t) on the capacitor if i(0) =0.4. Determine the charge and
current at � = 0.005� . Determine the charge as t → ∞
An electromotive force
120, 0 ≤ � ≥ 20
�(�) =
0, � < 20
is applied to an LR-series circuit in which the inductance is 20 henries and the resistance is 2 ohms. Find the
current i(t) if i(0) = 0.
An LR-series circuit has a variable inductor with the inductance defined by
1 0 ≤ � ≥ 10
�(�) = 1 − 10 �,
� < 10
0,
Find the current i(t) if the resistance is 0.2 ohm, the impressed voltage is E(t) =4, and i(0)=0
Applications of Differential Equations: Nonlinear
Models
Population Dynamics and Logistic Growth
The number N(t) of people in a community who are exposed to a particular advertisement is governed by the
logistic equation. Initially N(0)=500, and it is observed that N(1)=1000. Solve for N(t) if it is predicted that the
limiting number of people in the community who will see the advertisement is 50,000.
LECTURE 7
Systems of Linear Differential Equations
�
Here ��� (�) are polynomials in the differential operator � =
��
Restricting the study to first-order systems, the equations take the form:
��1
= �1 (�, �1 , �2 , . . . , �� )
��
��2
= �2 (�, �1 , �2 , . . . , �� )
��
: :
. .
���
= �� (�, �1 , �2 , . . . , �� )
��
Theory of Linear Systems
Definition
A linear system arises when each �� (�, �1 , �2 , . . . , �� ) is linear in the dependent variables �1 , �2 , . . . , �� . This
results in:
��1
= �11 (�)�1 + �12 (�)�2 + . . . + �1� (�)�� + �1 (�)
��
��2
= �11 (�)�1 + �22 (�)�2 + . . . + �2� (�)�� + �2 (�)
��
: :
. .
���
= ��1 (�)�1 + ��2 (�)�2 + . . . + ��� (�)�� + �� (�)
��
Homogeneous: If �� (�) = 0 for all � = 1,2, . . . , �.
Non-Homogeneous: If �� ≠ 0 for any �
Cont. Theory of Linear Systems
Matrix Representation of Linear Systems
�=
�� (�)
, �(t) = �21 (�) �22 (�) … �2� (�) , �(�) =
�� (�)
⋮ ⋮
⋮ ⋮ ⋮
�� (�) �31 (�) �32 (�) �3� (�) �� (�)
Homogeneous Case
If F(t)=0, the system reduces to:
�′ = ��
Theory of Linear Systems
Systems of n linear differential equations in n unknowns can be written in a normal form as first-order
differential equations. The equations are linear when the functions involved are linear in the dependent
variables.
Linear Systems and Notation
A first-order linear system is written as:
��
= �(�)� + �(�)
��
X′ = �� + �
Where:
• � is the column vector of unknown functions,
• �(�) is a matrix of coefficients,
• �(�) is a vector of nonhomogeneous terms.
For homogeneous systems, �(�) = � and the system simplifies to:
X′ = ��
Class Work and Practice problems set 3
� �
If � = ,or � then write the matrix form of the homogeneous linear system:
�
�
�� ��
�� �� a. = 4� − 7�, = 5�
�� ��
1)
��
= 3� + 4�,
��
= 5� − 7� A) = 3� − 5�, = 4� + 8�
�� ��
�� ��
��
�� b. =�−�
��
2)
��
= 6� + � + � + � B). ��
=− 3� + 4� − 9� ��
�� �� = � + 2�
��
��
= 8� + 7� − � + 10� = 6� − � ��
��
�� =− � + �
�� ��
= 10� + 4� + 3� ��
= 2� + 9� − � + 6� ��
��
��
�� c. =− 3� + 4� + �−� sin 2�
��
1)
��
= 6� + � + � + � C). =�−�+�+�−1 ��
��
��
�� = 5� + 9� + 4�−� cos 2�
��
��
= 8� + 7� − � + 10� = 2� + � − � − 3�2 ��
��
��
�� = � + 6� − �−�
�� = � + � + � + �2 − � + 2 ��
= 2� + 9� − � + 6� ��
��
Class Work and Practice problems set 3
Verify that the vector X is a solution of the given homogeneous linear system
��
1. ��
= 3� − 4�
2 1 1 4
��
= 4� − 7�; � =
1
�−3� 3. �′ = �; � = �� + ���
�� 2 −1 0 3 4
�� 1 2 1 1
2. =− 2� + 5� 4. �′ = 6 −1 0 �; � = 6
��
�� 5cos � −13
=− 2� + 4�; �= �� −1 −2 −1
�� 3cos �−sin �
����
1 1 0 1 1 1
3. �′ = −1 4 �; � =
−1 − sin �− cos �
2
�−3� 2
5. �′ = 1 1 0 �; � = 2 2
−����+����
1 −1 −2 0 −1
4.
Homogeneous Linear Systems
A homogeneous linear system is of the form:
�′ = ��
where A is a square matrix of constants, X is a vector of unknowns, and �′ is its derivative with respect to time
t
The goal is to determine solutions for X(t).
�(t) = K���
Characteristic Equation
For a nontrivial solution (� ≠ 0)the determinant of (A−λI) must be zero:
���(� − ��)=0
The resulting equation (polynomial) is the characteristic equation, and its roots are the eigenvalues � of �
Each eigenvalue λ has a corresponding eigenvector K, satisfies:
(A−λI)K=0.
Case 1: Distinct Real Eigenvalues
When A has n distinct real eigenvalues (�1 , �2 , . . . , �� ), the solution is:
Example
��
Solve 1. = 2� + 3�
��
�� ��
= 2� + �; 2. =− 4� + � + �
�� ��
��
= � + 5� − �
��
1 1 4 1 ��
3. �′ = 0 2 0 �; �(�) = 3 = � − 3�
��
0
1 1 1
Class Work and Practice problems set 3
Find the general solution of the given system
��
1. = � + 2�
�� ��
4. =�+�−� 1 0 4
��
6. �′ = 0 1 0 �
�� �� 1 0 1
= 4� + 3� = 2�
�� ��
�
�� � 3
=�−� 7. ′
� = �
�, �(�) =
�� 5 �� � 5
2. ��
=− � + 2�
2
� −
�
��
5. ��
= 2� − 7�
��
=
3
� − 3� −� �
�� 4 8. �′ = �,
�� −� �
= 5� + 10� + 4�
�� ��
3. =− 4� + 2�
��
��
= 5� + 2�
�� 5 ��
=− � + 2�
�� 2
Case 2: Repeated Eigenvalues
Not all eigenvalues (�1 , �2 , . . . , �� ) of an � � matrix A are distinct; some may be repeated.
if m is a positive integer and (� − �� )� is a factor of the characteristic equation while (� − �� )�+1 is not a
factor, then �� is an eigenvalue with multiplicity �
�� = K11 ��1�
�� = �21 ���1� + K22 ��1�
⋮
��−1 � �−2
�� = ��1 ��1� + ��2 ��1� + ⋯ + ��� ��1�
(� − 1)! (� − 2)!
Assume �1 is an eigenvalue of multiplicity 2 and only one eigenvector K can be associated with
it . Then, a second solution of �′ =AX will take the form
�� = ����1� + ���1�
Eigenvalue of Multiplicity Two
Assume �1 is an eigenvalue of multiplicity 2 and only one eigenvector K can be associated with it . Then,
a second solution of �′ =AX will take the form
�� = ����1� + ���1�
Substituting into the system yields
Now, assume �1 is an eigenvalue of multiplicity 3 and only one eigenvector can be associated with it .
Using a similar argument as above, a second solution would take the same form as �� (given above)
and a third solution will take the form
�� � �
�� = � � 1 + ����1� + ���1�
�
Where
(� − �1 �)� = 0
(� − �1 �)� = �
(� − �1 �)� = �
4 1 0
Solve �′ = 0 4 1 X
0 0 4
Class Work and Practice problems set 3
Find the general solution of the given system
��
1. = 3� − �
�� ��
4. = 3� − � − � 1 0 0
��
6. �′ = 2 2 −1 �
�� �� 0 1 0
= 9� − 3� =�+�−�
�� ��
� � −1
�� 7. �′ = �, �(�) =
=�−�+� −� � 6
�� ��
2. ��
=− 6� + 5�
−� �
5.
��
= 3� + 2� + 4� 8. �′ = �,
�� �� −� ��
=− 5� + 4�
��
��
= 2� + 2�
��
3.
��
= 4� + 2� + 3�
��
Case3: Complex Eigenvalues in Systems
For a coefficient matrix A, if the eigenvalues �1 = � + �� and �2 = � − �� (where β ≠ �) are complex, their
corresponding eigenvectors generally have complex entries.
Solution Construction:
A solution for the system can be expressed using eigenvalues and eigenvectors:
� = �1 �1 ��1� + �2 �2 ��2�
where �1 and �2 are eigenvectors corresponding to �1 and �2 , respectively.
Euler's Formula:
Using Euler’s formula (���� = cos (��) + �sin (��)), the exponential terms with complex eigenvalues can be
expressed as:
��1� = ��� (cos (��) + �sin (��))
Real Solutions:
The general solution is rewritten in terms of real functions cos (��) ��� sin (��), ensuring linear independence
and real entries:
�1 = �1 cos (��) + �2 sin (��)
� = �1 �1 + �2 �2
with arbitrary constants �1 ��� �2
Class Work and Practice problems set 3
find the general solution of the given system
� � 2
1. �′ = �, �(0) =
−� −� −1 ��
6. =�
�� ��
� −� 9. = 5� + �
2. �′ = �, ��
� −� �� ��
=− � =− 2� + 3�
�� ��
� −�
3. �′ = �, ��
� −� =� 10.
��
= 4� + 5�
�� ��
� ��
4. �′ = �, �� ��
−� −� 8. ��
= 2� + � + 2� =− 2� + 6�
��
−� � 0
5. �′ = �, �(0) = ��
= 3� + 6�
−� −� 3
��
� � ��
6. �′ = � =− 4� − 3�
−� � ��
LECTURE 8
Nonhomogeneous Linear Systems
vUndetermined Coefficients
vVariation of Parameters
Nonhomogeneous Linear Systems
Overview:
v Methods of undetermined coefficients and variation of parameters apply to nonhomogeneous linear systems.
v Undetermined Coefficients: Quick for constant, polynomial, exponential, sine, cosine functions.
v Variation of Parameters: More powerful, applied when undetermined coefficients aren’t sufficient.
General Solution:
� = �� + ��
Apply only when F(t) consists of constants, polynomials, exponentials, sines, cosines.
Steps:
1. Solve homogeneous system �� :
v Use characteristic equation of matrix A to find eigenvalues and eigenvectors.
Ex
−� � −�
1. �′ = �+
−� � �
� � ��
2. �′ = �+
� � −���+�
Method 2 - Variation of Parameters
Source:Zill (2017)
Method 2 - Variation of Parameters
Source:Zill (2017)
Method 2 - Variation of Parameters
Source:Zill (2017)
Class Work and Practice problems set 3
In Problems 1–6,use variation of parameters and in 7-11 , , use the method of undetermined coefficients
to solve the given system
−� � ����
9. �′ = �+
−� � ��
6.
��
= 2� − � −� � −�����
1. �′ = �+ �−�
� −� ��
��
� −� � = 3� − 2� + 4� � � � �
′ � �
2. � = � �+ � �� ′
10. � = � � � � + −��
−� −�
�
�� � � � ��
7. = 2� + 3� − 7
� −� ����� ��
3. �′ = �+ ��� ��
=− � − 2� + 5
� � ������
�� � � � �
� � �� ��
11. �′ = � � � � + −� ���
4. �′ = �+ � 8. = 5� + 9� + 2 � � � �
� −� �� ��
��
5. ��
= 3� − 3� + 4
��
�� =− � + 11� + 6
= 2� − 2� − 1 ��
��
LECTURE 8
The Laplace Transforms
Exponential Order
A function f(t) f(t) is of exponential order c if constants � > �, � > � , ��� � > �exist such that:
Existence Theorem:
If f (t) is continuous and exponential order with constant c, then
�(s) = L{f(t)}
is defined for all � > �
Class Work and Practice problems set 3
Use the Definition to find {�(�))}
−1, 0≤�<1
1) �(�) =
1, �≥1
4, 0≤�<2
2) �(�) =
0, �≥2
2� + 1, 0≤�<1
3) �(�) =
0, �≥1
−1, 0≤�<1
4) �(�) =
1, �≥1
����, 0≤�<�
5) �(�) =
0, �≥�
����, 0≤�<� 2
6) �(�) =
0, �≥� 2
The Inverse Laplace Transforms
if L{f(t)}=F(s), then the inverse Laplace transform is denoted as:
1 1
1. �{��� } = �−� ↔ �−1 �−�
= ���
1 1
2. �{1} = � ↔ �−1 �
=1
� �! −1 1 ��−1
3. �{� } = ��+1 ↔ � ��
=
(�−1)!
� �
4. �{cos �� } = �2+�2 ↔ �−1 �2 +�2
= cos ��
� 1 1
5. �{sin �� } = �2+�2 ↔ �−1 �2 +�2
= sin ��
�
Linearity of the Inverse Laplace Transform
The transform is linear, meaning:
−2�+6 �+1 �
3 10. �−1 2 15. �−1 �2 +2
20. �−1 (�−2)(�−3)(�−6)
−1 (�+1) � +4
5. �
�4
Transforms of Derivatives
Laplace transform converts derivatives into algebraic terms.
The first translation theorem allows for easy evaluation of Laplace transforms for exponential multiples of
functions.
�� ∞ −�� ��
where �{� �(�)} = 0
� � �(�) ��
∞ −(�−�)�
= 0
� �(�) �� == �(� − �)
Inverse Transformation:
To find �−1 {�(� − �)} identify F(s), take �−1 to find f(t), and multiply by ���
Class Work and Practice problems set 3
find either F (s) or f (t), as indicated
� 5�
14. �−1 15. �−1
{��−10� } 8. �−1
1 (�+1)2 (�−2)2
(�+2)3
2�−1
{�5� �3 } 1 16. �−1
9. �−1 �2 (�+1)3
(�−4)4
{�−2� cos 4� }
1
17. �′′ − 6�′ + 9� = �2 �3� , �(0) = 2, �′ (0) = 17
10. �−1
{�(�� + �2� )2 } �2 −6�+10
18. �′′ + 4�′ + 6� = 1 + �−� , �(0) = 0, �′ (0) = 0
1
{�2� (� − 1)2 } 11. �−1 �2 +2�+5 19. �′′ + 2�′ + � = 0, �(0) = 1, �′ (0) = 0
{�� sin 3� } 12. �−1
�
�2 +4�+5 20. �′ − � = 1 + ��� , �(0) = 0,
�
{�3� (9 − 4� + 10sin )}
2
13. �−1
2�+5 21. �′′ − 2�′ + 5� = 1 + �, �(0) = 0, �′ (0) = 4
�2 +6�+34
Unit Step Function
Unit Step Function
Translation on the t-Axis (Unit Step Function)
The unit step function �(� − �) is defined to be
0, 0≤�<�
�(� − �) = �(�) =
1, �≥�
Find the Laplace transform of the function f whose graph is given in Figure
In problem 1-12,ind either F(s) or f (t), In Problems 13–16, use the Laplace transform to solve the
as indicated given initial-value problem.
−��
−1 �
{cos � �(� − �)} 8. �
�2 +1
13. �′ + � = �(�), �(0) = 5, where
0, 0≤�<�
�(�) =
�−2� 3cos � , �≥�
{(� − 1)�(� − 1)} 9. � −1
�2 (�−1)
14. �′ + � = �(�), �(0) = 0 �ℎ��
{cos � �(� − �)} −�� 2 0, 0≤�<1
−1 �� �(�) =
10. �
�2 +4 5, �≥1
{cos 2� �(� − �)}
11. � −1 (1+�
−2� )2
15. �′ + � = �(�), �(0) = 0 �ℎ��
{(3� + 1)�(� − 1)} �+2 1, 0≤�<1
�(�) =
−2�
−1, �≥1
2−� −1 �
{� �(� − 1)} 12. �
�−4
16. �′′ + 4� = �(�), �(0) = 0, �′(0) =− 1, �ℎ��
−1 �
−2� 1, 0≤�<1
7. � �(�) =
�3 0, �≥1
1. The differential equation for the instantaneous
charge q(t) on the capacitor in an LRC-series circuit is write each function in terms of unit step functions.
Find the Laplace transform of the given function
2, 0≤�<3
a) �(�) =
−2, �≥3
�2� �� 1 1, 0≤�<4
� 2 + � + � = �(�). b) �(�) = 0, 4≤�<5
�� �� �
1, �≥5
Use the Laplace transform to find q(t) when � =
1ℎ, � = 20�ℎ��, � = 0.005�, �(�) = 150�, � > 0, 0≤�<1 Graph for Problem g)
c) �(�) =
0, �(0) = 0 ��� �(0) = 0. What is the current i(t)? �2 , �≥1
0, 0 ≤ � < 3� 2
2. Use the Laplace transform to find the charge q(t) d) �(�) =
����, � ≥ 3� 2
in an RC-series when q(0) = 0 and E(t) = �0 �−��
� > 0. Consider two cases: k ≠ 1/RC and k =1/RC.
Graph for Problem h)
In Problems 49–54, match the given graph with
one of the given functions in (a)–(f ). The
graph of f (t) is given in FIGURE 4.3.11.
a) �(�) − �(�)�(� − �)
b) �(� − �)�(� − �)
c) �(�)�(� − �)
d) �(�) − �(�)�(� − �)
e) �(�)�(� − �) − �(�)�(� − �)
f) �(� − �)�(� − �) − �(� − �)�(� − �)
Derivatives of Transforms
Multiplying a Function by ��
The Laplace transform of �� �(�) can be computed using derivatives of the Laplace transform of �(�)
If �(�) = �{�(�)} and if we assume that interchanging of differentiation and integration is possible, then
∞ ∞ ∞
� � � −��
�(�) = �−�� �(�)�� = [� �(�)�� =− �−�� ��(�)�� =− �{��(�)};
�� �� ��
0 0 0
�
that is �{��(�)} =− �{�(�)}
��
similarly
2
�
�{� �(�)} = �{� ∙ ��(�)} =− �{��(�)}
��
� 2
=− 2 �{�(�)}
��
�
• General result for �{� �(�)}, If �(�) = �{�(�)} ��� � = 1,2,3. . . , then
� �
�
�
�{� �(�)} = (−1) � �(�)
��
EXAMPLES
1. Evaluate �{������}. 2. Solve . �′′ + 16� = ���4�, �(0) = 0, �′ (0) = 1
Transform of an Integral
Convolution
The convolution of two functions f(t) and g(t), denoted (� ∗ �)(�), is defined as
�
(� ∗ �)(�) = 0
�(�)�(� − �) ��
The variable � (Greek letter “tau”) is a dummy variable used for integration.
The result of the convolution is a new function of t, emphasizing its role in the definition.
Key Properties of Convolution
Commutative Property:
Linearity: The convolution of sums or scalar multiples follows the rules of linearity
A function f(t) is periodic if �(� + �) = �(�) for all , where T is the period.
Fig 1: Graph for P 1 Fig 2: Graph for P 2 Fig 3: Graph for P 3 Fig 4: Graph for P 4 Fig 5: Graph for P 5
the Dirac Delta Function
The Dirac Delta Function, δ(t−t ), is a mathematical construct introduced to model forces of large magnitude
acting over very short periods, such as a hammer strike or lightning on a wing.