\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2017 (2017), No. 170, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2017 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2017/170\hfil Solvability of a nonlocal problem]
{Solvability of a nonlocal problem for a hyperbolic
 equation with integral conditions}

\author[A. T. Assanova \hfil EJDE-2017/170\hfilneg]
{Anar T. Assanova}

\address{Anar T. Assanova \newline
Department of Differential Equations,
Institute of Mathematics and Mathematical Modeling,
125, Pushkin str., 050010 Almaty,  Kazakhstan}
\email{assanova@math.kz, anarasanova@list.ru}

\dedicatory{Communicated by Ludmila S. Pulkina}

\thanks{Submitted March 17, 2017. Published July 6, 2017.}
\subjclass[2010]{35L51, 35L53, 35R30, 34B10}
\keywords{Hyperbolic equation; nonlocal problem;
integral condition; algorithm;
\hfill\break\indent  approximate solution}

\begin{abstract}
 We study a nonlocal problem with integral conditions for a hyperbolic equation
 two independent variables. By introducing additional functional parameters,
 we investigated the solvability and construction of approximate solutions.
 The original problem is reduced to an equivalent problem consisting of
 the Goursat problems for a hyperbolic equation with parameters and the
 boundary value problem with integral condition for the ordinary differential
 equations with respect to the  parameters. Based on the algorithms for finding
 solutions to the equivalent problem,  we propose algorithms for finding
 the approximate solutions, and prove their convergence.
 Coefficient criteria for the unique solvability of nonlocal problem with
 integral conditions for hyperbolic equation with mixed derivative are
 also established.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\allowdisplaybreaks

 \section{Introduction}

On the domain $ \Omega =[0,T]\times [0,\omega]$, we consider the nonlocal 
problem for hyperbolic equation with integral conditions
\begin{gather}
\frac{\partial ^2 u}{\partial t\partial x} 
= A(t,x)\frac{\partial u} {\partial x} + B(t,x)
 \frac{\partial u} {\partial t} + C(t,x)u + f(t,x),  \label{e1} \\
 \int ^a_0 K(t,\xi)u(t,\xi)d\xi = \psi (t), \quad t \in [0, T], \label{e2}\\
 \int ^b_0 M(\tau,x) u(\tau, x)d\tau =\varphi (x), \quad x \in [0, \omega], 
\label{e3}
\end{gather}
where  $u (t,x)$ is unknown function, the functions $ A(t,x)$, $B(t,x)$,
 $C(t,x)$, and $f(t,x)$ are continuous on $\Omega$, the functions  
$K(t,x)$ and $ \psi (t) $ are continuously differentiable by $t$ on
$\Omega$ and $[0, T]$, respectively, the functions $M(t,x)$ and $
\varphi (x) $ are continuously differentiable by  $x$ on $\Omega$
and $[0, \omega]$, respectively, $0 < a \leq \omega$, 
$ 0 < b \leq T$. The compatibility condition is given below.

Let $C(\Omega,\mathbb{R})$ be a space of functions  $u : \Omega \to R$
continuous on  $\Omega$ with norm  $ \|u\|_0 = \max_{(t,x)\in \Omega}|u(t,x)|$.

A function  $u(t,x) \in C(\Omega,\mathbb{R})$, having the partial
derivatives $\frac{\partial u(t,x)}{\partial x} \in C(\Omega,\mathbb{R})$,
$\frac{\partial u(t,x)}{\partial t} \in C(\Omega,\mathbb{R})$, and
$\frac{\partial ^2 u(t,x)}{\partial t \partial x} \in C(\Omega,\mathbb{R})$
is called a classical solution to problem \eqref{e1}--\eqref{e3}, if it
satisfies equation \eqref{e1} and integral conditions 
 \eqref{e2}, \eqref{e3}.

Mathematical modelling of various physical processes often leads
to the nonlocal problems for hyperbolic equations. Problems with
integral conditions arise while researching the processes of
heat distribution, plasma physics, clean technology of silicon
ores, moisture transfer in capillary-porous media, etc. 
\cite{b1,b2,g1,k1,k2,k3,n1,n2,n3,p1,s1,s2,s3,t1}.
Some classes of nonlocal boundary value problems
with integral conditions for hyperbolic equations are studied in
[6, 8-28]. Solvability conditions for the considered problems are
obtained in the different terms. Problem \eqref{e1}--\eqref{e3} 
for $K(t,x)=M(t,x)$ and $K(t,x)=M(t,x)=1$ are studied in 
\cite{g1,k1,p2,p3,p4,p5,p6}.
Under the assumptions of continuous differentiability of
the equation coefficients, the conditions for the unique
solvability of that problem have been obtained. For $K(t,x)=
K(x)$, $M(t,x)= M(t)$, the problem \eqref{e1}--\eqref{e3} for the system of
hyperbolic equations is studied in \cite{z1} by contractive mapping
principle.

For $K(t,x)= M(t,x)$, problem \eqref{e1}--\eqref{e3} for the system of
hyperbolic equations is studied in \cite{a7}. The  unique solvability
conditions for this problem are established in the terms of
initial data. Nonlocal problems with general integral conditions
for hyperbolic equations belong to the field of little studied
problems of mathematical physics. This formulation of problem is
considered for the first time.

The aim of this work  is  to construct  algorithms for finding a
solution to  problem \eqref{e1}--\eqref{e3} and establish  conditions for the
existence and uniqueness of classical solution to problem
\eqref{e1}--\eqref{e3}.

In  Section 2,  a scheme of  method used \cite{a4,a5,a6} is provided. By
introduction of new unknown functions being a linear combination
of solutions values on  characteristics,   the problem \eqref{e1}--\eqref{e3} is
reduced to an equivalent problem consisting of the Goursat problem
for hyperbolic equations with functional parameters and the
boundary value problems with integral conditi\-ons for ordinary
differential equations with respect to the  parameters entered.
Algorithm for finding the approximate solution to the investigated
problem is constructed. The algorithm consists of two parts: in
the first part, we solve two boundary value problems with integral
condition for ordinary differential equations, and in the second
part, we solve the Goursat problem for hyperbolic equation with
parameters. Boundary value problems with integral condition for
ordinary differential equations are intensively studied in recent
years, and they find  numerous applications in the applied
problems \cite{a1,a2,a3}. In  Section 3,  the conditions for the existence
of unique solution to the boundary value problems with integral
condition for ordinary differential equations are presented. In
 Section 4, the convergence of  algorithm is proved, and the
conditions for the unique solvability of  problem \eqref{e1}--\eqref{e3} are
given in the terms of initial data.

