Method of Systematic Inspection

TITLE: METHOD OF SYSTEMATIC INSPECTION FOR SOLVING DIFFERENTIAL EQUATIONS

AUTHOR: YOUSIF TAWFIQ NEMER SAMMOUR

---

Copyright©2009, Yousif Tawfiq Sammour. Publication of any part of this document “in any media format” must be approved by the author.

---

ABSTRACT

The method of systematic inspection solves or helps in discovering the behavior of differential equations. This method solves differential equations by creating functions of the independent variable(s) from two opposite small functions called seeds. Starting with a seed you can build a bigger function. Adding the two opposite functions together will give a function that satisfies the differential equation.

2000 mathematics subject classification. 34A30, 35G15, 34A34,35E99

1. Introduction

no previous research was carried out regarding systematic inspection. As a matter of fact inspecting a solution for a differential equation or for a prticular integral was a matter of trial and error. In this paper I tried to make inspection a systematic method. It is applied to linear and non-linear differential equations. Using this method doesn’t require much knowledge of differential equations. I adopted examples more than theory in this paper to clear up the concept. I solved simple examples although the method applies to complex ones.

2. How to inspect the solution

The procedure followed to inspect functions that satisfy an ordinary differential equation is the same as that for inspecting functions for a partial differential equation. Looking at the following example you will understand how to inspect a solution for an ordinary differential equation.

2.1 Example

Inspect the complementary function for the following differential equation :

steps

1-put the two terms in two rows

2- since then put opposite to andopposite tolike this

this is valid for linear differential equations only.

is called positive seed and,

is called negative seed.

They are called seeds because using them you can build two larger opposite functions as you will see soon

3- to simplify work , let us break step (2) into two cases and work them out separately.

First case

Second case

of both cases

4- now we start the process of what I call calculate and balance for the first case.

A – calculate

by differentiating once

B – balance 1

move it with an opposite sign

C – calculate

by differentiating once

D – balance 2

move with an opposite sign

continuing the process of calculate and balance we get a series like this,

… (1)

now we process the second case

A – calculate

by integrating once

B – balance 1

move with an opposite sign

C – calculate

by integrating once

D – balance 2

move with an opposite sign

repeat the process of calculate and balance till you get a satisfactory number of terms for as you see below,

… (2)

so the result is :

…

…

this is the inspected function form from which we can get a function that satisfies the given differential equation , and it is convergent for .

By substituting different values of n in the above form, we note that only one function satisfies the differential equation and the others will be similar to it as we’ll see soon, but if we have a second order differential equation then we will get two functions and the others will be similar to any of them, and three for third order and so on. (what I mean by similar is that the new function will be a previous inspected function multiplied by a constant). Now if we substitute into the inspected function form above we get,

… (this is Maclaurin’s series for .)

all other inspected functions that

we can get for will be similar to the above inspected function. For example if you get the inspected function for it will be:

…

which is equal to the inspected function for multiplied by

if you try other values of you will get similar functions.

2.2 Finding the particular integral

Suppose that our differential equation looks like this,

then the particular integral can be obtained by putting opposite to as follows,

where ca means calculate and b means balance .

in the case of every calculation we differentiate once with respect to and with balance we move with an opposite sign.

so particular integral is

we can get another particular integral by putting opposite to , but since it is divergent so it is neglected,

in the case of every calculation we integrate once with respect to .

then the general solution for the differential equation :

2.3 solve

solution

1- positive seed

put opposite to

in the case of calculation from to divide by then calculate by integrating once with respect to then multiply the result by .

in the case of calculation from to divide by then calculate by integrating twice with respect to .

if we proceed more, the power of will remain unchanged, that is .

2- negative seed

put opposite to

in the case of calculation from to differentiate once and multiply by .

in the case of calculation from to differentiate twice and multiply by .

since the resulting inspected function is a function of , then

substituting in the differential equation we get

minimizing

which gives us two inspected functions :

the complementary function for the above differential equation is then

in the above differential equation the method of systematic inspection could not find a direct solution but it could discover the behavior of the differential equation from which we derived the solution. We note here that when the dimension of the independent variable minus the dimension of the dependent variable is equal in any two terms in the differential equation then, if a seed is put opposite to any of these terms then the half of the inspected function created by this seed will be misleading by the process of calculate and balance, and you need to solve this half of the inspected function the same way I have done in the previous example. The following example gives more light on this problem.

