• Ешқандай Нәтиже Табылған Жоқ

PDF BULLETIN - rmebrk.kz

N/A
N/A
Protected

Academic year: 2024

Share "PDF BULLETIN - rmebrk.kz"

Copied!
141
0
0

Толық мәтін

(1)

ISSN 2518-7929 (Print) ISSN 2663–5011(Online) Индексi 74618

Индекс 74618

MATHEMATICS Series

№ 3(99)/2020

July-August-September September, 30th, 2020

Founded in 1996 Published 4 times a year

Karaganda, 2020

BULLETIN

OF THE KARAGANDA

UNIVERSITY

(2)

Main Editor

Candidate of Physics and Mathematics sciences

N.T. Orumbayeva

Responsible secretary PhD

M.T. Kosmakova

Editorial board

A. Ashyralyev, M.A. Sadybekov, M. Otelbayev, B.R. Rakishev, U.U. Umirbaev, T. Bekjan, B. Poizat, A.A. Shkalikov, A.S. Morozov, G. Akishev, A.T. Asanova, N.A. Bokaev, M.T. Jenaliyev, A.R. Yeshkeyev, K.T. Iskakov, L.K. Kusainova, E.D. Nursultanov, M.I. Ramazanov, E.S. Smailov,

Guest editor, Professor of Mathematics, Dr. of phys.-math. sciences, Near East University, Nicosia, TRNC, Mersin 10 (Turkey);

Guest editor, Corresponding member of NAS RK, Dr. of phys.-math. sciences, IMMM, Almaty (Kazakhstan);

Academician of NAS RK, Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Academician of NAS RK, Dr. of techn. sciences, Turysov IGOM, Almaty (Kazakhstan);

Corresponding member of NAS RK, Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Professor, Xinjiang University (China);

Professor, Universite Claude Bernard Lyon-1, Villeurbanne, (France);

Corresponding member of RAS RF, Dr. of phys.-math. sciences, Lomonosov Moscow State University (Russia);

Dr. of phys.-math. sciences, Sobolev Institute of Mathematics (Russia);

Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Dr. of phys.-math. sciences, IMMM, Almaty (Kazakhstan);

Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Dr. of phys.-math. sciences, IMMM, Almaty (Kazakhstan);

Dr. of phys.-math. sciences, Buketov KU (Kazakhstan);

Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Dr. of phys.-math. sciences, Gumilyov ENU, Nur-Sultan (Kazakhstan);

Dr. of phys.-math. sciences, KB Lomonosov MSU, Nur-Sultan (Kazakhstan);

Dr. of phys.-math. sciences, Buketov KU (Kazakhstan);

Dr. of phys.-math. sciences, IAM, Karaganda (Kazakhstan)

Postal address: 28, University Str., 100024, Kazakhstan, Karaganda Теl.: (7212) 77-04-38 (add. 1026); fax: (7212) 35-63-98.

E-mail: vestnikku@gmail.com. Web-site: mathematics-vestnik.ksu.kz

Editors

Zh.Т. Nurmukhanova, S.S. Balkeyeva, Т. Kokhanover Computer layout

G.K. Kalel

Bulletin of the Karaganda University. «Mathematics» series.

ISSN 2518-7929 (Print). ISSN 2663–5011 (Online).

Proprietary: NLC «Karagandy University of the name of academician E.A. Buketov».

Registered by the Ministry of Information and Social Development of the Republic of Kazakhstan.

Rediscount certificate No. KZ43VPY00027385 dated 30.09.2020.

Signed in print 29.09.2020. Format 60×84 1/8. Offset paper. Volume 17,62 p.sh. Circulation 200 copies.

Price upon request. Order № 58.

Printed in the Publishing house of NLC «Karagandy University of the name of acad. E.A. Buketov».

38, Gogol Str., 100012, Kazakhstan, Karaganda. Tеl.: (7212) 51-38-20. E-mail: izdkargu@mail.ru

c Karagandy University of the name of acad. E.A. Buketov, 2020

(3)

CONTENTS

MATHEMATICS