\section{Method's scheme and algorithm}

Notation: $ \mu (t) = u(t,0)- \frac{1}{2}u(0,0)$,
 $\lambda (x) = u(0,x)- \frac{1}{2}u(0,0)$, $\widetilde{u}(t,x)$,
 where  $\widetilde{u}(t,x)$ is a new unknown function.
We make a following replacement of desired function $u(t,x)$ in  problem  
\eqref{e1}--\eqref{e3}:
$ u(t,x) = \widetilde{u}(t,x) + \mu (t) + \lambda (x)$ and transit to the problem
\begin{gather}
\begin{aligned}
\frac{\partial ^2 \widetilde{u}}{\partial t \partial x} 
&= A(t,x) \frac{\partial \widetilde{u}} {\partial x} +
B(t,x)\frac{\partial \widetilde{u}} {\partial t} +
C(t,x)\widetilde{u} +   A(t,x)\dot{\lambda} (x) \\
&\quad  +B(t,x)
\dot{\mu}(t)+  C(t,x)\lambda (x) + C(t,x) \mu(t) + f(t,x),
\end{aligned}  \label{e1.1} \\
 \widetilde{u}(t,0) = 0, \quad t \in [0,T], \label{e1.2} \\
 \widetilde{u}(0,x) = 0, \quad  x \in [0, \omega], \label{e1.3} \\
\begin{gathered}
\int ^{a}_0 K(t,\xi)d\xi  \mu(t) 
+\int ^{a}_0 K(t,\xi) \widetilde{u}(t,\xi)d\xi  
+ \int ^{a}_0 K(t,\xi)\lambda(\xi)d\xi  = \psi (t), \\
 t\in [0,T], 
\end{gathered}\label{e1.4} \\
\begin{gathered}
\int ^b_0 M(\tau,x)d\tau  \lambda(x) 
+ \int ^b_0 M(\tau,x) \widetilde{u}(\tau,x)d\tau
+\int ^b_0 M(\tau,x)  \mu(\tau)d\tau =\varphi(x), \\
 x \in [0, \omega].
\end{gathered} \label{e1.5}
\end{gather}
A triplet of functions $(\widetilde{u}(t,x), \mu (t), \lambda
(x))$, satisfying the hyperbolic equation \eqref{e1.1}, the conditions
on characteristics  \eqref{e1.2}, \eqref{e1.3}, and the functional relations
\eqref{e1.4} and \eqref{e1.5} at $\mu(0)=\lambda(0)$, will be called a solution
to problem \eqref{e1.1}--\eqref{e1.5} if the function $\widetilde{u}(t,x) \in
C(\Omega, R)$ has the partial derivatives 
\[
\frac{\partial \widetilde{u}(t,x)}{\partial x} \in C(\Omega, R),\quad
\frac{\partial \widetilde{u}(t,x)}{\partial t} \in C(\Omega, R), \quad
\frac{\partial^2 \widetilde{u}(t,x)}{\partial t \partial x} \in C(\Omega, R),
\]
the functions  $\mu (t)$ and $\lambda (x)$
are continuously differentiable on $[0,T]$ and $[0,\omega]$,
respectively.

Relation $\mu(0)=\lambda(0)$ is a compatibility condition of data.

Problem \eqref{e1.1}--\eqref{e1.5} is equivalent to problem 
\eqref{e1}--\eqref{e3}. If the function
 $ u^{\ast}(t,x)$ is a solution to problem \eqref{e1}--\eqref{e3}, 
then the triplet of functions
 $(\widetilde{u}^{\ast}(t,x)$,  $ \mu^{\ast} (t)$, 
$\lambda^{\ast} (x))$, where 
$\widetilde{u}^{\ast}(t,x) =u^{\ast}(t,x) - \mu^{\ast} (t)- \lambda^{\ast} (x)$, 
$\mu^{\ast}(t) = u^{\ast}(t,0)- \frac{1}{2}u^{\ast}(0,0)$, 
$\lambda^{\ast}(x) = u^{\ast}(0,x) - \frac{1}{2}u^{\ast}(0,0)$, is a solution to
problem \eqref{e1.1}--\eqref{e1.5}. The converse is also true. If the triplet of
functions  $( \widetilde{u}^{\ast\ast}(t,x), \mu^{\ast \ast} (t)$,
$ \lambda^{\ast\ast} (x))$ is a solution to problem \eqref{e1.1}--\eqref{e1.5},
then the function $ u^{\ast\ast}(t,x)$ defined by the equality
\[
u^{\ast\ast}(t,x) = \widetilde{u}^{\ast\ast}(t,x)  +
\mu^{\ast\ast} (t)+ \lambda^{\ast\ast} (x),\]
where $u^{\ast\ast}(t,0) - \frac{1}{2}u^{\ast\ast}(0,0) = \mu^{\ast\ast}
(t)$,  $u^{\ast\ast}(0,x) - \frac{1}{2}u^{\ast\ast}(0,0)=
\lambda^{\ast\ast} (x)$ is a solution to problem \eqref{e1}--\eqref{e3}.

At fixed $ \mu (t)$, $\lambda (x)$, problem  \eqref{e1.1}--\eqref{e1.3} is the
Goursat  problem  with respect to the function
$\widetilde{u}(t,x)$ on the domain $\Omega$. Relations  \eqref{e1.4} and
\eqref{e1.5} allow us to determine the unknown parameters $ \mu (t)$,
$\lambda (x)$, where the functions  $ \mu (t)$, $\lambda (x)$
satisfy condition  $ \mu (0)= \lambda (0)$.

By conditions  \eqref{e1.2}, \eqref{e1.3}, relations  \eqref{e1.4} at $t=0$
and \eqref{e1.5} at $x=0$ yield
\begin{gather}
\int ^{a}_0 K(0,\xi) d\xi  \mu (0) +  \int ^{a}_0 K(0,\xi)\lambda(\xi)d\xi
 = \psi (0), \label{e1.6} \\
\int ^b_0 M(\tau,0)d\tau \lambda(0) +\int ^b_0 M(\tau,0)\mu(\tau)d\tau 
 = \varphi(0).
\label{e1.7}
\end{gather}
Taking into account  $ \mu (0)= \lambda (0)$, we obtain
\begin{gather} 
\int ^{a}_0 K(0,\xi) d\xi \lambda(0)  
+  \int ^{a}_0 K(0,\xi)\lambda(\xi)d\xi
= \psi (0),    \label{e1.8} \\
\int ^b_0 M(\tau,0)d\tau \mu (0) +\int ^b_0 M(\tau,0)\mu(\tau)d\tau  
= \varphi(0). \label{e1.9}
\end{gather}

Let us define the following condition.
\smallskip

\noindent\textbf{Condition (i)}.  Assume 
\begin{gather*}
B_1(t) = \int ^{a}_0 K(t,\xi)d\xi \neq 0\quad\text{for all } t \in [0,T], \\
B_2(x) = \int ^{b}_0 M(\tau,x) \,d\tau \neq 0\quad \text{for all } x \in [0,\omega].
\end{gather*}
From relation \eqref{e1.4} we determine the parameter 
\begin{equation}
\mu(t) = - \frac{1}{B_1(t)}\Big\{ \int ^{a}_0 K(t,\xi) \widetilde{u}(t,\xi)d\xi
+  \int ^{a}_0 K(t,\xi)\lambda(\xi)d\xi - \psi (t)\Big\}, \label{e1.10}
\end{equation}
 for $t\in [0,T]$.
Similarly, from relation \eqref{e1.5} we determine parameter
\begin{equation}
\lambda(x) = - \frac{1}{B_2(x)} \Big \{ \int ^b_0 M(\tau,x) \widetilde{u}(\tau,x)d\tau
+ \int ^b_0 M(\tau,x)  \mu(\tau)d\tau  -  \varphi(x)\Big\}, \label{e1.11}
\end{equation}
for $ x \in [0,\omega]$.
Assumptions on the data of problem \eqref{e1}--\eqref{e3} allow us to
differentiate \eqref{e1.10} and \eqref{e1.11} by $t$ and $x$, respectively. Then
we obtain
\begin{gather}
\begin{aligned}
\dot{\mu}(t) &= \frac{\dot{B}_1(t)}{B^2_1(t)} \int ^{a}_0 K(t,\xi)
 \widetilde{u}(t,\xi)d\xi -\frac{1}{B_1(t)} \Big\{ \int ^{a}_0 \frac{\partial
K(t,\xi)}{\partial t} \widetilde{u}(t,\xi)d\xi  \\
&\quad + \int ^{a}_0 K(t,\xi) \frac{\partial\widetilde{u}(t,\xi)}{\partial t} d\xi
 \Big\} +\frac{\dot{B}_1(t)}{B^2_1(t)}\int ^{a}_0 K(t,\xi)\lambda(\xi)d\xi  \\
&\quad - \frac{1}{B_1(t)}\int ^{a}_0 \frac{\partial K(t,\xi)}{\partial t}
\lambda(\xi)d\xi - \frac{\dot{B}_1(t)}{B^2_1(t)}\psi (t)
+ \frac{1}{B_1(t)} \dot{\psi} (t),
\quad t\in [0,T],
\end{aligned} \label{e1.12}\\
\begin{aligned}
\dot{\lambda} (x)
&= \frac{\dot{B}_2(x)}{B^2_2(x)} \int ^b_0 M(\tau,x)
\widetilde{u}(\tau,x)d\tau - \frac{1}{B_2(x)}
\Big\{ \int ^b_0 \frac{\partial M(\tau,x)}{\partial x} \widetilde{u}(\tau,x)d\tau \\
&\quad +  \int ^b_0 M(\tau,x) \frac{\partial \widetilde{u}(\tau,x)}{\partial x}
 d\tau \Big\} + \frac{\dot{B}_2(x)}{B^2_2(x)}\int ^b_0 M(\tau,x)  \mu(\tau)d\tau\\
&\quad - \frac{1}{B_2(x)}\int ^b_0 \frac{\partial M(\tau,x)}{\partial x}
\mu(\tau)d\tau  - \frac{\dot{B}_2(x)}{B^2_2(x)} \varphi(x)
+ \frac{1}{B_2(x)} \dot{\varphi} (x),
\end{aligned} \label{e1.13}
\end{gather}
for $x \in [0, \omega]$.
We introduce the new unknown functions
\[
{\widetilde{v}(t,x) = \frac{\partial \widetilde{u}(t,x)}{\partial x}}, \quad
{\widetilde{w}(t,x) = \frac{\partial \widetilde{u}(t,x)}{\partial t}},
\]
and the following notation
\begin{gather*}
\begin{aligned}
G_1(t,\widetilde{u}, \widetilde{w})
&= \frac{\dot{B}_1(t)}{B^2_1(t)} \int ^{a}_0 K(t,\xi) \widetilde{u}(t,\xi)d\xi \\
&\quad - \frac{1}{B_1(t)} \Big\{ \int ^{a}_0 \frac{\partial K(t,\xi)}{\partial t}
\widetilde{u}(t,\xi)d\xi  + \int ^{a}_0 K(t,\xi) \widetilde{w}(t,\xi) d\xi \Big\},
\end{aligned}\\
L_1(t,\lambda) =\frac{\dot{B}_1(t)}{B^2_1(t)}\int ^{a}_0 K(t,\xi)\lambda(\xi)d\xi
- \frac{1}{B_1(t)}\int ^{a}_0 \frac{\partial
K(t,\xi)}{\partial t} \lambda(\xi)d\xi, \\
\begin{aligned}
G_2(x,\widetilde{u}, \widetilde{v})
&= \frac{\dot{B}_2(x)}{B^2_2(x)} \int ^b_0 M(\tau,x) \widetilde{u}(\tau,x)d\tau\\
&\quad - \frac{1}{B_2(x)} \Big\{ \int ^b_0 \frac{\partial M(\tau,x)}
{\partial x} \widetilde{u}(\tau,x)d\tau + \int ^b_0 M(\tau,x)  \widetilde{v}(\tau,x)
 d\tau \Big\}.
\end{aligned}
\end{gather*}
Then equations \eqref{e1.12} and \eqref{e1.13} can be written in the form
\begin{gather}
 \dot{\mu}(t) = G_1(t,\widetilde{u}, \widetilde{w})
+ L_1(t,\lambda) + \frac{1}{B_1(t)} \dot{\psi} (t)
- \frac{\dot{B}_1(t)}{B^2_1(t)}\psi (t),  \quad t\in [0,T], \label{e1.14}\\
\dot{\lambda}(x)  = G_2(x,\widetilde{u}, \widetilde{v})
+ L_2(x,\mu) + \frac{1}{B_2(x)} \dot{\varphi} (x)- \frac{\dot{B}_2(x)}{B^2_2(x)}
\varphi(x), \quad x \in [0, \omega]. \label{e1.15}
\end{gather}
Thus, we have a closed system of equations \eqref{e1.1}--\eqref{e1.3}, \eqref{e1.14}
\eqref{e1.9}, \eqref{e1.15} \eqref{e1.8} for determining the  unknown functions
$\widetilde{v}(t,x)$, $\widetilde{w}(t,x)$, $\widetilde{u}(t,x)$,
  $ \dot{\lambda}(x)$, $ \lambda(x)$,  $ \dot{\mu}(t)$, $ \mu(t)$.

Relation  \eqref{e1.14} in conjunction with \eqref{e1.9} present a boundary
value problem with integral condition for a differential equation
with respect to $ \mu (t)$, and the relation \eqref{e1.15} in conjunction
with \eqref{e1.8} present a boundary value problem with integral
condition for a differential equation with respect to $ \lambda(x)$.

Boundary value problem with integral condition \eqref{e1.14}, \eqref{e1.9} is
equivalent to relation  \eqref{e1.4}, and boundary value problem with
integral condition  \eqref{e1.15}, \eqref{e1.8} is equivalent to 
relation \eqref{e1.5} at $ \mu (0)=\lambda (0)$.

If  $ \dot{\mu}(t)$,  $ \dot{\lambda}(x)$, $ \mu(t)$, 
$\lambda(x)$ are known, then we find  the functions
$\widetilde{v}(t,x)$, $ \widetilde{w}(t,x)$, $\widetilde{u}(t,x)$
from \eqref{e1.1}--\eqref{e1.3}. Conversely, if we know the functions
$\widetilde{v}(t,x)$,  $ \widetilde{w}(t,x)$,
$\widetilde{u}(t,x)$, then we can find $\dot{\mu}(t)$, $ \mu(t)$,
$ \dot{\lambda}(x)$, $ \lambda(x)$ from boundary value problems
\eqref{e1.14}, \eqref{e1.9} and \eqref{e1.15}, \eqref{e1.8}. 
The unknowns are both $\widetilde{v}(t,x)$, $ \widetilde{w}(t,x)$,
$\widetilde{u}(t,x)$, and  $ \dot{\mu}(t)$, $ \mu(t)$, 
$\dot{\lambda}(x)$,  $ \lambda(x)$. Therefore, to find solution of
problem \eqref{e1.1}--\eqref{e1.5}, we use an iterative method: determine the
triplet $(\widetilde{u}^{\ast}(t,x)$,
 $\mu^{\ast}(t)$, $\lambda^{\ast}(x))$  as a limit of sequence
$(\widetilde{u}^{(m)}(t,x)$, $\mu^{(m)}(t))$, $\lambda^{(m)}(x))$,
 $ m =0,1,2,\dots$, according to the following algorithm:
