Научная статья на тему 'Построение интерполяционного метода численного решения задачи Коши'

Построение интерполяционного метода численного решения задачи Коши Текст научной статьи по специальности «Математика»

CC BY
41
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МИНОРАНТ НЬЮТОНА / NEWTON'S MINORANT / ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ / DIFFERENTIAL EQUATIONS / ЗАДАЧА КОШИ / CAUCHY'S PROBLEM / ДИАГРАММА / ВЫПУКЛАЯ ФУНКЦИЯ / CONVEX FUNCTION / NEWTON'S DIAGRAM

Аннотация научной статьи по математике, автор научной работы — Bihun R., Tsehelyk G.

Построен интерполяционный численный метод решения задачи Коши для обыкновенных дифференциальных уравнений первого порядка с помощью аппарата неклассических минорант и диаграмм Ньютона функций, заданных таблично. Этот метод дает более точные результаты по сравнению с методом Эйлера в случае выпуклой функции. Доказана вычислительная устойчивость метода, то есть погрешность начальных данных не накапливается. Также показано, что метод имеет второй порядок точности

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Construction of interpolation method for numerical solution of the Cauchy''s problem

An interpolation numerical method is developed in order to solve the Cauchy's problem for ordinary first order differential equations using the apparatus of non-classical minorants and diagrams of Newton's functions, assigned in a tabular form. We have proven computational stability of the method, that is, an error of the initial data is not piled up. It is also shown that the method possesses a second order of accuracy and in the case of a convex function produces more accurate results than the Euler's method. The advantages also include simplicity and visual clarity of the method. Given this, it could gain widespread use in many areas, in particular mathematics, physics and mechanics. We also give an example of solving the Cauchy's problem applying the new method, the Euler's method, and the Runge-Kutta fourth order method, with the results compared. The proposed method does not require solving the systems of linear algebraic equations because we do not employ the Bernstein polynomials, and it is not required to superimpose additional conditions, in contrast to the method that applies the Haar functions

Текст научной работы на тему «Построение интерполяционного метода численного решения задачи Коши»

29. Pawlak, Z. Rough sets [Text] / Z. Pawlak // International Journal of Computer & Information Sciences. - 1982. - Vol. 11, Issue 5. - P. 341-356. doi: 10.1007/bf01001956

30. Raskin, L. Fuzzy models of rough mathematics [Text] / L. Raskin, O. Sira // Eastern-European Journal of Enterprise Technologies. -2016. - Vol. 6, Issue 4 (84). - P. 53-60. doi: 10.15587/1729-4061.2016.86739

31. Raskin, L. G. Prikladnoe kontinual'noe lineynoe programmirovanie [Text] / L. G. Raskin, I. O. Kirichenko, O. V. Seraya. - Kharkiv, 2013. - 293 p.

Побудовано ттерполяцшний чисельний метод розв'язування задачi Кошi для звичайних диференщ-альних рiвнянь першого порядку за допомогою апа-рату некласичних мторант та дiаграм Ньютона функцш, заданих таблично. Цей метод дае точт-шi результати вiд методу Ейлера у випадку опуклог функцш. Доведено обчислювальну стштсть методу, тобто похибка початкових даних не нагромаджу-еться. Також показано, що метод мае другий порядок точностi

Ключовi слова: мторанта Ньютона, диференщ-альш рiвняння, задача Кошi, дiаграма Ньютон, опукла функщя

Построен интерполяционный численный метод решения задачи Коши для обыкновенных дифференциальных уравнений первого порядка с помощью аппарата неклассических минорант и диаграмм Ньютона функций, заданных таблично. Этот метод дает более точные результаты по сравнению с методом Эйлера в случае выпуклой функции. Доказана вычислительная устойчивость метода, то есть погрешность начальных данных не накапливается. Также показано, что метод имеет второй порядок точности

Ключевые слова: минорант Ньютона, дифференциальные уравнения, задача Коши, диаграмма, выпуклая функция

-□ □-

UDC 519.6

|DOI: 10.15587/1729-4061.2017.108327|

CONSTRUCTION OF INTERPOLATION METHOD FOR NUMERICAL SOLUTION OF THE CAUCHY'S PROBLEM

R. B i h u n

Postgraduate student* E-mail: [email protected] G. Tsehelyk

Doctor of Physical and Mathematical Sciences, Professor, Head of Department* Е-mail: [email protected] *Department of Mathematical modeling of socio-economic processes Ivan Franko National University of Lviv Universytetska str., 1, Lviv, Ukraine, 79000

1. Introduction

Cauchy's problem is one of the main problems in the theory of differential equations, which comes down to finding a solution (integral) of the differential equation that satisfies initial conditions (original data).

Over many years, a numerical solution of the Cauchy's problem has been the focus of attention by scientists as it is widely used in different areas of science and technology. That is why there are a large number of developed methods for it. In spite of this, however, new methods are being devised, some of them with better properties than those preceding.

