Linear Differential Equations Research Articles

Description of Generalized Isomonodromic Deformations of Rank-Two Linear Differential Equations Using Apparent Singularities

In this paper, we consider the generalized isomonodromic deformations of rank-two irregular connections on the Riemann sphere. We introduce Darboux coordinates on the parameter space of a family of rank-two irregular connections by apparent singularities. Using the Darboux coordinates, we describe the generalized isomonodromic deformations as Hamiltonian systems.

NUMERICAL METHOD FOR SOLVING THE BOUNDARY VALUE PROBLEM FOR THE PARABOLIC EQUATION

A linear boundary value problem for a parabolic equation is considered in a closed domain. Based on the broken line method, the boundary value problem for a parabolic equation is replaced by a two-point boundary value problem for a system of linear ordinary differential equations by discretizing the unknown function u(t,x) with respect to the variable x. The obtained two-point boundary value problem is investigated by the parameterization method of Professor Dzhumabaev. Based on this method, an algorithm for finding a numerical solution to the two-point boundary value problem for a system of linear ordinary differential equations is constructed. The constructed algorithm is realized by applying known numerical methods. The constructiveness and efficiency of the parameterization method also allows us to construct a numerical solution of the considered linear boundary value problem for the parabolic equation. One numerical example is given to verify and illustrate the proposed algorithm.

Level-set based shape optimization for plane elastic structures using radial basis functions and Hilbertian descent direction

Structural optimization problems are often associated with the so-called shape functionals depending on a shape through its geometry and the state being a solution of given partial differential equation. In such a framework it is convenient to work with the gradient-like method based on a concept of a shape derivative and level set method. The key idea of level set method is to represent the structural boundary with zero level set of given function (level set function—LSF). Now, changing the shape of a structure under optimization is equivalent to transport the LSF in such a direction that ensures decreasing the value of the objective functional. To this end, we make use of coercive bilinear form taken from the weak formulation of elasticity problem to obtain descent direction at each iteration. This descent direction is a solution of an additional variational problem, involving the bilinear form mentioned above and the volumetric expression of the shape derivative plays the role of a linear form. In this paper, we combine level set method with radial basis functions (RBFs) used to approximate LSF. We focus on the so-called multiquadric RBFs, but other classes of RBFs are also briefly considered. This eventually leads to transformation of partial differential equation (linear transport equation governing the evolution of shapes) to a system of linear ordinary differential equations which admits analytical formula for the solution. We apply our method to compliance minimization of a cantilever problem as well as to total potential energy minimization of a structure with kinematic loading. To run all the numerical experiments, we wrote our own code in Wolfram Mathematica environment.

The operational matrices for Elliptic Partial Differential Equations with mixed boundary conditions

Abstract The purpose of this research is to implement the orthogonal polynomials associated with operational matrices to get the approximate solutions for solving two-dimensional elliptic partial differential equations (E-PDEs) with mixed boundary conditions. The orthogonal polynomials are based on the Standard polynomial (xi ), Legendre, Chebyshev, Bernoulli, Boubaker, and Genocchi polynomials. This study focuses on constructing quick and precise analytic approximations using a simple, elegant, and potent technique based on an orthogonal polynomial representation of the solution as a double power series. Consequently, a linear partial differential equation is transformed into a linear algebraic system which is solved by the Mathematica®12. Approximate solutions can be found if the answers are polynomials in and of itself. Three applications involving well-known linear problems Laplace, Poisson, and Helmholtz equations have been solved by using the proposed methods, and a comparison of the approaches has been provided. Furthermore, the computation of the error norm L∞ has been done to show the accuracy of the suggested approaches. The results clearly demonstrate how precise, efficient, and dependable the proposed methods are in obtaining rough solutions to the problem. Bernoulli was one of the best methods in most examples.

Tempered fractional differential equations on hyperbolic space