2.4 example

solve

provided

here we note that the second and the third terms are having the same difference between the dimension of and the dimension of which is and respectively.

solution

for positive seed we put opposite to . (the half of the inspected function created by this seed will be misleading).

…

for every calculation from to divide by then differentiate once then multiply by .

For every calculation from to divide by then differentiate twice.

We note here that if then will be equal to or , but for the repetition problem of the power of will appear so we use (because all other terms in the series of will be eliminated for ) and we do the same as we have done in the previous example.

substituting in the differential equation we get

or

then

then for

now we go for the negative seed

put opposite to (the half of the inspected function created by

this seed will be correct).

for every calculation from to we integrate once then multiply by .

for every calculation from to we integrate twice then multiply by.

for

…

for

…

final solution

for

…

for

…

3. solving partial differential equations

inspecting functions :

in the case of partial differential equations in terms of and let the starting seed to be .

3.1 example

solve the following differential equation (torsion problem)

where is the modulus of rigidity, is the angle of twist then is a constant

for the following boundary conditions

at&

at&

Solution

Positive seed

Put opposite to

in the case of calculation, integrate twice with respect to then calculate by differentiating twice with respect to .

to calculate integrate twice with respect to

is convergent for , any value

Negative seed

put opposite to

in the case of calculation integrate twice with respect to then calculate by differentiating twice with respect to .

to calculate integrate twice with respect to

is convergent for , any value

is convergent for , any value.

Radius of convergence is ,

3.2 Particular integral

then (the other way around is also ok i.e. )

now let us get enough functions by substituting different values of and into the above inspected function form :

here I have selected even values for & because of the symmetry of boundary conditions.

now multiplying any of these inspected functions with a constant will not affect it.

so to make them easier to handle multiply

by

by

by

by

by

by

by

As a result the inspected functions become

Now

or

now subjecting the above formula to the boundary conditions we get,

at and

now collecting the constants accompanied with , , , we get the following equations :

..(1)

..(2)

..(3)

..(4)

Subjecting the formula to the second boundary condition,

at and

now collecting the constants accompanied with , , , we get the following equations

..(5)

..(6)

..(7)

..(8)

Now solve these equations simultaneously to get the unknown constants C, C₀₀,

,,,,,.

In the above example I created 8 equations, but you can create more if you need more accuracy.

4. Solving differential equations that include functions

Find the particular integral for

solution

put opposite to

knowing that :

from this we find that :

in the case of calculation from to divide by -3 then integrate once w.r.t. x

in the case of calculation from to divide by –3 then differentiate once w.r.t. x

5. Solving non-linear differential equations

Solving a non-linear differential equation is like solving a polynomial, so at the beginning we will solve a polynomial. Let’s start by solving the following quadratic equation,

5.1 solve

put the equation in two rows

make a calculation

now subtract from the

now make a balance

now make the calculation by neglecting what you had previously for

make a balance

calculate by neglecting the previous value of

proceed for more accuracy until you get the difference (result) = 0, this will give us,

now we need to get the other root.

If you try to put opposite to as follows, you will get a divergent solution as you see down :

balance

the result is that is diverging as we proceed. So how to get the other root ?

our equation is

dividing the equation by will not affect it, so we get

or

now you can start the process of calculate and balance as follows,

now we make a calculation

subtract (or add 1)

make a balance

subtract (or add 1)

make a balance

subtract (or add 1)

if you proceed with the process of calculate and balance you will get :

, which is the second root.

This method works fine for real roots, but for imaginary roots it gives a series

that is a function of i ()

5.2 Solve

Solution

Calculate

result

balance

result

so the first root is

to find the second root try the following

calculate

result

Balance

Balance

proceed to get a more convergent solution, so

I did not check this solution. But if the solution is divergent or is identical to the first root, then you need to change the structure of the differential equation to get the second root as I have done with the quadratic equation.

6. Solving simultaneous ordinary linear differential equations

Solve the following set of ordinary differential equations :

…. (1)

…. (2)

…. (3)

Solution

Positive seed: put opposite to

make a calculation from in the first equation to in the third equation by integrating once then multiplying by . make a balance to then make a calculation to in the middle equation by integrating once and multiplying by .