Preface . . . . Ashyralyyev C., Cay A. Numerical solution to elliptic inverse problem with Neumann-type integral condition and overdetermination . . . . Karwan H.F. Jwamer and Rando R.Q. Rasul. A Comparison between the fourth order linear differential equation with its boundary value problem . . . . Mardanov M.J., Sharifov Y.A., Ismayilova K.E. Existence and uniqueness of solutions for the system ofintegro-differential equations with three-point and nonlinear integral boundary conditions . . . . Dovletov D.M.Nonlocal boundary value problem with Poissons operator on a rectangle and its difference interpretation . . . . Hincal E., Mohammed S., Kaymakamzade B.Stability analysis of an ecoepidemiological model consisting of a prey and tow competing predators with Si-disease in prey and toxicant . . . . Akat M., Kosker R., Sirma A.On the numerical schemes for Langevin-type equations . . . . Ashyralyev A., Sozen Y., Hezenci F.A remark on elliptic differential equations on manifold . Hincal E., Kaymakamzade B., Gokbulut N. Basic reproduction number and effective reproduction number for North Cyprus for fighting covid-19 . . . . Ashyralyev A., Ashyralyyev C., Zvyagin V.G.A note on well-posedness of source identification elliptic problem in a Banach space . . . . Ashyralyev A., Turk K., Agirseven D. On the stable difference scheme for the time delay telegraph equation . . . . Ashyraliyev M., Ashyralyyeva M., Ashyralyev A. A note on the hyperbolic-parabolic identification problem with involution and Dirichlet boundary condition . . . . Ashyralyev A., Erdogan A.S., Sarsenbi A.A note on the parabolic identification problem with involution and Dirichlet condition . . . . АВТОРЛАР ТУРАЛЫ МӘЛIМЕТТЕР — СВЕДЕНИЯ ОБ АВТОРАХ — INFORMATION ABOUT AUTHORS . . . .

4 5 18

26 38 55 62 75 86 96 105 120 130

140

(4)

Preface

This issue is a collection of 12 selected papers. These papers are presented at the Fifth International Conference on Analysis and Applied Mathematics (ICAAM 2020) organized by Near East University, Lefkosa (Nicosia), Mersin 10, Turkey.

The meeting was held on September 23–30, 2020 in North Cyprus, Turkey. The main organizer of the conference is Near East University, Nicosia (Lefkosa), Mersin 10, Turkey. The conference was also supported by Institute of Mathematics and Mathematical Modeling, Almaty, Kazakhstan and Analysis

&PDE Center, Ghent University, Belgium.

The conference is organized biannually. Previous conferences were held in Gumushane, Turkey in 2012; in Shymkent, Kazakhstan in 2014; in Almaty, Kazakhstan in 2016; in 2018 Lefkosa,Mersin 10, Turkey. The proceedings of ICAAM 2012, ICAAM 2014, ICAAM 2016, and ICAAM 2018 were published in AIP Conference Proceedings (American Institute of Physics) and in some rating scientific journals.

Near East University was pleased to host the fifth conference which was focused on various topics of analysis and its applications, applied mathematics and modeling. The main aim of the International Conferences on Analysis and Applied Mathematics (ICAAM) is to bring mathematicians working in the area of analysis and applied mathematics together to share new trends of applications of mathematics.

In mathematics, the developments in the field of applied mathematics open new research areas in analysis and vice versa. That is why, we planned to find the conference series to provide a forum for researches and scientists to communicate their recent developments and to present their original results in various fields of analysis and applied mathematics. This issue presents papers by authors from different countries: Azerbaijan, Iraq, Russia, Turkey, Turkmenistan, USA, Kazakhstan. Especially we are pleased with the fact that many articles are written by co-authors who work in different countries.

We are confident that such international integration provides an opportunity for a significant increase in the quality and quantity of scientific publications.

Finally, but not least, we would like to thank the Editorial board of the "Bulletin of the Karaganda University - Mathematics", who kindly provided an opportunity for the formation of this special issue.

July 2020

GUEST EDITORS:

Allaberen Ashyralyev

Department of Mathematics, Near East University, Nicosia, TRNC, Mersin 10, Turkey;

Institute of Mathematics and Mathematical Modeling, Almaty, Kazakhstan;

Peoples’ Friendship University of Russia (RUDN University), Moscow, Russia;

allaberen.ashyralyev@neu.edu.tr

Makhmud A. Sadybekov

Institute of Mathematics and Mathematical Modeling, Almaty, Kazakhstan;

sadybekov@math.kz

(5)

МАТЕМАТИКА MATHEMATICS

DOI 10.31489/2020M3/5-17 MSC 35N25, 65J22, 39A14

C. Ashyralyyev1,2, A. Cay3

1Department of Mathematical Engineering, Gumushane University, Gumushane, 29100, Turkey

2Department of Computer Technology, TAU, Ashgabat, 744000, Turkmenistan

3Netas Head Office Yenisehir Mah. Osmanli Bulvari 11 34912 Kurtkoy-Pendik, Istanbul (E-mail: charyar@gmail.com, ayselc@netrd.com.tr)

Numerical solution to elliptic inverse problem

with Neumann-type integral condition and overdetermination

In modeling various real processes, an important role is played by methods of solution source identification problem for partial differential equation. The current paper is devoted to approximate of elliptic over determined problem with integral condition for derivatives. In the beginning, inverse problem is reduced to some auxiliary nonlocal boundary value problem with integral boundary condition for derivatives. The parameter of equation is defined after solving that auxiliary nonlocal problem. The second order of accuracy difference scheme for approximately solving abstract elliptic overdetermined problem is proposed. By using operator approach existence of solution difference problem is proved. For solution of constructed difference scheme stability and coercive stability estimates are established. Later, obtained abstract results are applied to get stability estimates for solution Neumann-type overdetermined elliptic multidimensional difference problems with integral conditions. Finally, by using MATLAB program, we present numerical results for two dimensional and three dimensional test examples with short explanation on realization on computer.

