\documentclass[reqno]{amsart} \usepackage{graphicx} \AtBeginDocument{{\noindent\small {\em Electronic Journal of Differential Equations}, Vol. 2004(2004), No. 04, pp. 1--23.\newline ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu \newline ftp ejde.math.txstate.edu (login: ftp)} \thanks{\copyright 2004 Texas State University - San Marcos.} \vspace{9mm}} \begin{document} \title[\hfilneg EJDE-2004/04\hfil Oleck-Opial-Beesack-Troy inequalities] {On Oleck-Opial-Beesack-Troy integro-differential inequalities} \author[Evgeniy I. Bravyi \& Sergey S. Gusarenko\hfil EJDE-2004/04\hfilneg] {Evgeniy I. Bravyi \& Sergey S. Gusarenko} % in alphabetical order \address{Evgeniy I. Bravyi \hfill\break Research Centre ``Functional Differential Equations'', Perm State Technical University, Perm 614099, Russia} \email{bravyi@pi.ccl.ru} \address{Sergey S. Gusarenko \hfill\break Department of Mathematical Analysis, Perm State University, Perm 614099, Russia} \email{gusarenko@psu.ru} \date{} \thanks{Submitted October 27, 2003. Published January 2, 2004.} \thanks{Supported by grant 03-01-00255 from the Russian Foundation for Basic Research and by \hfill\break\indent the Program ``Universities of Russia" UR.04.01.001.} \subjclass[2000]{34K10, 34B30, 34B05, 41A44, 49J40, 58E35} \keywords{Integral inequalities, integro-differential inequalities, \hfill\break\indent functional differential equations, variational problems} \begin{abstract} We obtain necessary and sufficient conditions for the integro-differential inequality $$ \int_a^b\dot x^2(t)\,dt\geq\gamma\int_a^bq(t)\,|\dot x(t)x(t)|\,dt $$ to be valid with one of the three boundary conditions: $x(a)=0$, or $x(b)=0$, or $x(a)=x(b)=0$. For a power functions $q$, the best constants $\gamma$ are found. \end{abstract} \maketitle \numberwithin{equation}{section} \newtheorem{lem}{Lemma}[section] \newtheorem{thm}{Theorem}[section] \newtheorem{cor}{Corollary}[section] \section*{Introduction} Beesack \cite{Beesack} found that the best constant is $\overline{\gamma}=\frac2{b-a}$, such that the inequality \begin{equation}\label{uB} \int_a^b\dot x^2(t)\,dt\ge\gamma\int_a^b|\dot x(t)x(t)|\,dt \end{equation} holds for all $\gamma\le\overline{\gamma}$ and for all continuously differentiable functions $x:[a,b]\to\mathbb{R}$ satisfied the condition $x(a)=0$ (or $x(b)=0$). For inequality \eqref{uB} with two constrictions $x(a)=x(b)=0$ the best constant is $\overline{\gamma}=\frac4{b-a}$ which was found by Olech and Opial in \cite{Oleck, Opial}. Troy \cite{WTroy} obtained some estimates of the best constant $\gamma$ for the inequality \begin{equation}\label{uT} \int_a^b\dot x^2(t)\,dt\ge\gamma\int_a^b (t-a)^p\,|\dot x(t)x(t)|\,dt \end{equation} with the condition $x(a)=0$ (and also with $x(a)=x(b)=0$). But the question on the sharpness of these estimates remained open. The aim of the paper is to answer this question and find the best constants for extended inequalities with different kinds of the boundary conditions. The Perm Seminar on Functional Differential Equations has developed methods for the investigation nonclassical variational problems (see, for example, \cite{AMR2, Gusarenko1}). These methods allow to prove some known integral and integro-differential inequalities more effectively and obtain new inequalities, for example, for integrals with deviating arguments \cite{Grusia, Vestnik}. Here we apply these methods to the integro-differential inequality \begin{equation}\label{uG} \int_a^b\dot x^2(t)\,dt\ge\gamma\int_a^bq(t)\,|\dot x(t)x(t)|\,dt \end{equation} which generalizes the inequalities from \cite{Beesack,Oleck,Opial,WTroy}. We obtain estimates of such constants $\gamma$ for which \eqref{uG} is valid. We also compute the best constants in some cases including inequality \eqref{uT}. In particular, Troy have proved that \eqref{uG} holds for all continuously differentiable functions $x:[a,b]\to\mathbb{R}$ such that $x(a)=0$ if \begin{equation}\label{etroy} \gamma\le\frac{2\sqrt{p + 1}}{(b-a)^{p+1}},\quad p>-1. \end{equation} We prove that \eqref{uG} holds for all absolutely continuous functions $x:[a,b]\to\mathbb{R}$ such that $x(a)=0$ if and only if \begin{equation*} \gamma\le\overline{\gamma}_p=\frac{(p+1)^2\eta_p^2}{2|p|(b-a)^{p+1}},\quad p\ne0,\quad p>-1, \end{equation*} where $\eta_p$ is the smallest positive root of the Bessel function $J_{-\frac{1+2p}{1+p}}$ for $p<0$ and the smallest positive root of the modified Bessel function $I_{-\frac{1+2p}{1+p}}$ for $p>0$. In section 1.3.a some estimates for $\overline{\gamma}_p$, which are better than \eqref{etroy}, are obtained. Moreover, in this section the best constants $\overline{\gamma}_p$ are evaluated explicitly for some cases. The singular case $p=-1$ is considered in section 1.3.b. The method is based on the reduction of the integro-differential inequality to the problem of the minimization for the quadratic functional \begin{equation}\label{e} \mathcal{J}(z)=\int_a^b\bigl(z(t)-(Kz)(t)\bigr)z(t)\,dt\to\min \end{equation} in the space $\mathbf{L}_2$ of square integrable functions $z:[a,b]\to\mathbb{R}$ with the norm $\|z\|=\sqrt{\int_a^bz^2(t)\,dt}$. The linear operator $K:\mathbf{L}_2\to\mathbf{L}_2$ in \eqref{e} is supposed to be self-adjoint and bounded. It is clear that either a function $z\equiv0$ is a point of the a minimum of $\mathcal{J}$ and problem \eqref{e} is solvable ($\min_{z\in\mathbf{L}_2}\mathcal{J}(z)=0$) or $\mathcal{J}$ is unbounded below on $\mathbf{L}_2$. It's well known that problem \eqref{e} is solvable if and only if all points of the spectrum of the operator $K$ is not greater than one. Moreover, the following statement is valid \cite{Gusarenko2}. \begin{thm}\label{th} Suppose $\int_a^b(Kz)(t)z(t)\,dt\geq 0$ for all nonnegative functions $z\in\mathbf{L}_2$. Then problem \eqref{e} is solvable if and only if the norm (the spectral radius) of the operator $K$ is less than or equal to one. \end{thm} So if the conditions of Theorem \ref{th} are fulfilled, the question on the solvability of \eqref{e} is reduced to the computation or the estimation of $\|K\|$. The norm is equal to the spectral radius for self-adjoin operators in $\mathbf{L}_2$. If $K$ is completely continuous, then the norm is equal to the largest eigenvalue of $K$. Note that we can compute the norm (or the spectral radius) of the integral operator $Q:\mathbf{L}_2\to\mathbf{L}_2$ exactly only in rare, extraordinary cases. We can estimate its norm, in particular, with the help of the known inequality: \begin{equation}\label{HA} \|Q\|\le\Big(\int_a^b \int_a^bQ_n^2(t,s)\,ds\,dt\Big)^{1/(2n)}, \end{equation} $n=1,2,\ldots$, where $Q_n(t,s)$ is the kernel of the integral operator $Q^n$. In addition, we use the following helpful result (see, for example, \cite[c.406]{FAN}). \begin{lem}\label{l01} Let $Q:\mathbf{L}_2\to\mathbf{L}_2$ be a completely continuous integral operator with non-negative kernel. If there exists a positive function $v\in{\bf L}_2$ such that \begin{equation}\label{nnnnn} (Qv)(t)\le rv(t)\mbox{ almost everywhere on }[a,b], \end{equation} then $\|Q\|\le r$. \end{lem} \section{Integro-differential inequalities with the condition $x(a)=0$} By $\mathbf{W}_a$ denote the space of absolutely continuous functions $x:[a,b]\to\mathbb{R}$ such that $\dot x\in\mathbf{L}_2$ and $x(a)=0$. Suppose the function $q:[a,b]\to\mathbb{R}$ is measurable, nonnegative on $[a,b]$, and $\int_a^b(t-a)q^2(t)\,dt<\infty$. Our task is to find all $\gamma$ such that for any function $x\in\mathbf{W}_a$ inequality \eqref{uG} holds. Obviously, here and below we must consider only $\gamma>0$, since for all $\gamma\le0$ inequality \eqref{uG} is valid. Note that the functional on $\mathbf{W}_a$ \begin{equation*} x\longmapsto\int_a^b\Bigl(\dot x^2(t)-\gamma q(t)|\dot x(t)x(t)|\Bigr)\,dt \end{equation*} defined by \eqref{uG} is neither quadratic nor differentiable at zero. Therefore at first we reduce the question on the solvability of \eqref{uG} to the minimization problem for a functional of form \eqref{e}. \begin{lem}\label{l1} Inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$ if and only if the zero function is a solution of the variational problem \begin{equation}\label{v} \begin{gathered}J(x)=\int_a^b\Bigl(\dot x^2(t)- \gamma q(t)\dot x(t)x(t)\Bigr)\,dt\to\min,\\ x(a)=0. \end{gathered} \end{equation} \end{lem} \begin{proof} If $x\equiv0$ is not a solution to \eqref{v}, then the functional $J$ is negative for some function $x_0\in\mathbf{W}_a$. Hence, \begin{equation*} \int_a^b(\dot x_0^2(t)-\gamma q(t)|\dot x_0(t)x_0(t)|)\,dt\le\int_a^b(\dot x_0^2(t)-\gamma q(t)\dot x_0(t) x_0(t))\,dt < 0. \end{equation*} So inequality \eqref{uG} does not hold for the function $x_0$. If for a function $x_1\in\mathbf{W}_a$ inequality \eqref{uG} does not hold, then $J(x_2)< 0$ for the function $x_2(t)=\int_a^t |\dot x_1(s)|\,ds$. Indeed, we have $\dot x_2(t)= |\dot x_1(t)|$, $x_2(t) \ge|x_1(t)|$ on $[a,b]$. Therefore, $$ J(x_2)=\int_a^b\Bigl( \dot x_2^2(t)-\gamma q(t)|\dot x_2(t)x_2(t)| \Bigr)\,dt\le \int_a^b\Bigl( \dot x^2_1(t)-\gamma q(t)|\dot x_1(t)x_1(t)| \Bigr)\,dt < 0. $$ Thus inequality \eqref{uG} holds if and only if variational problem \eqref{v} is solvable. \end{proof} By $Q_b:\mathbf{L}_2\to\mathbf{L}_2$ denote the linear integral operator $$ (Q_bz)(t)=\int_a^b Q_b(t,s) z(s) \,ds := q(t)\int_a^tz(s)\,ds+\int_t^bq(s)z(s)\,ds, $$ where \begin{equation*} Q_b(t,s)= \begin{cases} q(t)\mbox{ if } a\le s< t\le b,\\ q(s)\mbox{ if } a\le t\le s\le b. \end{cases} \end{equation*} Since $$ \|Q_b\|^2\le\int_a^b\!\!\int_a^bQ_b^2(t,s)\,dt\,ds=2\int_a^b(t-a)q^2(t)\,dt<\infty, $$ the operator $Q_b$ is completely continuous. Our main result in this section is the following. \begin{thm}\label{thQ1} Inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$ if and only if \begin{equation*} \gamma\le\overline{\gamma}:=2/\|Q_b\|. \end{equation*} \end{thm} \begin{proof} The substitution $x(t)=\int_a^tz(s)\,ds$ reduces problem \eqref{v} to the problem of the form \eqref{e}: \begin{equation*} \int_a^b\Bigl(z^2(t)-\gamma q(t)z(t)\int_a^tz(s)\,ds\Bigr)dt= \int_a^b\bigl(z(t)-(Kz)(t)\bigr)z(t)\,dt\to\min, \end{equation*} where the integral operator $(Kz)(t)=\int_a^bK(t,s)z(s)\,ds$ has the symmetric kernel $K(t,s)=\dfrac\gamma2Q_b(t,s)$. By Theorem \ref{th}, problem \eqref{v} is solvable if and only if $\|K\|=\dfrac\gamma2\|Q_b\|\le1$. \end{proof} Estimating the norm of the operator $Q_b$ by inequality \eqref{HA} for $n=1$, we obtain \begin{cor} \label{s1} If \begin{equation}\label{app} \gamma\le\frac{\sqrt{2}}{\sqrt{\int_a^b(t-a)q^2(t)\,dt}}, \end{equation} then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$. \end{cor} Letting $v=1$ in inequality \eqref{nnnnn}, we obtain \begin{cor}\label{ss2} If \begin{equation}\label{s2} \gamma\le\frac2{\mathop{\rm ess\,sup}\limits_{t\in[a,b]}(q(t)(t-a)+\int_t^bq(s)\,ds)}, \end{equation} then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$. \end{cor} \begin{cor} If the function $q$ is non-increasing and \begin{equation*} \gamma\le\frac{2}{\int_a^bq(t)\,dt}, \end{equation*} then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$. \end{cor} \begin{proof} If the function $q$ is non-increasing, then $Q_b(t,s)\le q(s)$ for all $t\in[a,b]$. Therefore, $\|Q_b\|\le\int_a^bq(s)\,ds$. \end{proof} From Lemma \ref{l01}, we obtain \begin{cor}\label{ns1} If \begin{equation*} \gamma\le\frac{2}{(b-a)\mathop{\rm ess\,sup}\limits_{t\in[a,b]}q(t)}, \end{equation*} then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$. \end{cor} The corollaries given above yield lower estimates for the best constant $\overline{\gamma}$. It is easy to obtain upper estimates. For example, putting $x(t)=t-a$ in $\eqref{uG}$ gives $$ \overline{\gamma}\le\frac{b-a}{\int_a^b(t-a)q(t)\,dt}. $$ Another way for finding the best constant $\overline{\gamma}$ is presented by the following result. \begin{thm}\label{th2} Inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$ if and only if $\gamma\le\overline{\gamma}$, where $\overline{\gamma}$ is smallest $\gamma>0$ such that the Cuachy problem \begin{equation}\label{ee} \begin{gathered} \dot x(t)=\frac{\gamma}2\Bigl(q(t)x(t)+\int_t^bq(s)\dot x(s)\,ds\Bigr),\\ x(a)=0 \end{gathered} \end{equation} has a nonzero solution. \end{thm} \begin{proof} Substituting $x(t)=\int_a^tz(s)\,ds$ in \eqref{ee}, we obtain $\gamma^{-1}z=\frac12Q_bz$. So, $\overline{\gamma}^{-1}$ is the largest eigenvalue of the operator $\frac12Q_b$. Therefore, $\gamma\overline{\gamma}^{-1}=\|K\|$ is the largest eigenvalue of the operator $K$. Then it follows from Theorem \ref{th} that problem \eqref{e} is solvable if and only if $\gamma\le\overline{\gamma}$. \end{proof} If $q$ is absolutely continuous on $[a,b]$, problem \eqref{ee} can be written in the form \begin{equation}\label{eee} \begin{gathered} \dot x(t)=\frac{\gamma}2\Bigl(q(b)x(b)-\int_t^b\dot q(s)x(s)\,ds\Bigr),\\ x(a)=0. \end{gathered} \end{equation} \begin{cor} \label{s15} If the function $q$ is absolutely continuous and \begin{equation}\label{ss3} \gamma\le\frac{2}{(b-a)q(b)+\int_a^b(s-a)|\dot q(s)|\,ds}, \end{equation} then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$. \end{cor} \begin{proof} Problem \eqref{eee} is equivalent to the equation \begin{equation*} x(t)=\frac{\gamma}2\Bigl(q(b)x(b)(t-a)-\int_a^t\int_s^b\dot q(\tau)x(\tau)\,d\tau ds\Bigr). \end{equation*} Let \begin{equation*} \gamma<\gamma^*:=\frac{2}{(b-a)q(b)+\int_a^b(s-a)|\dot q(s)|\,ds}. \end{equation*} If $x$ is a solution to problem \eqref{eee}, then \begin{multline*} \max_{t\in[a,b]}|x(t)|\le \frac{\gamma}2\Bigl(q(b)|x(b)|(b-a)+\int_a^b\int_s^b\dot|q(\tau)|\,d\tau ds\max_{t\in[a,b]}|x(t)|\Bigr)\\ \le \frac{\gamma}2\Bigl(q(b)(b-a)+\int_a^b(s-a)|\dot q(s)|\,ds\Bigr)\max_{t\in[a,b]}|x(t)|\le\frac\gamma{\gamma^*} \max_{t\in[a,b]}|x(t)|. \end{multline*} Hence, $x=0$. Thus we get $\gamma*\le\overline{\gamma}$. \end{proof} From representation \eqref{eee}, it is easy to obtain the following theorem. \begin{thm} Let the function $q$ be absolutely continuous on $[a,b]$. Then inequality \eqref{uG} holds for all $x\in\mathbf{W}_a$ if and only if $\gamma\le\overline{\gamma}$, where $\overline{\gamma}$ is the smallest $\gamma>0$ such that the boundary value problem for the second order ordinary differential equation \begin{gather} \label{eE1}\ddot x(t)=\frac{\gamma}2\dot q(t)x(t),\\ \label{rr} x(a)=0,\\ \label{rdd} \dot x(b)=\frac\gamma2q(b)x(b) \end{gather} has a nonzero solution. \end{thm} If we can integrate equation \eqref{eE1}, then it is possible to obtain the exact value of the best constant. In this paper, we confine ourselves to the case of power functions $q$. First of all, we will gather some information on differential equations and Bessel functions. The general solution to the equation \begin{equation}\label{odm} \ddot x(t)=-\frac12\gamma|p|t^{p-1}x(t), \end{equation} can be written in the form \begin{equation}\label{p1} x(t)=c_1u^-(t)+c_2v^-(t), \end{equation} where $c_1$, $c_2$ are arbitrary constants, and the functions $u^-$, $v^-$ are defined with the help of the Bessel functions of the first and second kind respectively: \begin{equation*} u^-(t)=\sqrt{t}J_{\frac1{p+1}}(\frac{\sqrt{2\gamma|p|}}{p+1}t^{\frac{p+1}2}),\quad v^-(t)=\sqrt{t}Y_{\frac1{p+1}}(\frac{\sqrt{2\gamma|p|}}{p+1}t^{\frac{p+1}2}). \end{equation*} Note that there exist the finite non-zero limits $\lim_{t\to0}\dfrac{u^-(t)}t$ and $\lim_{t\to0}v^-(t)=v(0)$. The solution to \eqref{odm} satisfies the equation \begin{equation}\label{p2} \dot x(t)+\frac{\mathop{\rm sgn}(p)}2\gamma t^px(t)=c_1u^-_1(t)+c_2v^-_1(t), \end{equation} where $$ \mathop{\rm sgn}(p)=\begin{cases} 1&\text{if }p>0,\\ 0&\text{if }p=0,\\ -1&\text{if }p<0,\end{cases} $$ $$u^-_1(t)=-\frac{\mathop{\rm sgn}(p)}2\gamma t^{p+\frac12}J_{-\frac {2p+1}{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2}),$$ $$v^-_1(t)=-\frac{\mathop{\rm sgn}(p)}2\gamma t^{p+\frac12}Y_{-\frac {2p+1}{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2})$$ (here we use the recurrent formulae for Bessel functions $\Theta'_\nu(z)=\Theta_{\nu-1}(z)-\frac\nu z\Theta_\nu(z)$) and $\Theta_{\nu-1}(z)=2\frac\nu z\Theta_\nu(z)-\Theta_{\nu+1}(z)$). There exists the finite non-zero limits $$ \lim_{t\to0}u^-_1(t)=u^-_1(0)\text{\quad and\quad} \lim_{t\to0}v^-_1(t)=v^-(0). $$ If $\frac1{p+1}$ is half-integer (that is, for $p=\dots,-\frac75,-\frac53,-3,1,-\frac13,-\frac35,-\frac57,\dots$), then Bessel functions and solutions to equation \eqref{odm} can be expressed by elementary functions (see Appendix 1). Similarly, the general solution to the equation \begin{equation}\label{odp} \ddot x(t)=\frac12\gamma|p|t^{p-1}x(t) \end{equation} can be written in the form \begin{equation}\label{P1} x(t)=c_1u^+(t)+c_2v^+(t), \end{equation} where $c_1$, $c_2$ are arbitrary constants, and the functions $u^+$, $v^+$ are defined with the help of the modified Bessel functions of the first and second kind respectively: \begin{equation*} u^+(t)=\sqrt{t}I_{\frac1{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2}),\quad v^+(t)=\sqrt{t}K_{\frac1{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2}) \end{equation*} There exist the finite non-zero limits $\lim_{t\to0}\dfrac{u^+(t)}t$ and $\lim_{t\to0}v^+(t)=v^+(0)$. The solution to \eqref{odp} satisfies the equation \begin{equation}\label{P2} \dot x(t)-\frac{\mathop{\rm sgn}(p)}2\gamma t^px(t)=c_1u^+_1(t)+c_2v^+_1(t), \end{equation} where \begin{equation*} u^+_1(t)=-\frac{\mathop{\rm sgn}(p)}2\gamma t^{p+\frac12}I_{-\frac {2p+1}{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2}), \end{equation*} \begin{equation*} v^+_1(t)=-\frac{\mathop{\rm sgn}(p)}2\gamma t^{p+\frac12}K_{-\frac {2p+1}{p+1}}(\frac{\sqrt{2|p|\gamma}}{p+1}t^{\frac{p+1}2}) \end{equation*} (here we use the recurrent formulae for modified Bessel functions $\Theta'_\nu(z)=\Theta_{\nu-1}(z)-\frac\nu z\Theta_\nu(z)$) and $\Theta_{\nu-1}(z)=2\frac\nu z\Theta_\nu(z)+\Theta_{\nu+1}(z)$). There exists the finite non-zero limits $\lim_{t\to0}u^+_1(t)=u^+_1(0)$ and $\lim_{t\to0}v^+_1(t)=v^+_1(0)$ for $-\frac12
0$, and
$x(t)=c\sin\sqrt{-\frac{k\gamma}2}(t-a)$ if $k<0$.
Thus
$$
\overline{\gamma}=\frac{2\upsilon^2}{k(b-a)^2},\quad
\coth(\upsilon)=\upsilon(1+\frac r{k(b-a)})
$$
if $k>0$; and
$$
\overline{\gamma}=-\frac{2\upsilon^2}{k(b-a)^2},\quad
\cot(\upsilon)=\upsilon(1+\frac r{k(b-a)})
$$
if $k<0$, where $\upsilon$ is the smallest positive solution to
the respective equation.
\subsection*{1.3.a. The power function $q(t)=(t-a)^p$, $p>-1$}
Perform the change of variables $t-a=s$, $x(t)=y(s)$. Then
equation \eqref{eE1} has the form \eqref{odm} for $p<0$ and the
form \eqref{odp} for $p>0$. The boundary conditions
(\ref{rr}), (\ref{rdd}) have the form $y(0)=0$, $\dot
y(b-a)-\frac12\gamma (b-a)^py(b-a)=0$. From \eqref{p1} and
\eqref{P1}, it follows that $c_2=0$. From \eqref{p2} and
\eqref{P2}, it follows that $u^-_1(b-a)=0$ and $u^+_1(b-a)=0$,
respectively.
Denote by $\eta_p$ the smallest positive root of the Bessel
function $J_{-\frac{1+2p}{1+p}}$ for $p<0$ and the smallest
positive root of the modified Bessel function
$I_{-\frac{1+2p}{1+p}}$ for $p>0$. The the best constant for every
$p>-1$, $p\ne0$ has the representation
\begin{equation*}
\overline{\gamma}_p=\frac{(p+1)^2\eta_p^2}{2|p|(b-a)^{p+1}}.
\end{equation*}
Whence, in particular, we have \label{pagees}
$\overline{\gamma}_1=\frac{2\eta_1^2}{(b-a)^2}\approx\frac{2.87845768}{(b-a)^2}$,
where $\coth\eta_1=\eta_1$;
$\overline{\gamma}_{-1/3}=\dfrac{\pi^2}{6(b-a)^{2/3}}\approx\dfrac{1.64493407}{(b-a)^{2/3}}$;
$\overline{\gamma}_{-3/5}=\dfrac{2\pi^2}{15(b-a)^{2/5}}\approx\dfrac{1.315947254}{(b-a)^{2/5}}$;
$\overline{\gamma}_{-5/7}=\dfrac{2\eta^2_{-5/7}}{35(b-a)^{2/7}}\approx\dfrac{1.15375592}{(b-a)^{2/7}}$,
where $\eta_{-5/7}$ is the smallest positive solution to the
equation $\tan\eta=\eta$,
$\overline{\gamma}_{-7/9}=\dfrac{2\eta^2_{-7/9}}{63(b-a)^{2/9}}\approx\dfrac{1.05452260}{(b-a)^{2/9}}$,
where $\eta_{-7/9}$ is the smallest positive solution to the
equation $\tan\eta=\dfrac{3\eta}{3-\eta^2}$.
Some approximate values of $\overline{\gamma}_p$ for $b-a=1$ are
given in Appendix 2, and the graph is shown in Figure \ref{eta}.
\begin{figure}[ht]\label{eta}
% \includegraphics[bb=-60 180 550 330]{eta.bmp}
\includegraphics[width=0.6\textwidth]{eta.eps}
\caption{The best constant $\overline{\gamma}_p$ in case 1.3.a.}
\end{figure}
If the Bessel functions are not elementary, then we can obtain
some lower estimates for the best constant $\overline{\gamma}_p$.
From \eqref{HA} for $n=1,2,3,4$, we get respectively more and more
exact estimates:
\noindent %1)
$\overline{\gamma}_p\ge\gamma_p^{(1)}:=\frac{2\sqrt{p
+ 1}}{(b-a)^{p+1}}$ (it was obtained in \cite{WTroy} by W.C.
Troy);
\noindent %2)
$\overline{\gamma}_p\ge\gamma_p^{(2)}:=\frac{2}{(b-a)^{p+1}}\,\left({
\frac {(p + 1)(p+2)(2p + 3)}{p + 6}}\right)^{1/4}$;
\noindent %3)
$\overline{\gamma}_p\ge\gamma_p^{(3)}:=\frac{2}{(b-a)^{p+1}}\,\left({\frac
{(4+3p)(5+4p)(2p+3)(p+2)^{2}(p+1)}{6p^3+ 100p^{2}
+292p+240}}\right)^{1/6}$;
\noindent %4)
{\small $\overline{\gamma}_p\ge\gamma_p^{(4)}:=
\dfrac{2}{(b-a)^{p+1}}\left({\frac{(7+6p)(6 +
5p)(p+1)(p+2)^{2}(2p+3)^2(5+4p)(4+3p)}{180p^5+5431p^4+31882p^3+74652p^2+77832p
+30240}}\right)^{1/8}$.}
Estimating $\|K\|$ with the help of Lemma \ref{l01}, where
$v(t)=\left(1-\frac{\mu}{p+2}\right)
\frac{1}{1-\mu}+(\frac{t-a}{b-a})^{p+1}$ and $$\mu
=\frac{4}{(b-a)^{p+1}}{\frac{(p +
1)(p+2)}{\sqrt{8(p+1)^3+9(p+1)^2-2p-1}+3p+4}},$$ we also get
\noindent %5)
$\overline{\gamma}_p\ge\gamma_p ^{(5)}:=2\mu$ for
$p\ge0$;
\noindent estimating $\|K\|$ with the help of Lemma \ref{l01},
where $v(t)=(\frac{t-a}{b-a})^{-\frac12 +\varepsilon}$ and
$\varepsilon
> 0$ is small enough, we obtain
\noindent %6)
\begin{equation}\label{ps}
\overline{\gamma}_p\ge\gamma_p^{(6)}:={\frac{1}{2}}\,{\frac{1}
{|p|(b-a)^{p+1}}}\quad \text{for }p\in(-1, -1/2].
\end{equation}
Note that
$\gamma_p^{(1)}<\gamma_p^{(2)}<\gamma_p^{(3)}<\gamma_p^{(4)}$ for
all $p>-1$, $p\ne0$. Thus estimate 1) from \cite{WTroy} is not
sharp besides the known case $p=0$ \cite{Beesack}.
In particular, for $p=1$ and $b-a=1$, we have $\gamma_1 ^{(1)}
=2.8284271$, $\gamma_1 ^{(2)} =2.8776356$, $\gamma_1 ^{(3)}
=2.8784392$, $\gamma_1 ^{(4)} =2.8784572$ (the sharp constant
$\overline{\gamma}_1$ equals $2.8784577$ to the last decimal
place).
We have $\gamma_p^{(4)}<\gamma_p^{(5)}$ for $p > 50.1$, therefore
the estimate $\gamma_p^{(5)}$ yields the best result for large
values of $p$. The estimate $\gamma_p^{(6)}$ is best for small
$p$.
Setting $x(t)=(t-a)^{\frac{1+\sqrt{1+p}}2}$ in \eqref{uG}, we
obtain the upper estimate for $\overline{\gamma}_p$:
$$
\overline{\gamma}_p\le\gamma_p^{(0)}=\frac{(1+\sqrt{1+p})^2}2.
$$
The graphics of $\gamma^{(k)}_p$ for $b-a=1$ are shown in Figure
2.
%\ref{Gamma2}.
\begin{figure}[ht]\label{Gamma2}
% \includegraphics[bb=0 -20 500 300]{Gamma2.bmp}
\includegraphics[width=0.7\textwidth]{Gamma2.eps}
\caption{The estimates $\gamma^{(k)}_p$}
\end{figure}
\subsection*{1.3.b. The singular case $q(t)=\dfrac1{t-a}$.}
From \eqref{ps}, the passage to the limit $p\to-1$ yields
$\overline{\gamma} _{-1}\ge\frac{1}{2}$. On the other
hand, setting $x_n(t)=(t-a)^{\frac1n-\frac12}$, $n=2,3,\dots$, we
obtain $\overline{\gamma}_{-1}\le\frac12+\frac1n$. Therefore,
$$\overline{\gamma}_{-1}=\frac12.$$ Moreover, it can be shown that
$$
\begin{array}{l}
\int_a^b\dot x^2(t)\,dt>\frac{1}{2}\int_a^b\frac{|\dot x(t)x(t)|}{t-a}\,\,dt
\end{array}
$$
for all $x\in\mathbf{W}_a$, $x\not\equiv0$.
\subsection*{1.3.c. The power function $q(t)=(t-\tau)^p$, $\tau0$. The boundary conditions
(\ref{rr}),~(\ref{rdd}) have the form $y(a-\tau)=0$, $\dot
y(b-\tau)+\frac12\gamma(b-\tau)^py(b-\tau)=0$.
For $p<0$ from \eqref{p1} and \eqref{p2} we obtain
\begin{gather*}
c_1u^-(a-\tau)+c_2v^-(a-\tau)=0,\\
c_2u^-_1(b-\tau)+c_2v^-_1(b-\tau)=0.
\end{gather*}
Therefore, problem (\ref{eE1} - \ref{rdd}) has a nonzero solution
if and only if
$$
u^-(a-\tau)v^-_1(b-\tau)-v^-(a-\tau)u^-_1(b-\tau)=0,
$$
that is,
\begin{align*}
J_{\frac1{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(a-\tau)^{\frac{p+1}2})
Y_{-\frac{2p+1}{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(b-\tau)^{\frac{p+1}2})&\\
- J_{-\frac{2p+1}{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(b-\tau)^{\frac{p+1}2})
Y_{\frac1{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(a-\tau)^{\frac{p+1}2})&=0.
\end{align*}
For $p>0$ from \eqref{P1} and \eqref{P2} we obtain
\begin{gather*}
c_1u^+(a-\tau)+c_2v^+(a-\tau)=0,\\
c_2u^+_1(b-\tau)+c_2v^+_1(b-\tau)=0.
\end{gather*}
In this case problem (\ref{eE1} - \ref{rdd}) has a nonzero
solution if and only if
$$
u^+(a-\tau)v^+_1(b-\tau)-v^+(a-\tau)u^+_1(b-\tau)=0,
$$
that is,
\begin{align*}
I_{\frac1{p+1}}(\frac{\sqrt{2\gamma
p}}{p+1}(a-\tau)^{\frac{p+1}2})K_{-\frac{2p+1}{p+1}}(\frac{\sqrt{2\gamma
p}}{p+1}(b-\tau)^{\frac{p+1}2})&\\
- I_{-\frac{2p+1}{p+1}}(\frac{\sqrt{2\gamma
p}}{p+1}(b-\tau)^{\frac{p+1}2})K_{\frac1{p+1}}(\frac{\sqrt{2\gamma
p}}{p+1}(a-\tau)^{\frac{p+1}2})&=0.
\end{align*}
Denote by $\theta_p(k)$ the smallest positive root of the
function
\begin{equation*}
J_{\frac1{p+1}}(\theta k) Y_{-\frac{2p+1}{p+1}}(\theta)-
J_{-\frac{2p+1}{p+1}}(\theta) Y_{\frac1{p+1}}(\theta k)
\end{equation*}
for $p<0$ and the smallest positive root of the function
\begin{equation*}
I_{\frac1{p+1}}(\theta k)K_{-\frac{2p+1}{p+1}}(\theta)-I_{-\frac{2p+1}{p+1}}(\theta)
K_{\frac1{p+1}}(\theta k).
\end{equation*}
for $p>0$. Then the best constant is defined by the equality
\begin{equation*}
\overline{\gamma}=\frac{(p+1)^2\theta_p^2((\frac{a-\tau}{b-\tau})^{\frac{p+1}2})}{2|p|(b-\tau)^{p+1}}.
\end{equation*}
Note that $\theta_{-\frac13}(k)$ is the smallest positive solution
to the equation $\cot(1-k)\theta=k\theta$; $\theta_{-\frac35}(k)$
is the smallest positive solution to the equation
$\cot(1-k)\theta=\dfrac{k^2\theta^2-3}{3k\theta}$; $\theta_1(k)$
satisfies the equation $\coth(1-k)\theta=\theta$.
The graphics of $\theta_p(k)$ for some $k$ are shown in Figure 3.
\begin{figure}[ht]\label{theta}
% \includegraphics[bb=0 80 550 330]{theta.bmp}
\includegraphics[width=0.8\textwidth]{theta.eps}
\caption{The graphics of $\theta_p(k)$ and $\eta_p$}
\end{figure}
We have $\lim_{k\to0}\theta_p(k)=\eta_p$.
Now we give some estimate for the best constant
$\overline{\gamma}$.
By \eqref{app}, it follows that
\begin{equation*}
\overline{\gamma}\ge\frac{2\sqrt{2p+1}}
{\sqrt{2p+1-2(p+1)\frac{a-\tau}{b-\tau}+(\frac{a-\tau}{b-\tau})^{2p+2}}}\frac{\sqrt{p+1}}{(b-\tau)^{p+1}}.
\end{equation*}
By \eqref{s2}, it follows that
\begin{equation*}
\overline{\gamma}\ge\frac{2(p+1)}{(b-\tau)^{p+1}-(a-\tau)^{p+1}}\text{
for }p\le0,
\end{equation*}
and
\begin{equation*}
\overline{\gamma}\ge\frac2{(b-\tau)^p(b-a)}\text{ for }p\ge0.
\end{equation*}
\subsection*{1.4.a. The power function $q(t)=(b-t)^p$, $p>0$}
Perform the change of variables $b-t=s$, $x(t)=y(s)$. Then
\eqref{eE1} has the form \eqref{odm}. Boundary condition
\eqref{rdd} has the form $y(0)=0$. From \eqref{p2} it follows
that $c_1u_1^-(0)+c_2v_1^-(0)=0$. From boundary condition
\eqref{rr} and representation \eqref{p1}, it follows that
$c_1u^-(b-a)+c_2v(b-a)=0$. Therefore
$v_1^-(0)u^-(b-a)-u_1^-(0)v^-(b-a)=0$, that is,
$$
\cos(\frac{2p+1}{p+1}\pi)J_{\frac1{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(b-a)^{\frac{p+1}2})-
\sin(\frac{2p+1}{p+1}\pi)Y_{\frac1{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(b-a)^{\frac{p+1}2})=0.
$$
Denote by $\zeta_p$ the smallest positive root of the function
$$
\cos(\frac{2p+1}{p+1}\pi)J_{\frac1{p+1}}(\zeta)-
\sin(\frac{2p+1}{p+1}\pi)Y_{\frac1{p+1}}(\zeta)
$$
(the graphic of $\zeta_p$ is shown in Figure 4).
%\ref{zetaroot}).
\begin{figure}[ht]\label{zetaroot}
% \includegraphics[bb=-30 180 500 330]{zetaroot.bmp}
\includegraphics[width=0.7\textwidth]{zetaroot.eps}
\caption{The graphic of $\zeta_p$}
\end{figure}
In this case the best constant is
\begin{equation*}
\overline{\gamma}_p=\frac{(p+1)^2\zeta_p^2}{2p(b-a)^{p+1}}.
\end{equation*}
In particular,
\begin{equation*}
\overline{\gamma}_1=\frac{\pi^2}{2(b-a)^2}.
\end{equation*}
Some approximate values of $\overline{\gamma}_p$ for $b-a=1$ are
given in Appendix 3, and the graphic is shown in Figure 5.
%\ref{zeta}.
\begin{figure}[ht]\label{zeta}
% \includegraphics[bb=-50 180 550 330]{zeta.bmp}
\includegraphics[width=0.6\textwidth]{zeta.eps}
\caption{The graphic of $\overline{\gamma}_p$ in the case 1.4.a
for $b-a=1$}
\end{figure}
Note that inequality \eqref{app} implies the lower estimate
$$
\overline{\gamma}\ge\frac{2\sqrt{(2p+1)(p+1)}}{(b-a)^{p+1}}.
$$
\subsection*{1.4.b. The power function $q(t)=(\tau-t)^p$, $b<\tau$, $p\ne0$}
First let $p<0$. Perform the change of variables
$\tau-t=s$, $x(t)=y(s)$. Then equation \eqref{eE1} has the form
\eqref{odp}. And from condition \eqref{rr}, \eqref{rdd} and
representation \eqref{P1}, \eqref{P2}, we obtain
\begin{gather*}
c_1u^+(\tau-a)+c_2v^+(\tau-a)=0,\\
c_2u^+_1(\tau-b)+c_2v^+_1(\tau-b)=0.
\end{gather*}
Thus problem \eqref{eE1} - \eqref{rdd} has a nonzero solution if and
only if
$$
u^+(\tau-a)v^+_1(\tau-b)-v^+(\tau-b)u^+_1(\tau-b)=0,
$$
that is,
\begin{align*}
I_{\frac1{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(\tau-a)^{\frac{p+1}2})
K_{-\frac{2p+1}{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(\tau-b)^{\frac{p+1}2})&\\
-I_{-\frac{2p+1}{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(\tau-a)^{\frac{p+1}2})
K_{\frac1{p+1}}(\frac{\sqrt{-2\gamma p}}{p+1}(\tau-a)^{\frac{p+1}2})&=0.
\end{align*}
For $p>0$ the change of variables $t-\tau=s$, $x(t)=y(s)$ reduce
equation \eqref{eE1} to \eqref{odm}. From conditions \eqref{rr},
\eqref{rdd} and representation \eqref{p1}, \eqref{p2} we obtain
the conditions
\begin{gather*}
c_1u^-(\tau-a)+c_2v^-(\tau-a)=0,\\
c_2u^-_1(\tau-b)+c_2v^-_1(\tau-b)=0.
\end{gather*}
Thus problem \eqref{eE1}-\eqref{rdd} has a nonzero solution if and
only if
$$
y^-(\tau-a)v^-_1(\tau-b)-v^-(\tau-a)u^-_1(\tau-b)=0,
$$
that is,
\begin{align*}
J_{\frac1{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(\tau-a)^{\frac{p+1}2})
Y_{-\frac{2p+1}{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(\tau-b)^{\frac{p+1}2})&\\
-J_{-\frac{2p+1}{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(\tau-b)^{\frac{p+1}2})
Y_{\frac1{p+1}}(\frac{\sqrt{2\gamma p}}{p+1}(\tau-a)^{\frac{p+1}2})&=0.
\end{align*}
Denote by $\vartheta_p(k)$ the smallest positive root of the
function
\begin{equation*}
I_{\frac1{p+1}}(\vartheta
k)K_{-\frac{2p+1}{p+1}}(\vartheta)-I_{-\frac{2p+1}{p+1}}(\vartheta)
K_{\frac1{p+1}}(\vartheta k)
\end{equation*}
for $p<0$, and the smallest positive root of the function
\begin{equation*}
J_{\frac1{p+1}}(\vartheta k) Y_{-\frac{2p+1}{p+1}}(\vartheta)-
J_{-\frac{2p+1}{p+1}}(\vartheta) Y_{\frac1{p+1}}(\vartheta k)
\end{equation*}
for $p>0$.
Then the best constant is
\begin{equation*}
\overline{\gamma}=\frac{(p+1)^2\vartheta_p^2((\frac{\tau-a}{\tau-b})^{\frac{p+1}2})}{2|p|(\tau-b)^{p+1}}.
\end{equation*}
In particular, $\vartheta_{-\frac13}(k)$ is the smallest positive
solution to the equation $\coth(k-1)\vartheta=k\vartheta$,
$\vartheta_{-\frac35}$ is the smallest positive solution to the
equation
$\coth(k-1)\vartheta=\frac{k^2\vartheta^2+3}{3k\vartheta}$.
The graphics of $\vartheta_p(k)$ for some $k$ are shown in Figure 6.
\begin{figure}[ht]\label{vartheta}
% \includegraphics[bb=0 150 430 350]{vartheta.bmp}
\includegraphics[width=0.8\textwidth]{vartheta.eps}
\caption{The graphic of $\vartheta_p(k)$}
\end{figure}
Let us get lower estimates for $\overline{\gamma}_p$.
By \eqref{app}, it follows that
\begin{equation*}
\overline{\gamma}\ge\frac{2\sqrt{p+1}\sqrt{2p+1}}
{\sqrt{(2(p+1)a-(2p+1)b-\tau)(b-\tau)^{2p+1}+(\tau-a)^{2p+2}}}.
\end{equation*}
By \eqref{s2}, it follows that
\begin{equation}\label{new25}
\overline{\gamma}\ge\frac2{(\tau-b)^p(b-a)}\text{ for }p\le0
\end{equation}
and
\begin{equation*}
\overline{\gamma}\ge\frac{2(p+1)}{(\tau-a)^{p+1}-(\tau-b)^{p+1}}\text{
for }p\ge0.
\end{equation*}
\section{Integro-differential inequalities with the condition $x(b)=0$}
Denote by $\mathbf{W}^b$ the space of absolutely continuous
functions $x:[a,b]\to\mathbb{R}$ such that $\dot x\in\mathbf{L}_2$
and $x(b)=0$.
Suppose the function $q:[a,b]\to\mathbb{R}$ is measurable,
nonnegative on $[a,b]$, and $\int_a^b(b-t)q^2(t)\,dt<\infty$.
The proofs of all statements in this section are similar to the
proofs of the respective statements in section 1.
\begin{lem}\label{l2}
Inequality \eqref{uG} holds for all $x\in\mathbf{W}^b$ if and only
if the zero function is a solution to the variational problem
\begin{equation}
\begin{gathered}\label{p9}
\int_a^b\left(\dot x^2(t)+\gamma q(t)\dot x(t)x(t)\right)\,dt\to\min,\\
x(b)=0.
\end{gathered}
\end{equation}
\end{lem}
The substitution $x(t)=-\int_t^bz(s)\,ds$ reduces problem
\eqref{p9} to problem \eqref{e}, where $K=\frac\gamma2 Q^a$, and
$Q^a$ is the integral completely continuous operator with the
kernel
$$
Q^a(t,s)=\begin{cases}q(s)&\text{if } a\le s -1$. From here we can conclude
that
$$
\int_a^b\dot x^2(t)\,dt\ge \frac12\int_a^b\frac1{t-a}\,|\dot x(t)x(t)|\,dt
$$
for all $x\in\mathbf{W}_a^b$. The constant $\frac{1}{2}$ is exact.
Note that the best constants in the cases 1.3.b and 3.3.b coincide
and do not depend on $b-a$.
\subsection*{3.4.a. The power function $q(t)=(b-t)^p$, $p>-1$}
It follows from 1.4.b that
$\overline{\gamma}_c=\frac{(p+1)^2\vartheta_p^2((\frac{b-a}{b-c})^{\frac{p+1}2})}{2|p|(b-c)^{p+1}}$.
Using 2.3.a, we get
$\overline{\gamma}^c=\frac{(p+1)^2\eta_p^2}{2|p|(b-c)^{p+1}}$.
Then the exact constant equals
$$
\overline{\gamma}=\frac{(p+1)^2(\eta_p\vartheta_p^{-1}(\eta_p))^2}{2|p|(b-a)^{p+1}},
$$
that is, $\overline{\gamma}$ coincides with the constant from
3.3.a.
\subsection*{3.4.b. The singular case $q(t)=\frac1{b-t}$}
Similarly 3.3.b we have
$$
\int_a^b\dot x^2(t)\,dt\ge \frac12\int_a^b\frac1{b-t}\,|\dot x(t)x(t)|\,dt
$$
for all $x\in\mathbf{W}_a^b$.
\section{Appendix}
\subsection{Bessel functions for some half-integer indexes} \quad\\
$J_{-\frac32}(z)=-\sqrt{\frac2{\pi z}}\frac1z(z\sin z+\cos z)$, \quad
$Y_{-\frac32}(z)=\sqrt{\frac2{\pi z}}\frac1z(z\cos z-\sin z)$,\\
$J_{-\frac12}(z)=\sqrt{\frac2{\pi z}}\cos z$, \quad
$Y_{-\frac12}(z)=\sqrt{\frac2{\pi z}}\sin z$, \quad
$J_{\frac12}(z)=\sqrt{\frac2{\pi z}}\sin z$,\\
$Y_{\frac12}(z)=-\sqrt{\frac2{\pi z}}\cos z$, \quad
$J_{\frac32}(z)=\sqrt{\frac2{\pi z}}\frac1z(-z\cos z+\sin z)$, \\
$Y_{\frac32}(z)=-\sqrt{\frac2{\pi z}}\frac1z(z\sin z+\cos z)$, \quad
$J_{\frac52}(z)=\sqrt{\frac2{\pi z}}\frac1{z^2}(-z^2\sin z-3z\cos z+3\sin z)$, \quad
$Y_{\frac52}(z)=\sqrt{\frac2{\pi z}}\frac1{z^2}(z^2\cos z-3z\sin z-3\cos z)$.\\
$I_{-\frac32}(z)=\sqrt{\frac2{\pi z}}\frac1z(z\sinh z-\cosh z)$, \quad
$K_{-\frac32}(z)=\sqrt{\frac\pi{2z}}\frac1ze^{-z}(z+1)$,\\
$I_{-\frac12}(z)=\sqrt{\frac2{\pi z}}\cosh z$, \quad
$K_{-\frac12}(z)=\sqrt{\frac\pi{2z}}e^{-z}$,\\
$I_{\frac12}(z)=\sqrt{\frac2{\pi z}}\sinh z$, \quad
$K_{\frac12}(z)=\sqrt{\frac\pi{2z}}e^{-z}$,\\
$I_{\frac32}(z)=\sqrt{\frac2{\pi z}}\frac1z(z\cosh z-\sinh z)$,\quad
$K_{\frac32}(z)=\sqrt{\frac\pi {2z}}\frac1ze^{-z}(z+1)$, \\
$I_{\frac52}(z)=\sqrt{\frac2{\pi z}}\frac1{z^2}(z^2\sinh z-3z\cosh z+3\sinh z)$,
$K_{\frac52}(z)=\sqrt{\frac\pi{2z}}\frac1{z^2}e^{-z}(z^2+3z+3)$.
\subsection{Approximate best constants for some $p$ in the case 1.3.a for $b-a=1$}
\begin{equation*}
\overline{\gamma}_p=\frac{(p+1)^2\eta_p^2}{2|p|}, \quad p\ne0.
\end{equation*}
{\footnotesize
\begin{tabbing}
$\overline{\gamma}_p$\quad\= $0.500$\quad\= $0.576$\quad\=
$0.618$\quad\= $0.653$\quad\= $0.712$\quad \=
$0.830$\quad\=$0.929$\quad\= $1.018$\quad\=$1.099$\quad\=
$1.175$\quad\= $1.316$\kill $p$\>$-1.0$\> $-0.99$\>$-0.98$\>
$-0.97$\>$-0.95$\> $-0.90$\>
$-0.85$\> $-0.8$\> $-0.75$\>$-0.7$\> $-0.6$\\
$\overline{\gamma}_p$\> $0.500$\> $0.576$\> $0.618$\> $0.653$\>
$0.712$\> $0.830$\>
$0.929$\>$1.018$\>$1.099$\>$1.175$\> $1.316$\\
\end{tabbing}}
\noindent {\footnotesize
\begin{tabbing}
$\overline{\gamma}_p$\quad\=$1.446$\quad\=$1.508$\quad\=$1.567$\quad\=$1.626$\quad\=$1.683$\quad\=$1.738$\quad\=$1.793$\quad\=$1.846$
\quad\=$1.898$\quad\=$1.950$\quad\=$2.000$\kill
$p$\>$-0.5$\>$-0.45$\>
$-0.4$\>$-0.35$\>$-0.3$\>$-0.25$\>$-0.2$\>$-0.15$\>$-0.1$\>$-0.05$\>$0$\\
$\overline{\gamma}_p$\>$1.446$\>$1.508$\>$1.567$\>$1.626$\>$1.683$\>$1.738$\> $1.793$\>$1.846$\>$1.898$\>$1.950$\>$2.000$\\
\end{tabbing}}
\noindent{\footnotesize
\begin{tabbing}
$\overline{\gamma}_p$\quad\=$2$\quad\=$2.098$\quad\=$2.194$\quad\=$2.286$\quad\=$2.465$\quad\=
$2.636$\quad\=$2.879$\quad\=$3.033$\quad\=$3.256$\quad\=$3.606$\quad\=$3.935$\quad\=$4.123$\kill
$p$\> $0$\> $0.1$\> $0.2$\> $ 0.3$\> $0.5$\>
$0.7$ \>$1$\> $1.2$ \>$1.5$ \>$2$\> $2.5$\>$2.8$\\
$\overline{\gamma}_p$\> $2$\> $2.098$\> $2.194$\> $2.286$\>
$2.465$\>
$2.636$\> $2.879$\> $3.033$\> $3.256$\> $3.606$\>$3.935$\>$4.123$\\
\end{tabbing}}
\noindent{\footnotesize
\begin{tabbing}
$\overline{\gamma}_p$\quad\=$4.245$\quad\=$4.541$\quad\=$4.824$\quad\=$5.357$\quad\=$5.854$\quad\=$6.766$\quad\=$7.593$\quad\=$10.976$\quad\=$18.000$\quad\=$26.103$\quad\=$86.938$\kill
$p$\>$3$\>$3.5$\>$4$\>$5$\>$6$\>$8$\>$10$\>$20$\>$50$\>$100$\>$1000$\\
$\overline{\gamma}_p$\>$4.245$\>$4.541$\>$4.824$\>$5.357$\>$5.854$\>$6.766$\>$7.593$\>$10.976$\>$18.000$\>$26.103$\>$86.938$\\
\end{tabbing}}
\subsection{Approximate best constants for some $p$ in the case 1.4.a for $b-a=1$}
\begin{equation*}
\overline{\gamma}=\frac{(p+1)^2\zeta_p^2}{2p},\quad p\ne0.
\end{equation*}
{\footnotesize
\begin{tabbing}
$\overline{\gamma}$\quad\= $2.000$\quad\=$2.299$\quad\=
$2.595$\quad\= $3.477$\quad\= $4.935$\quad\=
$7.837$\quad\=$10.734$\quad\=$13.628$\quad\=$16.522$\quad\=
$22.307$\quad\=$30.984$\kill
$p$\>$0$\>$0.1$\>
$0.2$\>$0.5$\> $1.0$\>
$2$\>$3$\>$4$\>$5$\>$7$\>$10$\\
$\overline{\gamma}$\> $2.000$\>$2.299$\> $2.595$\> $3.477$\>
$4.935$\>
$7.837$\>$10.734$\>$13.628$\>$16.522$\> $22.307$\>$30.984$\\
\end{tabbing}}
\subsection{Approximate best constants for some
$p$ in the case 3.3.a for $b-a=1$}
$$
\overline{\gamma}=\frac{(p+1)^2(\eta_p\vartheta_p^{-1}(\eta_p))^2}{2|p|},\quad
p\ne0.
$$
{\footnotesize
\begin{tabbing}
$\overline{\gamma}$\quad\= $0.500$\quad\= $0.5890$\quad\=
$0.6943$\quad\= $0.7833$\quad\= $0.9826$\quad\=
$1.69271$\quad\=$2.3657$\quad\=$3.0252$\quad\=$3.6766$\quad\=
$4.0000$\kill $p$\>$-1.0$\> $-0.99$\>$-0.97$\>$-0.95$\> $-0.90$\>
$-0.7$\> $-0.5$\>$-0.3$\>$-0.1$\>$0$\\
$\overline{\gamma}$\> $0.500$\> $0.5890$\> $0.6943$\> $0.7833$\>
$0.9826$\>
$1.69271$\>$2.3657$\>$3.0252$\>$3.6766$\> $4.0000$\\
\end{tabbing}}
\noindent{\footnotesize
\begin{tabbing}
$\overline{\gamma}$\quad\= $0.500$\quad\= $0.5890$\quad\=
$0.6943$\quad\= $0.7833$\quad\= $0.9826$\quad\=
$1.69271$\quad\=$2.3657$\quad\=$3.0252$\quad\=$3.6766$\quad\=
$4.0000$\kill $p$\>$0.1$\> $0.3$\>$0.5$\> $1.0$\>
$1.5$\> $2$\>$3$\>$4$\>$5$\>$10$\\
$\overline{\gamma}$\> $4.3211$\> $4.9630$\> $5.6003$\> $7.7863$\>
$8.74181$\>$10.293$\>$13.368$\>$16.415$\> $19.443$\>$34.396$\\
\end{tabbing}}
\subsection*{Acknowledgment}
The authors would like to thank professor N. V. Azbelev for
his helpful discussions and comments.
\begin{thebibliography}{00}
\bibitem{AMR2} \textsc{N. V. Azbelev and L. F. Rakhmatullina,} Theory of linear
abstract functional differential equations and applications. Mem.
Differential Equations Math. Phys. {\bf 8} (1996), 1-102.
\bibitem{Beesack} \textsc{P. R. Beesack,} On an integral inequality of
Z. Opial. Trans. Amer. Math. Soc. {\bf 104}
(1962), 470-475.
\bibitem{Grusia}\textsc{E. I. Bravyi, S. A. Gusarenko,}
On a class of generalized variational inequalities. Memoirs on
Mem. Differential Equations Math.Phys. {\bf 26} (2002), 43-54.
\bibitem{Vestnik} \textsc{E. I. Bravyi, S. A. Gusarenko,} On a class of
integral inequalities with deviating argumemt. (Russian) Perm.
Functional Differential Equations. J. of Perm State Technical
University, (2002), 74-86.
\bibitem{Gusarenko1} \textsc{S. A. Gusarenko,}
On variational problems with linear restrictions. (Russian) Izv.
vusov. Matematika. N 2, (1999), 30-44.
\bibitem{Gusarenko2} \textsc{S. A. Gusarenko,}
On the conditions of positife definiteness of operators in Hilbert
space. (Russian) Izv. vusov. Matematika. 2002, N 1, 1-4.
\bibitem{FAN} {\em Functional Analisys.} (Russian) Moscow,
"Nauka", 1972.
\bibitem{Oleck}\textsc{C. Olech,} A simple proof of a certain result of
Z. Opial, Ann. Polon. Math. {\bf 8} (1960),
61-63.
\bibitem{Opial}\textsc{Z. Opial,} Sur une inegalite, Ann. Polon. Math.
{\bf 8} (1960), 29-32.
\bibitem{WTroy}\textsc{W. C. Troy,} On the Opial-Olech-Beesack inequalities.
USA-Chile Workshop on Nonlinear Analysis,
Electron. J. Diff. Eqns., Conf. 06, 2001, 297--301. http://ejde.math.txstate.edu or
http://ejde.math.unt.edu
\end{thebibliography}
\end{document}