SlideShare a Scribd company logo
EDUCATION HOLE PRESENTS
Engineering Mathematics – II
Unit-I
Differential Equations.............................................................................................................. 2
Linear differential equations of nth order with constant coefficients..................................................................2
Complementary function and Particular integral..............................................................................4
To find the Complementary Function...................................................................................................................4
Particular Integral.............................................................................................................................................4
Simultaneous linear differential equations.......................................................................................5
Using the D Operator............................................................................................................................................5
Using Laplace Transform.......................................................................................................................................6
Solution of second order differential equations by changing dependent & independent variables....................7
Normal form of differential equation ...............................................................................................8
Types of Differential Equations.............................................................................................................................9
Methods................................................................................................................................................................9
Ordinary Differential Equations........................................................................................................................9
Linear Differential Equation Solutions............................................................................................................10
First Order Linear Equations...........................................................................................................................10
Method of variation of parameters....................................................................................................................10
Applications to engineering problems ............................................................................................ 13
Differential Equations
Linear differential equations of nth order with constant coefficients
The linear homogeneous differential equation of the n-th order with constant coefficients can be
written as
where a1, a2,..., an constants which may be real or complex.
Using the linear differential operator L(D), this equation can be represented as
where
For each differential operator with constant coefficients, we can introduce the characteristic
polynomial
The algebraic equation
is called the characteristic equation of the differential equation.
According to the fundamental theorem of algebra, a polynomial of degree n has exactly n roots,
counting multiplicity. In this case the roots can be both real and complex (even if all the
coefficients of a1, a2,..., an are real).
Let us consider in more detail the different cases of the roots of the characteristic equation and
the corresponding formulas for the general solution of differential equations.
Case 1. All Roots of the Characteristic Equation are Real and Distinct
Assume that the characteristic equation L(λ) = 0 has n roots λ1, λ2,..., λn. In this case the general
solution of the differential equation is written in a simple form:
where C1, C2,..., Cn are constants depending on initial conditions.
Case 2. The Roots of the Characteristic Equation are Real and Multiple
Let the characteristic equation L(λ) = 0 of degree n have m roots λ1, λ2,..., λm, the multiplicity of
which, respectively, is equal to k1, k2,..., km. It is clear that the following condition is satisfied:
Then the general solution of the homogeneous differential equations with constant coefficients
has the form
It is seen that the formula of the general solution has exactly ki terms corresponding to each root
λi of multiplicity ki. These terms are formed by multiplying x to a certain degree by the
exponential function exp(λi x). The degree of x varies in the range from 0 to ki − 1, where ki is
the multiplicity of the root λi.
Case 3. The Roots of the Characteristic Equation are Complex and Distinct
If the coefficients of the differential equation are real numbers, the complex roots of the
characteristic equation will be presented in the form of conjugate pairs of complex numbers:
In this case the general solution is written as
Case 4. The Roots of the Characteristic Equation are Complex and Multiple
Here, each pair of complex conjugate roots α ± iβ of multiplicity k produces 2k particular
solutions:
Then the part of the general solution of the differential equation corresponding to a given pair of
complex conjugate roots is constructed as follows:
Complementary function and Particular integral
Standard Form
where are constants.
The general solution of such a differential equation contains two parts; the Complementary
Function and the Particular Integral.
To find the Complementary Function
Write the auxiliary equation in the form: where the differential equation
is in the standard form.
Solving the auxiliary equation, a quadratic, will yield two roots, say m1 and m2.
The Complementary Function (CF) is then written based on the roots.
• If m1 and m2 are real and different, then the CF is
• If m1 and m2 are real and equal, ie m1 = m2 = m, then the CF is
• If m1 and m2 are complex roots of the form and , then the CF is
where and are arbitrary constants.
Particular Integral
The Particular Integral is similar to the function of x, on the RHS of the differential
equation and it is determined by the method of undetermined coefficients.
If the RHS of the equation is a constant, then we use a trial PI of K (a different constant). A trial
PI of a general first degree form and a general second degree form
is used for functions of and respectively. If the function is of the form ,
then we use , the same exponential function. Similarly, trigonometric functions will have
trial PIs of the form .
are undetermined coefficients. By finding and and substituting and
into the given differential equation, the values of the undetermined coefficients can be
obtained. Hence the Particular Integral is found.
Combine the Complementary Function and the Particular Integral and we will get the General
Solution of the Second Order Differential Equation with constant coefficients
Simultaneous linear differential equations
A simultaneous differential equation is one of the mathematical equations for an indefinite
function of one or more than one variables that relate the values of the function. Differentiation
equation in various orders, Differential equations plays an important function in engineering,
physics, economics, and other disciplines. This analysis concentrates on linear equations with
Constant Coefficients.
Using the D Operator
The D operator is a linear operator applied to functions and which is defined as .
For example
Example:
Example - First example
Problem
Workings
The equations may be written as:
(1)
(2)
Eliminating y from equations (1) and (2)
i.e.
Solution
From equation (1)
Therefore
Using Laplace Transform
Laplace transform or the Laplace operator is a linear operator applied to functions and which is
defined as where
Problem
Apply the Laplace Transform and find and :
Give that at
Workings
The equations may be written as:
Hence
Therefore
(3)
(4)
Then we multiply the equation (3) with :
(5)
We subtract equation (4) from equation (5) :
Solution
Therefore
From equation (4) we get :
Solution of second order differential equations by changing dependent &
independent variables
We suppose that we have a second order differential equation with dependent variable u and independent
variable t. The harmonic oscillator or pendulums are good examples. We suppose, further, that we are
given initial values of u and u', as well as a formula for u" in terms of u, u' and t
u" = f(u, u', t)
If f does not involve u or u' we could integrate both sides of the equation once to find u', and again to find
u, with a linear function ct + d that must be determined from the initial values
Thus the problem we face is not unlike that of performing a double integral. However, it is surprisingly
easy to do, on a spreadsheet.
We will use an approximation technique that has the order of accuracy of the trapezoid rule;
improvements can be obtained by extrapolating as with most numerical methods.
We begin by describing the basic approach.
We will use a column to represent each the variables t, u, u' and u". The first three of these will start at the
given initial values of these variables; and the initial value of u" can be computed from them.
In each successive row, t will increase by a constant amount, d. (You can later choose to make d smaller
than your initial choice when variables change too much in one d interval.)
We will set
(Which are generic statements, independent of the equation itself.)
Finally we provide an expression for u"(t + d). In doing so we cannot use u(t + d) or u'(t + d) or else our
definitions would be circular. On the other hand we would like to use something which averages the
second derivative at the ends of the time interval between t and t + d (at least to some order), just as the is
done above for u and u'. For the value of u" at t we use f(t, u(t), u'(t)). For the value of u" at t + d we use
u"(t + d) = f(t + d, u(t) + d * u'(t), u'(t) + d * f(t, u(t), u'(t)))
Thus we use the linear approximation to u and u' in f(t + d, u, u') defined at t and evaluated at t + d. We
could do slightly better by using a quadratic approximation
Which would be correct to second order, but the trapezoid rule is already in error in second order, so this
will not do much good, usually.
You can observe the behavior of u and u' as a function of t by plotting the first three columns (using the
"charting" capability of the spreadsheet program with xy scatter plots). You can observe the "phase plane"
behavior of the solution by plotting the second and third columns, namely u' vs. u.
If you do things right you can change d or parameters in f in one keystroke, and watch how the solutions
change as you change them.
You can handle equations with several dependent variables, like those of planetary motion, by the same
approach; with motion in the xy plane you can have a column for t, x, y, x', y', x" and y", and can observe
trajectories and look at behavior in the xy plane as you change parameters.
Normal form of differential equation
A differential equation is an equation that relates a function to one or more of its derivatives. Differential
equations are especially applicable when the tools of algebra, which are ideally suited for static systems,
are not enough. Many physical systems are modeled by solving differential equations, although their
usefulness extends well into other fields of science such as chemistry and economics.
Types of Differential Equations
There are two main types of differential equations: Ordinary Differential Equations and Partial
Differential Equations.[1]
The former is simpler of the two, as it can be written in the normal form
for a simple function y = g(x)
The function F consists of the function y and its derivatives up to the nth order. Notice that y is comprised
of only one independent variable x. A differential equation is considered ordinary if the function y in F is
dependent on only one variable. It is important to note that, while most ordinary differential equations can
be written in the normal form (isolating the highest derivative on one side of the equation and moving all
other variables to the other), there are equations in which this cannot be done.
If y were a function of multiple variables, for example
y = F(u,v)
then the derivatives of y in the ordinary equation may be partial derivatives with respect to either u or v. In that
case, any differential equation that has partial derivatives is called a partial differential equation. For
example, the 1-dimensional wave equation[2]
:
Obviously, partial differential equations are much more complicated to solve.
Methods
There are many ways to find solutions to differential equations.
Ordinary Differential Equations
The simplest differential equations to solve are separable differential equations. A differential equation is
separable if it can be written in the form
Then we can separate the two variables, collect the x's on one side and the y's on the other side, then
integrate to get the (n-1) derivative, and integrating again to get the (n-2) derivative, until we have found
the function y. For example, for the derivative n = 1:
The solution is then given implicitly by the expression:
where C is an arbitrary constant.
Linear Differential Equation Solutions
If a differential equation can be written in the normal form
It is considered a linear differential equation.
First Order Linear Equations
The order of a differential equation is equal to the degree of the highest derivative in the equation. For
example, the above equations are order n equations. A first order linear equation appears in the form:
Method of variation of parameters
This method has no prior conditions to be satisfied. Therefore, it may sound more general than
the previous method. We will see that this method depends on integration while the previous one
is purely algebraic which, for some at least, is an advantage.
Consider the equation
In order to use the method of variation of parameters we need to know that is a set of
fundamental solutions of the associated homogeneous equation y'' + p(x)y' + q(x)y = 0. We know
that, in this case, the general solution of the associated homogeneous equation is
. The idea behind the method of variation of parameters is to look for a
particular solution such as
where and are functions. From this, the method got its name.
The functions and are solutions to the system
,
which implies
,
where is the wronskian of and . Therefore, we have
Summary:Let us summarize the steps to follow in applying this method:
(1)
Find a set of fundamental solutions of the associated homogeneous equation
y'' + p(x)y' + q(x)y = 0.
;
(2)
Write down the form of the particular solution
;
(3)
Write down the system
;
(4)
Solve it. That is, find and ;
(5)
Plug and into the equation giving the particular solution.
Example: Find the particular solution to
Solution: Let us follow the steps:
(1)
A set of fundamental solutions of the equation y'' + y = 0 is ;
(2)
The particular solution is given as
(3)
We have the system
;
(4)
We solve for and , and get
Using techniques of integration, we get
;
(5)
The particular solution is:
,
or
Remark: Note that since the equation is linear, we may still split if necessary. For example, we
may split the equation
,
into the two equations
then, find the particular solutions for (1) and for (2), to generate a particular solution for the
original equation by
Applications to engineering problems
• Uniquely covers mathematical modeling, the numerical solution, ordinary and partial
differential equations, and both finite difference and finite element methods
• Introduces readers to scientific computing, which is the integrated science of solving
problems with mathematical, numerical, and programming tools
• Originally written for the course 'Numerical Solution of Differential Equations' offered
at KTH, Stockholm and has been classroom tested over several years
• Features a 'Five M' approach: Modeling, Mathematics, Methods, MATLAB®, and Multi
physics.
• Supplemented with a related Web Site that features solutions to problems, demos,
MATLAB programs, etc.

