Method of Systematic Inspection

2009 October 22
by admin

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, 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.


252 Responses leave one →
  1. engerrish permalink
    December 19, 2009

    Потрясающий материал. Много контрастов и ценного материала. Особая благодарность автору!

  2. PerryXV permalink
    December 22, 2009

    Здравсивуйте. поддерживаю мнение предыдущего автора!

  3. December 30, 2009

    Привет всем, подскажет кто нибудь такой же блог по тематике рыбалка, а то везде один бред.

  4. Dima.Gulin permalink
    December 31, 2009

    И что вы все тут нашли такого, непонятно.

  5. JessieHO permalink
    January 3, 2010

    С праздником вас, спасибо за статью!

  6. January 9, 2010

    Нашел этот блог и очень рад этому. как и многие другие.

  7. January 10, 2010

    Неплохо сделали блог, молодцы.

  8. January 10, 2010

    Внимательно изучил, в целом материал интересный.Спасибо

  9. January 11, 2010

    Спасибо за материал, хороший блог, только вот не могу рсс ленту найти, она тут есть?

  10. admin permalink*
    January 14, 2010

    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.

  11. Aristofanis Halivopoulos permalink
    January 15, 2010

    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!)

  12. Slavik.Zero permalink
    January 20, 2010

    Хорош блог, выполнен на отлично.

  13. admin permalink*
    January 20, 2010

    Dear Mr. Aristofanis Halivopoulos,
    Dear Mr. Slavik.Zero,
    I appreciate your beautiful comments,
    Best Regards.

  14. January 24, 2010

    Может кто нибудь все таки что то умное напишет, как роботы одно и тоже все пишите.

  15. January 25, 2010

    А почему Вы сдесь об этом пишите, можно в другом месте обсудить.

  16. February 4, 2010

    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

  17. admin permalink*
    February 5, 2010

    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.

     

  18. admin permalink*
    February 6, 2010

    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…

     

  19. February 7, 2010

    winsome answers i like it

  20. February 7, 2010

    Здравствуйте,прелесный материал. очень глубоко захвотили данную тему. большое вам спасибо.

  21. MartyTZ permalink
    February 7, 2010

    Здравствуйте, очеьнь понравился ваш блог, отлично подбираете и доносите для читатиеля материал.

  22. itmagz permalink
    February 9, 2010

    Приветствую, где вы заказывали дизайн для сайта? очень очень понравился.. хочу такоже

  23. February 9, 2010

    Так же поддерживаю выше сказанное о блоге, хорош

  24. February 9, 2010

    Об этом можно и не расписывать на каждом углу, и так все понятно.

  25. webmagz permalink
    February 9, 2010

    Странно, я не когда об этом не задумывался…надо будет хорошо подумать.

  26. February 9, 2010

    Отличный блог, тематика та что я искал и сделано неплохо.

  27. February 9, 2010

    Согласен, тематика подобрана неплохо, и по качеству тут тоже не поспоришь.

  28. businessmagz permalink
    February 9, 2010

    Привет, вы очень доходчиго пишите , нам – читателям – это очень нравиться – спасибо!

  29. February 10, 2010

    Да комменты надо подчищать, чтоб ничто не портило впечатление от блога.

  30. February 10, 2010

    А чему вы тут учится многие собрались, все по обычным стандартам сделано.

  31. February 10, 2010

    This is a superior web page.

    I’ve bookmarked this web page and also I will inform my pals about it.

    Thankfully

  32. February 10, 2010

    Согласен, блог действительно сделан для людей.

  33. February 10, 2010

    Пособие то нормальное, только можно это ни здесь расписывать, а в другом месте.

  34. February 10, 2010

    This is a highly rated website.

    I have bookmarked this particular site and I am going to tell my friends over it.

    Thank you

  35. MounkVodotuct permalink
    February 10, 2010

    А то! :) согласен. там реально полезные мысли

  36. Timur.Bahtin permalink
    February 11, 2010

    Отличный блог никто не спорит, еще бы почаще обновления были.

  37. February 11, 2010

    Тонко и иронично… согласен. там реально полезные мысли

  38. February 11, 2010

    Все понравится конечно не может, но в основном положительное отношение к блогу.

  39. Nikolay.Shipenkov permalink
    February 11, 2010

    Усть что то сдесь такое привлекательное не спорю.

  40. February 11, 2010

    А то! :) согласен. достаточно познавательная

  41. February 11, 2010

    Да уж. А еще говорят, что мировой кризис по миру гялет…

  42. admin permalink*
    February 11, 2010

    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.

  43. February 12, 2010

    Блог то отличный, только что об этом расписывать на каждом углу, и так всем понятно.

  44. February 12, 2010

    Вообще да, мне это напоминает ситуацию с выборами призидента Украины. Выборы прошли, а не факт, кто призедентом будет. Да, блин, демократия. Хех, простите за флуд небольшой, навеяло просто.

  45. February 12, 2010

    Отличный блог, все без исключения понравилось. так держать.

  46. February 12, 2010

    давно искал такой материал, перерыл весь рунет! рад что наконец нащёл

  47. February 13, 2010

    блог вполне нормальный, но дизайн поменять стоит.

  48. February 13, 2010

    Согласен с аппонентом, слишком много лишнего, надо почистить.

  49. February 13, 2010

    Хватит засирать блог, панацея прям какая то.

  50. February 13, 2010

    В общем блог не плохой, но и соглашусь с тем что почистить не мешало бы.

  51. February 13, 2010

    давно искал подобный материал, перерыл весь инет, а тут на тебе. благодарю автора

  52. February 13, 2010

    так же думаю, что это лучше в другой теме обсуждать и не засирать комментами блог.

  53. February 13, 2010

    Поддерживаю большенство и обязательно добавляю в избранное

  54. February 13, 2010

    материал заслуживает повышенного внимания. дабы изучить более подробно, добавлю в избранное

  55. February 13, 2010

    А мне все понравилось, качество есть, а зсрано во всех блогах почти.

  56. February 13, 2010

    Отличный блог, мне тут все по душе.

  57. February 13, 2010

    Неплохо сделали, тематика нужная для меня лично.

  58. February 13, 2010

    Это точно, зато все лезут со своей колокольни.

  59. February 13, 2010

    Да уж, тут смотрю и немногие в этом разбераются.

  60. February 13, 2010

    Автор жжет, спасибо – интересная новость.

Leave a Reply

Note: You can use basic XHTML in your comments. Your email address will never be published.

Subscribe to this comment feed via RSS