Keywords: difference scheme, inverse elliptic problem, overdetermination, source identification problem, stability, coercive stability, estimate.

Introduction

Methods of solutions and theory nonlocal boundary value problems (BVPs) for differential equations have been studied by numerous authors (see [1–5, 7–12, 14–16, 18, 19] and references herein).

Let us I is identity operator andA is a selfadjoint and positive definite operator (SAPDO) in an arbitrary Hilbert space H. It is known thatA > δI for some positive numberδ, and the operator C= τ2A+

q

A+τ24A2) is also SAPDO.

Assume that given function f ∈ C1([0, T], H), elements φ, η, ζ ∈ H, number λ0 ∈ [0,1]. Denote by [0,1]τ ={ti =iτ, i= 1,· · · , N, τ N =T} the uniform grid space with step size τ >0,whereN is a fixed integer number. Let β be known scalar continuous function satisfying condition

N

X

j=1

β

tj−1 2

τ <1. (1)

(6)

C.Ashyralyyev, A. Cay

In the study [10] established well-possedness of elliptic inverse problem with Neumann-type over- determination and integral condition for obtaining a functionu∈C2([0, T], H)∩C([0, T], D(A))and an element p∈H such that

−u00(t) + Au(t) = f(t) + p, t∈(0, T), u0(0) =φ, u0(T) =

T

R

0

β(λ)u0(λ)dλ+η, u(λ0) =ζ. (2) Moreover, in [10], the stability inequalities for solution of inverse problem (2) were applied to investigate the following source identificating problem (SIP) for multi dimensional elliptic partial differential equation









−utt(t, x) −

n

P

r=1

(ar(x)uxr(t, x))xr+σu(t, x) = f(t, x) + p(x), (t, x)∈(0, T)×Ω, ut(0, x) =φ(x), ut(T, x) =

T

R

0

β(γ)uγ(γ, x)d γ+η(x), u(λ0, x) =ζ(x), x∈Ω, u(t, x) = 0, (t, x) ∈[0, T]×S.

(3)

HereΩ = (0, T)n is open cube inRnwith boundaryS, Ω = Ω∪S;ar, ζ, φ, η , f are given sufficiently smooth functions;∀x∈Ω, ar(x)≥a0 >0;σ >0,0< λ0 < T are known numbers.

We denote byR, P, and D,the corresponding operatorsR= (I+τ C)−1, P = (I−R2N)−1, D= (I+τ C)(2I +τ C)−1C−1.

Now, let us to give some lemmas that will be used in further.

Lemma 1. [8] The following estimates hold:

kRkkH→H≤M(δ) (1 +δ12τ)−k, kCRkkH→H≤ 1

kτM(δ), k≥1,kP kH→H≤M(δ), δ >0. (4) Lemma 2.

Suppose that inequality (1) is satisfied, then the operator

G2 =

−3(I −R2N) + 4 R−R2N−1

− R2−R2N−2h 3−τ β

tN−3

2

(I−R2N) +

−4−τ β tN−5

2

R−R2N−1 +

1−τ β tN−7

2

+τ β tN−3

2

R2−R2N−2 +τ β

t3 2

RN−1−RN+1 +

N−3

P

i=2

τh β

ti+1 2

−β ti−3

2

i

RN−i−RN+i

RN−1−RN+1−RN−2+RN+2h

4 +τ β tN−5

2

RN−1−RN+1 +

1−τ β tN7

2

+τ β tN−3

2

RN−2−RN+2 +

N−3

P

i=2

τh β

ti+1 2

−β ti−3

2

i

Ri−R2N−i +τ β

t3

2

R−R2N−1 +τ β

t1

2

(I−R2N) i

(5) has an inverse G−12 and its norm is bounded, i.e.

kG−12 kH→H≤M(δ). (6) In the paper [8] , for givenv0 andvN,the solution of difference scheme

−τ−2(vi+1−2vi+vi−1) +Avi=fi, 1≤i≤N −1 (7) was represented by formula

vi =P

Ri−R2N−i

v0+ RN−i−RN+i vN

−P RN−i−RN+i D

×

N−1

P

j=1

RN−j−RN+j

fjτ +D

N−1

P

j=1

R|i−j|−Ri+j

fjτ , 1≤i≤N −1. (8)

(7)

Numerical solution to elliptic ...

Letα∈(0,1)is a given number. Introduce notations forCτ(H), Cτα(H),andCτα,α(H),the Banach spaces ofH-valued grid functions wτ ={wk}N−1k=1 with the corresponding norms,

kwτkC

τ(H) = max

1≤k≤N−1kwkkH, kwτkCα

τ(H) = sup

1≤k<k+n≤N−1

(nτ)−αkwk+n−wkkH+kwτkC