More Related Content

PDF
Numerical Solution of Nth - Order Fuzzy Initial Value Problems by Fourth Orde...
PDF
Nonlinear perturbed difference equations
PDF
Polya recurrence
PPT
Partial
PDF
International Journal of Mathematics and Statistics Invention (IJMSI)
PDF
umerical algorithm for solving second order nonlinear fuzzy initial value pro...
PDF
NPDE-TCA
PDF
Derivation and Application of Multistep Methods to a Class of First-order Ord...
Numerical Solution of Nth - Order Fuzzy Initial Value Problems by Fourth Orde...
Nonlinear perturbed difference equations
Polya recurrence
Partial
International Journal of Mathematics and Statistics Invention (IJMSI)
umerical algorithm for solving second order nonlinear fuzzy initial value pro...
NPDE-TCA
Derivation and Application of Multistep Methods to a Class of First-order Ord...

What's hot (19)

PDF
Engineering Mathematics-IV_B.Tech_Semester-IV_Unit-I
PDF
Conformal field theories and three point functions
DOCX
COMPLEX NUMBER
PDF
progfree++
PDF
Formal expansion method for solving an electrical circuit model
PDF
Solution of second kind volterra integro equations using linear
PPT
Presentation03 ss
PDF
A semi analytic method for solving two-dimensional fractional dispersion equa...
PDF
Fuzzy Group Ideals and Rings
PDF
Linear regression
PDF
11.a family of implicit higher order methods for the numerical integration of...
DOCX
Indefinite Integral
PDF
Coueete project
PDF
A Two Grid Discretization Method For Decoupling Time-Harmonic Maxwell’s Equat...
PDF
Complex Numbers
PDF
ON COMPUTATIONS OF THPD MATRICES
PPT
Constant strain triangular
PDF
On the construction and comparison of an explicit iterative
PDF
Numerical solution of linear volterra fredholm integro-
Engineering Mathematics-IV_B.Tech_Semester-IV_Unit-I
Conformal field theories and three point functions
COMPLEX NUMBER
progfree++
Formal expansion method for solving an electrical circuit model
Solution of second kind volterra integro equations using linear
Presentation03 ss
A semi analytic method for solving two-dimensional fractional dispersion equa...
Fuzzy Group Ideals and Rings
Linear regression
11.a family of implicit higher order methods for the numerical integration of...
Indefinite Integral
Coueete project
A Two Grid Discretization Method For Decoupling Time-Harmonic Maxwell’s Equat...
Complex Numbers
ON COMPUTATIONS OF THPD MATRICES
Constant strain triangular
On the construction and comparison of an explicit iterative
Numerical solution of linear volterra fredholm integro-
Ad

