\documentclass[12pt,draft,a4paper]{amsart}
\usepackage{amsmath,amsthm}
\usepackage[T2A]{fontenc}
\usepackage[cp1251]{inputenc}
\usepackage[english]{babel}
\usepackage{amsfonts}
\usepackage{latexsym}
\usepackage{cite}

\makeatletter
\renewcommand{\@biblabel}[1]{#1.}
\makeatother
\tolerance 650

\newtheorem{theorem}{Theorem}
\newtheorem{lemma}{Lemma}
\newtheorem{corollary}{Corollary}
\newtheorem{proposition}{Proposition}

\newcommand{\No}{No}
\newcommand*{\ov}{\overline}
\newcommand*{\ei}{e^{i\theta}}
\newcommand{\cd}{(\cdot)}
\newcommand*{\la}{\langle}
\newcommand*{\ra}{\rangle}

\DeclareMathOperator*{\infp}{inf\vphantom p}
\DeclareMathOperator*{\codim}{codim}
\DeclareMathOperator*{\spa}{span}
\DeclareMathOperator*{\Ker}{Ker}
\DeclareMathOperator*{\supp}{supp}
\DeclareMathOperator*{\IM}{Im}
\DeclareMathOperator*{\RE}{Re}
\DeclareMathOperator{\sn}{sn}
\DeclareMathOperator{\cn}{cn}
\DeclareMathOperator{\dn}{dn}
\DeclareMathOperator*{\sign}{sign}
\DeclareMathOperator*{\card}{card}
%\DeclareMathOperator*{\liminf}{lim\,inf}
\DeclareMathOperator*{\vrai}{vrai\,sup}

\begin{document}
\thispagestyle{empty}
\begin{center}
USSR ACADEMY OF SCIENCE\\
SCIENTIFIC COUNCIL FOR CYBERNETICS
\end{center}
\vspace{50pt}
\vspace{60pt}
\begin{center}
{\bf ON SOME PROBLEMS OF OPTIMAL RECOVERY OF 
ANALYTIC AND HARMONIC FUNCTIONS FROM INACCURATE DATA}\\

\vspace{50pt}
K.Yu.~Osipenko and M.I.~Stessin
\vspace{50pt}

\vspace{100pt}

Preprint\\

\vspace{200pt}

Moscow 1990
\end{center}

\newpage


\underline{Introduction}. Let $X$ be a linear space, $Y$, $Z$ --- normed spaces, $W\subset X$ --- subset of $X$. Consider the problem of optimal recovery of the operator $L\colon W\to Z$ using the values of the information operator $I\colon W\to Y$ in the case those values are inaccurate ones. More precisely, let us consider the extremal problem
\begin{equation}\label{1}
E(L,I,\delta)=\infp_S\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-Sy\|,
\end{equation}
where $S\colon Y\to Z$ is some mapping (algorithm). $E(L,I,\delta)$ is called the intrinsic error of recovery. An algorithm $S_0$ is called an optimal one if
$$\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=E(L,I,\delta).$$
If $S_0$ is an optimal algorithm, $x_0\in W$, and
$$\sup_{\|Ix_0-y\|\le\delta}\|Lx_0-S_0y\|=E(L,I,\delta),$$
then $x_0$ is called a worst element.

The investigations of the problem \eqref{1} were initiated in \cite{1} for the case $\dim Y<\infty$. The case $\dim Y=\infty$ was worked out in \cite{2} (see also \cite{3}). In this paper we consider some problems of optimal recovery of analytic functions from the Hardy space $H_p$ and the Bergman space $A_p$, $1\le p\le\infty$. We also consider the same problems for harmonic functions from similar classes $h_p$ and $a_p$. Some results related to $H_\infty$, can be found in \cite{4,5,6,7}.

In Section 1 we prove some general theorems on optimal recovery from inaccurate data, which closely relate to results obtained in \cite{2,3,8,9}. In Section 2 we apply these theorems to finding optimal recovery algorithms for functions from $H_p$, $A_p$, $h_p$, and $a_2$ in some point of the unit disc of $\mathbb C$, when the disposed data is the inaccurate value of these functions in some other point. In particular we obtain some generalizations of Schwartz’s
Lemma. In the last section we find optimal algorithms of recovery of $f'(0)$ from inaccurate data $f(-h)$ and $f(h)$, $h\in(0,1)$ in $H_p$ spaces. In $H_\infty$ we also find the optimal value of $h$ for which the intrinsic error is minimal.

1. \underline{Some general theorems on recovery from inaccurate data}. Now our aim is proving the sufficiency of some conditions for the algorithm $S_0$ to be an optimal one and $x_0$ to be a worst element. These conditions were originally found by Micchelli and Rivlin \cite{2}. Though it is closely connected with Micchelli and Rivlin’s result the theorem we need slightly differs from it.

\begin{theorem}\label{T1}
Let $x_0\in W$, $-x_0\in W$, $L(-x_0)=-Lx_0$, $\|Ix_0\|\le\delta$, $\|I(-x_0)\|\le\delta$, and
$$\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=\|Lx_0\|.$$
Then
\begin{itemize}
\item[i)] $S_0$ is an optimal algorithm,
\item[ii)] $x_0$ is a worst element,
\item[iii)] the intrinsic error is $E(L,I,\delta)=\|Lx_0\|$.
\end{itemize}
\end{theorem}

\begin{proof}
It follows from \eqref{1} that
$$E(L,I,\delta)\le\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=\|Lx_0\|.$$
On the other hand, for any algorithm $S$ we have
\begin{equation}\label{2}
\|Lx_0-S(0)\|+\|L(-x_0)-S(0)\|\ge2\|Lx_0\|
\end{equation}
and therefore
$$\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-Sy\|\ge\max\{\|Lx_0-S(0)\|,\|L(-x_0)-S(0)\|\}\ge\|Lx_0\|.$$
Thus $E(L,I,\delta)=\|Lx_0\|$ and $S_0$ is an optimal algorithm. Now suppose that $x_0$ is not a worst element, i.e.,
$$\sup_{\|Ix_0-y\|\le\delta}\|Lx_0-S_0y\|<\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=\|Lx_0\|.$$
Then
$$\|Lx_0-S_0(0)\|<\|Lx_0\|,\quad\|L(-x_0)-S_0(0)\|\le\|Lx_0\|,$$
which contradicts \eqref{2}. Theorem~\ref{T1} is proved.
\end{proof}

\begin{corollary}\label{C1}
Let $x_0\in W$, $-x_0\in W$, $L(-x_0)=-Lx_0$, $\|Ix_0\|\le\delta$, $\|I(-x_0)\|\le\delta$, $S_0$ is a linear operator and
$$\sup_{x\in W}\|Lx-S_0Ix\|=\|Lx_0\|-\delta\|S_0\|.$$
Then
\begin{itemize}
\item[i)] $S_0$ is an optimal algorithm,
\item[ii)] $x_0$ is a worst element,
\item[iii)] the intrinsic error is $\displaystyle E(L,I,\delta)=\sup_{\substack{x\in W\\\|Ix\|\le\delta}}\|Lx\|=\|Lx_0\|$.
\end{itemize}
\end{corollary}

\begin{proof}
Note that
\begin{multline*}
\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0Ix+S_0(Ix-y)\|\\
\le\sup_{x\in W}\|Lx-S_0Ix\|+\delta\|S_0\|=\|Lx_0\|.
\end{multline*}
Since $\|Ix_0\|\le\delta$ we have
$$\sup_{\substack{x\in W\\\|Ix\|\le\delta}}\|Lx\|\ge\|Lx_0\|\ge\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|\ge\sup_{\substack{x\in W\\\|Ix\|\le\delta}}\|Lx\|.$$
Thus
$$\sup_{\substack{x\in W\\\|Ix-y\|\le\delta}}\|Lx-S_0y\|=\|Lx_0\|=\sup_{\substack{x\in W\\\|Ix\|\le\delta}}\|Lx\|.$$
Now the corollary follows from Theorem~\ref{T1}.
\end{proof}