make a balance from to in the middle equation then make a calculation from in the middle equation to in the first equation by integrating once and multiplying by then make a balance from to in the first equation. Repeat the process till you get a satisfactory number of terms.

this gives

…

…

…

negative seed

put opposite to

make a calculation from in the first equation to in the middle equation by dividing by then differentiating once then make a balance in the middle equation form to . make a calculation from in the middle equation to in the third equation by dividing by then differentiating once.

Make a balance in the third equation from to in the same equation then make a calculation from in the third equation to in the first equation by dividing by and differentiating once.

Make a balance in the first equation from to in the same equation then repeat the process mentioned above.

this gives

…

For

for

for

if you try you will get similar functions .

final solution

in the previous example we had two terms per equation so the process was very simple. But if we have more than two terms in one or some of the simultaneous differential equations then the process will be more difficult. The following example shows this and shows also how to deal with simultaneous partial differential equations.

7. solving simultaneous linear partial differential equations

Find one variation of the particular integral for the following set of linear partial differential equations:

the method of solution is as follows

1 – put opposite to , opposite to , and opposite to

2 – calculate from and from

3 – calculate from and from

4 – calculate from and from

5a – in the second equation balance 2x from ∂u/∂x into ∂v/∂y and balance z³/3 from ∂w/∂z into ∂u/∂x. (not necessarily into ∂u/∂x, you can balance ∂w/∂z into ∂v/∂y).

5b – in the third equation balance 2 from ∂²u/∂y² into ∂²v/∂y² and balance 2y from ∂²v/∂y² into ∂²u/∂x².

6 – now make a calculation from the balances mentioned in (5a)&(5b) as follows,

6a – in the second equation make a calculation form –2x of ∂v/∂y into v in the first equation, this gives –2xy and make a calculation from –z³/3 of ∂u/∂x into u in the first equation, this gives –z³x/3 6b – in the third equation, make a calculation from –2 of ∂²v/∂y² into ∂v/∂y in the second equation, this gives –2y, and make a calculation from –2 of ∂²v/∂y² into v in the first equation, this gives –2y².

7 – now make the first balance in the first equation by balancing v into u. that is to balance (y³/3 – 2xy –y²) from v into u, and balance u and w into v, that is to balance (yx² - z³x) from u and z⁴/12 from w into v.

by finishing this process we finish the first run of calculate and balance process.

Next we will start the second run.

8a –make a calculation from u into ∂u/∂x, this gives 2y in ∂u/∂x, and make a calculation from u into ∂²u/∂x², this gives 0.

8b – make a calculation from v into ∂v/∂y, this gives x² in ∂v/∂y, and make a

calculation from v into ∂²v/∂y², this gives 0.

9 – now make balances in every equation as follows,

9a – in the second equation, we need to balance 2y form ∂u/∂x into ∂v/∂y, but since the second equation is already balanced, that is –2y already exists, so the balance is not necessary.

(note) : after every calculation run for an equation, we should check the balance of the equation so that no disturbance in the balance occurs.

9b – in the third equation, no balance required, since the balance results are zeros.

10 – in the second equation, make a calculation for –x² from ∂u/∂x into u, this gives -x³/3, and make a calculation for –x² from ∂u/∂x into ∂²u/∂x²,

this gives –2x.

11 – in the third equation make a balance for –2x from ∂²u/∂x² into ∂²v/∂y².

12 – in the third equation make a calculation for +2x from ∂²v/∂y² into

∂v/∂y in the second equation, this gives +2xy, and make a calculation for +2x from ∂²v/∂y² into v in the first equation, this gives xy².

13 – now make the second run balance in the first equation, that is to balance from u into v which balances –x³/3 from u into v, and make a balance from v into z which balances xy² from v into z.

by this we terminate the process of calculation, because any calculation process gives zero in the second and the third equation.

Note that since the number of arrows pointing from one term to the other is much, so I did not put all the arrows.

Solution

Conclusion

I want to conclude this work by stating that this method gives us a new way of thinking of differential equations. that’s why I used simple ones. Just to clear up the concept. Every differential equation has its own behavior and should be treated and analyzed differently. Differential equations with more than two terms need computer to be solved, since the number of the processes of calculate and balance becomes huge.