The aim of this work is to propose a fast and reliable algorithm for computing integrals of the type $$\int_{-\infty}^{\infty} f(x) e^{\scriptstyle -x^2 -\frac{\scriptstyle 1}{\scriptstyle x^2}} dx,$$ where $f(x)$ is a sufficiently smooth function, in floating point arithmetic. The algorithm is based on a product integration rule, whose rate of convergence depends only on the regularity of $f$, since the coefficients of the rule are ``exactly'' computed by means of suitable recurrence relations here derived. We prove stability and convergence in the space of locally continuous functions

We numerically study the dynamics of quasi-two dimensional cholesteric liquid crystal droplets in the presence of a time-dependent electric field, rotating at constant angular velocity. A surfactant sitting at the droplet interface is also introduced to prevent droplet coalescence. The dynamics is modeled following a hybrid numerical approach, where a standard lattice Boltzmann technique solves the Navier-Stokes equation and a finite difference scheme integrates the evolution equations of liquid crystal and surfactant.

In this paper we consider a non-standard discretization to a Volterra integro-dierential system which includes a number of age-of-infection models in the literature. The aim is to provide a general framework to analyze the proposed scheme for the numerical solution of a class of problems whose continuous dynamic is well known in the literature and allow a deeper analysis in cases where the theory lacks

We propose a numerical method for a general integro-differential system of equations which includes a number of age-of-infection epidemic models in the literature [1, 2]. The numerical solution is obtained by a non-standard discretization of the nonlinear terms in the system, and agrees with the analytical solution in many important qualitative aspects. Both the behaviour at finite time and the asymptotic properties of the solution are preserved for any value of the discretization parameter.

This paper deals with the limit cases for s-fractional heat flows in a cylindrical domain, with homogeneous Dirichlet boundary conditions, as s-+ 0+ and s-+ 1-. We describe the fractional heat flows as minimizing move-ments of the corresponding Gagliardo seminorms, with re-spect to the L2 metric. To this end, we first provide a Gamma-convergence analysis for the s-Gagliardo seminorms as s-+ 0+ and s-+ 1-; then, we exploit an abstract stability result for minimizing movements in Hilbert spaces, with respect to a sequence of Gamma-converging uniformly lambda-convex energy function-als.

In this manuscript we propose a numerical method for non-linear integro-differential systems arising in age-of-infection models in a heterogeneously mixed population. The discrete scheme is based on direct quadrature methods and provides an unconditionally positive and bounded solution. Furthermore, we prove the existence of the numerical final size of the epidemic and show that it tends to its continuous equivalent as the discretization steplength vanishes.

In this paper we analyze the stability of the problem of performing a rational QZ$step with a shift that is an eigenvalue of a given regular pencil H-lambda K in unreduced Hessenberg-Hessenberg form. In exact arithmetic, the backward rational QZ step moves the eigenvalue to the top of the pencil, while the rest of the pencil is maintained in Hessenberg-Hessenberg form, which then yields a deflation of the given shift. But in finite-precision the rational QZ step gets ``blurred'' and precludes the deflation of the given shift at the top of the pencil.

Background: Major Depressive Disorder (MDD) is a psychiatric illness that is often associated with potentially life -threatening physiological changes and increased risk for suicidal behavior. Electroencephalography (EEG) research suggests an association between depression and specific frequency imbalances in the frontal brain re-gion.

In this paper, we extend the result on the probability of (falsely) connecting two distinct components when learning a GGM (Gaussian Graphical Model) by the joint regression based technique. While the classical method of regression based technique learns the neighbours of each node one at a time through a Lasso penalized regression, its joint modification, considered here, learns the neighbours of each node simultaneously through a group Lasso penalized regression.

Epidemic models structured by the age of infection can be formulated in terms of a system of renewal equations and represent a very general mathematical framework for the analysis of infectious diseases ([1, 2]). Here, we propose a formulation of renawal equations that takes into account of the behavioral response of individuals to infection. We use the so called "information index", which is a distributed delay that summarizes the information available on current and past disease trend, and extend some results regarding compartmental behavioral models [3, 4, 5].