Let $\Omega$ be a subset of $\mathbb C^n$ and $\mu$ be a nonnegative measure on $\Omega$. Denote
by $L_p(\Omega,\mu)$ the Lebesgue space of complex (or real) -valued functions with the usual norm
\begin{gather*}
\|f\|_p=\biggl(\int_\Omega|f(z)|^p\,d\mu(z)\biggr)^{1/p},\quad1\le p<\infty,\\
\|f\|_\infty=\vrai_{z\in\Omega}|f(z)|.
\end{gather*}

Let $X_p$ be some linear subspace of $L_p(\Omega,\mu)$ and $BX_p=\{f\in X_p:\|f\|_p\le1\}$. Consider the problem \eqref{1} for $X=X_p$, $W=BX_p$ and $Z=\mathbb C(\mathbb R)$.

The following theorem is a generalization of the appropriate results from \cite{8,9} obtained for the case $\delta=0$.

\begin{theorem}\label{T2}
Let $g\in X_p$, $\|g\|_p\ne0$, $g_0=g/\|g\|_p$, $L$ be a functional on $X_p$, $L(-g_0)=-Lg_0$, $\|Ig_0\|\le\delta$, $\|I(-g_0)\|\le\delta$, $S_0$ is a linear functional, $S_0Ig_0=\delta\|S_0\|$, and for every $f\in BX_p$
\begin{equation}\label{3}
Lf-S_0If=\begin{cases}\displaystyle\alpha\int_\Omega\ov{g(z)}|g(z)|^{p-2}f(z)\,d\mu(z),\quad1\le p<\infty,\\[10pt]
\displaystyle\int_\Omega\ov{g(z)}|\varphi(z)|f(z)\,d\mu(z),\quad p=\infty,\end{cases}
\end{equation}
where $\alpha>0$, $\varphi\in L_1(\Omega,\mu)$ and if $p=\infty$ \ $|g(z)|=1$ almost everywhere on $\Omega$ with respect to measure $\mu$. Then
\begin{itemize}
\item[i)] $S_0$ is an optimal algorithm,
\item[ii)] $g_0$ is a worst element,
\item[iii)] the intrinsic error is
$$E(L,I,\delta)=\sup_{\substack{f\in BX_p\\\|If\|\le\delta}}|Lf|=Lg_0=\begin{cases}\alpha\|g\|_p^{p-1}+\delta
\|S_0\|,&1\le p<\infty,\\
\|\varphi\|_1+\delta\|S_0\|,&p=\infty.\end{cases}$$
\end{itemize}
\end{theorem}

\begin{proof}
For every $f\in BX_p$ from \eqref{3} and the H\"older inequality we have
$$|Lf-S_0If|\le\begin{cases}\alpha\|g\|_p^{p-1},&1\le p<\infty,\\
\|\varphi\|_1,&p=\infty.\end{cases}$$
On the other hand we obtain
$$Lg_0-\delta\|S_0\|=Lg_0-S_0Ig_0=\begin{cases}\alpha\|g\|_p^{p-1},&1\le p<\infty,\\
\|\varphi\|_1,&p=\infty.\end{cases}$$
Hence
$$\sup_{f\in BX_p}|Lf-S_0If|=Lg_0-\delta\|S_0\|.$$
Now the theorem follows from Corollary~\ref{C1}.
\end{proof}

Let $l_q^m$ be the space $\mathbb C^m$ supplied with the norm
$$\|a\|_q=\|(a_1,\ldots,a_m)\|_q=\begin{cases}\displaystyle\biggl(\sum_{j=1}^m|a_j|^q\biggr)
^{1/q},&1\le q<\infty,\\
\displaystyle\max_{1\le j\le m}|a_j|,&q=\infty,\end{cases}$$
and by $(a,b)$ denote the Hermitian inner product
$$(a,b)=\sum_{j=1}^ma_j\ov b_j.$$
Let $a\ne0$, $a^*\in Bl_q^m$ and
$$(a,a^*)=\|a\|_{q'},\quad1/q+1/q'=1.$$
It is easy to see that
$$a_j^*=\frac{a_j|a_j|^{q'-2}}{\|a\|_{q'}^{q'-1}},\quad1\le q'<\infty,$$
and for $q'=\infty$
$$a_j^*=\begin{cases}0,&j\ne j_0,\\
a_{j_0}/|a_{j_0}|,&j=j_0,\end{cases}$$
where $j_0$, such that $|a_{j_0}|=\displaystyle\max_{1\le j\le m}|a_j|$.

\begin{corollary}\label{C2}
Let $I\colon BX_p\to l_q^m$, $S_0y=(y,a)$, $a\in\mathbb C^m$, $g\in X_p$, $\|g\|_p\ne0$, $g_0=g/\|g\|_p$, $L$ be a functional, $L(-g_0)=-Lg_0$, $\|I(-g_0)\|_q\le\delta$.
Suppose that for every $f\in BX_p$ the equality \eqref{3} holds and $Ig_0=\delta a^*$ if $a\ne0$, or $\|Ig_0\|_q\le\delta$ if $a=0$. Then
\begin{itemize}
\item[i)]$S_0$ is an optimal algorithm,
\item[ii)]$g_0$ is a worst element,
\item[iii)]the intriilsic error is
$$E(L,I,\delta)=\sup_{\substack{f\in BX_p\\\|If\|_q\le\delta}}|Lf|=Lg_0=\begin{cases}\alpha\|g\|_p^{p-1}+\delta
\|a\|_{q'},&1\le p<\infty,\\
\|\varphi\|_1+\delta\|a\|_{q'},&p=\infty.\end{cases}$$
\end{itemize}
\end{corollary}

2. \underline{Optimal rcovery of analytic and harmonic functions}. Let $D=\{z\in\mathbb C:|z|<1\}$ and $H_p$ is the Hardy space, i.e., the space of functions which are analytic in $D$ and 
\begin{equation}\label{4}
\begin{aligned}
\|f\|_{H_p}&=\sup_{0<r<1}\left(\frac1{2\pi}\int_0^{2\pi}|f(r\ei)|^p\,d
\theta\right)^{1/p}<\infty,\quad1\le p<\infty,\\
\|f\|_{H_\infty}&=\sup_{z\in D}|f(z)|<\infty.
\end{aligned}
\end{equation}
It is well known that the functions from $H_p$ have boundary values almost everywhere and therefore $H_p$ can be considered as a subspace of $L_p(\Omega,\mu)$ for $\Omega=\{z\in\mathbb C:|z|=1\}$ and $d\mu(\ei)=\dfrac1{2\pi}d\theta$.

Recall that the Bergman space $A_p$ is the space of analytic functions which satisfy the inequality
\begin{equation}\label{5}
\|f\|_{A_p}=\left(\frac1\pi\int_D|f(z)|^p\,d\sigma(z)\right)^{1/p}<\infty,\quad1\le p<\infty,
\end{equation}
where $\sigma(z)$ is the Lebesgue measure on $D$ (for $p=\infty$  \ $A_\infty=H_\infty$). Thus the space $A_p$ is the subspace of $L_p(D,\mu)$ for $d\mu(z)=\dfrac1\pi d\sigma(z)$.

Denote by $h_p$ and $a_p$ the spaces of harmonic functions in $D$ which satisfy \eqref{4} and \eqref{5}, respectively.

Consider the problem \eqref{1} when $X$ is one of the spaces $H_p$, $A_p$, $h_p$, or $a_p$,
$W=BX$, $Lf=f(\xi)$, $If=f(z_1)$, $\xi$ and $z_1$ are distinct points in $D$. The relative intrinsic error will be denoted by $E(\xi,z_1,\delta,X)$.

Put
\begin{equation}\label{6}
\rho=\left|\frac{\xi-z_1}{1-\ov z_1\xi}\right|,\quad W(z)=e^{i\varphi}\frac{z-z_1}{1-\ov z_1z},
\end{equation}
where $\varphi$ is determined from the condition $W(\xi)=\rho$,
\begin{gather*}
\delta_1=\left(\frac{1-\rho}{2(1-|z_1|^2)}\right)^{1/p},\quad\delta_2=
\left(\frac{1-\rho^2}{1-|z_1|^2}\right)^{1/p},\\
h(z)=\begin{cases}1,&0\le\delta\le\delta_1,\\
(W(z)+a)/(1+aW(z)),&\delta_1<\delta<\delta_2,\\
W(z),&\delta\ge\delta_2,\end{cases}
\end{gather*}
where $a\in[0,1]$ and satisfies the equation
\begin{equation}\label{7}
\frac{a\delta_2}{h(z_1)(1+a\rho+a^2)^{1/p}}=\delta,\quad0\le\delta<\delta_2
\end{equation}
(the existence of a solution follows from the continuity of the function from the left hand side of \eqref{7}). Put $a=0$ for $\delta\ge\delta_2$.

\begin{proposition}\label{P1}
Let $X=H_p$. Then for every $1\le p<\infty$ and $\delta\ge0$
\begin{itemize}
\item[i)]
$$S_0y=\frac{h(z_1)(1-\rho^2)}{h(\xi)(1+a\rho)^{2(p-1)/p}}\left(\frac{1-\ov\xi z_1}{1-|\xi|^2}\right)^{2/p}y$$
is an optimal algorithm,
\item[ii)]
\begin{equation}\label{8}
g_0(z)=\left(\frac{1-|\xi|^2}{1+2a\rho+a^2}\right)^{1/p}\frac{(W(z)+a)(1+aW(z))^{(2-p)/p}}
{h(z)(1-\ov\xi z)^{2/p}}
\end{equation}
is a worst function,
\item[iii)]the intrinsic error is
$$E(\xi,z_1,\delta,H_p)=\frac{(\rho+a)(1+a\rho)^{(2-p)/p}}{h(\xi)(1+2a\rho+a^2)^{1/p}(1-|\xi|^2)
^{1/p}}.$$
\end{itemize}
\end{proposition}