\smallskip

\noindent\textbf{Step 0.} (1)
Assuming $\widetilde{u}(t,x) =0$, $\widetilde{w}(t,x) =0$,
$\lambda(x)=0$, in the right-hand side of equation \eqref{e1.14}, we find
initial approximations $\dot{\mu}^{(0)}(t)$, $ \mu^{(0)}(t)$,
$t\in [0,T]$ from the boundary value problem  with  integral
condition \eqref{e1.14}, \eqref{e1.9}. Assuming $\widetilde{u}(t,x) =0$,
$\widetilde{v}(t,x) =0$, $\mu (t)=0$ in the right-hand side of
equation \eqref{e1.15},    we find initial approximations
$\dot{\lambda}^{(0)}(x)$, $ \lambda^{(0)}(x)$, $x\in [0,\omega]$
from the boundary value problem  with integral condition \eqref{e1.15},
\eqref{e1.8}.

(2) Find $\widetilde{v}^{(0)}(t,x)$, $\widetilde{w}^{(0)}(t,x)$, 
$\widetilde{u}^{(0)}(t,x)$, $(t,x) \in \Omega$ from the Goursat
problem \eqref{e1.1}--\eqref{e1.3} for 
$ \dot{\lambda}(x)=\dot{\lambda}^{(0)}(x)$, 
$\dot{\mu}(t)= \dot{\mu}^{(0)}(t)$,
$\lambda (x) = \lambda^{(0)}(x)$,  
$ \mu (t) = \mu^{(0)}(t)$.
\smallskip

