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
and
opposite to
like 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, C00,
,
,
,
,
,
.
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 z3/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 ∂2u/∂y2 into ∂2v/∂y2 and balance 2y from ∂2v/∂y2 into ∂2u/∂x2.
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 –z3/3 of ∂u/∂x into u in the first equation, this gives –z3x/3 6b – in the third equation, make a calculation from –2 of ∂2v/∂y2 into ∂v/∂y in the second equation, this gives –2y, and make a calculation from –2 of ∂2v/∂y2 into v in the first equation, this gives –2y2.
7 – now make the first balance in the first equation by balancing v into u. that is to balance (y3/3 – 2xy –y2) from v into u, and balance u and w into v, that is to balance (yx2 - z3x) from u and z4/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 ∂2u/∂x2, this gives 0.
8b – make a calculation from v into ∂v/∂y, this gives x2 in ∂v/∂y, and make a
calculation from v into ∂2v/∂y2, 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 –x2 from ∂u/∂x into u, this gives -x3/3, and make a calculation for –x2 from ∂u/∂x into ∂2u/∂x2,
this gives –2x.
11 – in the third equation make a balance for –2x from ∂2u/∂x2 into ∂2v/∂y2.
12 – in the third equation make a calculation for +2x from ∂2v/∂y2 into
∂v/∂y in the second equation, this gives +2xy, and make a calculation for +2x from ∂2v/∂y2 into v in the first equation, this gives xy2.
13 – now make the second run balance in the first equation, that is to balance from u into v which balances –x3/3 from u into v, and make a balance from v into z which balances xy2 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.
Потрясающий материал. Много контрастов и ценного материала. Особая благодарность автору!
Здравсивуйте. поддерживаю мнение предыдущего автора!
Привет всем, подскажет кто нибудь такой же блог по тематике рыбалка, а то везде один бред.
И что вы все тут нашли такого, непонятно.
С праздником вас, спасибо за статью!
Нашел этот блог и очень рад этому. как и многие другие.
Неплохо сделали блог, молодцы.
Внимательно изучил, в целом материал интересный.Спасибо
Спасибо за материал, хороший блог, только вот не могу рсс ленту найти, она тут есть?
Dear Mr. RHStephen,
Dear Mr. engerrish,
Dear Mr. PerryXV,
Dear Mr. Ivan.Kulikov,
Dear Mr. Dima.Gulin,
Dear Mr. JessieHO,
Dear Mr. ecommerce,
Dear Mr. psccp,
I appreciate your beautiful comments,
Best Regards.
Dear Mr. belvuhotel,
It’s down there before the “Submit Comment” button.
Thanks and Best Regards.
Thank you very much for showing and verifying that complex thought is composed by small simple parts (an ancient way of thinking but after thousand of years seems yet the rightest!)
Хорош блог, выполнен на отлично.
Dear Mr. Aristofanis Halivopoulos,
Dear Mr. Slavik.Zero,
I appreciate your beautiful comments,
Best Regards.
Может кто нибудь все таки что то умное напишет, как роботы одно и тоже все пишите.
А почему Вы сдесь об этом пишите, можно в другом месте обсудить.
Very interesting.
Can you propose general solution for:
Y’- YY=f(x) yy= y to the second power
This is very important Riccati Equation. Please advise
Kindest Regards,
Nicos H. Andreas
Riccati Equation is no exception but you must define
f(x), so let’s assume that f(x)=x. then the equation to solve is the following:
y’ – y2 = x
so lets start
for the first root try
y’ = x
-y2 = -x4/4
Sum = x – x4/4
-x
= -x4/4
Balance
y’ = x + x4/4
-y2 = -(x2/2 + x5/20)2
Sum = x + x4/4 – (x2/2 + x5/20)2
Or
y’ = x + x4/4
-y2 = – x4/4 – x7/20
– x10/400
Sum = x – x7/20 – x10/400
-x
= – x7/20 – x10/400
Balance
y’ = x + x4/4 + x7/20 + x10/400
-y2 = -x4/4
– x7/20 – 7x10/800 – 3x13/3520 – 87x16/1408000
– x19/352000 – x22/19360000
Sum = x – x10/160 – 3x13/3520 –
87x16/1408000 – x19/352000 – x22/19360000
-x
Balance
This gives us y’ which is
y’ = x + x4/4 + x7/20 + 7x10/800
+ 3x13/3520 + 87x16/1408000 + x19/352000 + x22/19360000
I will stop here to get y … for you, you can continue
for more accuracy.
y = x2/2 + x5/20 + x8/160
+ 7x11/8800 + 3x14/49280 + 87x17/23936000 + x20/7040000
+ x23/445280000
for the second root try
-y2 = x
y’ =
Calculate
-y2 = x
y’ = ix-1/2/2
sum = x + ix-1/2/2
-x
= ix-1/2/2
-y2 = x – ix-1/2/2
y’ = (-x + ix-1/2/2)-1/2.(-1 –
ix-3/2/4)/2
proceed this way until you get a satisfactory number of
terms then calculate y from y’. this will be the second root.
Trying to get general solution
First Root
y’ = f(x)
-y2 = -(∫f(x)dx)2
y = ∫f(x)dx
sum = f(x) – (∫f(x)dx)2
-f(x)
= -(∫f(x)dx)2
Balance
y’ = f(x) + (∫f(x)dx)2
-y2 = -(∫(f(x) + (∫f(x)dx)2)dx)2
y = ∫(f(x) + (∫f(x)dx)2)dx
sum = f(x) + (∫f(x)dx)2 – (∫(f(x) + (∫f(x)dx)2)dx)2
-f(x)
= (∫f(x)dx)2 – (∫(f(x) + (∫f(x)dx)2)dx)2
Balance
y’ = f(x) + (∫(f(x) + (∫f(x)dx)2)dx)2
-y2 = -(∫(f(x) + (∫(f(x) + (∫f(x)dx)2)dx)2)dx)2
y = ∫(f(x) + (∫(f(x) + (∫f(x)dx)2)dx)2)dx
proceed for more accuracy…
Second Root
-y2 = f(x)
y’ = d/dx(-f(x))1/2
y = (-f(x))1/2
sum = f(x) + d/dx(-f(x))1/2
-f(x)
= d/dx(-f(x))1/2
Balance
-y2 = f(x) – d/dx(-f(x))1/2
y’ = d/dx(-(f(x)-d/dx(-f(x))1/2))1/2
y = (-(f(x) – d/dx(-f(x))1/2))1/2
sum = f(x) – d/dx(-f(x))1/2 + d/dx(-(f(x) – d/dx(-f(x))1/2))1/2
-f(x)
= -d/dx(-f(x))1/2 + d/dx(-(f(x) – d/dx(-f(x))1/2))1/2
Balance
-y2 = f(x) – d/dx(-(f(x) – d/dx(-f(x))1/2))1/2
y’ = d/dx(-(f(x) – d/dx(-f(x) – d/dx(-f(x))1/2))1/2))1/2
y = (-(f(x) – d/dx(-(f(x) – d/dx(-f(x))1/2))1/2))1/2
proceed for more accuracy…
winsome answers i like it
Здравствуйте,прелесный материал. очень глубоко захвотили данную тему. большое вам спасибо.
Здравствуйте, очеьнь понравился ваш блог, отлично подбираете и доносите для читатиеля материал.
Приветствую, где вы заказывали дизайн для сайта? очень очень понравился.. хочу такоже
Так же поддерживаю выше сказанное о блоге, хорош
Об этом можно и не расписывать на каждом углу, и так все понятно.
Странно, я не когда об этом не задумывался…надо будет хорошо подумать.
Отличный блог, тематика та что я искал и сделано неплохо.
Согласен, тематика подобрана неплохо, и по качеству тут тоже не поспоришь.
Привет, вы очень доходчиго пишите , нам – читателям – это очень нравиться – спасибо!
Да комменты надо подчищать, чтоб ничто не портило впечатление от блога.
А чему вы тут учится многие собрались, все по обычным стандартам сделано.
This is a superior web page.
I’ve bookmarked this web page and also I will inform my pals about it.
Thankfully
Согласен, блог действительно сделан для людей.
Пособие то нормальное, только можно это ни здесь расписывать, а в другом месте.
This is a highly rated website.
I have bookmarked this particular site and I am going to tell my friends over it.
Thank you
А то!
согласен. там реально полезные мысли
Отличный блог никто не спорит, еще бы почаще обновления были.
Тонко и иронично… согласен. там реально полезные мысли
Все понравится конечно не может, но в основном положительное отношение к блогу.
Усть что то сдесь такое привлекательное не спорю.
А то!
согласен. достаточно познавательная
Да уж. А еще говорят, что мировой кризис по миру гялет…
Dear Mr. DenueHakevack,
Dear Mr. seoshnikov,
Dear Mr. MartyTZ,
Dear Mr. itmagz,
Dear Mr. Sergey.Grinev,
Dear Mr. Oleg.Trubinev,
Dear Mr. webmagz,
Dear Mr. Kolyan.Dubev,
Dear Mr. Oleg.Proninko,
Dear Mr. businessmagz,
Dear Mr. Slava.Nikolaiko,
Dear Mr. Artem.Torgashin,
Dear Mr. credit repairfanceexpern,
Dear Mr. Vlad.Braginev,
Dear Mr. Grigoriy.Bunich,
Dear Mr. improve credit scorefanceexpern,
Dear Mr. MounkVodotuct,
Dear Mr. Timur.Bahtin,
Dear Mr. Mupduegendrag,
Dear Mr. Viktor.Senin,
Dear Mr. Nikolay.Shipenkov,
Dear Mr. Dorruinoben,
Dear Mr. Martinis,
I appreciate your encouraging comments.
Thanks and Best Regards.
Yousif.
Блог то отличный, только что об этом расписывать на каждом углу, и так всем понятно.
Вообще да, мне это напоминает ситуацию с выборами призидента Украины. Выборы прошли, а не факт, кто призедентом будет. Да, блин, демократия. Хех, простите за флуд небольшой, навеяло просто.
Отличный блог, все без исключения понравилось. так держать.
давно искал такой материал, перерыл весь рунет! рад что наконец нащёл
блог вполне нормальный, но дизайн поменять стоит.
Согласен с аппонентом, слишком много лишнего, надо почистить.
Хватит засирать блог, панацея прям какая то.
В общем блог не плохой, но и соглашусь с тем что почистить не мешало бы.
давно искал подобный материал, перерыл весь инет, а тут на тебе. благодарю автора
так же думаю, что это лучше в другой теме обсуждать и не засирать комментами блог.
Поддерживаю большенство и обязательно добавляю в избранное
материал заслуживает повышенного внимания. дабы изучить более подробно, добавлю в избранное
А мне все понравилось, качество есть, а зсрано во всех блогах почти.
Отличный блог, мне тут все по душе.
Неплохо сделали, тематика нужная для меня лично.
Это точно, зато все лезут со своей колокольни.
Да уж, тут смотрю и немногие в этом разбераются.
Автор жжет, спасибо – интересная новость.