\begin{proof}
Put
$$g(z)=\frac{(W(z)+a)(1+aW(z))^{(2-p)/p}}{h(z)(1-\ov\xi z)^{2/p}},\quad\alpha=\frac{\rho(1-|\xi|^2)^{(p-2)/p}}{h(\xi)(1+a\rho)^{2(p-1)/p}}.$$
By the residue theorem we have for every $f\in H_p$
\begin{multline*}
\alpha\frac1{2\pi}\int_0^{2\pi}\ov{g(\ei)}|g(\ei)|^{p-2}f(\ei)\,d\theta\\
=\alpha\frac1{2\pi i}\int_{|z|=1}\frac{(1+aW(z))^{2(p-1)/p}h(z)}{W(z)(z-\xi)(1-\ov\xi z)^{(p-2)/p}}f(z)\,dz=f(\xi)-S_0f(z_1).
\end{multline*}
Besides that
\begin{multline*}
\|g\|_{H_p}^p=\frac1{2\pi}\int_0^{2\pi}\left|\frac{1+aW(\ei)}{1-\ov\xi\ei}\right|^2\,d\theta\\
=\frac1{2\pi i}\int_{|z|=1}\frac{(1+aW(z))(W(z)+a)}{W(z)(1-\ov\xi z)(z-\xi)}\,dz=\frac{1+2a\rho+a^2}{1-|\xi|^2}.
\end{multline*}
It follows from \eqref{7} that for $0\le\delta<\delta_2$ \ $|g_0(z_1)|=\delta$ and for $\delta\ge\delta_2$ \ $|g_0(z_1)|=\delta_2\le\delta$, $S_0y\equiv0$. Since $S_0g_0(z_1)\ge0$, we have $S_0g_0(z_1)=\delta\|S_0\|$ for every $\delta\ge0$. Now the proposition follows from Theorem~\ref{T2}.
\end{proof}

Note that in virtue of Theorem~\ref{T2} the following generalization of the Schwartz Lemma can be obtained from Proposition~\ref{P1}:
\begin{multline}\label{9}
\sup_{\substack{f\in BH_p\\|f(0)|\le\delta}}|f(z)|\\
=\begin{cases}\dfrac{(|z|+a)(1+a|z|)^{(2-p)/p}}{(1-|z|^2)^{1/p}(1+2a|z|+a^2)^{1/p}},
&0\le\delta\le\left(\dfrac{1-|z|}2\right)^{1/p},\\
\dfrac{(1+a|z|)^{2/p}}{(1-|z|^2)^{1/p}(1+2a|z|+a^2)^{1/p}},&\left(\dfrac{1-|z|}2\right)^{1/p}\le
(1-|z|^2)^{1/p},\\
(1-|z|^2)^{-1/p},&\delta\ge(1-|z|^2)^{1/p},\end{cases}
\end{multline}
where $a$ is defined by \eqref{7} for $z_1=0$.

Now consider the same problem for $X=A_p$. Put
\begin{gather*}
\delta_1=\frac{(2+\rho)^{2/p}(1-\rho)^{2/p}}{2^{1/p}(3-\rho^2)^{1/p}(1-|z_1|^2)^{2/p}},\quad
\delta_2=\left(\frac{1-\rho^2}{1-|z_1|^2}\right)^{2/p},\\
b=\begin{cases}\dfrac1{1+a\rho},&0\le\delta<\delta_1,\\[10pt]
\dfrac a{a+\rho},&\delta\ge\delta_1,\end{cases}
\end{gather*}
where $a\in[0,1]$ and satisfies the equations
\begin{equation}\label{10}
\frac{a\rho^{2/p}((p/2-1)(1-a^2)b^2+b+b^2)^{2/p}(1-|z_1|^2)^{-2/p}}
{\left[\dfrac{(p/2-1)(1-a^2)(1+2a\rho+a^2)\rho^2b^4}{1-\rho^2}+\dfrac{\rho^2(2-\rho^2)}{
(1-\rho^2)^2}+(1-b^2)^2\right]^{1/p}}
=\delta
\end{equation}
for $0\le\delta<\delta_1$ and
$$\frac{(1-b^2)^{2/p}(1-\rho^2)^{2/p}}{(1-2(1-\rho^2)^2b^2+(1-\rho^2)^2b^4)^{1/p}
(1-|z_1|^2)^{2/p}}=\delta$$
for $\delta_1\le\delta<\delta_2$ (the solution of the last equation may be given in direct
form and the existence of solution \eqref{10} will be shown below). Put $a=0$ for $\delta\ge\delta_2$. Let
\begin{gather*}
\varphi(z)=\begin{cases}(p/2-1)(1-a^2)(1-\rho W(z))&\\
\hspace{75pt}+(1+aW(z))(2+a\rho-\rho W(z)),&0\le\delta<\delta_1,\\
(1+aW(z))(2a+\rho-a\rho W(z)),&\delta\ge\delta_1,\end{cases}\\
g(z)=\begin{cases}\dfrac{W(z)+a}{1+aW(z)}\dfrac{(\varphi(z))^{2/p}}{(1-\ov\xi z)^{4/p}},&0\le\delta<\delta_1,\\
\dfrac{(\varphi(z))^{2/p}}{(1-\ov\xi z)^{4/p}},&\delta\ge\delta_1.\end{cases}
\end{gather*}

\begin{proposition}\label{P2}
Let $X=A_p$. Then for every $1\le p<\infty$ and $\delta\ge0$
\begin{itemize}
\item[i)]
$$S_0y=b^2(1-\rho^2)^2\left(\frac{1-\ov\xi z_1}{1-|\xi|^2}\right)^{4/p}\left(\frac{\varphi(z_1)}{\varphi(\xi)}\right)^{(p-2)/p}y$$
is an optimal algorithm,
\item[ii)]$g_0=g/\|g\|_{A_p}$ is a worst function,
\item[iii)]the intrinsic error is
$$E(\xi,z_1,\delta,A_p)=\begin{cases}\dfrac{\rho((p/2)(1-\rho^2)+1)^{1/p}}{(1-|\xi|^2)^{2/p}}
,&\hspace{-9pt}\delta=0,\\
\delta b\dfrac{a+\rho}a\left(\dfrac{1-|z_1|^2}{(1-|\xi|^2)(1-\rho^2)}\right)^{2/p}
\left(\dfrac{\varphi(\xi)}{\varphi(z_1)}\right)^{2/p},&\\
&\hspace{-34pt}0<\delta<\delta_2,\\
\dfrac1{(1-|\xi|^2)^{2/p}},&\hspace{-13pt}\delta\ge\delta_2.\end{cases}$$
\end{itemize}
\end{proposition}