Abstract In this paper we study linear fractional differential equations involving tempered Caputo-type derivatives in the hyperbolic space. We consider in detail the three-dimensional case for its simple and useful structure. We also discuss the probabilistic meaning of our results in relation to the distribution of an hyperbolic Brownian motion time-changed with the inverse of a tempered stable subordinator. The generalization to an arbitrary dimension n can be easily obtained. We also show that it is possible to construct a particular solution for the non-linear porous-medium type tempered equation by using elementary functions.

The boundary value problem for an ordinary linear half-order differential equation

This study is devoted to the study of the solution of a boundary value problem for an ordinary linear differential equation of half order with constant coefficients. Using of the fundamental solution of the main part of the considered equation, we obtained the principal relations, from which we obtain the necessary conditions for the Fredholm property of the original problem. Further, using the Mittag-Leffler function, a general solution of the homogeneous equation is obtained. Finally, the problem under consideration is reduced to an integral Fredholm equation of the second kind with a non-singular kernel, i.e., the Fredholm property of the stated problem is proved.

ЗАДАЧА ОПТИМАЛЬНОГО СТОХАСТИЧЕСКОГО УПРАВЛЕНИЯ И ОЦЕНКИ СТОИМОСТИ КОРПОРАТИВНЫХ ОБЛИГАЦИЙ, ЗАВИСЯЩИХ ОТ ВРЕМЕНИ И СЛУЧАЙНЫХ СОСТОЯНИЙ ПАРАМЕТРОВ

In the optimal stochastic control problem, one estimates, using the utility indifference method, the bond price and the premium of a default swap contract (Credit default swap) when the model parameters (market interest rate, drift coef-ficient, and volatility of risky underlying assets) are random functions of time and state. Namely, the interest rate of the risk-free asset depends on time, and the prices of risky assets are described by linear homogeneous stochastic dif-ferential equations (SDEs) with multiplicative noise. To do this, the authors consider a portfolio with a risk-free asset and a risky asset with no default risk. For each such portfolio, the authors determine the amount of risky assets that maximizes the expected utility of its final wealth. This quantity allows one to solve the parabolic partial differential equations arising from the Hamilton-Jacobi-Bellman (HJB) equation by transforming them into ordinary differential equations using the method of variable separation to obtain the instantaneous value function of each portfolio. The authors derive the bond price and the default swap-contract premium (CDS) rate, which are the amounts that provide the same level of the expected utility by investing all of one’s wealth in a portfolio that does not contain these credit instruments, or by investing these amounts in credit instruments and the remainder of one’s wealth in the portfolio.

On separability of semigroups: application to models of physics and discrete optimization

Nilpotent semigroups are used in various fields of physics to model processes like system decay in quantum mechanics, state transitions in statistical mechanics, and simplifying dynamical systems. They also assist in optimizing control processes. The set of nilpotent matrices forms a semigroup and plays a key role in Jordan decomposition, which has many physical applications. In quantum mechanics, the Jordan form helps analyze linear operators on Hilbert spaces, especially with systems involving eigenvalue multiplicities. In vibration analysis, it aids in studying normal modes of mechanical systems with degenerate frequencies. Additionally, Jordan decomposition simplifies control theory for linear time-invariant systems, stability analysis in dynamical systems, and the behavior of coupled oscillators, as well as solving systems of linear differential equations. Certainly, the semigroup theory has even more applications in theoretical computer science: suffice it to say, for example, that some authors identify semigroups and finite automata. In this paper, we address the problem of P-separability of semigroups with respect to three key predicates: equality-separability, divisibility-separability, and subsemigroup-separability, achieved through homomorphisms into a nilpotent semigroup. We present some significant results that advance the understanding of these concepts.

Deep learning for solving initial path optimization of mean-field systems with memory