Cauchy's problem usually emerges during analysis of the processes predetermined by the differential law and original state. Mathematical notation of such equations is an equation and the initial condition.

The difference between the boundary-value problems and the Cauchy's problem is that the region over which the desired solution should be determined is not specified in the latter in advance. However, the Cauchy's problem can be considered as one of the boundary-value problems.

2. Literature review and problem statement

Numerical methods of the Cauchy's problem solution are divided into 3 groups [1]:

- one-point;

- multipoint (methods of prediction and correction);

- methods with automatic choice of integration step.

The one-point methods include methods that have certain

common features, such as:

1. Underlying all one-point methods is the function decomposition into Taylor's series, which preserves members that have h in a power to k inclusive. An integer k is called the order of the method. Error on a step has an order of k+1.

2. All one-point methods do not require a valid computation of derivatives, because only the function itself is calculated, however, one may require its values in some intermediate points. This entails, of course, additional cost of time and effort.

3. In order to receive information in a new point, it is necessary to have data only from the previous point. This property can be called "self-starting". A capability to "self-start" makes it possible to easily change the magnitude of step h.

©

4. Compared with the one-point methods, the methods of prediction and correction possess a number of special features [2].

1) To implement the methods of prediction and correction, it is necessary to have information about several of the previous points (they do not belong to the "self-starting" methods), which is why, in order to obtain additional information, it is necessary to apply the one-point method. If in the process of solving differential equations by the method of prediction and correction the step changes, then one has typically to switch over temporarily to the one-point method.

2) One-point methods and methods of prediction and correction provide approximately the same accuracy of the results. However, the latter, in contrast to the former, make it possible to estimate only an error in a step. For this reason, when employing the one-point methods, the magnitude of step h is typically chosen slightly smaller than it is required, which is why the methods of prediction and correction prove to be the most effective.

3) When applying the Runge-Kutta method [2] of the fourth-order accuracy, at each step one has to calculate four values of the functions, but for the convergence of the method of prediction and correction of the same order of accuracy, it is often sufficient to have two values of the function. That is why methods of prediction and correction require almost twice less computing time than the Runge-Kutta methods of comparable accuracy.

To solve differential equation y'= f (x, y) by a numerical method means to find for the assigned sequence of arguments x0,x1v..,xn and y0 such values of y0,y1v..,yn so that yi = F(xi) (i = 1,2,...,n) and F(x0) = y0. Thus, the numerical methods make it possible, instead of deriving function y = F (x), to receive a table of the values of the given function for the assigned sequence of arguments. The magnitude h = xk - xk1 is called a step of integration.

Graphically, numerical solution [3] is a sequence of short straight-line segments, by which analytical solution y = F(x) of the equation is approximated (a piecewise-linear approximation).

There are methods of differential transformation (MDT) of solution to the Cauchy's problem. Basic definitions and fundamental theorems of a one-dimensional MDT and its suitability for different types of differential and integraldifferential equations are given in [4].

A reliable, yet very simple, numerical method to solve different cases of a singular Cauchy-type integral equation is developed in [5]. For this purpose, first Bernstein polynomials are derived, which are used to approximate a solution of the given singular integral equation. This, however, leads to solving the system of linear algebraic equations (SLAE), which sometimes is difficult to resolve.

Article [6] examines numerical solution of the class of systems of singular integral Cauchy's equations with constant coefficients. The proposed procedure consists of two main stages: the first is to consider a modified problem, equivalent to the original under appropriate conditions, the second is to bring its solution using a vector of polynomial functions. But the solution comes down to solving the linear systems.

By applying the Haar functions [7], it is possible to receive a solution with a very small error, more accurate in some cases than the solution derived by the second order Runge-Kutta method. But the function must be superimposed with certain conditions.

3. The aim and objectives of the study

The goal of present work is to construct a numerical method to solve the Cauchy's problem for ordinary first order differential equations, which would yield more accurate results than the classical methods. The new method should not require solving a system of linear algebraic equations and should not require superimposing of conditions on the function.

To accomplish the goal, the following tasks have been set:

- to develop a new interpolation numerical method, employing the apparatus of non-classical Newton's minorants in order to solve the Cauchy's problem for ordinary first order differential equations;

- to prove the computational stability and convergence of the method;

- to assess accuracy of the solution and error of the method.

4. Materials and methods for examining a solution to the Cauchy's problem

4. 1. Formulae of the minorant type for the approximated calculation of definite integrals

There are a number of approaches in order to construct formulae for the approximated calculation of definite integrals: replacing an integrand function with interpolation polynomial [8], the use of Bernoulli numbers and polynomials [9]. In practice, however, most often used are quadrature formulae of the interpolation type. A question then arose: is it at all possible to employ the apparatus of non-classical minorants of Newton's functions, assigned in a tabular form, to solve the Cauchy's problem as well? In this case, the new method would not have to solve SLAE since this is a cumbersome process. In [10], the apparatus of non-classical minorants of Newton's functions, assigned in a tabular form, was used for finding zeros of the function.