\begin{proof}
Note that the functions
$$(p/2-1)(1-a^2)(1-\rho w)+(1+aw)(2+a\rho-\rho w)$$
and
$$(1+aw)(2a+\rho-a\rho w),$$
as the functions of $w$, have real zeros which are outside the interval $(-1,1)$. Therefore $\varphi(z)$ is zero free in $D$. Let $0\le\delta<\delta_1$. For $f\in H_\infty$ denote
$$Jf=\frac1{2\pi i}\int_{|z|=1}\frac{(1+aW(z))^2(\varphi(z))^{(p-2)/p}}{W(z)(W(z)-\rho)(1-\ov\xi z)^{2(p-2)/p}}f(z)\,dz.$$
Since 
$$W(z)-\rho=e^{i\varphi}\frac{(z-\xi)(1-|z_1|^2)}{(1-\ov z_1z)(1-\ov z_1\xi)}$$
we obtain by the residue theorem
$$e^{i\varphi}\frac{\rho(1-|z_1|^2)(1-|\xi|^2)^{2(p-2)/p}}{(1+a\rho)^2(1-\ov z_1\xi)(\varphi(\xi))^{(p-2)/p}}Jf=f(\xi)-S_0f(z_1).$$
On the other hand, in view of the equality $\ov{W(\ei)}=W^{-1}(\ei)$, we obtain by
Stokes’ formula
\begin{multline*}
Jf=\frac1{2\pi i}\int_{|z|=1}\left(\frac{\ov{W(z)}+a}{1+a\ov{W(z)}}\right)^2\frac{(1+a\ov{W(z)})^2
(\varphi(z))^{(p-2)/p}}{(1-\rho\ov{W(z)})(1-\ov\xi z)^{2(p-2)/p}}f(z)\,dz\\
=\frac1{2\pi i}\int_{|z|=1}\left(\frac{\ov{W(z)}+a}{1+a\ov{W(z)}}\right)^{p/2+1}\frac{(1+a\ov{W(z)})^2}
{1-\rho\ov{W(z)}}\left(\frac{W(z)+a}{1+aW(z)}\right)^{p/2-1}\\
\times\frac{(\varphi(z))^{(p-2)/p}}{(1-\ov\xi z)^{2(p-2)/p}}f(z)\,dz=e^{-i\varphi}\frac{(1-\ov z_1\xi)^2}{1-|z_1|^2}\frac1\pi\int_D\left(\frac{\ov{W(z)}+a}{1+a\ov{W(z)}}\right)^{p/2}\\
\times\frac{\ov{\varphi(z)}}{(1-\xi\ov z)^2}\left(\frac{W(z)+a}{1+aW(z)}\right)^{p/2-1}\frac{(\varphi(z))^{(p-2)/p}}{(1-\ov\xi z)^{2(p-2)/p}}f(z)\,d\sigma(z)\\
=e^{-i\varphi}\frac{(1-\ov z_1\xi)^2}{1-|z_1|^2}\frac1\pi\int_D\ov{g(z)}|g(z)|^{p-2}f(z)\,d\sigma(z).
\end{multline*}
Thus for every $f\in H_\infty$ we have
\begin{equation}\label{11}
\frac{\rho(1-|\xi|^2)^{2(p-2)/p}}{(1+a\rho)^2(\varphi(\xi))^{(p-2)/p}}\frac1\pi\int_D
\ov{g(z)}|g(z)|^{p-2}f(z)\,d\sigma(z)=f(\xi)-S_0f(z_1).
\end{equation}

As functions from $H_\infty$ are dense in $A_p$ for every $1\le p<\infty$, the equality \eqref{11} holds for every function from $A_p$. It is easily seen that $S_0g(z_1)\ge0$.
Therefore it follows from Theorem~\ref{T2} that if $a\in[0,1]$ satisfies the condition $|g(z_1)|/\|g\|_{A_p}=\delta$, then $S_0$ is an optimal algorithm. For $f=g$, from \eqref{11} we
have
\begin{multline*}
\frac{\rho(1-|\xi|^2)^{2(p-2)/p}}{(1+a\rho)^2(\varphi(\xi))^{(p-2)/p}}\|g\|_{A_p}^p\\
=\frac{(\rho+a)(\varphi(\xi))^{2/p}}{(1+a\rho)(1-|\xi|^2)^{4/p}}-\left(\frac{1-\rho^2}{1+a\rho}
\right)^2\frac{a\varphi(z_1)}{(1-|\xi|^2)^{4/p}(\varphi(\xi))^{(p-2)/p}}.
\end{multline*}
Hence
$$\|g\|_{A_p}^p=\frac{(\rho+a)(1+a\rho)}{\rho(1-|\xi|^2)^2}\varphi(\xi)-
\frac{a(1-\rho^2)^2}{\rho(1-|\xi|^2)^2}\varphi(z_1).$$
We find by direct calculation that
\begin{multline*}
\|g\|_{A_p}^p=\frac1{b^4\rho^2(1-|\xi|^2)^2}\biggl(\biggl(\frac p2-1\biggr)(1-a^2)(1-\rho^2)(1+2a\rho+a^2)\rho^2b^4\\
+1-2(1-\rho^2)^2b^2+(1-\rho^2)^2b^4\biggr).
\end{multline*}
Since $\|g\|_{A_p}>0$ for every $a\in[0,1]$, the function in the left-hand side
of \eqref{10} is continuous as a function of $a$ ($a\in[0,1]$), and therefore the
equation \eqref{10} has a solution for every $\delta\in[0,\delta_1)$. We have
\begin{multline*}
|g(z_1)|=a\frac{|\varphi(z_1)|^{2/p}}{|1-\ov\xi z_1|^{4/p}}=a\frac{((p/2-1)(1-a^2)+2+a\rho)^{2/p}}{|1-\ov\xi z_1|^{4/p}}\\
=\frac{a((p/2-1)(1-a^2)+1+b^{-1})^{2/p}(1-\rho^2)^{2/p}}{(1-|\xi|^2)^{2/p}(1-|z_1|^2)^{2/p}}
\end{multline*}
and so the equation \eqref{10} means that $|g(z_1)|/\|g\|_{A_p}=\delta$.

The case $\delta\in[\delta_1,\delta_2)$ can be considered in the same way if we set
\begin{multline*}
Jf=\frac1{2\pi i}\int_{|z|=1}\frac{(W(z)+a)^2(\varphi(z))^{(p-2)/p}}{W(z)(W(z)-\rho)(1-\ov\xi z)^{2(p-2)/p}}f(z)\,dz\\
=\frac1{2\pi i}\int_{|z|=1}\frac{(1+a\ov{W(z)})^2(\varphi(z))^{(p-2)/p}}{(1-\rho\ov{W(z)})(1-\ov\xi z)^{2(p-2)/p}}f(z)\,dz.
\end{multline*}
Let now $\delta\ge\delta_2$. Then $a=0$, $g(z)=\rho^{2/p}(1-\ov\xi z)^{-4/p}$, and
\begin{multline*}
\frac1\pi\int_D\ov{g(z)}|g(z)|^{p-2}f(z)\,d\sigma(z)
=\frac{\rho^{2(p-1)/p}}\pi\int_D\frac{f(z)\,d\sigma(z)}{(1-\xi\ov z)^2(1-\ov\xi z)^{2(p-2)/p}}\\
=\rho^{2(p-1)/p}\frac{f(\xi)}{(1-|\xi|^2)^{2(p-2)/p}}
\end{multline*}
(here we use the fact that the Bergman kernel $(1-\ov\xi z)^{-2}$ is the reproducing kernel on $A_p$). Thus
\begin{equation}\label{12}
\frac{(1-|\xi|^2)^{2(p-2)/p}}{\rho^{2(p-1)/p}}\frac1\pi\int_D\ov{g(z)}|g(z)|^{p-2}f(z)\,d\sigma(z)
=f(\xi)-S_0f(z_1).
\end{equation}