\noindent\textbf{Step 1.} (1) Assuming   
$\widetilde{u}(t,x)=\widetilde{u}^{(0)}(t,x)$,
$\widetilde{w}(t,x)=\widetilde{w}^{(0)}(t,x)$, $\lambda(x)=\lambda^{(0)}(x)$ in the
right-hand side of equation \eqref{e1.14}, we find $\dot{\mu}^{(1)}(t)$,
$ \mu^{(1)}(t)$, $t\in [0,T]$ from the boundary value problem with
integral condition \eqref{e1.14}, \eqref{e1.9}.  Assuming $\widetilde{u}(t,x)
=\widetilde{u}^{(0)}(t,x)$, $\widetilde{v}(t,x)
=\widetilde{v}^{(0)}(t,x)$, $\mu (t)=\mu^{(0)}(t)$  in the
right-hand side of equation \eqref{e1.15}, we find
$\dot{\lambda}^{(1)}(x)$, $ \lambda^{(1)}(x)$, $x\in [0,\omega]$
from the boundary value problem  with integral condition \eqref{e1.15},
\eqref{e1.8}.

(2) Find  $\widetilde{v}^{(1)}(t,x), \widetilde{w}^{(1)}(t,x),
\widetilde{u}^{(1)}(t,x), $ $(t,x) \in \Omega$ from the Goursat
problem \eqref{e1.1}--\eqref{e1.3} for $ \dot{\lambda}(x)=
\dot{\lambda}^{(1)}(x)$, $\dot{\mu}(t)= \dot{\mu}^{(1)}(t)$,
$\lambda (x) = \lambda^{(1)}(x)$, $ \mu (t) = \mu^{(1)}(t)$.  \\
And so on.
\smallskip

\noindent\textbf{Step $m$.} (1) Assuming  $\widetilde{u}(t,x)
=\widetilde{u}^{(m-1)}(t,x)$, $\widetilde{w}(t,x)
=\widetilde{w}^{(m-1)}(t,x)$, $\lambda(x)=\lambda^{(m-1)}(x)$ in
the right-hand side of equation \eqref{e1.14}, we find
$\dot{\mu}^{(m)}(t)$, $ \mu^{(m)}(t)$, $t\in [0,T]$ from the
boundary value problem  with  integral condition  \eqref{e1.14}, \eqref{e1.9}.
Assuming
 $\widetilde{u}(t,x)
=\widetilde{u}^{(m-1)}(t,x)$, $\widetilde{v}(t,x)
=\widetilde{v}^{(m-1)}(t,x)$, $\mu (t)=\mu^{(m-1)}(t)$  in the
right-hand side of equation \eqref{e1.15}, we find
$\dot{\lambda}^{(m)}(x)$, $ \lambda^{(m)}(x)$, $x\in [0,\omega]$
from the boundary value problem  with integral condition  \eqref{e1.15},
\eqref{e1.8}.

(2) Find  $\widetilde{v}^{(m)}(t,x)$, $ \widetilde{w}^{(m)}(t,x)$,
$ \widetilde{u}^{(m)}(t,x), $ $(t,x) \in \Omega$, from the Goursat
problem \eqref{e1.1}--\eqref{e1.3} for  $ \dot{\lambda}(x)=
\dot{\lambda}^{(m)}(x)$, $\dot{\mu}(t)= \dot{\mu}^{(m)}(t)$,
$\lambda (x) = \lambda^{(m)}(x)$, $ \mu (t) = \mu^{(m)}(t)$, 
$m =1,2,\dots $.

The constructed algorithm consists of two parts: 
we solve the boundary value problems with  integral condition for the ordinary
differential equations \eqref{e1.14}, \eqref{e1.9} and \eqref{e1.15}, 
\eqref{e1.8}  in the first part, and we solve the Goursat problem for hyperbolic
equations with functional parameters in the second part.



\section{Boundary value problems with integral condition for the
differential equations}

Consider the boundary value problem with integral condition for
the ordinary differential equations
\begin{gather}
\dot{\mu}(t) =  \dot{g}_1(t), \quad t\in [0, T], \label{e2.1}\\
\int ^b_0 M(\tau,0)d\tau \mu (0) +\int ^b_0 M(\tau,0)\mu(\tau)d\tau  = \varphi(0),
\label{e2.2}
\end{gather}
where the function $g_1(t)$  is continuously differentiable on
$[0, T]$, the function $M(t,x)$ is continuous on $\Omega$,
$\varphi(0)$ is a constant, $ 0 < b \leq T$.

The function  $ \mu(t) \in C([0,T],\mathbb{R})$ having the derivative
$\dot{\mu}(t) \in C([0,T],\mathbb{R})$, is called a solution to problem
\eqref{e2.1}, \eqref{e2.2}, if it satisfies ordinary differential
 equation \eqref{e2.1} and boundary condition \eqref{e2.2}.

We also consider the boundary value problem with integral
condition for the ordinary differential equation of the 
type
\begin{gather}
\dot{\lambda}(x) = \dot{g}_2(x), \quad x\in [0,\omega], \label{e2.3}\\
\int ^{a}_0 K(0,\xi) d\xi \lambda(0)  +  \int ^{a}_0 K(0,\xi)\lambda(\xi)d\xi
= \psi (0),    \label{e2.4}
\end{gather}
where the function $g_2(x)$ is continuously differentiable on $[0,\omega]$, 
the function  $K(t,x)$ is continuous on $ \Omega$,
 $\psi (0)$ is a constant, $ 0 < a \leq \omega$.

The function $ \lambda (x) \in C([0,\omega],\mathbb{R})$ having the
derivative $\dot{\lambda} (x) \in C([0,\omega],\mathbb{R})$, is called a
solution to problem  \eqref{e2.3}, \eqref{e2.4}, if it satisfies  ordinary
differential equation \eqref{e2.3} and boundary condition \eqref{e2.4}.

General solution to equation \eqref{e2.1} has the form
$$
\mu(t) = g_1(t) + C_1, \quad t\in [0, T], 
$$
where $C_1$ is a constant.

Since  Condition (i)  holds,  the constant $C_1$ is uniquely  determined
 by \eqref{e2.2}:
$$ 
C_1 = \frac{1}{2 B_2(0)}  \varphi(0) -  \frac{1}{2 B_2(0)} 
\int ^b_0 M(\tau,0)\bigl\{ g_1(\tau) + g_1(0) \bigr\}d\tau. 
$$
Then the unique solution to problem \eqref{e2.1}, \eqref{e2.2} has the form
\begin{equation}
\mu(t) = g_1(t) + \frac{1}{2 B_2(0)}  \varphi(0)
-  \frac{1}{2 B_2(0)} \int ^b_0 M(\tau,0)\bigl\{ g_1(\tau)
+ g_1(0) \bigr\}d\tau,  \label{e2.5}
\end{equation}
for $t\in [0, T]$.
Analogously, the general solution to equation \eqref{e2.3} has the form
$$
\lambda (x) = g_2(x) + C_2, \quad x\in [0, \omega],
 $$
where $C_2$ is a constant.