τ(H), kwτkCα,α

τ (H)=kwτkC

τ(H)+ sup

1≤k<k+n≤N−1

(1−kτ)α(nτ)−α(kτ +nτ)αkwk+n−wkkH.

In the current study, we construct the second order accuracy difference scheme (ADS) for approximately solution of inverse problem (2) and study well-posedness of difference problem. Then, we discuss the second order ADS for SIP (3).

The second order of ADS for SIP (3) Now, we study second order of ADS









−τ−2(uk+1−2uk+uk−1) +Auk =fk+p, fk=f(tk),1≤k≤N−1,

−3u0+ 4u1−u2 = 2τ φ,3uN−4uN−1+uN−2 =

N−1

P

i=1

τ β ti−1

2

(ui+1−ui−1) + 2τ η, ul+µ(ul+1−ul) =ζ

µ= λτ0 −l

(9)

for approximate solution inverse problem (2).

Theorem 1. Let us φ, η, ζ ∈D(A),and fτ ∈Cτ(H) and inequality (1) is satisfied. Then, solution

{uk}Nk=1−1, p

of difference problem (9) exists inCτ(H)×Hand the next stability estimates for solution

{uk}N−1k=1

Cτ(H) ≤ M(δ)

kφkH+kζkH +kηkH+kfτkC

τ(H)

, (10)

A−1p

H ≤ M(δ)

kφkH+kζkH +kηkH+kfτkC

τ(H)

(11) are fulfilled.

Proof. Firstly, by using

uk=vk+A−1p, (12)

we get auxiliary difference problem for unknowns {vk}Nk=0 :













−τ−2(vk+1−2vk+vk−1) +Avk=fk , 1≤k≤N −1,

−3v0+ 4v1−v2 = 2τ φ,

3−τ β

tN−3

2

vN +

−4−τ β

tN−5

2

vN−1

+

1−τ β

tN7

2

+τ β

tN−3

2

vN−2+

N−3

P

i=2

τ h

β

ti+1

2

−β

ti−3

2

i vi

+τ β t3

2

v1 +τ β t1

2

v0 = 2τ η.

(13)

We seek solution of (13) by (8). By using (8), from first condition of difference problem (13), we get equation

−3(I−R2N) + 4 R−R2N−1

− R2−R2N−2 v0

+

4 RN−1−RN+1

− RN−2−RN+2

vN =F1, (14)

for unknowns v0 and vN,where

F1= 2τ(I−R2N)φ+ 4 RN−1−RN+1 D

N−1

P

j=1

RN−j−RN+j

fjτ −4(I−R2N)D

×

N−1

P

j=1

R|1−j|−R1+j

fjτ − RN−2−RN+2 D

N−1

P

j=1

RN−j−RN+j fjτ +(I−R2N)D

N−1

P

j=1

R|2−j|−R2+j fjτ.

(8)

C.Ashyralyyev, A. Cay

From integral condition follows the next equation

3−τ β tN−3

2

(I −R2N)vN +

−4−τ β tN−5

2

RN−1−RN+1

v0+ R−R2N−1 vN +

1−τ β tN7

2

+τ β tN−3

2

RN−2−RN+2

v0+ R2−R2N−2 vN +

N−3

P

i=2

τh β

ti+1 2

−β ti−3

2

i

Ri−R2N−i

v0+ RN−i−RN+i vN +τ β

t3

2

R−R2N−1

v0+ RN−1−RN+1 vN

+τ β

t1

2

(I−R2N)v0 =F2

(15) for unknowns v0 and vN,where

F2 =

−4−τ β tN−5

2

"

R−R2N−1 D

N−1

P

j=1

RN−j −RN+j

fjτ −(I −R2N)D

×N

−1

P

j=1

R|N−1−j|−RN−1+j fjτ

# +

1−τ β tN−7

2

+τ β tN3

2

×

"

R2−R2N−2 D

N−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|N−2−j|−RN−2+j fjτ

#

N

−3

P

i=2

τh β

ti+1 2

−β ti−3

2

i

"

RN−i−RN+i D

N−1

P

j=1

RN−j−RN+j fjτ

−(I−R2N)D

N−1

P

j=1

R|i−j|−Ri+j fjτ

#

−τ β t3

2

RN−1−RN+1 D

×

N−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|1−j|−R1+j