Now let us verify that $|g(z_1)|/\|g\|_{A_p}\le\delta$. Substituting $f=g$ in \eqref{12}, we
obtain
$$(1-|\xi|^2)^{2(p-2)/p}\|g\|_{A_p}^p=\frac{\rho^{2/p}}{(1-|\xi|^2)^{4/p}}$$
which yields
$$\frac{|g(z_1)|}{\|g\|_{A_p}}=\left(\frac{1-\rho^2}{1-|z_1|^2}\right)^{2/p}=\delta_2\le\delta.$$
The proposition is proved.
\end{proof}

Now we consider the same problem for $X=h_p$, $p>1$. Put
$$\alpha(\lambda)=\frac{\displaystyle\dfrac1{2\pi}\int_0^{2\pi}P(z_1,\ei)(P(\xi,\ei)-\lambda P(z_1,\ei))_{(q)}\,d\theta}{\|P(\xi,\cdot)-\lambda P(z_1,\cdot)\|_q^{q-1}},\quad\delta_1=\alpha(0),$$
where $P(\xi,z)=\dfrac{1-|\xi|^2}{|1-\ov\xi z|^2}$ is the Poisson kernel, $\dfrac1p+\dfrac1q=1$, $(x)_{(q)}=|x|^{q-1}\sign x$, and
$$\|f\|_q=\left(\frac1{2\pi}\int_0^{2\pi}|f(\ei)|^q\,d\theta\right)^{1/q}.$$
Let us show that for every $0\le\delta\le\delta_1$ the equation
\begin{equation}\label{13}
\alpha(\lambda)=\delta
\end{equation}
has a solution $\lambda\in\left[0,\dfrac{1+\rho}{1-\rho}\right]$. For $z=\ei$ and $\zeta=W(z)=e^{i\varphi}\dfrac{z-z_1}{1-\ov z_1z}$ we have
\begin{equation}\label{14}
\frac{P(\xi,z)}{P(z_1,z)}=P(\rho,\zeta)=\frac{1-\rho^2}{|1-\rho\zeta|^2}\le\frac{1-\rho}{1+\rho}
\end{equation}
($\rho$ and $\varphi$ are the same as in \eqref{6}). Thus for every $z=\ei$
$$P(\xi,z)-\frac{1+\rho}{1-\rho}P(z_1,z)\le0.$$
Therefore
$$\alpha\left(\frac{1+\rho}{1-\rho}\right)<0.$$
Since $\alpha(\lambda)$ is continuous for $\lambda\in\left[0,\dfrac{1+\rho}{1-\rho}\right]$ and $\alpha(0)=\delta_1$, the equation \eqref{13} has a solution in this interval for every $0\le\delta\le\delta_1$. We denote this solution $C_p(\xi,z_1,\delta)$. For $\delta>\delta_1$ we put $C_p(\xi,z_1,\delta)=0$.

\begin{proposition}\label{P3}
For $X=h_p$, $p>1$,
\begin{itemize}
\item[i)]$S_0y=C_p(\xi,z_1,\delta)y$ is an optimal algorithm,
\item[ii)]
$$u_0(\zeta)=\frac{\displaystyle\dfrac1{2\pi}\int_0^{2\pi}P(\zeta,\ei)(P(\xi,\ei)-
C_p(\xi,z_1,\delta)P(z_1,\ei))_{(q)}\,d\theta}{\|P(\xi,\cdot)-C_p(\xi,z_1,\delta)
P(z_1,\cdot)\|_q^{q-1}}$$
is a worst function,
\item[iii)]the intrinsic error is
$$E(\xi,z_1,\delta,h_p)=u_0(\xi)=\|P(\xi,\cdot)-C_p(\xi,z_1,\delta)P(z_1,\cdot)\|_q+\delta
C_p(\xi,z_1,\delta).$$
\end{itemize}
\end{proposition}

\begin{proof}
It is known (see \cite{10}) that every function from $h_p$, $p>1$, has boundary values almost everywhere. It is also known that boundary values reconstruct this function by Poisson transformation. So for every $u\in Bh_p$, $p>1$, we have by the H\"older inequality
\begin{multline*}
|u(\xi)-C_p(\xi,z_1,\delta)u(z_1)|\\=\biggl|\frac1{2\pi}\int_0^{2\pi}(P(\xi,\ei)-
C_p(\xi,z_1,\delta)P(z_1,\ei))u(\ei)\,d\theta\biggr|\\
\le\|P(\xi,\cdot)-C_p(\xi,z_1,\delta)P(z_1,\cdot)\|_q.
\end{multline*}
On the other hand, the function
$$f(\theta)=\frac{(P(\xi,\ei)-C_p(\xi,z_1,\delta)P(z_1,\ei))_{(q)}}
{\|P(\xi,\cdot)-C_p(\xi,z_1,\delta)P(z_1,\cdot)\|_q^{q-1}}\in BL_p(0,2\pi)$$
and therefore (see \cite{10}) $u_0\in Bh_p$, and has almost everywhere boundary values which coincide with $f(\theta)$. Thus we obtain that
$$u_0(\xi)-C_p(\xi,z_1,\delta)u_0(z_1)=\|P(\xi,\cdot)-C_p(\xi,z_1,\delta)P(z_1,\cdot)\|_q.$$
Hence
\begin{equation}\label{15}
\sup_{u\in Bh_p}|u(\xi)-C_p(\xi,z_1,\delta)u(z_1)|=u_0(\xi)-C_p(\xi,z_1,\delta)u_0(z_1).
\end{equation}
Let $0\le\delta\le\delta_1$. It follows from the definition of $C_p(\xi,z_1,\delta)$ that $u_0(z_1)=\delta$. Since $C_p(\xi,z_1,\delta)\ge0$, we have from \eqref{15} that $u_0(\xi)\ge0$ and
$$\sup_{u\in Bh_p}|u(\xi)-C_p(\xi,z_1,\delta)u(z_1)|=|u_0(\xi)|-\delta C_p(\xi,z_1,\delta).$$
For $\delta>\delta_1$ the same equality holds because $C_p(\xi,z_1,\delta)=0$. Now the proposition follows from Corollary~\ref{C1}.
\end{proof}

We can easily find $C_\infty(\xi,z_1,\delta)$. In this case $q=1$ and
$$\alpha(\lambda)=\frac1{2\pi}\int_0^{2\pi}P(z_1,\ei)\sign(P(\xi,\ei)-\lambda P(z_1,\ei))\,d\theta.$$
In view of \eqref{14} substitution 
$$z=\frac{e^{-i\varphi}\zeta+z_1}{1+\ov z_1e^{-i\varphi}\zeta}$$ 
yields
\begin{multline*}
\alpha(\lambda)=\frac1{2\pi}\int_0^{2\pi}\sign(P(\rho,\ei)-\lambda)\,d\theta\\
=
\frac1\pi\int_0^\pi\sign\left(\frac{1-\rho^2}{1-2\rho\cos\theta+\rho^2}-\lambda\right)\,d\theta\\
=\begin{cases}1,&\lambda\le\dfrac{1-\rho}{1+\rho},\\
\dfrac2\pi\arccos\dfrac{\lambda(1+\rho^2)-(1-\rho^2)}{2\rho\lambda}-1,&\dfrac{1-\rho}{1+\rho}
\le\lambda\le\dfrac{1+\rho}{1-\rho},\\
-1,&\lambda\ge\dfrac{1+\rho}{1-\rho}.\end{cases}
\end{multline*}
Hence for $0\le\delta<1$ the solution of \eqref{13} is
$$C_\infty(\xi,z_1,\delta)=\frac{1-\rho^2}{1+2\rho\sin(\pi\delta/2)+\rho^2}.$$
If $\delta=1$, every $\lambda\in\left[0,\dfrac{1+\rho}{1-\rho}\right]$ is a solution of \eqref{13}.