Viewers also liked (7)

PDF
N.P. BALI ENGINEERING MATH REAL
DOC
Chapter 2 (maths 3)
DOC
Chapter 4 (maths 3)
DOC
Chapter 3 (maths 3)
DOC
Chapter 1 (maths 3)
DOC
Chapter 5 (maths 3)
PPT
Free Download Powerpoint Slides
N.P. BALI ENGINEERING MATH REAL
Chapter 2 (maths 3)
Chapter 4 (maths 3)
Chapter 3 (maths 3)
Chapter 1 (maths 3)
Chapter 5 (maths 3)
Free Download Powerpoint Slides
Ad

Similar to Free Ebooks Download ! Edhole (20)

PPTX
formulation of first order linear and nonlinear 2nd order differential equation
PPTX
Solution to second order pde
PDF
Higher Differential Equation
PPTX
Unit_3_Differential_Equations_PPT a final .pptx
PDF
Maths Notes - Differential Equations
PPTX
Network Analysis Introduction class.pptx
PDF
Higher Order Differential Equation
PPTX
3.1 Ordinary Differential equation Higher Order
PPT
Ordinary differential equation Presentation
PDF
3c3 nonhomgen lineqns-stuhuhuhuhuhh
PPTX
Unit_3_Lecture_Notes_PPT.pptxon differential Equation
PPTX
3.1 Ordinary Differential equation with Higher Order
PDF
Diffy Q Paper
PDF
Solving Linear Differential Equations with Constant Coefficients
PDF
UNIT-III.pdf
DOCX
DIFFERENTIAL EQUATION
PPTX
M - II (R23) UNIT - 2 - Linear D.E's of Higher Order-1.pptx
PPTX
Higher order differential equation
PPTX
odes1.pptx
PDF
formulation of first order linear and nonlinear 2nd order differential equation
Solution to second order pde
Higher Differential Equation
Unit_3_Differential_Equations_PPT a final .pptx
Maths Notes - Differential Equations
Network Analysis Introduction class.pptx
Higher Order Differential Equation
3.1 Ordinary Differential equation Higher Order
Ordinary differential equation Presentation
3c3 nonhomgen lineqns-stuhuhuhuhuhh
Unit_3_Lecture_Notes_PPT.pptxon differential Equation
3.1 Ordinary Differential equation with Higher Order
Diffy Q Paper
Solving Linear Differential Equations with Constant Coefficients
UNIT-III.pdf
DIFFERENTIAL EQUATION
M - II (R23) UNIT - 2 - Linear D.E's of Higher Order-1.pptx
Higher order differential equation
odes1.pptx