fjτ + 2τ(I −R2N

# .

Thus, determinant operatorG2of linear system equatıon (14), (15) has bounded inverseG−12 .Therefore solution of linear system equatıon (14), (15) is defined by

v0 =G−12 nh

3−τ β

tN−3

2

(I−R2N) + −4−τ β tN−2τ2

R−R2N−1 +

1−τ β tN−7

2

+τ β tN3

2

R2−R2N−2 +

N−3

P

i=2

τh β

ti+1

2

−β ti−3

2

i

RN−i−RN+i

+τ β t3

2

RN−1−RN+1

×

"

2τ(I−R2N)φ+ 4 RN−1−RN+1 D

N−1

P

j=1

RN−j−RN+j

fjτ−4(I−R2N)D

×

N−1

P

j=1

R|1−j|−R1+j

fjτ− RN−2−RN+2 D

N−1

P

j=1

RN−j−RN+j fjτ +(I−R2N)D

N−1

P

j=1

R|2−j|−R2+j fjτ

#

− RN−1−RN+1−RN−2+RN+2

× (

2τ(I−R2N)η+

−4−τ β tN−5

2

"

R−R2N−1 D

N−1

P

j=1

RN−j−RN+j fjτ

× −(I−R2N)D

N−1

P

j=1

R|N−1−j|−RN−1+j fjτ

# +

1−τ β tN−7

2

+τ β tN−3

2

×

"

R2−R2N−2 D

N−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|N−2−j|−RN−2+j fjτ

#

(9)

Numerical solution to elliptic ...

N−3

P

i=2

τ

α ti+1τ2

−α ti−1τ2

"

RN−i−RN+i D

N−1

P

j=1

RN−j−RN+j fjτ

−(I −R2N)D

N−1

P

j=1

R|i−j|−Ri+j fjτ

#

−τ α t2τ2 RN−1−RN+1 D

×N

−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|1−j|−R1+j fjτ

##)) ,

(16)

and

vN =G−12

−3(I−R2N) + 4 R−R2N−1

− R2−R2N−2

2τ(I−R2N)η + −4−τ β tN−2τ2

"

R−R2N−1 D

N−1

P

j=1

RN−j−RN+j fjτ

−(I−R2N)D

N−1

P

j=1

R|N−1−j|−RN−1+j fjτ

#

+ 1−τ β tN−3τ2

+τ β tN−1τ2

×

"

R2−R2N−2 D

N−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|N−2−j|−RN−2+j fjτ

#

N−3

P

i=2

τ h

β

ti+1

2

−β

ti−3

2

i

"

RN−i−RN+i D

N−1

P

j=1

RN−j−RN+j

fjτ−(I−R2N)D

N−1

P

j=1

R|i−j|−Ri+j fjτ

#

−τ β t3

2

"

RN−1−RN+1 D

N−1

P

j=1

RN−j −RN+j

fjτ −(I−R2N)D

N−1

P

j=1

R|1−j|−R1+j fjτ

#

−h

4 +τ β tN−5

2

RN−1−RN+1 +

1−τ β

tN7

2

+τ β

tN−3

2

RN−2−RN+2 +τ β

t1

2

(I−R2N) +

N−3

P

i=2

τ h

β

ti+1

2

−β

ti−3

2

i

Ri−R2N−i +τ β

t3

2

R−R2N−1

×

"

2τ(I−R2N)φ+ 4 RN−1−RN+1 D

N−1

P

j=1

RN−j −RN+j fjτ

−4(I−R2N)D

N−1

P

j=1

R|1−j|−R1+j fjτ

− RN−2−RN+2 D

N−1

P

j=1

RN−j−RN+j

fjτ + (I−R2N)D

N−1

P

j=1

R|2−j|−R2+j fjτ

#) .

(17) Thus solution of difference problem (13) exists and it is defined by (8) with the corresponding v0 and vN via (16) and (17). From (8), (16),(17), estimates (4), (6), it follows that for solution of difference problem (13) stability estimates

{vk}N−1k=1 Cτ(H)

≤M(δ)

kφkH +kζkH+kηkH +kfτkC

τ(H)

, (18)

{Avk}Nk=1−1

Cα,α τ (H)

+

nv

k+1−2vk+vk−1

τ2

oN−1 k=1

Cα,α τ (H)

≤M(δ) 1

α(1−α)kfτk

Cα,α

τ (H)+kAζkH+kAφkH +kAηkH .

(19)

are fulfilled. (12) and estimates (18) permit us to get estimates estimates (11) (10) and (19).

(10)

C.Ashyralyyev, A. Cay

Theorem 2. Let us fτ ∈ Cτα,α(H), and φ, ζ, η ∈ D(A) and inequality (1) is satisfied. Then, for solution

{uk}N−1k=1 , p

of difference problem (9) the coercive stability inequality

nu

k+1−2uk+uk−1

τ2

oN−1 k=1

Cτα,α(H)

+

{Auk}N−1k=1

Cα,α τ (H)

+kpkH

≤M(δ) 1

α(1−α)kfτk

Cα,α

τ (H) +kAζkH +kAφkH +kAηkH (20)

is valid.

The proof of inequality (20) is based on formulas (8), (12), (16), (17), and (19).

Approximation of (3)

Denote by Ωeh =

x= (h1m1, ..., hnmn);m= (m1, ..., mn), mi= 0, Mi, hiMi = 1, i= 1, n , Ωh=Ωeh∩Ω, Sh=Ωeh∩S

and byAxh difference operator

Axhuh(x) =−

n

X

i=1

ai(x)uhxi(x)

xi,ji

+σuh(x)

acting in the space of grid functionsuh(x),satisfying boundary conditionuh(x) = 0 for all x∈Sh. In the beginning, by using approximation in variable x and later by approximation in variablet, one can get the following difference scheme for approximately solution of SIP (3):

−τ−2 uhk+1(x)−2uhk(x) +uhk−1(x)

+Auhk(x) =fkh(x) +ph(x), 1≤k≤N −1, x∈Ωh

−3uh0(x) + 4uh1(x)−uh2(x) =τ φh(x), uhl(x) +µ(uhl+1(x)−uhl(x)) =ζh(x) 3uhN(x)−4uhN−1(x) +uhN−2(x) =

N−1

P

i=1

τ α tiτ2

uhi+1(x)−uhi(x)

+ 2τ ηh(x), x∈Ωeh. .

(21)

Let L2h =L2(Ωeh) and W2h2 =W22(Ωeh), the Banach spaces of the grid functions uh(x) ={u(h1m1,· · ·, hnmn)} defined on Ωeh,equipped with the corresponding norms

uh

L2h= (P

x∈eh|uh(x)|2h1· · ·hn)1/2,

uh

W2h2 = uh

L2h+ +(P

x∈eh

Pn i=1

(uh(x))xixi, mi

2h1· · ·hn)1/2.

Theorem 3. Assume that (1) is valid, fτ ∈ Cτα,α(L2h), and φh, ηh, ζh ∈ D(Axh)∩L2h. Then, the solution of difference problem (21) exists and for solution the stability estimates hold:

uhk N1 −1 C

τ(L2h)≤M(δ) [

φh

L2h+ ηh

L2h+ ζh

L2h+kfτkC

τ(L2h)

,

ph

L2h≤M(δ)

ζh

W2h2 + ηh

W2h2 + φh

W2h2 +α(1−α)1 kfτkCα,α

τ (L2h)

.

Theorem 4.Assume that (1) is true, fτ ∈Cτα,α(W2h2 ), andφh, ηh, ζh ∈D(Axh)∩W2h2 . Then, for the solution of difference problem (21) the coercive stability estimate obeys

uhk+1−2uhk+ukk−1

τ2 )

N−1 1

Cτ(L2h)

+

uhk N−11

Cτ(W2h2 )+ ph

L2h

≤M(δ) [

ζh

W2h2 + ηh

W2h2 + φh

W2h2 +α(1−α)1 kfτkCα,α

τ (W2h2 )

.

(11)

Numerical solution to elliptic ...

The proofs of Theorems 3 and 4 are based on the symmetry property of the operatorAxh in the Hilbert spaceL2h and the corresponding theorem in [20] on the coercivity stability inequality for the solution of the elliptic difference problem inL2h with first kind boundary condition.

Test examples

In the present section, we illustrate computed results for twodimensional and threedimentional examples of inverse elliptic problem with Neumann-type overdetermination and integral condition. All computed results are carried out by using MATLAB.

2D example Notice that pair functions (p(x), u(t, x)) = ( π2+ 1

sin(πx), e−t +t+ 1

sin(πx)) is exact solution of the following 2D overdetermined elliptic problem with integral boundary condition:





−utt(t, x)−uxx (t, x) +u(t, x) =f(t, x) +p(x), t, x∈(0,1), ut(0, x) = 0, u(0.3, x) =ζ(x), ut(1, x) =

R1 0

e−λuλ(λ)dλ +η(x), x∈[0,1], u(t,0) = 0, u(t,1) = 0, t∈[0,1],

(22)

where

f(t, x) =

−e−t+ π2+ 1

e−t+t sin(πx), ζ(x) = e−0.3 + 1.3

sin(πx) η(x) =1

212e−2

sin(πx).

The notation [0,1]τ×[0,1]h means the set of grid points

[0,1]τ×[0,1]h ={(ti, xn) :ti=iτ, i= 0, N , xn=nh, n= 0, M}, which depends on the small parameters τ and h such thatN τ = 1, M h= 1.Let us

l0 =

0.3τ−1

, µ0 = 0.3τ−1−l0n= 0, ηn=η(xn), ζn=ζ(xn), n= 0, M; fnk=f(tk, xn), k= 0, N , n= 0, M .

To approximately soving (22), we use algorithm which contains three stages. Firstly, we find approximately solution of auxiliary NBVP









τ−2 vnk+1−2vnk+vnk−1

+h−2 vkn+1−2vkn+vn−1k

−vnk=−f(tk, xn), k= 1, N −1, n= 1, M −1,

vk0 =vkM = 0, k= 0, N , −3v0n+ 4vn1−vn2 = 0, 3vnN −4vnN−1+vnN−2=

N−1

P

j=1 τ

2e(tjτ

2)

vj+1n −vnj−1+vnj −vj−2n

+ 2τ ηn, n= 0, M .

(23)

Secondly, we find pn.It is caried out by

pn=−h12[(ζn+1−(µ0vn+1l0+1−(µ0−1)vln+10 ))−2(ζn−(µ0vln0+1−(µ0−1)vln0))

+ (ζn−1−(µ0vn−1l0+1−(µ0−1)vln−10 ))] +ζn−(µ0vln0+1−(µ0−1)vnl0), , n= 1, M −1.

Difference problem (23) can be rewritten in the matrix form

Avn+1+Bvn+Cvn−1=Ig(n), n= 1, M −1, v0=−→

0, vM =−→

0. (24)

(12)

C.Ashyralyyev, A. Cay

Here,A, B, C, Iare(N+1)×(N+1)square matrices, andI is identity matrix,vs, s=n−1, n, n+1, g(n) are column matrices with(N+ 1) rows,vs=

vs0 . . . vsN t

.Denote by a= 1

h2, c= 1

h2, q=− 2 h2 − 2

τ2 −1, r= 1 τ2. Then,

An=diag(0, a, a, . . . , a,0), Cn=An, gk(n)=−f(tk, xn), k= 1, N −1, n= 1, M −1,

bi,i =q, bi−1,i=r, bi,i−1 =r, i= 2, N , b1,1 =−3, b1,2 = 4, b1,3 =−1, bN+1,N+1= 2τ

e

tN−3 2

4 +e−tN−32

−3, bN+1,N = 2τ

etN−5

4 +e

−tN−3 2

4 −e−tN−12

+ 4, bN+1,N−1 = τ e

−tN−7 2

2 +τ e

−tN−5 2

2τ e

−tN−3 2

2 −1, bN+1,1 = 2τ

e

−t3 2

4 −e−t12

, bN+1,2 = 2τ

e

−t3 2

4e

−t5 2

4 +e−t12

,

bN+1,3 = τ e

−t3 2

2τ e

−t5 2

2τ e

−t7 2

2 , bN+1,j = τ2

e−tj−32 +e−tj−12 −e−tj+ 12 −e−tj+ 32

, j= 4, ..., N −2;

bi,j = 0, for otheriand j;g0n= 2τ φn, gNn = 2τ ηn, n= 1, M−1.

To solve (24), we use modified Gauss elimination method.

Thirdly, we define

ukn by ukn=vnkn− µ0vln0+1−(µ0−1)vln0 .

Errors are presented in Tables 1-3 for second order ADS in case N=M=10,20, 40, 80, 160 and 320.

It can be seen from Tables 1-3 when N, M are increased two times that errors are decreased with approximately ratio 14.

T a b l e 1 Test example (22) - errorv

DS\(N, M) (10,10) (20,20) (40,40) (80,80) (160,160) (320,320) 2nd order of ADS 6.29×10−3 1.57×10−3 3.93×10−4 9.84×10−5 2.46×10−5 6.15×10−6

T a b l e 2 Test example (22) - erroru

DS \ (N, M) (10,10) (20,20) (40,40) (80,80) (160,160) (320,320) 2nd order of ADS 3.13×10−4 7.95×10−5 2.02×10−5 5.10×10−6 1.28×10−6 3.22×10−7

T a b l e 3 Test example (22) - errorp

Appr. \ (N, M) (10,10) (20,20) (40,40) (80,80) (160,160) (320,320) 2nd order 5.03×10−3 1.28×10−3 3.21×10−4 8.06×10−5 2.02×10−5 5.05×10−6

(13)

Numerical solution to elliptic ...

3D example

Now, consider the three dimensional inverse elliptic problem with integral condition









−utt(t, x, y)−uxx(t, x, y)−uyy(t, x, y) +u(t, x, y) =f(t, x, y) +p(x, y), x, y, t∈(0,1), u(t,0, y) =u(t,1, y) = 0, y, t∈[0,1], u(t, x,0) =u(t, x,1) = 0, x, t∈[0,1],

ut(0, x, y) =φ(x, y), u(0.6, x, y) =ζ(x, y), ut(1, x, y)−

1

R

0

e−λuλ(λ, x, y)dλ =η(x, y), x, y∈[0,1],

(25)

where

f(t, x, y) = 2π2e−tq(x, y), φ(x, y) =−q(x, y), η(x, y) =

−e−1+13 e−0.6+e−1.2

q(x, y), ζ(x, y) =

e35 + 1

q(x, y),q(x, y) = sin(πx) sin(πy) It is clear that pair funcions p(x, y) = 2π2+ 1

q(x, y) and u(t, x, y) = e−t+ 1

q(x, y) is exact solution of (25).

Denote by [0,1]τ×[0,1]h ×[0,1]h set of grid points depending on the small parameters τ andh [0,1]τ×[0,1]2h ={(ti, xn, ym) :ti =iτ, i= 0, N , xn=nh, n= 0, M ,

ym=mh, m= 0, M , τ N = 1, hM = 1}.

Let us l0 =

0.3τ−1

, µ0= 0.3τ−1−l0 , φm,n =φ(xn, ym), ηm,n =η(xn, ym), ζm,n =ξ(xn, ym), n= 0, M , m= 0, M;fm,ni =f(ti, xn, ym), i= 0, N , n= 0, M, m= 0, M .

Firstly, difference scheme for approximate solution of NBVP can be written in the following form:

























−τ−2 vk+1m,n −2vm,nk +vk−1m,n

−h−2 vm,n+1k −2vkm,n+vkm,n−1

−h−2 vm+1,nk −2vkm,n+vm−1,nk

+vkm,n =fm,nk , k= 1, N −1, n= 1, M −1, m= 1, M −1,

v0,nk =vM,nk =vkm,n=vm,Mk = 0, k = 0,· · · , N, n= 1, M−1, m= 1, M −1,

−3v0m,n+ 4vm,n1 −vm,n2 = 2τ φm,n, 3vNm,n−4vm,nN−1+vm,nN−2

=

N−1

P

j=1 τ

2e(tjτ

2)

vj+1m,n−vm,nj−1+vm,nj −vm,nj−2

+ 2τ ηm,n, n= 1, M−1, n= 1, M −1.

(26)

Secondly, calculation ofpn n= 1, M −1, m= 1, M −1

is caried out by pm,n =−h12

nh

ζm,n+1

µ0vm,n+1l0+1 −(µ0−1)vm,n+1l0 i

−2

ζm,n− µ0vm,nl0+1−(µ0−1)vm,nl0 +h

ζm,n−1

µ0vm,n−1l0+1 −(µ0−1)vm,n−1l0 io

h12 nh

ζm+1,n

µ0vm+1,nl0+1 −(µ0−1)vm+1,nl0 i

−2

ζm,n− µ0vm,nl0+1−(µ0−1)vm,nl0 +h

ζm−1,n

µ0vlm−1,n0+1 −(µ0−1)vlm−1,n0 io . Thirdly, we calculate

ukn by

ukm,n=vkm,nm,n

µ0vlm,n0+1−(µ0−1)vm,nl0 .

Difference problem (26) can be rewritten in the matrix form (24). In this case, gn is a column matrix with(N + 1)(M+ 1) elements,A, B, C, I are square matrices with(N+ 1)(M+ 1)rows and columns, and I is the identity matrix, vs is column matrix with (N + 1)(M+ 1)elements such that

vs =

v0,s0 · · · v0,sN v01,s · · · vN1,s · · · v0M,s · · · vM,sN t

, s=n−1, n, n+ 1.

(14)

C.Ashyralyyev, A. Cay

Denote by

a= 1

h2, q= 1 + 2 τ2 + 4

h2, r= 1 τ2. Then,

A=C=

O O · · · O O O E · · · O O

· · · . .. · · · · O O · · · E O O O · · · O O

, B =

Q O · · · O O O D · · · O O

· · · . .. · · · · O O · · · D O O O · · · O Q

 ,

E=diag(0, a, a, . . . , a,0), Q=I(N+1)×(N+1), O=O(N+1)×(N+1), gkm,n=−f(tk, xn, ym), k= 1, N −1, n= 1, M−1, m= 1, M −1, di,i=q, di−1,i=r, di,i−1 =r, i= 2, N , d1,1 =−3, d<

Сурет

Figure 1. The convergence of expected value of the squared position and velocity functions
Figure 2. Error measured as the difference between the exact solution and the numerical solution at particular points in time between 0 and 1
Figure 1. The basic reproduction numbers for pandemics comparing to SARS-Cov-2 (Covid-19)
Figure 2. R 0 and R t values from March 11, 2020 to May 15, 2020
+2

Ақпарат көздері

СӘЙКЕС КЕЛЕТІН ҚҰЖАТТАР

Loaded differential equations arise naturally in the study of nonlinear equations, particles transport equations and optimal control problems, in the numerical solution

Loaded differential equations arise naturally in the study of nonlinear equations, particles transport equations and optimal control problems, in the numerical solution

Keywords: nonlinear Fredholm integro-differential equation, boundary value problem, special Cauchy problem, iterative process, isolated solution, algorithm, Dzhumabaev

In the present paper we study of a questions the existence and uniqueness of solutions to family of multipoint- integral boundary value problems for the third-order

Key words: nonlinear boundary value problem for the Fredholm integro-differential equation, special Cauchy problem, Dzhumabaev parameterization

Asymptotical representation of singularly perturbed boundary value problems for integro-differential equations// News of the National Academy of Sciences of the Republic

We prove the existence and uniqueness of the solution of mixed value problem for one-dimensional nonlinear nonstationary Boltzmann moment system of equations in first and

Dzhumabaev, the necessary and sufficient conditions for the unique solvability of boundary value problems for integro-differential equations are obtained using the basic