For $0\le\delta<1$ and $z=\ei$ we have
\begin{multline*}
\sign(P(\xi,z)-C_\infty(\xi,z_1,\delta)P(z_1,z))\\=\sign\left(\frac{1-\rho^2}{|1-\rho W(z)|^2}-
C_\infty(\xi,z_1,\delta)\right)
=\sign\left(\RE W(z)+\sin\frac\pi2\delta\right)\\=\sign\RE\frac{W(z)+\tan(\pi\delta/4)}{1+
W(z)\tan(\pi\delta/4)}
=\frac4\pi\RE\arctan\frac{W(z)+\tan(\pi\delta/4)}{1+
W(z)\tan(\pi\delta/4)}.
\end{multline*}
In the case $p=\infty$, $0\le\delta<1$, we obtain
$$u_0(\zeta)=\frac4\pi\RE\arctan\frac{W(\zeta)+\tan(\pi\delta/4)}{1+
W(\zeta)\tan(\pi\delta/4)}.$$
Thus the next corollary follows from Proposition~\ref{P3}.

\begin{corollary}\label{C3}
For $X=h_\infty$
\begin{itemize}
\item[i)]
$$S_0y=\begin{cases}\dfrac{1-\rho^2}{1+2\rho\sin(\pi\delta/2)+\rho^2}y,&0\le\delta<1,\\
c\dfrac{1-\rho}{1+\rho}y,&\delta=1,\ c\in[0,1],\\
0,&\delta>1,\end{cases}$$
is an optimal algorithm $($in the case $\delta=1$, $c$ is an arbitrary value in $[0,1]$$)$.
\item[ii)]
$$u_0(z)=\frac4\pi\RE\arctan\frac{W(z)+\Delta}{1+\Delta W(z)},$$
where
$$\Delta=\begin{cases}\tan(\pi\delta/4),&0\le\delta<1,\\
1,&\delta\ge1,\end{cases}$$
is a worst function,
\item[iii)]the intrinsic error is
$$E(\xi,z_1,\delta,h_\infty)=u_0(\xi)=\frac4\pi\arctan\frac{\rho+\Delta}{1+\Delta\rho}.$$
\end{itemize}
\end{corollary}

The solution of the equation \eqref{13} may also be obtained in direct form for $p=2$. Nevertheless, we shall prove a more general result for the Hilbert space.

Let $X$ be a complex (or real) Hilbert space. Consider the problem \eqref{1} in the case $W=BX$, $Y=Z=\mathbb C\ (\mathbb R)$, $Lx=(x,x_1)$, $Ix=(x,x_2)$, $x_1,x_2\in X$. The intrinsic error will be denoted by $E(x_1,x_2,\delta,X)$.

\begin{proposition}\label{P4}
Let $x_1$ and $x_2$ be linear independent elements from the Hilbert space $X$. Put
$$\varepsilon=\min\left\{\delta,\frac{|(x_1,x_2)|}{\|x_1\|}\right\}.$$
Then
\begin{itemize}
\item[i)] $S_0y=\lambda\dfrac{(x_2,x_1)}{\|x_2\|^2}y$, where
$$\lambda=1-\frac\varepsilon{|(x_1,x_2)|}\sqrt{\frac{\|x_1\|^2\|x_2\|^2-|(x_1,x_2)|^2}{\|x_2\|^2
-\varepsilon^2}},$$
is an optimal algorithm,
\item[ii)]
$$x_0=\sqrt{\frac{\|x_2\|^2-\varepsilon^2}{\|x_1\|^2\|x_2\|^2-|(x_1,x_2)|^2}}\left(x_1-\lambda
\frac{(x_1,x_2)}{\|x_2\|^2}x_2\right)$$
is a worst element,
\item[iii)]the intrinsic error is
$$E(x_1,x_2,\delta,X)=\sqrt{1-\frac{\varepsilon^2}{\|x_2\|^2}}\sqrt{\|x_1\|^2-
\frac{|(x_1,x_2)|^2}{\|x_2\|^2}}+\varepsilon\frac{|(x_1,x_2)|}{\|x_2\|^2}.$$
\end{itemize}
\end{proposition}

\begin{proof}
We have
\begin{multline*}
\sup_{\|x\|\le1}|(x,x_1)-S_0(x,x_2)|=\sup_{\|x\|\le1}\left|\left(x,x_1-\lambda
\frac{(x_1,x_2)}{\|x_2\|^2}x_2\right)\right|\\
=\left\|x_1-\lambda\frac{(x_1,x_2)}{\|x_2\|^2}x_2\right\|=\sqrt{\frac{\|x_1\|^2\|x_2\|^2-|(x_1,x_2)|^2}{\|x_2\|^2
-\varepsilon^2}}.
\end{multline*}
Moreover
$$(x_0,x_1)-S_0(x_0,x_2)=\left(x_0,x_1-\lambda
\frac{(x_1,x_2)}{\|x_2\|^2}x_2\right)=\left\|x_1-\lambda\frac{(x_1,x_2)}{\|x_2\|^2}x_2\right\|.$$
Thus $\|x_0\|=1$ and 
$$\sup_{\|x\|\le1}|(x,x_1)-S_0(x,x_2)|=(x_0,x_1)-S_0(x_0,x_2).$$ 
Since
$$S_0(x_0,x_2)=\frac{\lambda(1-\lambda)|(x_1,x_2)|^2}{\|x_2\|^2
\left\|x_1-\lambda\dfrac{(x_1,x_2)}{\|x_2\|^2}x_2\right\|}\ge0$$
and
$$|(x_0,x_2)|=\sqrt{\frac{\|x_2\|^2-\varepsilon^2}{\|x_1\|^2\|x_2\|^2-|(x_1,x_2)|^2}}|(x_1,x_2)|
(1-\lambda)=\varepsilon,$$
we obtain $S_0(x_0,x_2)=\delta\|S_0\|$. To finish the proof of the proposition we should only apply Corollary~\ref{C1}.
\end{proof}

The problems of optimal recovery in Hilbert spaces from inaccurate data were investigated in \cite{11} (see also \cite{2} for a more general situation).

Let $X$ be a Hilbert space of functions $f\colon\Omega\to\mathbb C\ (\mathbb R)$ with the reproducing kernel $K\colon\Omega\times\Omega\to\mathbb C\ (\mathbb R)$, i.e.,
$$f(z)=(f\cd,K(\cdot,z))$$
for every $f\in X$ and $z\in\Omega$. Consider the problem \eqref{1} for $W=BX$, $Lf=f(\xi)$, $If=f(z_1)$. If $K(\cdot,\xi)$ and $K(\cdot,z_1)$ are linearly independent (i.e., the class $BX$
distinguishes the points $\xi$ and $z_1$), then from Proposition~\ref{P4} we get 

\begin{corollary}\label{C4}
Put
$$\varepsilon=\min\left\{\delta,\frac{|K(z_1,\xi)|}{\sqrt{K(\xi,\xi)}}\right\}.$$
Then
\begin{itemize}
\item[i)] $S_0y=\lambda\dfrac{K(\xi,z_1)}{K(z_1,z_1)}y$, where
$$\lambda=1-\frac\varepsilon{|K(z_1,\xi)|}\sqrt{\frac{K(\xi,\xi)K(z_1,z_1)-|K(z_1,\xi)|^2}
{K(z_1,z_1)-\varepsilon^2}}$$
is an optimal algorithm,
\item[ii)]
\begin{multline*}
f_0(z)=\sqrt{\frac{K(z_1,z_1)-\varepsilon^2}{K(\xi,\xi)K(z_1,z_1)-|K(z_1,\xi)|^2}}\\
\times\left(K(z,\xi)-\lambda\frac{K(z_1,\xi)}{K(z_1,z_1)}K(z,z_1)\right)
\end{multline*}
is a worst function, and
\item[iii)]the intrinsic error is
$$E(\xi,z_1,\delta,X)=\sqrt{1-\frac{\varepsilon^2}{K(z_1,z_1)}}\sqrt{K(\xi,\xi)-
\frac{|K(z_1,\xi)|^2}{K(z_1,z_1)}}+\varepsilon\frac{|K(z_1,\xi)|}{K(z_1,z_1)}.$$
\end{itemize}
\end{corollary}