More from Edhole.com (20)

PPT
Ca in patna
PPT
Chartered accountant in dwarka
PPT
Ca in dwarka
PPT
Ca firm in dwarka
PPT
Website development company surat
PPTX
Website designing company in surat
PPTX
Website dsigning company in india
PPT
Website designing company in delhi
PPT
Ca in patna
PPT
Chartered accountant in dwarka
PPT
Ca firm in dwarka
PPTX
Ca in dwarka
PPTX
Website development company surat
PPT
Website designing company in surat
PPT
Website designing company in india
PPT
Website designing company in delhi
PPT
Website designing company in mumbai
PPT
Website development company surat
PPT
Website desinging company in surat
PPT
Website designing company in india
Ca in patna
Chartered accountant in dwarka
Ca in dwarka
Ca firm in dwarka
Website development company surat
Website designing company in surat
Website dsigning company in india
Website designing company in delhi
Ca in patna
Chartered accountant in dwarka
Ca firm in dwarka
Ca in dwarka
Website development company surat
Website designing company in surat
Website designing company in india
Website designing company in delhi
Website designing company in mumbai
Website development company surat
Website desinging company in surat
Website designing company in india

Recently uploaded (20)

PDF
Supply Chain Operations Speaking Notes -ICLT Program
PDF
VCE English Exam - Section C Student Revision Booklet
PDF
Abdominal Access Techniques with Prof. Dr. R K Mishra
PDF
Insiders guide to clinical Medicine.pdf
PPTX
Introduction to Child Health Nursing – Unit I | Child Health Nursing I | B.Sc...
PPTX
Week 4 Term 3 Study Techniques revisited.pptx
PDF
Chapter 2 Heredity, Prenatal Development, and Birth.pdf
PPTX
Renaissance Architecture: A Journey from Faith to Humanism
PDF
Business Ethics Teaching Materials for college
PDF
Microbial disease of the cardiovascular and lymphatic systems
PPTX
human mycosis Human fungal infections are called human mycosis..pptx
PDF
Classroom Observation Tools for Teachers
PPTX
Introduction_to_Human_Anatomy_and_Physiology_for_B.Pharm.pptx
PDF
ANTIBIOTICS.pptx.pdf………………… xxxxxxxxxxxxx
PDF
Saundersa Comprehensive Review for the NCLEX-RN Examination.pdf
PPTX
Pharmacology of Heart Failure /Pharmacotherapy of CHF
PDF
Physiotherapy_for_Respiratory_and_Cardiac_Problems WEBBER.pdf
PDF
O7-L3 Supply Chain Operations - ICLT Program
PPTX
Pharma ospi slides which help in ospi learning
PDF
RMMM.pdf make it easy to upload and study
Supply Chain Operations Speaking Notes -ICLT Program
VCE English Exam - Section C Student Revision Booklet
Abdominal Access Techniques with Prof. Dr. R K Mishra
Insiders guide to clinical Medicine.pdf
Introduction to Child Health Nursing – Unit I | Child Health Nursing I | B.Sc...
Week 4 Term 3 Study Techniques revisited.pptx
Chapter 2 Heredity, Prenatal Development, and Birth.pdf
Renaissance Architecture: A Journey from Faith to Humanism
Business Ethics Teaching Materials for college
Microbial disease of the cardiovascular and lymphatic systems
human mycosis Human fungal infections are called human mycosis..pptx
Classroom Observation Tools for Teachers
Introduction_to_Human_Anatomy_and_Physiology_for_B.Pharm.pptx
ANTIBIOTICS.pptx.pdf………………… xxxxxxxxxxxxx
Saundersa Comprehensive Review for the NCLEX-RN Examination.pdf
Pharmacology of Heart Failure /Pharmacotherapy of CHF
Physiotherapy_for_Respiratory_and_Cardiac_Problems WEBBER.pdf
O7-L3 Supply Chain Operations - ICLT Program
Pharma ospi slides which help in ospi learning
RMMM.pdf make it easy to upload and study