Since Condition (i)  holds,  the constant $C_2$ is uniquely  determined 
by \eqref{e2.4}:
$$ 
C_2 = \frac{1}{2 B_1(0)} \psi(0) -  \frac{1}{2 B_1(0)} \int ^a_0 K(0,\xi) 
\bigl\{ g_2(\xi) + g_2(0) \bigr\} d\xi. 
$$
Then unique solution to problem \eqref{e2.3}, \eqref{e2.4} has the form
\begin{equation}
\lambda (x) = g_2(x) + \frac{1}{2 B_1(0)}  \psi(0)
-  \frac{1}{2 B_1(0)} \int ^a_0 K(0,\xi) \bigl\{ g_2(\xi)
+ g_2(0) \bigr\} d\xi, \label{e2.6}
\end{equation}
for  $x\in [0, \omega]$.
Below we give conditions for the unique solvability of boundary
value problems with integral condition \eqref{e2.1}, \eqref{e2.2} 
and \eqref{e2.3}, \eqref{e2.4}.

\begin{theorem} \label{thm2.1} 
 Suppose Condition (i) is holds.
Then  problem  \eqref{e2.1}, \eqref{e2.2} has a unique solution  
$\mu^{\ast}(t)\in C([0,T],\mathbb{R})$  representable in the form \eqref{e2.5}, and
\begin{equation}
\max  _{t\in [0,T]}|\mu^{\ast}(t)| \leq \mathcal{K}_1 \max \Bigl( \max  _{t\in
[0,T]}|g_1(t)|, |\varphi(0)|\Bigr), \label{e2.7}
\end{equation}
 where
\[
\mathcal{K}_1  = 1 + \frac{1}{2|B_2(0)|} \Bigl [ 1 + 2b \max _{t\in
[0,b]} |M(t,0)|\Bigr].
\]
\end{theorem}

\begin{theorem} \label{thm2.2} 
Suppose Condition (i) holds.
Then  problem \eqref{e2.3}, \eqref{e2.4} has a unique solution 
$\lambda^{\ast}(x)\in C([0,\omega],\mathbb{R})$ representable in the form
\eqref{e2.6},  and 
\begin{equation}
\max  _{x\in [0,\omega]}|\lambda^{\ast}(x)|
\leq \mathcal{K}_2 \max \Bigl( \max  _{x\in
[0,\omega]}|g_2(x)|, |\psi(0)|\Bigr), \label{e2.8}
\end{equation}
 where
\[
\mathcal{K}_2  =1 + \frac{1}{2|B_1(0)|}\Bigl [ 1 +  2a \max _{x\in
[0,a]} |K(0,x)|\Bigr].
\]
\end{theorem}



