This paper is devoted to the investigation of mixed Fourier-Bessel transformation (※Equations, See Full-text) We apply the method of [6] which provides the estimate for weighted L(sub)$infty$-norm of the spherical mean of │f│$^2$ via its weighted L$_1$-norm (generally it is wrong without the requirement of the non-negativity of f). We prove that in the case of Fourier-Bessel transformatin the mentioned method provides (in dependence on the relation between the dimension of the space of non-special variables n and the length of multiindex ν) similar estimates for weighted spherical means of │f│$^2$, the allowed powers of weights are also defined by multiindex ν. Further those estimates are applied to partial differential equations with singular Bessel operators with respect to y$_1$, …, y(sub)m and we obtain the corresponding estimates for solutions of the mentioned equations.