The automated modeling of multi-physical dynamical systems is usually realized by coupling different subsystems together via certain interface or coupling conditions. This approach results in large-scale high-index differential-algebraic equations (DAEs). Since the direct numerical simulation of these kind of systems leads to instabilities and possibly non-convergence of the numerical methods a...

In many mathematical models of physical phenomenons and engineering fields, such as electrical circuits or mechanical multibody systems, which generate the differential algebraic equations (DAEs) systems naturally. In general, the feature of DAEs is a sparse large scale system of fully nonlinear and high index. To make use of its sparsity, this paper provides a simple and efficient algorithm fo...

In this paper, by using majorization inequalities, upper bounds on summations of eigenvalues (including the trace) of the solution for the Lyapunov matrix differential equation are obtained. In the limiting cases, the results reduce to bounds of the algebraic Lyapunov matrix equation. The effectiveness of the results are illustrated by numerical examples.

In this paper, by using majorization inequalities, upper bounds on summations of eigenvalues (including the trace) of the solution for the Lyapunov matrix differential equation are obtained. In the limiting cases, the results reduce to bounds of the algebraic Lyapunov matrix equation. The effectiveness of the results are illustrated by numerical examples.

Nonsmooth equation-solving and optimization algorithms which require local sensitivity information are extended to systems with nonsmooth parametric differential-algebraic equations embedded. Nonsmooth differential-algebraic equations refers here to semi-explicit differential-algebraic equations with algebraic equations satisfying local Lipschitz continuity and differential right-hand side func...

We consider a nonsymmetric algebraic matrix Riccati equation arising from transport theory. The nonnegative solutions of the equation can be explicitly constructed via the inversion formula of a Cauchy matrix. An error analysis and numerical results are given. We also show a comparison theorem of the nonnegative solutions.

In this paper, a new numerical method for solving the fractional Riccati differential  equation is presented. The fractional derivatives are described in the Caputo sense. The method is based upon  fractional-order Bernoulli functions approximations. First, the  fractional-order Bernoulli functions and  their properties are  presented. Then, an operational matrix of fractional order integration...

In this paper, we develop a numerical resolution of the space-time fractional advection-dispersion equation. We utilize spectral-collocation method combining with a product integration technique in order to discretize the terms involving spatial fractional order derivatives that leads to a simple evaluation of the related terms. By using Bernstein polynomial basis, the problem is transformed in...

a novel and eective method based on haar wavelets and block pulse functions(bpfs) is proposed to solve nonlinear fredholm integro-dierential equations of fractional order.the operational matrix of haar wavelets via bpfs is derived and together with haar waveletoperational matrix of fractional integration are used to transform the mentioned equation to asystem of algebraic equations. our new m...