\section{ Algorithm's convergence and main result}

In Section 2, an algorithm for finding a solution to problem
\eqref{e1.1}--\eqref{e1.5}, which is equivalent to problem \eqref{e1}--\eqref{e3}, 
is constructed. To formulate the main result,  we let us give  few
assumptions and notation.
Let  Condition (i) hold, and introduce the notation:
\begin{gather*}
\alpha = \max  _{(t,x) \in \Omega} |A(t,x)|, \quad
\beta =\max  _{(t,x) \in \Omega} |B(t,x)|, \quad
\gamma = \max_{(t,x) \in \Omega} |C(t,x)|, \\
 H = \alpha + \beta + \gamma, \quad
\kappa_1 = \max  _{(t,x) \in \Omega} |K(t,x)|,\quad
\kappa_2 = \max  _{(t,x) \in \Omega}\Bigl|
\frac{\partial K(t,x)}{\partial t}\Bigr|, \\
\sigma_1 = \max  _{(t,x) \in \Omega} |M(t,x)|,\quad
\sigma_2 = \max  _{(t,x) \in \Omega}\Bigl| \frac{\partial
M(t,x)}{\partial x}\Bigr|, \\
\beta_1 = \max _{t \in [0,T]} |[B_1(t)]^{-1}|, \quad
\beta_2 = \max _{x \in [0,\omega]} |[B_2(x)]^{-1}|, \\
\delta_1 = \max _{t \in [0,T]} |\dot{B}_1(t)|, \quad
\delta_2 = \max _{x \in [0,\omega]} |\dot{B}_2(x)|, \\
l_{11} (a) = a \beta_1 \kappa_1 
 \Bigl\{ 1 +  \max (T,\omega)H e^{H(T + \omega)}\Bigr\},  \\
l_{21} (b) = b \beta_2 \sigma_1 \Bigl\{
1 + \max (T, \omega)H e^{H(T +\omega)} \Bigr\}, \\
l_{12} (a) = a \beta_1 \Bigl\{ \delta_1 \beta_1\kappa_1 +  \kappa_2  
 +  (\delta_1 \beta_1\kappa_1 +  \kappa_1 +\kappa_2)\max (T,
 \omega)H e^{H(T + \omega)} \Bigr\}, \\
l_{22} (b) = b \beta_2 \Bigl\{
\delta_2\beta_2 \sigma_1 + \sigma_2  +  (\delta_2\beta_2 \sigma_1 
+ \sigma_1+\sigma_2) \max (T, \omega)H e^{H(T +\omega)} \Bigr\}.
\end{gather*}

In Section  3, the conditions for unique solvability of boundary
value problems with integral condition \eqref{e2.1}, \eqref{e2.2} and \eqref{e2.3},
\eqref{e2.4} are established. At fixed $\widetilde{v}(t,x)$,
$\widetilde{w}(t,x)$, $\widetilde{u}(t,x)$ on each step of the
algorithm we solve the  boundary value problems with integral
condition  \eqref{e1.12}, \eqref{e1.9}  and \eqref{e1.13}, \eqref{e1.8}. 
Besides, in problem
\eqref{e1.12}, \eqref{e1.9}, we consider $ \lambda(x)$ is also known. Same for
 $\mu(t) $ in problem \eqref{e1.13}, \eqref{e1.8}. At fixed $ \dot{\lambda}(x)$, 
$\dot{\mu}(t)$, $ \lambda(x)$, $ \mu(t)$, we solved the Goursat
problem  \eqref{e1.1}--\eqref{e1.3}.

The following statement gives the conditions for the convergence
of proposed algorithm  and the existence of unique solution to
problem \eqref{e1.1}--\eqref{e1.5}.

\begin{theorem} \label{thm3.1} 
Let:
\begin{itemize}
\item[(1)] the functions $ A(t,x)$, $B(t,x)$, $C(t,x)$,  and $f(t,x)$ be
continuous on $\Omega$;

\item[(2)] the functions $K(t,x)$ and  $ \psi (t) $ be continuously
differentiable by $t$ on $\Omega$ and $[0, T]$, respectively; the
functions $M(t,x)$ and $ \varphi (x) $ be continuously
differentiable by  $x$ on $\Omega$ and $[0, \omega]$,
respectively;

\item[(3)]  Condition (i) hold;

\item[(4)] the inequality
$$ 
q = \max \Bigl(\mathcal{K}_1 l_{11}(a) 
+ \mathcal{K}_2 l_{21}(b), l_{12}(a), l_{22}(b)\Bigl) <1
$$
be fulfilled.
\end{itemize}
Then  problem \eqref{e1.1}--\eqref{e1.5} has a unique solution.
\end{theorem}

\begin{proof} Let  conditions  (1)--(3)  be fulfilled.
We use step 0 of algorithm  and consider the
 boundary value problem with integral condition
\begin{gather}
 \dot{\mu}(t) =  \frac{1}{B_1(t)} \dot{\psi} (t)
- \frac{\dot{B}_1(t)}{B^2_1(t)}\psi (t), \quad  t\in [0,T],
\label{e3.1}\\
 \int ^b_0 M(\tau,0)d\tau \mu (0) +\int ^b_0 M(\tau,0)\mu(\tau)d\tau  = \varphi(0).
\label{e3.2} \\
\dot{\lambda}(x)  =\frac{1}{B_2(x)} \dot{\varphi} (x)
- \frac{\dot{B}_2(x)}{B^2_2(x)} \varphi(x), \quad x \in [0, \omega], \label{e3.3}\\
 \int ^{a}_0 K(0,\xi) d\xi \lambda(0)  +  \int ^{a}_0 K(0,\xi)\lambda(\xi)d\xi
= \psi (0).    \label{e3.4}
\end{gather}
Condition (3) and the conditions
of Theorems \ref{thm2.1} and \ref{thm2.2} yield the unique solvability of problems
\eqref{e3.1}, \eqref{e3.2} and \eqref{e3.3}, \eqref{e3.4}. 
We find initial approximations
$\mu^{(0)}(t)$ and $\lambda^{(0)}(x)$   from the boundary value
problems \eqref{e3.1}, \eqref{e3.2} and \eqref{e3.3}, \eqref{e3.4}. 
Then, similar to the estimates \eqref{e2.7} and \eqref{e2.8}, 
for the functions $\mu^{(0)}(t)$,
$\lambda^{(0)}(x)$ and their derivatives 
$  \dot{\mu}^{(0)}(t)$, $ \dot{\lambda}^{(0)}(x)$
the estimates hold:
\begin{gather} 
\max  _{t\in [0, T]}|\mu^{(0)}(t)| \leq \mathcal{K}_1
\max \Bigl( \beta_1 \max  _{t\in [0,
T]}|\psi(t)|, |\varphi(0)|\Bigr), \label{e3.5} \\
 \max  _{t\in [0, T]}|\dot{\mu}^{(0)}(t)| \leq \beta_1 \max  _{t\in [0,
T]}|\dot{\psi}(t)| + \delta_1\beta^2_1\max  _{t\in [0,
T]}|\psi(t)|. \label{e3.6} \\
 \max  _{x\in [0,\omega]}|\lambda^{(0)}(x)| \leq \mathcal{K}_2
\max \Bigl( \beta_2 \max  _{x\in
[0,\omega]}|\varphi(x)|, |\psi(0)|\Bigr), \label{e3.7} \\
 \max  _{x\in [0,\omega]}|\dot{\lambda}^{(0)}(x)|\leq
\beta_2 \max  _{x\in
[0,\omega]}|\dot{\varphi}(x)| + \delta_2 \beta^2_2\max  _{x\in
[0,\omega]}|\varphi(x)|.
\label{e3.8}
\end{gather}
Solving the Goursat problem  \eqref{e1.1}--\eqref{e1.3} for the found values of
parameters,  we find $\widetilde{v}^{(0)}(t,x)$,
$\widetilde{w}^{(0)}(t,x)$, $ \widetilde{u}^{(0)}(t,x)$ for all
$(t,x)\in \Omega$.

The following inequalities are valid:
\begin{gather*}
 |\widetilde{v}^{(0)}(t,x)| \leq \max (T, \omega)
 e^{H(T + \omega)} \max  _{(t,x) \in \Omega}
|\widetilde{f}(t,x)|, \\
|\widetilde{w}^{(0)}(t,x)| \leq \max (T, \omega)
e^{H(T + \omega)} \max  _{(t,x) \in \Omega}
|\widetilde{f}(t,x)|, \\
 |\widetilde{u}^{(0)}(t,x)| \leq \max (T, \omega)
 e^{H(T + \omega)}\max  _{(t,x) \in \Omega}
|\widetilde{f}(t,x)|,
\end{gather*}
where 
\[
 \widetilde{f}(t,x) = A(t,x) \dot{\lambda}^{(0)}(x) +
B(t,x)\dot{\mu}^{(0)}(t) + C(t,x)\Bigl [\lambda^{(0)}(x) +
\mu^{(0)}(t)\Bigr] + f(t,x).
\]
Then, we determine the functions
 $\mu^{(m)}(t)$,  $\lambda^{(m)}(x)$, 
 $\dot{\mu}^{(m)}(t)$, $ \dot{\lambda}^{ (m)}(x)$, 
$\widetilde{v}^{(m)}(t,x)$, $ \widetilde{w}^{(m)}(t,x)$, 
$\widetilde{u}^{(m)}(t,x) $ from the $m$th step of the algorithm,
and we obtain $\mu^{(m+1)}(t)$, $\lambda^{(m+1)}(x)$, 
$\dot{\mu}^{(m+1)}(t)$, $ \dot{\lambda}^{(m+1)}(x)$,   
$\widetilde{v}^{(m+1)}(t,x)$, $ \widetilde{w}^{(m+1)}(t,x)$, 
$\widetilde{u}^{(m+1)}(t,x) $, from step $(m+1)$, $m=1, 2,
\dots$.

 Evaluating the corresponding differences of successive
approximations, we obtain
\begin{gather}
\begin{aligned} 
&\max  _{t\in [0, T]}|\mu^{(m+1)}(t) - \mu^{(m)}(t) | \\
&\leq  \mathcal{K}_1 \beta_1 \max  _{t\in [0, T]}\Big
[ \int ^a_0 |K(t,\xi) (\widetilde{u}^{(m)}(t,\xi) 
 - \widetilde{u}^{(m-1)}(t,\xi)) |d\xi \\
&\quad + \int ^a_0 |K(t,\xi) (\lambda^{(m)}(\xi)
  - \lambda^{(m-1)}(\xi))|d\xi \Big], 
\end{aligned} \label{e3.9} \\
\begin{aligned}
&\max  _{x\in [0,\omega]}|\lambda^{(m+1)}(x) - \lambda^{(m)}(x) | \\
&\leq \mathcal{K}_2  \beta_2 \max  _{x\in [0,\omega]}
\Big [\int ^b_0 |M(\tau,x)(\widetilde{u}^{(m)}(\tau,x)
 - \widetilde{u}^{(m-1)}(\tau,x))|d\tau \\
&\quad  + \int ^b_0 |M(\tau,x)(\mu^{(m)}(\tau) - \mu^{(m-1)}(\tau))|d\tau
\Big], 
\end{aligned} \label{e3.10} \\
\begin{aligned}
&\max  _{t\in [0, T]}|\dot{\mu}^{(m+1)}(t) - \dot{\mu}^{(m)}(t) | \\
&\leq \max  _{t\in [0,T]}\Bigl [|L_1(t,\lambda^{(m)} - \lambda^{(m-1)})|
+ |G_1(t,\widetilde{u}^{(m)}- \widetilde{u}^{(m-1)},
\widetilde{w}^{(m)}- \widetilde{w}^{(m-1)})| \Bigr],
\end{aligned} \label{e3.11} \\
\begin{aligned}
&\max  _{x\in [0,\omega]}|\dot{\lambda}^{(m+1)}(x) - \dot{\lambda}^{(m)}(x) | \\
&\leq \max  _{x\in [0,\omega]}\Bigl [|L_2(x,\mu^{(m)} - \mu^{(m-1)})|
+ |G_2(x,\widetilde{u}^{(m)}- \widetilde{u}^{(m-1)},
\widetilde{v}^{(m)}- \widetilde{v}^{(m-1)})|\Bigr], 
\end{aligned} \label{e3.12} \\
\begin{aligned}
&|\widetilde{v}^{(m+1)}(t,x) - \widetilde{v}^{(m)}(t,x)|\\
& \leq \max (T, \omega) e^{H(T + \omega)} \Big \{
\alpha \max _{x\in [0,\omega]}|\dot{\lambda}^{(m+1)}(x) -
\dot{\lambda}^{(m)}(x) | \\
&\quad + \beta \max _{t\in
[0,T]}|\dot{\mu}^{(m+1)}(t) - \dot{\mu}^{(m)}(t)| 
 + \gamma \Bigl[\max _{x\in [0,\omega]}|\lambda^{(m+1)}(x) -
\lambda^{(m)}(x) | \\
&\quad  + \max _{t\in [0,T]}|\mu^{(m+1)}(t) -
\mu^{(m)}(t)|\Bigr ] \Big\},
\end{aligned}  \label{e3.13} \\
\begin{aligned}
&|\widetilde{w}^{(m+1)}(t,x) - \widetilde{w}^{(m)}(t,x)|\\
& \leq \max (T, \omega) e^{H(T + \omega)} \Big \{
\alpha \max _{x\in [0,\omega]}|\dot{\lambda}^{(m+1)}(x) 
-\dot{\lambda}^{(m)}(x) | \\
&\quad +  \beta \max _{t\in [0,T]}|\dot{\mu}^{(m+1)}(t) 
 - \dot{\mu}^{(m)}(t)|   +\gamma \Bigl[\max _{x\in [0,\omega]}|\lambda^{(m+1)}(x) 
-\lambda^{(m)}(x) | \\
&\quad  + \max _{t\in [0,T]}|\mu^{(m+1)}(t) 
 -\mu^{(m)}(t)|\Bigr ] \Big\}, 
\end{aligned}\label{e3.14} \\
\begin{aligned}
&|\widetilde{u}^{(m+1)}(t,x) - \widetilde{u}^{(m)}(t,x)| \\
&\leq \max (T, \omega) e^{H(T + \omega)}
 \Big \{ \alpha \max_{x\in [0,\omega]}|\dot{\lambda}^{(m+1)}(x)
  -\dot{\lambda}^{(m)}(x) | \\
&\quad + \beta \max _{t\in [0,T]}|\dot{\mu}^{(m+1)}(t) - \dot{\mu}^{(m)}(t)| 
 + \gamma \Bigl[\max _{x\in [0,\omega]}|\lambda^{(m+1)}(x) 
 -\lambda^{(m)}(x) |  \\
&\quad + \max _{t\in [0,T]}|\mu^{(m+1)}(t)  -\mu^{(m)}(t)|\Bigr ] \Big\}. 
\end{aligned} \label{e3.15} 
\end{gather}
Suppose that
\begin{align*}
 \Delta _{m+1} 
&= \max \Bigl( \max_{x\in [0,\omega]}|\lambda^{(m+1)}(x) - \lambda^{(m)}(x)| 
  +\max _{t\in [0,T]}|\mu^{(m+1)}(t) - \mu^{(m)}(t)|, \\
&\quad  \max_{x\in [0,\omega]}|\dot{\lambda}^{(m+1)}(x) -
\dot{\lambda}^{(m)}(x) |, \max _{t\in
[0,T]}|\dot{\mu}^{(m+1)}(t) - \dot{\mu}^{(m)}(t)|\Bigr).
\end{align*}

Then, from relations \eqref{e3.9}--\eqref{e3.12}, taking into account the
notation introduced and estima\-ti\-ons \eqref{e3.13}--\eqref{e3.15}, 
we obtain the main inequality
\begin{equation}
 \Delta _{m+1} \leq q \Delta_m.  \label{e3.16}
\end{equation}
 Condition  (4) of the theorem leads to the convergence of sequence
$ \Delta _m \to 0 $  as $m \to \infty $, i.e., $\Delta _{\ast} =0$.
This gives the uniform convergence of sequences  $\lambda^{(m)}(x)$,
$ \dot{\lambda}^{(m)}(x)$, $\mu^{(m)}(t)$, $ \dot{\mu}^{(m)}(t)$, to
  $\lambda^{\ast}(x)$,  $\dot{\lambda}^{\ast}(x)$, $\mu^{\ast}(t)$,
$ \dot{\mu}^{\ast}(t)$,respectively, as $m \to \infty $. Functions
$\lambda^{\ast}(x)$ and $\mu^{\ast}(t)$ are continuous and continuously
differentiable on $[0,\omega]$ and $[0,T]$, respectively. Based on estimates
\eqref{e3.13}--\eqref{e3.15}, we establish the uniform convergence of sequences
$\widetilde {v}^{(m)}(t,x)$, $ \widetilde {w}^{(m)}(t,x)$,
$\widetilde {u}^{(m)}(t,x)$  to the functions
 $\widetilde {v}^{\ast}(t,x)$,
$ \widetilde{w}^{\ast}(t,x)$, $\widetilde {u}^{\ast}(t,x)$, respectively, with
respect to $(t,x)\in \Omega$. Obviously, the functions
$\widetilde{u}^{\ast}(t,x)$, $\widetilde {v}^{\ast}(t,x)$, and
$ \widetilde {w}^{\ast}(t,x)$ are continuous on $\Omega$. Solving the problems
on the $(m + 1)$th step of the algorithm and passing to the limit
as $m \to \infty $, we obtain that the functions
$\widetilde{u}^{\ast}(t,x), \lambda^{\ast} (x)$,
$ \mu ^{\ast}(t)$ together with their derivatives  satisfy the Goursat
problem \eqref{e1.1}--\eqref{e1.3} and boundary value problems with integral
condition  \eqref{e1.14}, \eqref{e1.9} and \eqref{e1.15}, \eqref{e1.8}.

We carry out the inverse transition from  problem \eqref{e1.14}, \eqref{e1.9}
 to  relation \eqref{e1.4}, and pass from problem \eqref{e1.15}, \eqref{e1.8}   to
relation \eqref{e1.5}. Then the triplet of functions
$(\widetilde{u}^{\ast}(t,x), \lambda^{\ast} (x), \mu ^{\ast}(t))$
is solution to problem \eqref{e1.1}--\eqref{e1.5}.

Next we prove the uniqueness of solution to \eqref{e1.1}--\eqref{e1.5}.  
Let  the triplet  of functions $(\widetilde
{u}^{\ast}(t,x), \lambda^{\ast}(x), \mu^{\ast}(t))$ and another
triplet  of functions $(\widetilde {u}^{\ast \ast}(t,x),
\lambda^{\ast \ast}(x), \mu^{\ast \ast}(t) )$ be two solutions to
the problem. We introduce the notation
\begin{align*}
 \widetilde{\Delta} 
&= \max \Bigl (\max _{x\in [0,\omega]}|\lambda^{\ast}(x) - \lambda^{\ast
\ast}(x) | + \max _{t\in [0,T]}|\mu^{\ast}(t) - \mu^{\ast
\ast}(t)|, \\
&\quad \max_{x\in [0,\omega]}|\dot{\lambda}^{\ast}(x) -
\dot{\lambda}^{\ast \ast}(x) |, 
\max _{t\in [0,T]}|\dot{\mu}^{\ast}(t) - \dot{\mu}^{\ast \ast}(t)|, \Bigr).
\end{align*}
After calculation,  analogous to \eqref{e3.9}--\eqref{e3.15}, we obtain
\begin{equation}
 \widetilde{\Delta} \leq q \widetilde{\Delta}.\label{e3.17}
\end{equation}
By condition (4) of the theorem, we have $ q < 1$. Then
 inequality  \eqref{e3.17} takes place only  for $ \widetilde{\Delta}
\equiv 0.$ This gives us  $ \lambda^{\ast}(x) =\lambda^{\ast
\ast}(x)$,
$ \mu^{\ast}(t) = \mu^{\ast \ast}(t)$ and
$ \widetilde {u}^{\ast}(t,x) = \widetilde {u}^{\ast \ast}(t,x)$. Therefore, the
solution to problem \eqref{e1.1}--\eqref{e1.5} is unique.
\end{proof}

The next assertion follows from the equivalence of problem  \eqref{e1}--\eqref{e3} and 
problem  \eqref{e1.1}--\eqref{e1.5}.

\begin{theorem} \label{thm3.2} 
Let  conditions  (1)--(4) of Theorem \ref{thm3.1} be fulfilled.
Then problem  \eqref{e1}--\eqref{e3} has a unique classical solution.
\end{theorem}

\begin{proof}  Conditions (1)--(4) of Theorem \ref{thm3.1} imply the existence
of a unique solution to  \eqref{e1.1}--\eqref{e1.5}, the triplet of
functions $(\widetilde {u}^{\ast}(t,x), \lambda^{\ast}(x),
\mu^{\ast}(t))$. According to the  algorithm presented above, for
each $m =0, 1,2, \dots $, this triplet  is determined as a limit of
sequence triplets $(\widetilde{u}^{(m)}(t,x)$, $\mu^{(m)}(t))$,
$\lambda^{(m)}(x))$ as $m \to \infty$.