We consider the problem of finding the optimal initial investment strategy for a system modelled by a linear McKean–Vlasov (mean-field) stochastic differential equation with delay, driven by Brownian motion and a pure jump Poisson random measure. The goal is to determine the optimal initial values for the system in the period [ − δ , 0 ] , where δ > 0 is a delay constant, before the system starts at t = 0. Due to the delay in the dynamics, the system will, after startup, be influenced by these initial investment values. It is known that linear stochastic delay differential equations are equivalent to stochastic Volterra integral equations. By utilizing this equivalence, we can find implicit expressions for the optimal investment. Moreover, we propose a deep neural network-based algorithm to solve the stochastic control problem with delay. Specifically, we employ a multi-layer feed-forward neural network for control modelling in the interval [ − δ , 0 ] , and use back-propagation to train the feed-forward neural network. The gradient of the loss function is computed using stochastic gradient descent (SGD) with respect to the weights of the network.

Analytical Solution of Some Non – Linear Delay Differential Equations Using Adomian Decomposition Method

Analytical Solution of Some Non – Linear Delay Differential Equations Using Adomian Decomposition Method

The Riemann–Hilbert Approach to the Higher-Order Gerdjikov–Ivanov Equation on the Half Line

The Fokas method exhibits remarkable versatility in solving boundary value problems associated with both linear and nonlinear partial differential equations, particularly when conventional approaches encounter challenges in handling intricate boundary conditions. The existing literature often lacks the incorporation of unconventional boundary conditions, and this study addresses this issue by extending the application of the Fokas method to the higher-order Gerdjikov-Ivanov equation on the half line (−∞,0]. We have demonstrated the exclusive representation of the potential function u(z,t) in the higher-order Gerdjikov–Ivanov equation through the solution of a Riemann–Hilbert problem. The characteristic function is partitioned on the complex plane, and we obtain the jump matrix between each partition based on the positive and negative values of the partition as well as the spectral matrix determined by the initial data and boundary value data. The findings suggest that the spectral functions are not mutually independent; instead, they conform to a global relationship. The novel aspect of this study is the application of the Fokas method to a previously unexplored case, contributing to the theoretical and practical understanding of complex partial differential equations and filling a gap in the treatment of boundary conditions.

A computational probabilistic procedure to quantify the time of breast cancer recurrence after chemotherapy administration

We propose a random linear differential equation with jumps to model the dynamics of breast tumor growth using real patients’ data. The model considers the effect of chemotherapy administration and tumor resection. Also, the inherent randomness of the data and the model are taken into account.The model is probabilistically solved by combining two main methods belonging to Uncertainty Quantification: the principle of maximum entropy (PME) and the random variable transformation technique (RVT). The PME is applied to assign proper probability distributions to model parameters. The RVT technique is used to determine the probability density function (PDF) of the solution of the random differential equation, which is a stochastic process.We apply computational optimization techniques to determine the PDF of model parameters so that the PDF of the solution matches the ones assigned to real patients’ data.Once random tumor dynamics has been modeled, we simulate, after surgery, different adjuvant chemotherapy strategies with the aim of delaying tumor recurrence as long as possible. Results are in concordance with medical literature, and tumor relapse is delayed as more cycles of chemotherapy are administered.Although the proposed model has a common general structure, it is shown how it can be customized to patients in order to construct projections about the tumor’s growth. To better illustrate the applicability of the proposed approach, we carefully show how our model can be tailored to two real patients treated at the Hospital Clínico de Valencia (Spain). The obtained results have been overseen by doctors from this hospital.

Application of Chelyshkov polynomials in solving stochastic model with fractional Brownian motion

Abstract This paper introduces a computational spectral collocation approach utilizing orthonormal Chelyshkov polynomials to estimate solutions for both linear and nonlinear stochastic differential equations containing fractional Brownian motion characterized by the Hurst parameter H. The method employs Newton–Cotes nodes to transform these integral equations into manageable linear and nonlinear algebraic forms. The reliability of the proposed method is established through theorems on error estimation and convergence analysis. Moreover, some examples are given to demonstrate the viability, credibility, coherence, and dependability of the approach. Also, a comparative evaluation is conducted to contrast the results obtained from the suggested approach with those produced by the spectral collocation method utilizing orthonormal Bernoulli polynomials. Furthermore, a 95 % confidence interval is constructed for the error mean, revealing that the approximations maintain high accuracy.