We list some examples of Hilbert spaces with reproducing kernels:
\begin{enumerate}
\renewcommand{\labelenumi}{\theenumi)}
\item $H_2$, $\quad K(\xi,z)=(1-\xi\ov z)^{-1}$, $\displaystyle\quad(f,g)=\dfrac1{2\pi}\int_0^{2\pi}f(\ei)\ov{g(\ei)}\,d\theta$,
\item $A_2$, $\quad K(\xi,z)=(1-\xi\ov z)^{-2}$, $\displaystyle\quad (f,g)=\dfrac1\pi\int_Df(z)\ov{g(z)}\,d\sigma(z)$,
\item $h_2$, $\quad K(\xi,z)=2\RE(1-\xi\ov z)^{-1}-1$,
\item[]\hspace{152pt}$\displaystyle\quad(u,v)=\dfrac1{2\pi}\int_0^{2\pi}u(\ei)v(\ei)\,d\theta$,
\item $a_2$, $\quad K(\xi,z)=2\RE(1-\xi\ov z)^{-2}-1$,
\item[]\hspace{163pt}$\displaystyle\quad (u,v)=\dfrac1\pi\int_Du(z)v(z)\,d\sigma(z)$.
\end{enumerate}

Note that we can obtain the generalization of Schwarz’s Lemma in the same way as \eqref{9}:
$$\sup_{\substack{f\in BX\\|f(z_1)|\le\delta}}|f(\xi)|=E(\xi,z_1,\delta,X),$$
where $E(\xi,z_1,\delta,X)$ can be found from the corresponding proposition for
$X=H_p,A_p,h_p$, and $a_2$.

Denote
$$D(\xi,z_1,\delta,X)=\{z\in D:|g_0(z)|\le\delta\},$$
where $X=H_p,A_p,h_p$, or $a_2$ and $g_0(z)$ is a worst function for the appropriate recovery problem. Consider the information operator $\widetilde If=f_{|D(\xi,z_1,\delta,X)}$ instead of $If=f(z_1)$ and the space $Y$ is a space of continuous in $D(\xi,z_1,\delta,X)$ functions with the norm
$$\|y\|=\sup_{z\in D(\xi,z_1,\delta,X)}|y(z)|.$$
It follows from Corollary~\ref{C1} that the optimal algorithm, the worst function, and the intrinsic error will be stay the same. Thus the additional information (with the same inaccuracy) about the behaviour of the function $f$ in $D(\xi,z_1,\delta,X)$ will not decrease the intrinsic error. In other words, the point $z_1$ forms some ``shadow'' set in which any additional information is useless.

3. \underline{Optimal recovery of the derivative from inaccurate data}. We turn now to the problem \eqref{1} for $X=H_p$, $Z=\mathbb C$, $Y=l_q^2$, $Lf=f'(0)$, $If=(f(-h),f(h))$, $h\in(0,1)$. The intrinsic error will be denoted by $E'_q(h,\delta,H_p)$.

There is the well-known algorithm
$$f'(0)\approx\frac{f(h)-f(-h)}{2h}$$
which is not optimal even in the case $\delta=0$ (see \cite{12}). It was shown in \cite{2} that
$$f'(0)\approx(1-h^4)\frac{f(h)-f(-h)}{2h}$$
is an optimal algorithm for $\delta=0$ and $p=\infty$. It follows from \cite{8} that this
algorithm is an optimal one for $\delta=0$ and every $1\le p\le\infty$. Moreover it is also optimal for $\delta=0$ and $X=h_\infty$ (see \cite{9}).

Now we consider the case when the value of functions in the points $-h$ and $h$ are known with an error $\le\delta$ in the norm of $l_q^2$, that is, we know $y_1$ and $y_2$ such that
\begin{gather*}
|f(-h)-y_1|^q+|f(h)-y_2|^q\le\delta^q,\quad1\le q<\infty,\\
\max\{|f(-h)-y_1|,|f(h)-y_2|\}\le\delta,\quad q=\infty.
\end{gather*}
Put
\begin{gather*}
\varepsilon_p=\begin{cases}1/p,&1\le p<\infty,\\
0,&p=\infty,\end{cases}\quad\delta_1=h2^{\varepsilon_q-\varepsilon_p}(1+h^2)^{-\varepsilon_p},
\quad\delta_2=h2^{\varepsilon_q},\\
\alpha(z)=\begin{cases}1,&0\le\delta<\delta_1,\\
\dfrac{a^2-z^2}{1-a^2z^2},&\delta\ge\delta_1\end{cases}
\end{gather*}
Let $a\in[h,1]$ be a solution of the equation
\begin{equation}\label{16}
\frac{h(a^2-h^2)(1-a^2h^2)^{2\varepsilon_p-1}}{\alpha(h)(1-h^4)^{\varepsilon_p}
(1-2a^2h^2+a^4)^{\varepsilon_p}}=\delta2^{-\varepsilon_q},
\end{equation}
where $0\le\delta\le\delta_2$ (the existence of the solution follows from the continuity
of the function in the left hand side of this equation). Put $a=h$ for $\delta>\delta_2$.

\begin{proposition}\label{P5}
For every $\delta\ge0$, $1\le p,q\le\infty$
\begin{itemize}
\item[i)]
$$f'(0)\approx S_0y=\frac{\alpha(h)(1-a^2h^2)^{2(1-\varepsilon_p)}}{\alpha(0)(1-h^4)^{1-2\varepsilon_p}}
\frac{y_2-y_1}{2h}$$
is an optimal algorithm,
\item[ii)]
$$g_0(z)=\left(\frac{1-h^4}{1-2a^2h^2+a^4}\right)^{\varepsilon_p}\frac{z(a^2-z^2)
(1-a^2z^2)^{2\varepsilon_p-1}}{\alpha(z)(1-h^2z^2)^{2\varepsilon_p}}$$
is a worst function,
\item[iii)]the intrinsic error is
$$E'_q(h,\delta,H_p)=\frac{a^2}{\alpha(0)}\left(\frac{1-h^4}{1-2a^2h^2+a^4}\right)^{\varepsilon_p}
.$$
\end{itemize}
\end{proposition}

\begin{proof}
Put
$$g(z)=\frac{z(a^2-z^2)(1-a^2z^2)^{2\varepsilon_p-1}}{\alpha(z)(1-h^2z^2)^{2\varepsilon_p}},
\quad\varphi(z)=\left(\frac{1-a^2z^2}{1-h^2z^2}\right)^2.$$
For every $f\in H_p$ and $z=\ei$ we obtain by the residue theorem
\begin{multline*}
f'(0)-S_0If=-\frac{h^2}{\alpha(0)}\frac1{2\pi i}\int_{|z|=1}\frac{\alpha(z)
(1-a^2z^2)^{2(1-\varepsilon_p)}}{z^2(z^2-h^2)(1-h^2z^2)^{1-2\varepsilon_p}}f(z)\,dz\\
=\begin{cases}\displaystyle\frac{h^2}{\alpha(0)}\frac1{2\pi}\int_0^{2\pi}\ov{g(z)}
|g(z)|^{p-2}f(z)\,d\theta,&1\le p<\infty,\\[10pt]
\displaystyle\frac{h^2}{\alpha(0)}\frac1{2\pi}\int_0^{2\pi}\ov{g(z)}
|\varphi(z)|f(z)\,d\theta,&p=\infty.\end{cases}
\end{multline*}
For $f=g$ we have from these equations
$$\|g\|_{H_p}^p=\|\varphi\|_{H_1}=\frac{1-2a^2h^2+a^4}{1-h^4}.$$
Note that if $p=\infty$, $|g(\ei)|\equiv1$. Now to use Corollary~\ref{C2} we must prove that
$$Ig_0=\delta a^*=\frac\delta{2^{\varepsilon_q}}(-1,1)$$
if $0\le\delta<\delta_2$ and $\|Ig_0\|_{l_q^2}\le\delta$ if $\delta\ge\delta_2$. Let $0\le\delta<\delta_2$. Since $g_0(-h)=-g_0(h)$ it is sufficient to prove the equation
$$g_0(h)=\delta2^{-\varepsilon_q}$$
which coincides with \eqref{16}. If $\delta\ge\delta_2$, $g_0(z)=z$ and $\|Ig_0\|_{l_q^2}=h2^{\varepsilon_q}=\delta_2\le\delta$. This completes the proof of the proposition.
\end{proof}