Then  solution to problem  \eqref{e1}--\eqref{e3}, the function
${u}^{\ast}(t,x)$, exists and is determined by the equality
${u}^{\ast}(t,x) = \widetilde {u}^{\ast}(t,x) + \lambda^{\ast}(x)
+ \mu^{\ast}(t)$.
\end{proof}


\subsection*{Acknowledgements} The author thanks the anonymous referee
for the careful reading of this article and  the useful
suggestions. This research was  partially
supported by Grant of Ministry of Education and Science of the
Republic of Kazakhstan, No 0822/$\Gamma \Phi$4.


\begin{thebibliography}{00}

\bibitem{a1} A. A. Abramov, L. F. Yukhno;
\emph{Nonlinear eigenvalue problem for a system of ordinary differential 
equations subject to a nonlocal condition}, 
Computational Mathematics and Mathematical
Physics,  52 (2012), No 2, pp. 213--218.

\bibitem{a2} A. A. Abramov, L. F. Yukhno;
\emph{ Solving a system of linear
ordinary differential equations with redundant conditions},
Computational Mathematics and Mathematical
Physics,  54 (2014), No 4, pp. 598--603.

\bibitem{a3} A. A. Abramov, L. F. Yukhno;
\emph{A solution method for a nonlocal
problem for a system of linear differential equations},
Computational Mathematics and Mathematical
Physics,  54 (2014), No 11, pp. 1686--1689.