Interval observers for hybrid continuous-time stationary systems

The problem of interval observer design for hybrid continuous-time stationary systems under external disturbances and measurement noise is studied. It is assumed that continuous-time dynamic of such systems is described by linear or nonlinear differential state equations with linear function of output. System parameters depend on system states and are switched based on the control system with finite set of states. The relations allowing designing hybrid interval observer of minimal dimension estimating the set of admissible values of the prescribed linear vector function of the system states are derived. To solve the problem, algebra of partitions and linear algebra are used. Theoretical results are illustrated by example.

On intermediate exceptional series

The Freudenthal–Tits magic square m(A1,A2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {m}(\\mathbb {A}_1,\\mathbb {A}_2)$$\\end{document} for A=R,C,H,O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb {A}=\\mathbb {R},\\mathbb {C},\\mathbb {H},\\mathbb {O}$$\\end{document} of semi-simple Lie algebras can be extended by including the sextonions S\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb {S}$$\\end{document}. A series of non-reductive Lie algebras naturally appear in the new row associated with the sextonions, which we will call the intermediate exceptional series, with the largest one as the intermediate Lie algebra E7+1/2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$E_{7+1/2}$$\\end{document} constructed by Landsberg–Manivel. We study various aspects of the intermediate vertex operator (super)algebras associated with the intermediate exceptional series, including rationality, coset constructions, irreducible modules, (super)characters and modular linear differential equations. For all gI\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}_I$$\\end{document} belonging to the intermediate exceptional series, the intermediate VOA L1(gI)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_1(\\mathfrak {g}_I)$$\\end{document} has characters of irreducible modules coinciding with those of the simple rational C2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C_2$$\\end{document}-cofinite W-algebra W-h∨/6(g,fθ)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$W_{-h^\\vee /6}(\\mathfrak {g},f_\ heta )$$\\end{document} studied by Kawasetsu, with g\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g} $$\\end{document} belonging to the Cvitanović–Deligne exceptional series. We propose some new intermediate VOA Lk(gI)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_k(\\mathfrak {g}_I)$$\\end{document} with integer level k and investigate their properties. For example, for the intermediate Lie algebra D6+1/2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D_{6+1/2}$$\\end{document} between D6\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D_6$$\\end{document} and E7\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$E_7$$\\end{document} in the subexceptional series and also in Vogel’s projective plane, we find that the intermediate VOA L2(D6+1/2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_2(D_{6+1/2})$$\\end{document} has a simple current extension to a SVOA with four irreducible Neveu–Schwarz modules, and the supercharacters can be realized by a fermionic Hecke operator on the N=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=1$$\\end{document} Virasoro minimal model L(c16,2,0)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L(c_{16,2},0)$$\\end{document}. We also provide some (super) coset constructions such as L2(E7)/L2(D6+1/2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_2(E_7)/L_2(D_{6+1/2})$$\\end{document} and L1(D6+1/2)⊗2/L2(D6+1/2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_1(D_{6+1/2})^{\\otimes 2}\\!/L_2(D_{6+1/2})$$\\end{document}. In the end, we find that the theta blocks associated with the intermediate exceptional series produce some new holomorphic Jacobi forms of critical weight and lattice index.

Best Ulam constants for damped linear oscillators with variable coefficients

An associated Riccati equation is used to study the Ulam stability of non-autonomous linear differential equations that model the damped linear oscillator. In particular, the best (minimal) Ulam constants for these equations are derived. These robust results apply to equations with solutions that blow up in finite time, as well as to equations with solutions that exist globally on (−∞,∞). Illustrative, non-trivial examples are presented, highlighting the main results.

Cole–Hopf linearization of the thermocapillary Marangoni dynamics of a two-dimensional bubble with insoluble surfactant

The Marangoni stress-induced dynamics of a two-dimensional inviscid bubble loaded with insoluble surfactant moving in a linear temperature gradient is studied in the limit of small Reynolds, capillary and thermal Péclet numbers. The bubble moves due to the combined effects of thermocapillary Marangoni stresses and those induced by the advective-diffusive spreading of an initial concentration of surfactant on the bubble boundary. It is shown that this nonlinear multiphysics initial value problem is linearizable at any finite non-zero value of a surface Péclet number Pes\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Pe_s$$\\end{document} governing the size of surface diffusion of the insoluble surfactant relative to its advective spreading. This is done by showing that the dynamics can be encoded in the evolution of a function, analytic and single-valued outside the bubble, that satisfies a complex partial differential equation of Burgers type. It is shown that this equation can be linearized by a complex generalization of the classical Cole–Hopf transformation. A numerical method is formulated to solve this linear partial differential equation and determine the Marangoni dynamics of a bubble with some initial surfactant concentration and illustrative calculations are carried out. Results are shown to be consistent with exact equilibrium solutions available from the formulation as Pes→0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Pe_s \\rightarrow 0$$\\end{document} and Pes→∞\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Pe_s \\rightarrow \\infty $$\\end{document} and a perturbative formula for the bubble migration velocity at large but finite Pes\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Pe_s$$\\end{document}.

Differential equations for Carrollian amplitudes

Differential equations are powerful tools in the study of correlation functions in conformal field theories (CFTs). Carrollian amplitudes behave as correlation functions of Carrollian CFT that holographically describes asymptotically flat spacetime. We derive linear differential equations satisfied by Carrollian MHV gluon and graviton amplitudes. We obtain non-distributional solutions for both the gluon and graviton cases. We perform various consistency checks for these differential equations, including compatibility with conformal Carrollian symmetries.

Development and Analysis of novel Integrable Nonlinear Dynamical Systems on Quasi-One-Dimensional Lattices. Parametrically Driven Nonlinear System of Pseudo-Excitations on a Two-Leg Ladder Lattice

Following the main principles of developing the evolutionary nonlinear integrable systems on quasi-one-dimensional lattices, we suggest a novel nonlinear integrable system of parametrically driven pseudo-excitations on a regular two-leg ladder lattice. The initial (prototype) form of the system is derivable in the framework of semi-discrete zero-curvature equation with the spectral and evolution operators specified by the properly organized 3 × 3 square matrices. Although the lowest conserved local densities found via the direct recursive method do not prompt us the algebraic structure of system’s Hamiltonian function, but the heuristically substantiated search for the suitable two-stage transformation of prototype field functions to the physically motivated ones has allowed to disclose the physically meaningful nonlinear integrable system with time-dependent longitudinal and transverse inter-site coupling parameters. The time dependencies of inter-site coupling parameters in the transformed system are consistently defined in terms of the accompanying parametric driver formalized by the set of four homogeneous ordinary linear differential equations with the time-dependent coefficients. The physically meaningful parametrically driven nonlinear system permits its concise Hamiltonian formulation with the two pairs of field functions serving as the two pairs of canonically conjugated field amplitudes. The explicit example of oscillatory parametric drive is described in full mathematical details.

A feasible numerical computation based on matrix operations and collocation points to solve linear system of partial differential equations

In this work, a system of linear partial differential equations with constant and variable coefficients via Cauchy conditions is handled by applying the numerical algorithm based on operational matrices and equally-spaced collocation points. To demonstrate the applicability and efficiency of the method, four illustrative examples are tested along with absolute error, maximum absolute error, RMS error, and CPU times. The approximate solutions are compared with the analytical solutions and other numerical results in literature. The obtained numerical results are scrutinized by means of tables and graphics. These comparisons show accuracy and productivity of our method for the linear systems of partial differential equations. Besides, an algorithm is described that summarizes the formulation of the presented method. This algorithm can be adapted to well-known computer programs.