Note that $S_0y\equiv0$ for $\delta\ge2^{\varepsilon_q}$. If $\delta<2^{\varepsilon_q}$ we can consider the problem of finding such a value $h_0$, that
$$E'_q(h_0,\delta,H_p)=\min_{h\in(0,1)}E'_q(h,\delta,H_p).$$
The value $h_0$ is called an optimal value of $h$. We shall give the solution of this problem for $p=\infty$.

\begin{proposition}\label{P6}
For $p=\infty$, $1\le q\le\infty$, and $0\le\delta<2^{\varepsilon_q}$ the optimal value $h_0$ satisfies the equation
\begin{equation}\label{17}
\delta h_0^4+2^{1+\varepsilon_q}h_0^3-\delta^22^{1-\varepsilon_q}h_0-\delta=0.
\end{equation}
The equality
$$\min_{h\in(0,1)}E'_q(h,\delta,H_\infty)=h_0^2$$
holds. The optimal value $h_0$ can also be found from the equality
$$h_0=\sqrt k\sn(K/3,k),$$
where $k$ is determined by the equation
\begin{equation}\label{18}
\sqrt k=2h_1^{1/4}\frac{\displaystyle\sum_{m=0}^\infty h_ 1^{m(m+1)}}{\displaystyle1+
2\sum_{m=1}^\infty h_1^{m^2}},\quad h_1=e^{\textstyle-\frac{\pi\Lambda'}{3\Lambda}}
\end{equation}
or
\begin{equation}\label{19}
\frac{K'}K=\frac{\Lambda'}{3\Lambda};
\end{equation}
here $K$, $\Lambda$ denote the complete elliptic integrals of the first kind with respective moduli $k$, $\lambda=\delta^24^{-\varepsilon_q}$ and $K'$, $\Lambda'$ denote the ones with complementary moduli.
\end{proposition}

\begin{proof}
From Proposition~\ref{P5} we have
$$E'_q(h_0,\delta,H_\infty)=\begin{cases}a^2,&0\le\delta<h2^{\varepsilon_q},\\
1,&\delta\ge h2^{\varepsilon_q},\end{cases}$$
where $a\in[h,1]$ and determines by the equation
\begin{equation}\label{20}
h\frac{a^2-h^2}{1-a^2h^2}=\delta2^{-\varepsilon_q}.
\end{equation}
Extracting $a^2$ from this equation and minimizing it as a function of $h\in(0,1)$ we obtain that the minimum $h_0$ is unique and satisfies the equation \eqref{17}. Taking a derivative from \eqref{20}, we have
$$\frac{a^2-h^2}{1-a^2h^2}-2h^2\frac{1-a^4}{(1-a^2h^2)^2}+2aa'h\frac{1-h^4}{(1-a^2h^2)^2}=0.$$

Thus if $h_0$ is minimum then $g_0'(h_0)=0$, where
\begin{equation}\label{21}
g_0(z)=z\frac{a^2-z^2}{1-a^2z^2}.
\end{equation}
Now it is sufficient to find a function $g_0(z)$ of kind \eqref{21} such that for some $h_0\in(0,1)$, $g_0(h_0)=\delta2^{-\varepsilon_q}$ and $g_0'(h_0)=0$. It follows from Lemma~2.2 of \cite{7} that this function is a Blaschke product of order $3$ with minimal norm
$$\|g_0\|=\max_{z\in[-\sqrt k,\sqrt k]}|g_0(z)|=\delta2^{-\varepsilon_q},$$
where $k$ is determined by the conditions $|g_0(-\sqrt k)|=|g_0(\sqrt k)|=\delta2^{-\varepsilon_q}$. From \cite{13} this function can be written in the form
$$g_0(z)=z\frac{k\sn^2(2K/3,k)-z^2}{1-k\sn^2(2K/3,k)z^2}.$$
This function can be rewritten by using the first fundamental transformation of degree $3$ (see \cite{14})
$$g_0(z)=\sqrt\lambda\sn\left(3\Lambda u/K,\lambda\right),\quad z=\sqrt k\sn(u,k),$$
where $\lambda=\delta^24^{-\varepsilon_q}$ and $k$ satisfies \eqref{18}, \eqref{19}. If we put $h_0=\sqrt k\sn(K/3,k)$ then $g_0(h_0)=\delta2^{-\varepsilon_q}$ and $g_0'(h_0)=0$. This completes the proof of the proposition.
\end{proof}

It is easily shown from \eqref{17} that
$$h_0=2^{-\textstyle\frac{1+\varepsilon_q}3}\delta^{1/3}+O\left(\delta^{5/3}\right),$$
and consequently
$$\min_{h\in(0,1)}E'_q(h_0,\delta,H_\infty)=4^{-\textstyle\frac{1+\varepsilon_q}3}\delta^{2/3}+
O\left(\delta^2\right).$$


\begin{thebibliography}{99}

\bibitem{1} A.G.Marchuk and K.Yu.Osipenko, Best approximation of functions specified with an error at a finite number of points, Mat. Zameki, 17 (1975), 359--368; English transl. in Math. Notes 17 (1975).

\bibitem{2} C.A.Micchelli and T.J.Rivlin, A survey of optimal recovery, Optimal Estimation in Approximation Theory, C.A.Micchelli, T.J.Rivlin, eds.), Plenum Press, 1977, pp.~l--53.

\bibitem{3} C.A.Micchelli and T.J.Rivlin, Lectures on Optimal Recovery, Lect. Notes Math., 1129 (1985), 21--93.

\bibitem{4} T.J.Rivlin, A survey of recent results in optimal recovery, Polynomial and Spline Approx. Proc., NATO Adv. Study Inst., Calgary, 1978, pp.~225--245 (Dordrecht e.a., 1979)  .

\bibitem{5} K.Yu.Osipenko, Best methods of approximating analytic functions given with an error, Mat. Sb., 118 (160) (1982), 350--370; English transl. in Math. USSR-Sb., 46 (1983).

\bibitem{6}	K.Yu.Osipenko, The Heins problem and optimal extrapolation of analytic functions given with an error, Mat. Sb., 126 (168) (1985), 566--575.

\bibitem{7}	K.Yu.Osipenko, On optimal extrapolation and interpolation of fuzzy analytic functions, Anal. Math., 13 (1987), 199--210.

\bibitem{8} K.Yu.Osipenko and M.I.Stessin, On recovery problems in the Hardy and Bergman spaces, Mat. Zametki, to appear.

\bibitem{9} K.Yu.Osipenko, Best and optimal methods of recovery on classes of harmonic functions, Mat. Sb., to appear.

\bibitem{10} G.M.Goluzin, Geometric theory of functions of a complex variable, Translations of Mathematical Monographs, 26 (American Mathematical Society, Providence, Rhode Island, 1969).

\bibitem{11} A.A.Melkman and C.A.Micchelli, Optimal estimation of linear operators in Hilbert spaces from inaccurate data, SIAM J. Numer. Anal., 16 (1979), 87--105.

\bibitem{12} T.J.Rivlin, The optimal recovery of functions, Comemp. Mafh., 9 (1982), 121--151.

\bibitem{13} K.Yu.Osipenko, Optimal interpolation of analytic functions, Mat. Zametki, 12 (1972), 465--476; English transl. in Math. Notes 12 (1972).

\bibitem{14} N.I.Akhiezer, ``Elements of the theory of elliptic functions, 2nd rev. ed., ``Nauka'', Moscow, 1970 (in Russian).
\end{thebibliography}
\end{document}