\bibitem{a4} A. T. Asanova, D. S. Dzhumabaev;
\emph{Unique Solvability of the Boundary Value Problem for Systems of 
Hyperbolic Equations with Data on the Characteristics}, 
Computational Mathematics and Mathematical Physics,  42 (2002), No 11, 
pp. 1609--1621.

\bibitem{a5} A. T. Asanova, D. S. Dzhumabaev;
\emph{Unique solvability ofnonlocal boundary value  problems for systems  
of hyperbolic equations}, Differential Equations,  39 (2003), No 10, 
pp. 1414--1427. 

\bibitem{a6} A. T. Asanova, D. S. Dzhumabaev;
\emph{Well-posedness of nonlocalboundary value problems with integral 
condition for the system of hyperbolic equations},Journal of Mathematical
Analysis and Applications, 402 (2013), No 1, pp.167--178. 

\bibitem{a7} A. T. Assanova;
\emph{Nonlocal problem with integral conditions
for the system of hyperbolic equations in the characteristic
rectangle}, Russian Math. (Iz. VUZ), 61 (2017), No 5, pp.  7--20.

\bibitem{b1} A. Bouziani;
\emph{Solution forte d'un probleme mixte avec
conditions non locales pour une classe   d'equations
hyperboliques}, Bull. CI. Sci. Acad. Roy. Belg. 1997. Vol. 8. P. 53--70. 

\bibitem{b2} L. Byszewski;
\emph{Existence and uniqueness of solutions of nonlocal problems for 
hyperbolic equation $u_{xt}=F(x,t,u,u_x)$},
Journal of Applied Mathematics and Stochastic Analysis, 3 (1990),
No 3, pp. 163--168.

\bibitem{g1}  N. D. Golubeva,  L. S. Pul'kina;
\emph{A nonlocal problem with integral conditions}, 
Math. Notes, 59(1996), No 3, pp. 326--328.

\bibitem{k1} O. M. Kechina;
\emph{Nonlocal problem for hyperbolic equation with
conditions given into characteristic rectangle}, Vestn. SamGU.
Estestvonauchn. ser. 72 (2009), No 6. pp. 50--56. (in Russian)

\bibitem{k2} T. Kiguradze;
\emph{Some boundary value problems for systems of
linear partial differential equations of hyperbolic type}, Mem.
Differential Equations and Math. Phys. 1994. Vol. 1.  P. 1--144.

\bibitem{k3} A. I. Kozhanov, L. S. Pul'kina;
\emph{On the solvability of boundary value problems with a nonlocal 
boundary condition of integral form for multidimensional hyperbolic equations}, 
Differ. Equ., 42 (2006), No 9, pp. 1233--1246.

\bibitem{m1} E. I. Moiseev, V. I. Korzyuk, I. S. Kozlovskaya;
\emph{Classical solution of a problem with an integral condition for the
one-dimensional wave equation},  Differ. Equ., 50 (2014),
No 10. C. 1364--1377. 

\bibitem{n1} A. M. Nakhushev;
\emph{Approximate method of solving boundary-value problems for differential 
equations and its application to the dynamics of soil moisture and groundwater},
Differ. Equ., 18 (1982), No 1,  pp. 60--67.

\bibitem{n2} A. M. Nakhushev;
\emph{Problems with shift for partial differential  equations},
 M. Nauka, 2006. (in Russian)

\bibitem{n3} Z. A. Nakhusheva;
\emph{On a nonlocal problem for partial differential equations},
 Differ. Equ., 22 (1986), No 1, pp. 171--174. (in Russian) 

\bibitem{p1} B. I. Ptashnyck;
\emph{Ill-posed boundary value problems for
partial differential equations}, Naukova Dumka, Kiev, Ukraine,
1984. (in Russian)

\bibitem{p2} L. S. Pulkina;
\emph{A nonlocal problem with integral conditions
for the quasilinear hyperbolic equation},
Electronic Journal of Differential Equations,  1999 (1999), No. 45, pp. 1--6.

\bibitem{p3} L. S. Pulkina;
\emph{The $L_2$ solvability of a nonlocal problem
with integral conditions for a hyperbolic equation}, Differ. Equ.,
36 (2000), No 2, pp. 316--318.

\bibitem{p4} L. S. Pulkina;
\emph{A nonlocal problem with integral conditions for hyperbolic equations},
 Math. Notes, 70 (2001), No 1, pp.79--85.

\bibitem{p5} L. S. Pul'kina;
\emph{A nonlocal problem for a hyperbolic equation with integral conditions
 of the 1st kind with time-dependent kernels}, 
Russian Math. (Iz. VUZ), 58 (2012), No 10, pp.  26--37.

\bibitem{p6} L. S. Pul'kina, O. M. Kechina;
\emph{A nonlocal problem with integral conditions for hyperbolic equations 
in  characteristic rectangle}, Vestn. SamGU. Estestvonauchn. ser. 68 (2009),
 No 2, pp. 80--88. (in Russian) 

\bibitem{s1} K. B. Sabitov;
\emph{Boundary value problem  for a parabolic-hyperbolic equation with 
a nonlocal integral condition}, Differ. Equ., 46 (2010), No 10, pp. 1472--1481.

\bibitem{s2} Yu. K. Sabitova;
\emph{Nonlocal initial-boundary-value problem for
a degenerate hyperbolic equation},  Russian Math. (Iz. VUZ), 53
(2009), No 12, pp. 49--58.

\bibitem{s3} Yu. K. Sabitova;
\emph{Boundary-value problem with nonlocal integral
condition for mixed type equations with degeneracy on the transition
line}, Math. Notes, 98 (2015), No 3, pp.454--465.

\bibitem{t1} B. P. Tkach,  L. B. Urmancheva;
\emph{Numerical-analytical method for finding solutions of systems with distributed
parameters and integral condition}, Nonlinear Oscillations,
12 (2009), No 1, pp. 110--119.

\bibitem{z1} S. V. Zhestkov;
\emph{The Goursat problem with integral boundary conditions}, 
Ukranian Mathematical Journal, 42 (1990), No 1, pp. 119--122. 

\end{thebibliography}

\end{document}