Assume that it is required to compute a definite integral

J = \f (x) dx.

(1)

Without loss of generality, we shall consider that f (x)> 0 for all x e[fl, b]. Let us build for function y = f (x) a Newton's minorant [11] by two points x1 = a, x2 = b. We obtain

) = ( Ab

(2)

where A = f (a), B = f (b). We shall replace at interval [a, b] function f (x) with Newton's minorant mf (x) [12]. We receive

J = }f (x) dx = \mf (x) dx + R (f),

where R (f) is the remainder. Calculate

b b 1 S = fmf (x) dx = J(Ab-xBx-a dx.

At A * B (if A = B, then S = A(b- a) = obtain:

(3)

(4)

- a)) we

( Ab

S =

Ba

VB',

VB,

A ^

Ba ,

Thus,

b

= J ( A a x \b-a dx

b - a

ln (B)- ln ( A)

b - a

ln (B)- ln ( A)

b - a

B

i

b-a

b - a

-i B

ln ( B)-ln ( A)l, A ) a

b-a b

B ) b-a ( B \ b-a

A) 1A,

B(B | b-

AIA J

B | b-A J

(5)

ln (B)-ln (A) I, A) A f (b)- f (a

b-aB - A ,, , B - A = (b - a)-

'ln ( B)- ln ( A)'

[f(x)dx = (b-a)—+ R(f). (6)

' 1 j ln (/ (b) / / (a)) f

The resulting formula is called a small formula of the minorant type [11] for approximated calculation of definite integrals.

Let us prove that

limffL = f (a )

fUiln ( f (b) / f (a )) f (a)-

t(a)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Actually,

llm f(b)-f(a) = llm

f^lnf (b) / f (a))

ln

M vf (a ),

f (b)- f (a.

llm ln

fbUi f (a)

m

J (ab

f (b)- f (a

llm ln

f(b)Ui

f (a)

f («

llm ln

f(bUi

f (a)

1 +

f (b)-f (a) f (a)

f (a) f (b)-f (a

1 + f (b)-f (a ) . f (a )

- = f (a ).

f (b)-f (a!

If function f (x)eC[a,b] on Interval [a,b] s satisfies the Lipschitz condition with constant L and within this interval f '(x) does not change a sign, then

( f HI f ( x )-mf ( x )| dx <

a

b b <J| f (b)-f (a) dx < Lf\(b - a) dx < L (b - a )2.

(9)

where x0 = a, xn = b, and at every interval

[xi,xi+1], i = 0,1,...,n-1,

we shall substitute function / (x) with Newton's minorant m/ (x). We receive

b n-1 b n-1 b

J = [/(x)dx = £[/(x)dx = Xjm/ (x)dx + Rn(/). (10)

Since at f (x, )* f ( x!+i ) xi+i

J Mf (x )dx =( x,+i- x, )

f ( x+i )-f ( x, )

ln

fix,

f ( x

(ii)

then

\f ( x ) dx = h|f ( x^" f(*) + Rn ( f ).

ln

f(*

f ( x

(i2)

The resulting formula is called a combinedformula of the (7) minorant type [11] for an approximate calculation of definite integrals.

We shall estimate the remainder Rn (/) for the constructed formula of the minorant type. Since

n-i ~ i+i

Rn ( f ) = U( f ( x )-Mf ( x )) dx,

then

n-1 i+i

R(f)|<IJ \f (x)-mf (x)|dx.

(i3)

(i4)

If we assume that function f (x) on interval [a, b] satisfies the Lipschitz condition with constant L and on each of (8) the intervals [ xt, x,+i ] is monotonous, then

n-i i+i

R ( f )|<lJ\f ( x,+i )-f ( x, )| dx <

i=0 x,

n-i xi+i

< LZ J ( x,+i- x, )dx.

i=0 x,

Hence

R ( f )|< L

(b - a

(i5)

(i6)

Thus, the following theorem holds. Theorem 1. If function / (x) on interval [a, b] satisfies In order to record the constructed formula of the the Lipschitz condition with constant L and on each of the minorant type for calculating definite integrals, we shall intervals [xi-1,xi], i = 1,2,...,n is monotonous, the formu-split the interval of integration [a, b] into n equal parts of la (12) holds and the evaluation of accuracy is performed

length h =

b - a

with points

xi = a + ih, i = 0,i,...,n,

(16). If function /(x) is convex on interval [a,b], then the inequality holds

a n-1 xi+1

[/ (x) dx <£[ mf (x )dx. (17)

b i=0 x,

Theorem2. If /(x)eC2[a,b] and /(x)>0,xe[a,b],

n

n

then

\R (f )| <— max f (x)max

I 8 y ' *E[ab]

& (.f (x I)

(18)

that is, a quadrature formula (12) will be of second-order accuracy.

The quadrature formula constructed has the same complexity and the same order of accuracy (in case the conditions of Theorem 1 or 2 are fulfilled) compared with the quadrature formula of trapezoids. At the same time, if function f (x) takes the form

f (x )= b exp (cix ),x e[ xiv xi ], i = 1,2,..., n,

then the quadrature formula of the minorant type produces the exact value of the desired integral. This is the advantage of the method.

5. Results of research into an interpolation numerical method to solve the Cauchy's problem for ordinary first order differential equations

5. 1. Algorithm of the method to solve the Cauchy's problem for ordinary first order differential equations

Consider the Cauchy's problem for ordinary first order differential equation

y'=f (x, y\ y (xo )=y>.

(19)

(20)

Integrate this identity at each of the intervals [xi ,xi+1 ] (i = 0,1,..., n -1). We receive

y (xi+1) = y (xi)+J f (x, y (x))dx. (22)

Without reducing generality, we shall assume that f (x, y (x ))> 0 for all

x e[ x0,x0 + a].

Construct for function f (x, y (x)) Newton's minorant mf (x) by two points

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(xi, f (xi,y(xi))) and (xi+l, f (xi+l,y(xi+1))).

We receive

mf (x) = (Af+1-xBf-x )h, (23)

where Ai = f (xi,y(xi)), Bi = f (xi+1,y(xi+1)). Next, we shall substitute integrand function f fx,y(x)) in (22) with a Newton's minorant. We obtain

x+1

y(xi+1 ) = y(xi)+ J mf (x)^ + Ri+l,

where is the remainder. Then, by computing integral at A- ^ Bi, we receive

Let us assume that the solution to this problem should be found on certain interval [x0,x0 + a], where a>0. In this case, we shall consider that in region D, which contains a rectangle

R = {x0 < x < x0 + a,|y - y„| < b}

function f (x, y) is continuous and satisfies the Lipschitz condition by y, that is,

|f (x, y)-f (x, y2 )|< L|y - y21,

where L is some constant while (x, y1) and (x, y2) are any two points of region D. Choose on interval [ x0, x0 + a], system of points x0, x1,..., xn, where

xk = x0 + kh (k = 0,1,...,n), h> 0, xn < x0 + a.

Then, by employing the apparatus of non-classical minorants and diagrams of the Newton's functions, assigned in a tabular form, we shall build interpolation numerical method for solving the problem (19), (20). That is, we shall develop a method for finding the approximated values y0,y1,.,yn of the exact solution y = y(x) to the problem (19), (20) in points x0, x1,., xn.

Let us clarify the issue of convergence, accuracy and computational stability of the techniques.

Let y = y (x) be the desired solution to the problem (19), (20). By substituting it into equation (19), we shall obtain the identity

y (xi+1 ) =

= y(x) + h f (xi+1,y(xi+1))-f (xi,y(xi)) n + R ' ln (f(xi+l, y(xi+1))/ f(xi, y(xi))

(24)

At A. = B

i i

y (xi+1 ) = y (xi) + hf (xi, y (xi)) + Ri+1.

Note that on the basis of superior boundary lim (1 + x) x = e

we receive

lim

h^0

f (xi+l, y (xi+1))-f (xi, y (xi))

ln (f (xi+l, y (xi+l))/ f (xi, y (xi)) = f (xi, y (xi)).

In fact,

lim f (y(xi+l))-f (xi,y(xi)).

h^0

(25)

ln

f (xi+l, y (xi+l)) f (xi, y (xi))

f (xi, y (xi))

y' (x ) = f (x, y (x)).

(21)

lim |ln

h^0

f (y (xi+l))- f (xi, y (xi)) f (xi, y (xi))

f (xi y(xi))

f (x'+l ,y(xj+1))- f ,y(xi))

Since at h ^ 0

f ( x,+v y ( x,+i ))-f (x,, y (x, )) f (x,, y (x, ))

Since

then it means that the boundary (25) holds.

Thus, in order to find the approximate values y0, y1,..., yn of solution

y=f (x)

to the problem (12), (16), we shall obtain formula

f (xl+v y+1)- f (x, )

y,+i = yt + h

ln (f ( x,+v V,+i )/ f ( x,, y, ))

(i — 0,1,..., n -1),

where y0 = y(x0). If

f+df f dx dy

<N <~,

~k+\

J f (x, y(x)) dx = ( x - Xk+i)f (x, y(x ))

xk xk+i

xk+i xk

- J (x - xk+i)d~ dx =

xk+1 1 r

=hf (xk,y(xk))- J(x-xk+i)"dxdx,

(32)

then

Aek — h

f (Xk+l,y.+i)-f (x>,y) - f (xk,y(xk))

(27)

ln

v v

f ( yk+i ) f ( xk, yk )

xk+1 df -J( x - xk+i)dxdx-

(33)

f ( x,+v y,+i ) = f ( xi, y ( xi ))

for certain i, then, as it follows from (25), yi+i will be searched for by formula

y,+i = yt+hf (x,, yt ). (28)

Theorem 3. If in region D, which contains a rectangle

R — {x0 < x < x0 + a,|y - y0| < b},

function f ( x, y) is continuous, satisfies the Lipschitz condition at variable y with constant L and

We can write based on (30)

f(xk+i,yk+i)-f(xk,yk) = f(x ) 8 ' f (xk+i, yk+i Y = f(xk,yk)-Sk

ln

f ( xk, yk )

where 8k (h) ^ 0 at h ^ 0 . That is why

f ( xk+i, yk+i)-f ( xk, yk) f ( x ) — ln f(x(+i,yk+i)1 -f(xk,yk) =

V f (xk, yk) ,

= f ( xk, yt )-f ( xt, y( xk ))-8k (h).

(34)

(29)

where N is some constant, then the approximated values y1, y2,..., yn, found from formula (7), at h ^ 0 uniformly, relative to x, converge to the exact solution

y=y (x)

to the problem (19), (20).

Theorem proving. First of all, note that, similar to (25), the following boundary holds

By employing the Lipschitz condition for function f (x, y) along y, we receive

\f (xk,Vk)-f (xk,y(xk))|<L\Vk-y(xk)| = L\ek|.

And based on condition (29), we obtain

Nh2

(35)

i A f

J (x - xk+i ) d dx

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

<N J

x x(+i | dx —

(30)

IT f ^-f (*k' yk\ — f (xk, yk ).

ln( f (xk+l, V(+i)/ f (xk, yk)

Let

ek = yk- y (xk )

be an error of the approximated value of solution y — y ( x ) to the problem (i9), (20) in point x — xk. Then the error gain at the (k+i)-th step will be equal to

Aek = ek+i -ek = ( V(+i - yt )-( y ( xk+i )- y ( xt )) =

=h .nwii,mC!r?(-(x(31)

Thus,

i

| Aek|<Lh\ek| + ^Nh2 + h| 8k(h)|. From inequality

-v. (36)

(37)

it follows that

e k+i-le d< Lh\ed+^ Nh+h\8 k (h)|

ek+i < (1 + Lh)\eh\ + - Nh2 + h\8 k ( h)|.

or

We have derived a recurrent formula for the estimation of error at the (k+1)-th step via the error at the k-th step. Substituting this formula with k=0.1,..., n-1, we receive

| ej<(i + Lh)|e0| +1Nh2 + h| S0(h)|

|ej < - Nh2 + h| S0(h)|,

since eo = yo - y(xo) = 0;

| e2| <(i + Lh)|ej +1Nh2 + h| Si(h)|<

< iNh2 ((i+Lh) +1) + h ((i + Lh) | S0 (h) | + Si (h) | )

(38)

(39)

where yi+1 is the chosen zero approximation, for example, by the Euler method:

y'+i = y + h/(x,,yt), yo = y(x0).

Accordingly, if for some i and k

f (^ y!1)=f ( x , )

is satisfied, then at this step yj++1) will be searched for by formula

y,(++i) = + hf ( xi, ).

(42)

| en\ <1 Nh2 ((1+Lh)n-1 + (1 + Lh)n-2 + ■■■ +1) + +h((1 + Lh)n-11 8o(h) | + (1 + Lh)n-21 81(h)) + ■■■ + 8„_1(h)|). Assume max I 8,(h) I = 8(h).

0 < k < n-11 '

Then

| en\ < ^ 1 Nh2 + h8(h) j ■ ((1 + Lh)n-1 + (1+Lh)n-2 + ■■■ +1)

or

|e„| < L ^ 2 Nh +8( h)j((1+Lh)n -1).

Because at u>0, the inequality eu > 1+u, holds, then

|e„| < L (2Nh +8(h) j ■ (enLh -1).

If one considers that nh = xn -x0 < a and 8(h) = Ch, where C is some constant, then we finally receive

We shall find at which condition this iterative process converges. An attribute of the end of the iterative process is the condition: |y(+Î!+i) — yi+Îil < e. After this, we accept

y,+i »y(++i) (i = 0,i,n-1).

Consider

ur - yj=

=h

=h

f ( x,+i, y£) - f (x,, yt )

f ( xi+v y,+i) - f ( x,, yt )

ln (f(x,+i, yl+))lf(xi, yt ) ) ln ( f (x!+l, y!+iV f (x,, y, )

f( x,, y,) - f ( xi+v y(k)) f ( x,, y,) - f ( xl+l, y,+i)

ln ( f (x,, y, )/ f (x,+i, yW ) ) ln ( f (x,, y, V f (x^ y,+i ) )

Since at h ^ 0

f (x,, yt)- f (x,+i, y^)

f (x,+i, y^) f ( x,, y,) - f ( x y,+i)

0,

f (xl+l, y,+i)

0,

then

f (x,, y, ) f ( x,+i, y(+i)) = f( x y(k))-S(k)(h) ln(f(x,,y,Vf(x^) ) = f (Xl+1,y,+i)-Sl+l(h),

f ( x,, y,) - f ( x y,+i)

= f (xi+V y,+1)-Sl+1(h),

I l^h ( 1

e < — I — ' n L

(2N + Cj (eaL -1).

(40)

It follows from here that at h ^ 0, regardless of x, we have |e J ^ 0. This means that the approximates values y1,y2,...,yn at h^0 uniformly relative to x converge to the exact solution y = y(x). The theorem is proven.

It follows from the proven theorem that the method possesses first order of accuracy relative to h.

Formula (24) is actually the equation for finding yi+1. That is why, in order to compute yi+1, we shall construct an iterative process

yâ+11 = y,+h

(, = 0,1,..., n -1),

f ( xl+l, yh )- f (x,, y, ) ln (f ( xl+l, yh )/ f (x,, y, )

(41)

ln (f (x, y, Vf (xl+l, y,+i) )

where S,+1)(h) ^ 0 and S,+1(h) ^ 0 at h ^ 0. That is why

|y(r - y,+1 =

= h\ (f(x,+i,y$)-f(x,+i,Уl+l))-(S(+кl)(h)-S,+i(h)) |<

<h| f(x,+i,y^i)-f(x,+i,y,+i) I + hS,+1)(h),

where

S,+1)(h) =S<+)(h)-S,+i(h)|.

By employing the Lipschitz condition at variable y for f (x, y), we receive

|y<r - y,+i| < hL\y(kl - y,+i| + h8,+1)(h). Hence, at k = 0,1, ...n -1, we obtain

or

y™ - y,+i < hL\yd - y,+i + h8 ,+1)(h),

Since, based on superior boundary

\y(i - y,+i\< hL^) - y,+i| + h8 ,+^h) < < (hL)21 - yj + h ( hLS,+1)(h) + Si+^h))

| y m - y,+i\< (hL)^y(+°i) - y,+1 + + h ((hL)n-i 8i+i)(h) + (hL)n-2 Si+;)(h)-

Assume

0 man, S,+m)(h) =8 ,+i(h).

0 < m < n-1

Then

- y^ < (hL)^y(+°i) - y,+1 + +hS i+1(h) ((hL)n-i + (hL)n-2 + ■■■ + hL +1 )

-hL8i+1-2)(h) + 8,+n-1)(h))

yC/ - yj < (hL)n yi+7 - yj + hS ,+i(h)

1 - (hL)n 1 - hL

lim (1 + x) x — e.

We obtain

lim f(x,+i, ^ - f(x,,y,) = f (x,, y ), hmln ( f (x,+i, ^i+iV f (x,, y, ) ) f ( i, ^ *

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

li^^^^i^W^,^ — f(x y )

^ln (f (x,+i, y+i)/f (x,,/,) ) f (xi,yi),

and then

f (xi+v yi+i)-f (xi, y:)

ln (f (xi+l, y i+iVf (x,, y:) f ( xi+l, y,+i) - f (xi, )

— f ( x,, y : )-8, (h),

— f ( x,, yt )-8, (h),

It follows from what we received that at hL < 1 the iterative process converges.

Because a Newton's minorant consists of convex arcs, this method produces more accurate results than known two-point methods, in the case when the function f (x, y(x)) is convex.

ln (f (x:+v y+O/f (x,, y: )',

where 8, (h) ^ 0 at h ^ 0. That is why

e'+i < e' + h ( f (x,, y, ) - f ( x,, y, )) + (8, (h) -8, (h)) |.

Function f (x, y) satisfies the Lipschitz condition at variable y with constant L, which is why

e'+i << + ^1 yi - y,l + h| 8,(h)-8: (h)|

6. Discussion of results: a study of computational stability of the method

Let us consider an issue of computational stability of the given method.

Let y0 be an approximated value of exact initial value y0 while e0 is the absolute error of the initial approximation, that is,

e„ —

e'+1 < (1 + Lh) e' + h8:(h), where

8 (h)—h| 8 (h)-8S (h)|.

But 8, (h) ^ 0 at h ^ 0. Then at i — 0,1,..., n -1 we receive By denoting

Then, instead of formula (22) to cal- e'<(1 +Lh) e° +^h); _ _

culate the approximated values of solution e^ < (1 + Lh) e' + h8,(h) < (1+Lh)2 e0 + h ((1 + Lh) 80 (h) + 8, (h)) ; y — y(x) in points x,,x2,...,xn, we shall receive formula

e'3 < (1 + Lh) e2 + h82 (h) < (1 + Lh)3 e0 + h ((1+Lh)2 80 (h) + (1 + Lh) 8,(h) + 82(h)) ;

, f( x-,, y.- ) - f(x , y )

ln (f ^ yi+iVf (x,, y:),

(i=0.1,...,n-1).

If we denote

e' H yi- y, I (i=0.1,•••, n),

then

e;+i H y:+i - y,+1 —

— ( yi- yt ) + h

e'n < (1 + Lh)n e0 + h ((1+Lh)n-i 80 ( h) + (1+Lh)n-2 8,( h) + ■■■ + ( h) ).

max 8 (h) —8(h).

0 < : < n-1 iW W

We shall obtain

e'n < (1 + Lh)n e0 +

+h8( h) ((1 + Lh)n-i + (1+Lh)n

■1)

f ( y+i) - f (xi, Hi) f ( xi+l, y,+i) - f (xi, y,)

<e' + h

ln (f (xi+l, y:+i Vf (x,, ) ) ln (f ^ y,+i Vf (x,, yi) ! f ( y :+i) - f (x,, y,) f ( xi+l, y,+i) - f ( x,, yi)

e'n < (1 + Lh)n e0 + - 8( h) ((1+Lh)n -1).

ln (f(xi+l, ^i+iVf(x,, ) ) ln (f(xi+l, y,+iVf(x,, y,) ^

Considering that at u > 0 the inequality eu > 1 + u, holds, then

or

or

or

e'„ < e^ + L 8(A) (eLnh -1)

y' = e2x + ex -2yex + y2, y(0) = 0.5

(42)

in interval [0;1]. Exact solution

£'„ < eLae0 + -8(h)(eLa -1).

We have on the basis of the derived inequality that the error of initial data is not piled up, that is, the method possesses computational stability.

Example.

It is required to find a numerical solution to the Cauchy's problem

y = e

1

x + 2

is obtained by introducing a new variable and by taking into account the original condition. Let us also compare the solution found by using the new method with the Euler method, and the Runge-Kutta method of fourth order. The solution to the problem is given in Table 1.

Table 1

Solution to problem (42) obtained by different methods

or

No. x ym yE yR-K y'

0 0 0.5 0.5 0.5 0.5

1 0.02 0.52515 0.525 0.52515 0.52515

2 0.04 0.55062 0.55031 0.55061 0.55061

3 0.06 0.5764 0.57594 0.5764 0.5764

4 0.08 0.60252 0.6019 0.60252 0.60252

5 0.1 0.62898 0.6282 0.62898 0.62898

6 0.12 0.6558 0.65485 0.6558 0.6558

7 0.14 0.68299 0.68187 0.68298 0.68298

8 0.16 0.71055 0.70926 0.71055 0.71055

9 0.18 0.73851 0.73704 0.7385 0.7385

10 0.2 0.76686 0.76522 0.76686 0.76686

11 0.22 0.79563 0.79381 0.79563 0.79563

12 0.24 0.82483 0.82282 0.82482 0.82482

13 0.26 0.85446 0.85227 0.85445 0.85445

14 0.28 0.88454 0.88216 0.88453 0.88453

15 0.3 0.91508 0.91251 0.91508 0.91508

16 0.32 0.9461 0.94334 0.94609 0.94609

17 0.34 0.9776 0.97464 0.9776 0.9776

18 0.36 1.00961 1.00645 1.0096 1.0096

19 0.38 1.04212 1.03876 1.04212 1.04212

20 0.4 1.07516 1.07159 1.07516 1.07516

21 0.42 1.10874 1.10496 1.10874 1.10874

22 0.44 1.14288 1.13888 1.14287 1.14287

23 0.46 1.17758 1.17335 1.17757 1.17757

24 0.48 1.21285 1.20841 1.21285 1.21285

25 0.5 1.24873 1.24406 1.24872 1.24872

26 0.52 1.28521 1.28031 1.2852 1.2852

27 0.54 1.32231 1.31717 1.32231 1.32231

28 0.56 1.36005 1.35467 1.36005 1.36005

29 0.58 1.39845 1.39282 1.39844 1.39844

30 0.6 1.43751 1.43164 1.4375 1.4375

31 0.62 1.47725 1.47113 1.47725 1.47725

32 0.64 1.5177 1.51132 1.51769 1.51769

33 0.66 1.55886 1.55221 1.55885 1.55885

34 0.68 1.60075 1.59384 1.60074 1.60074

35 0.7 1.64339 1.6362 1.64338 1.64338

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

36 0.72 1.68679 1.67933 1.68679 1.68679

37 0.74 1.73097 1.72323 1.73097 1.73097

38 0.76 1.77596 1.76793 1.77596 1.77596

39 0.78 1.82176 1.81344 1.82176 1.82176

40 0.8 1.8684 1.85977 1.8684 1.8684

41 0.82 1.91589 1.90696 1.91589 1.91589

42 0.84 1.96426 1.95501 1.96425 1.96425

43 0.86 2.01351 2.00395 2.01351 2.01351

44 0.88 2.06368 2.0538 2.06368 2.06368

45 0.9 2.11478 2.10457 2.11478 2.11478

46 0.92 2.16682 2.15628 2.16682 2.16682

47 0.94 2.21984 2.20896 2.21985 2.21985

48 0.96 2.27386 2.26262 2.27386 2.27386

49 0.98 2.32888 2.31729 2.32889 2.32889

50 1 2.38495 2.37299 2.38495 2.38495

In Table 1, y* is the exact solution, ym is the solution obtained using the new interpolation method of the minorant type (over 2 iterations), yE is the solution obtained applying the Euler's method, and yR-K is the solution obtained employing the Runge-Kutta method. Thus, the new interpolation method produces a solution that is very close to the exact solution.

As one can see, it is necessary to run a lot of iterations in order to solve the Cauchy's problem by the new method, which is the main shortcoming of the method.

6. Conclusions

1. We developed and constructed an interpolation numerical method to solve the Cauchy's problem for ordinary

first order differential equations. Underlying the method is the so-called apparatus of non-classical minorants and diagrams of Newton's functions, assigned in a tabular form. In the case of a convex function, this method produces more accurate results than the Euler's method. The method also does not require solving the systems of linear algebraic equations or superimposing additional conditions on the equations.

2. The order of accuracy is proven, as well as computational stability, convergence of the new method, and an error of approximated value is estimated by the new method. The method has a second order of accuracy, similar to the Euler's method, yet produces more accurate results in the case of a convex function.

References

1. Verzhbitskiy, V. M. Chislennyye metody (matematicheskiy analiz i obyknovennyye differentsialnyye uravneniya) [Text] / V. M. Verzhbitskiy. - Moscow: Vysshaya shkola, 2001. - 382 p.

2. Süli, E. An Introduction to Numerical Analysis [Text] / E. Süli, D. Mayers. - Cambridge University Press, 2003. - 435 p.

3. Zadachyn, V. M. Chyselni metody [Text] / V. M. Zadachyn. - Kharkiv: Vyd. KhNEU im. S. Kuznetsya, 2014. - 180 p.

4. Abdulkawi, M. Solution of Cauchy type singular integral equations of the first kind by using differential transform method [Text] / M. Abdulkawi // Applied Mathematical Modelling. - 2015. - Vol. 39 Issue 8. - P. 2107-2118. doi: 10.1016/j.apm.2014.10.003

5. Setia, A. Numerical solution of various cases of Cauchy type singular integral equation [Text] / A. Setia // Applied Mathematics and Computation. - 2014. - Vol. 230. - P. 200-207. doi: 10.1016/j.amc.2013.12.114

6. De Bonis, M. C. Numerical solution of systems of Cauchy singular integral equations with constant coefficients [Text] / M. C. De Bonis, C. Laurita // Applied Mathematics and Computation. - 2012. - Vol. 219 Issue 4. - P. 1391-1410. doi: 10.1016/j.amc.2012.08.022

7. Lukomskii, D. S. Solution of Cauchy Problem for Equation First Order Via Haar Functions [Text] / D. S. Lukomskii, S. F. Lukomskii, P. A. Terekhin // Izvestiya of Saratov University New Series, Series: Mathematics Mechanics Informatics. -2016. - Vol. 16, Issue 2. - P. 151-159. doi: 10.18500/1816-9791-2016-16-2-151-159

8. Lyashenko, M. Ya. Chysel'ni metody [Text] / M. Ya. Lyashenko, M. S. Holovan. - Kyiv: «Lybid'», 1996. - 285 p.

9. Varenych, I. I. Vyshcha matematyka: matematychnyy analiz, dyferentsialni rivnyannya [Text] / I. I. Varenych. - Kyiv: DiaSoft, 2008. - 267 p.

10. Bihun, R. R. Construction of a numerical method for finding the zeros of both smooth and nonsmooth functions [Text] / R. R. Bihun, G. G. Tsehelyk // Eastern-European Journal of Enterprise Technologies. - 2017. - Vol 2, Issue 4 (86). - P. 58-64. doi: 10.15587/1729-4061.2017.99273

11. Bihun, R. R. Formulas of minorant type to approximate computing definite integrals [Text] / R. R. Bihun, G. G. Tsehelyk // Scientific light. - 2017. - Vol. 1, Issue 7. - P. 76-79.

12. Bihun, R. R. Numerical Method for Finding All Points of Extremum of Random as Smooth and NonSmooth Functions of One Variable [Text] / R. R. Bihun, G. G. Tsehelyk // Global Journal of Science Frontier Research: F Mathematics & Decision Sciences. - 2015. -Vol. 15 Issue 2. - P. 87-93.

i Надоели баннеры? Вы всегда можете отключить рекламу.