Related papers: Alternative representation of Magnus series by exa…
We compare and discuss the respective efficiency of three methods (with two variants for each of them), based respectively on Taylor (Maclaurin) series, Pad\'{e} approximants and conformal mappings, for solving quasi-analytically a…
I compare the matrix representation of the basic statements of Special Relativity with the conventional vector space representation. It is shown, that the matrix form reproduces all equations in a very concise and elegant form, namely:…
Let v be a real polynomial of even degree, and let \rho be the equilibrium probability measure for v with support S; so that v(x)\geq 2\int \log |x-y| \rho (dy)+C_v for some constant C_v with support S. Then S is the union of finitely many…
Matrix seriation, the problem of permuting the rows and columns of a matrix to uncover latent structure, is a fundamental technique in data science, particularly in the visualization and analysis of relational data. Applications span…
For almost 75 years, the general solution for the Schr\"odinger equation was assumed to be generated by an exponential or a time-ordered exponential known as the Dyson series. We study the unitarity of a solution in the case of a singular…
Moving beyond the classical additive and multiplicative approaches, we present an "exponential" method for perturbative renormalization. Using Dyson's identity for Green's functions as well as the link between the Faa di Bruno Hopf algebra…
This paper investigates analytical approximate solutions for a system of multipantograph delay differential equations using the residual power series method (RPSM), which obtains a Taylor expansion of the solutions and produces the exact…
Banded matrices can be used as precision matrices in several models including linear state-space models, some Gaussian processes, and Gaussian Markov random fields. The aim of the paper is to make modern inference methods (such as…
Bernstein's theorem (also called Hausdorff--Bernstein--Widder theorem) enables the integral representation of a completely monotonic function. We introduce a finite completely monotonic function, which is a completely monotonic function…
Schr\"odinger equations with time-dependent potentials are of central importance in quantum physics and theoretical chemistry, where they aid in the simulation and design of systems and processes at atomic scales. Numerical approximation of…
Based on functional analysis, we propose an algorithm for finite-norm solutions of higher-order linear Fuchsian-type ordinary differential equations (ODEs) P(x,d/dx)f(x)=0 with P(x,d/dx):=[\sum_m p_m (x) (d/dx)^m] by using only the four…
We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real-time using higher-order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this…
We give an algorithm for reversion of formal power series, based on an efficient way to implement the Lagrange inversion formula. Our algorithm requires $O(n^{1/2}(M(n) + MM(n^{1/2})))$ operations where $M(n)$ and $MM(n)$ are the costs of…
The exponential of an operator or matrix is widely used in quantum theory, but it sometimes can be a challenge to evaluate. For non-commutative operators ${\bf X}$ and ${\bf Y}$, according to the Campbell-Baker-Hausdorff-Dynkin theorem,…
An algorithm for numerically computing the exponential of a matrix is presented. We have derived a polynomial expansion of $e^x$ by computing it as an initial value problem using a symbolic programming language. This algorithm is shown to…
A new representation for a regular solution of the perturbed Bessel equation of the form $Lu=-u"+\left( \frac{l(l+1)}{x^2}+q(x)\right)u=\omega^2u$ is obtained. The solution is represented as a Neumann series of Bessel functions uniformly…
We consider some non-linear non-homogeneous partial differential equations (PDEs) and derive their exact Green function solution as a functional Taylor expansion in powers of the source. The kind of PDEs we consider are dispersive ones…
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the…
We consider a finite family of invertible $2 \times 2$ real matrices and a transitive Markov shift on the index set. Let $\lambda$ be the top Lyapunov exponent for random matrix products driven by the Markov shift. We prove that, if the…
Burgers' equation with fixed Dirichlet boundary conditions is considered on generic bounded intervals. By using the Hopf-Cole transformation and the exact operational solution recently established for linear reaction-diffusion equations…