Free Ebooks Download ! Edhole

  • 1. EDUCATION HOLE PRESENTS Engineering Mathematics – II Unit-I
  • 2. Differential Equations.............................................................................................................. 2 Linear differential equations of nth order with constant coefficients..................................................................2 Complementary function and Particular integral..............................................................................4 To find the Complementary Function...................................................................................................................4 Particular Integral.............................................................................................................................................4 Simultaneous linear differential equations.......................................................................................5 Using the D Operator............................................................................................................................................5 Using Laplace Transform.......................................................................................................................................6 Solution of second order differential equations by changing dependent & independent variables....................7 Normal form of differential equation ...............................................................................................8 Types of Differential Equations.............................................................................................................................9 Methods................................................................................................................................................................9 Ordinary Differential Equations........................................................................................................................9 Linear Differential Equation Solutions............................................................................................................10 First Order Linear Equations...........................................................................................................................10 Method of variation of parameters....................................................................................................................10 Applications to engineering problems ............................................................................................ 13 Differential Equations Linear differential equations of nth order with constant coefficients The linear homogeneous differential equation of the n-th order with constant coefficients can be written as where a1, a2,..., an constants which may be real or complex. Using the linear differential operator L(D), this equation can be represented as where For each differential operator with constant coefficients, we can introduce the characteristic polynomial
  • 3. The algebraic equation is called the characteristic equation of the differential equation. According to the fundamental theorem of algebra, a polynomial of degree n has exactly n roots, counting multiplicity. In this case the roots can be both real and complex (even if all the coefficients of a1, a2,..., an are real). Let us consider in more detail the different cases of the roots of the characteristic equation and the corresponding formulas for the general solution of differential equations. Case 1. All Roots of the Characteristic Equation are Real and Distinct Assume that the characteristic equation L(λ) = 0 has n roots λ1, λ2,..., λn. In this case the general solution of the differential equation is written in a simple form: where C1, C2,..., Cn are constants depending on initial conditions. Case 2. The Roots of the Characteristic Equation are Real and Multiple Let the characteristic equation L(λ) = 0 of degree n have m roots λ1, λ2,..., λm, the multiplicity of which, respectively, is equal to k1, k2,..., km. It is clear that the following condition is satisfied: Then the general solution of the homogeneous differential equations with constant coefficients has the form It is seen that the formula of the general solution has exactly ki terms corresponding to each root λi of multiplicity ki. These terms are formed by multiplying x to a certain degree by the exponential function exp(λi x). The degree of x varies in the range from 0 to ki − 1, where ki is the multiplicity of the root λi. Case 3. The Roots of the Characteristic Equation are Complex and Distinct If the coefficients of the differential equation are real numbers, the complex roots of the characteristic equation will be presented in the form of conjugate pairs of complex numbers: In this case the general solution is written as
  • 4. Case 4. The Roots of the Characteristic Equation are Complex and Multiple Here, each pair of complex conjugate roots α ± iβ of multiplicity k produces 2k particular solutions: Then the part of the general solution of the differential equation corresponding to a given pair of complex conjugate roots is constructed as follows: Complementary function and Particular integral Standard Form where are constants. The general solution of such a differential equation contains two parts; the Complementary Function and the Particular Integral. To find the Complementary Function Write the auxiliary equation in the form: where the differential equation is in the standard form. Solving the auxiliary equation, a quadratic, will yield two roots, say m1 and m2. The Complementary Function (CF) is then written based on the roots. • If m1 and m2 are real and different, then the CF is • If m1 and m2 are real and equal, ie m1 = m2 = m, then the CF is • If m1 and m2 are complex roots of the form and , then the CF is where and are arbitrary constants. Particular Integral The Particular Integral is similar to the function of x, on the RHS of the differential equation and it is determined by the method of undetermined coefficients.
  • 5. If the RHS of the equation is a constant, then we use a trial PI of K (a different constant). A trial PI of a general first degree form and a general second degree form is used for functions of and respectively. If the function is of the form , then we use , the same exponential function. Similarly, trigonometric functions will have trial PIs of the form . are undetermined coefficients. By finding and and substituting and into the given differential equation, the values of the undetermined coefficients can be obtained. Hence the Particular Integral is found. Combine the Complementary Function and the Particular Integral and we will get the General Solution of the Second Order Differential Equation with constant coefficients Simultaneous linear differential equations A simultaneous differential equation is one of the mathematical equations for an indefinite function of one or more than one variables that relate the values of the function. Differentiation equation in various orders, Differential equations plays an important function in engineering, physics, economics, and other disciplines. This analysis concentrates on linear equations with Constant Coefficients. Using the D Operator The D operator is a linear operator applied to functions and which is defined as . For example Example: Example - First example Problem Workings The equations may be written as:
  • 6. (1) (2) Eliminating y from equations (1) and (2) i.e. Solution From equation (1) Therefore Using Laplace Transform Laplace transform or the Laplace operator is a linear operator applied to functions and which is defined as where Problem Apply the Laplace Transform and find and : Give that at Workings The equations may be written as: Hence Therefore
  • 7. (3) (4) Then we multiply the equation (3) with : (5) We subtract equation (4) from equation (5) : Solution Therefore From equation (4) we get : Solution of second order differential equations by changing dependent & independent variables We suppose that we have a second order differential equation with dependent variable u and independent variable t. The harmonic oscillator or pendulums are good examples. We suppose, further, that we are given initial values of u and u', as well as a formula for u" in terms of u, u' and t u" = f(u, u', t) If f does not involve u or u' we could integrate both sides of the equation once to find u', and again to find u, with a linear function ct + d that must be determined from the initial values Thus the problem we face is not unlike that of performing a double integral. However, it is surprisingly easy to do, on a spreadsheet. We will use an approximation technique that has the order of accuracy of the trapezoid rule; improvements can be obtained by extrapolating as with most numerical methods. We begin by describing the basic approach. We will use a column to represent each the variables t, u, u' and u". The first three of these will start at the given initial values of these variables; and the initial value of u" can be computed from them.
  • 8. In each successive row, t will increase by a constant amount, d. (You can later choose to make d smaller than your initial choice when variables change too much in one d interval.) We will set (Which are generic statements, independent of the equation itself.) Finally we provide an expression for u"(t + d). In doing so we cannot use u(t + d) or u'(t + d) or else our definitions would be circular. On the other hand we would like to use something which averages the second derivative at the ends of the time interval between t and t + d (at least to some order), just as the is done above for u and u'. For the value of u" at t we use f(t, u(t), u'(t)). For the value of u" at t + d we use u"(t + d) = f(t + d, u(t) + d * u'(t), u'(t) + d * f(t, u(t), u'(t))) Thus we use the linear approximation to u and u' in f(t + d, u, u') defined at t and evaluated at t + d. We could do slightly better by using a quadratic approximation Which would be correct to second order, but the trapezoid rule is already in error in second order, so this will not do much good, usually. You can observe the behavior of u and u' as a function of t by plotting the first three columns (using the "charting" capability of the spreadsheet program with xy scatter plots). You can observe the "phase plane" behavior of the solution by plotting the second and third columns, namely u' vs. u. If you do things right you can change d or parameters in f in one keystroke, and watch how the solutions change as you change them. You can handle equations with several dependent variables, like those of planetary motion, by the same approach; with motion in the xy plane you can have a column for t, x, y, x', y', x" and y", and can observe trajectories and look at behavior in the xy plane as you change parameters. Normal form of differential equation A differential equation is an equation that relates a function to one or more of its derivatives. Differential equations are especially applicable when the tools of algebra, which are ideally suited for static systems, are not enough. Many physical systems are modeled by solving differential equations, although their usefulness extends well into other fields of science such as chemistry and economics.
  • 9. Types of Differential Equations There are two main types of differential equations: Ordinary Differential Equations and Partial Differential Equations.[1] The former is simpler of the two, as it can be written in the normal form for a simple function y = g(x) The function F consists of the function y and its derivatives up to the nth order. Notice that y is comprised of only one independent variable x. A differential equation is considered ordinary if the function y in F is dependent on only one variable. It is important to note that, while most ordinary differential equations can be written in the normal form (isolating the highest derivative on one side of the equation and moving all other variables to the other), there are equations in which this cannot be done. If y were a function of multiple variables, for example y = F(u,v) then the derivatives of y in the ordinary equation may be partial derivatives with respect to either u or v. In that case, any differential equation that has partial derivatives is called a partial differential equation. For example, the 1-dimensional wave equation[2] : Obviously, partial differential equations are much more complicated to solve. Methods There are many ways to find solutions to differential equations. Ordinary Differential Equations The simplest differential equations to solve are separable differential equations. A differential equation is separable if it can be written in the form Then we can separate the two variables, collect the x's on one side and the y's on the other side, then integrate to get the (n-1) derivative, and integrating again to get the (n-2) derivative, until we have found the function y. For example, for the derivative n = 1:
  • 10. The solution is then given implicitly by the expression: where C is an arbitrary constant. Linear Differential Equation Solutions If a differential equation can be written in the normal form It is considered a linear differential equation. First Order Linear Equations The order of a differential equation is equal to the degree of the highest derivative in the equation. For example, the above equations are order n equations. A first order linear equation appears in the form: Method of variation of parameters This method has no prior conditions to be satisfied. Therefore, it may sound more general than the previous method. We will see that this method depends on integration while the previous one is purely algebraic which, for some at least, is an advantage. Consider the equation
  • 11. In order to use the method of variation of parameters we need to know that is a set of fundamental solutions of the associated homogeneous equation y'' + p(x)y' + q(x)y = 0. We know that, in this case, the general solution of the associated homogeneous equation is . The idea behind the method of variation of parameters is to look for a particular solution such as where and are functions. From this, the method got its name. The functions and are solutions to the system , which implies , where is the wronskian of and . Therefore, we have Summary:Let us summarize the steps to follow in applying this method: (1) Find a set of fundamental solutions of the associated homogeneous equation y'' + p(x)y' + q(x)y = 0. ; (2) Write down the form of the particular solution ;
  • 12. (3) Write down the system ; (4) Solve it. That is, find and ; (5) Plug and into the equation giving the particular solution. Example: Find the particular solution to Solution: Let us follow the steps: (1) A set of fundamental solutions of the equation y'' + y = 0 is ; (2) The particular solution is given as (3) We have the system ; (4) We solve for and , and get Using techniques of integration, we get
  • 13. ; (5) The particular solution is: , or Remark: Note that since the equation is linear, we may still split if necessary. For example, we may split the equation , into the two equations then, find the particular solutions for (1) and for (2), to generate a particular solution for the original equation by Applications to engineering problems • Uniquely covers mathematical modeling, the numerical solution, ordinary and partial differential equations, and both finite difference and finite element methods • Introduces readers to scientific computing, which is the integrated science of solving problems with mathematical, numerical, and programming tools • Originally written for the course 'Numerical Solution of Differential Equations' offered at KTH, Stockholm and has been classroom tested over several years • Features a 'Five M' approach: Modeling, Mathematics, Methods, MATLAB®, and Multi physics.
  • 14. • Supplemented with a related Web Site that features solutions to problems, demos, MATLAB programs, etc.