We provide methodology and numerical results for the Hattendorff differential equa- tion for the continuous time evolution of the variance of L(j)t , the loss at time t random variable for a multi-state process, given that the state at time t is j.

DGJ (Daftardar-Gejii-Jafaris) method is used to obtain numerical solution of the third order fractional differential equation. Providing the DGJ method converges, the approximate solution is a good and effective numerical result which is close to the exact solution or the exact solution. For this,the examples of the explaning the method are presented. The proposed method is implemented for the approximation solution of the third order nonlinear fractional partial differential equations. The method was shown to be unsuitable and inconsistent for an example of a nonlinear fractional partial differential equation depend on initial-boundary value conditions. The fact that these numerical results are not consistent can be explained by the fact that the method is not convergent.

The application of evolutionary optimization algorithms in problem solving is currently gaining wide popularity. Use of Differential Evolution (DE) algorithm in obtaining analytically approximate solution of unstable second order initial value Ordinary Differential Equation (ODE) is presented in this work. The methodology involves solving an associated problem of optimization with con-strains to get an analytically approximate solution for the ODE under consider-ation. Three test cases were used to demonstrate the efficiency of our method. In comparison with other methods discussed in the literature, our method gave significant improvement on the accuracy of the obtained results.

Homotopy perturbation and analysis methods have been widely used to obtain both approximate and exact so- lutions to nonlinear problems. In general, these two methods are based on the Taylor series with respect to an embedding parameter. Many researchers have compared the two methods and raised more concerns on the homo- topy perturbation method (HPM) because the homotopy analysis method (HAM) contains a convergence-control parameter ~: For this reason, in this article, a more general form of HPM is introduced as the -homotopy per- turbation method (-HPM), which contains a control parameter : The introduction of parameter  in this new modication gives a better way to adjust and control the convergence region and the rate of the series solution. We conrm through the given examples in this study that the HPM is a special case of the -HPM. The error and convergence analysis of this proposed method are also presented

In this work, the analytical treatment of Fourth order pd equations has been discussed with the variational iteration algorithm with an auxiliary parameter. The technique has the capability of reducing the size of computational work and easily overcomes the difficulty of the perturbation method or Adomian polynomials. Comparison with the classic variational iteration algorithm- is carried out, showing that the modification is more efficient and reliable.

In this short note we present a technique using which one attributes frequency and wavevector to (almost) arbitrary scalar fields. Our proposed definition is then applied to the classical wave equation to yield a novel nonlinear PDE.

In this paper, we investigate some kind of Dynkin game under $g$-expectation induced by backward stochastic differential equation (shortly for BSDE). We define the lower and upper value functions $\underline{V}_t=ess\sup\limits_{\tau\in{\mathcal{T}_t}} ess\inf\limits_{\sigma\in{\mathcal{T}_t}}\mathcal{E}^g_t[R(\tau,\sigma)]$ and $\overline{V}_t=ess\inf\limits_{\sigma\in{\mathcal{T}_t}}ess\sup\limits_{\tau\in{\mathcal{T}_t}}\mathcal{E}^g_t[R(\tau,\sigma)]$, respectively. Under some regular assumptions, a pair of saddle point is obtained and the value function of Dynkin game $V(t)=\underline{V}_t=\overline{V}_t$ follows. Furthermore, the constrained case of Dynkin game is also considered.

The main object of this work is to establish the existence and uniqueness of a solution to the 3D Navier-Stokes (NS) system for an incompressible fluid with viscosity. The nonlinearity of the NS system, as well as the need to estimate velocity and pressure for every value of the viscosity parameter [1] make them challenging to solve. In this regard, in the present work, we study the Navier-Stokes system, which describe the flow of a viscous incompressible fluid and the solution was obtained for velocity and pressure in an analytical form. In addition, the found pressure distribution law, which is described by a Poisson type equation and plays a fundamental role in the theory of Navier-Stokes systems in constructing analytic smooth (conditionally smooth) solutions.

A nonlinear singularly perturbed Cauchy problem with confluent fuchsian singularities is examined. This problem involves coefficients with polynomial dependence in time. A similar initial value problem with logarithmic reliance in time has been investigated by the author in a recent work, for which sets of holomorphic inner and outer solutions were built up and expressed as a Laplace transform with logarithmic kernel. Here, a family of holomorphic inner solutions are constructed by means of exponential transseries expansions containing infinitely many Laplace transforms with special kernel. Furthermore, asymptotic expansions of Gevrey type for these solutions relatively to the perturbation parameter are established.

We investigate a singularly perturbed q-difference differential Cauchy problem with polynomial coefficients in complex time and space with quadratic nonlinearity. We construct local holomorphic solutions on sectors in the complex plane with respect to the perturbation parameter with values in some Banach space of formal power series in space with analytic coefficients on shrinking domains in time. Two aspects of these solutions are addressed. One feature concerns asymptotic expansions in the parameter for which a Gevrey type structure is unveiled. The other fact deals with confluence properties as q tends to 1. In particular the built up Banach valued solutions are shown to merge in norm to a fully bounded holomorphic map in all its arguments that solves a nonlinear partial differential Cauchy problem.

We examine a family of nonlinear q-difference-differential Cauchy problems obtained as a coupling of linear Cauchy problems containing dilation q-difference operators, recently investigated by the author, and quasi-linear Kowalevski type problems that involve contraction q-difference operators. We build up local holomorphic solutions to these problems. Two aspects of these solutions are explored. One facet deals with asymptotic expansions in the complex time variable for which a mixed type Gevrey and q-Gevrey structure is exhibited. The other feature concerns the problem of confluence of these solutions as q tends to 1.

A linear Cauchy problem with polynomial coefficients wich combines q-difference operators for q>1 and differential operators of irregular type is examined. A finite set of sectorial holomorphic solutions w.r.t the complex time is constructed by means of classical Laplace transforms. These functions share a common asymptotic expansion in the time variable which turns out to carry a double layers structure which couples q-Gevrey and Gevrey bounds. In the last part of the work, the problem of confluence of these solutions as q tends to 1 is investigated.

Polynomial recurrence bounds for a class of stochastic differential equations with a gradient type drift and an additive Wiener process are studied without Lyapunov functions.

The differential transformation, an approach based on Taylor's theorem, is proposed as convenient for finding exact or approximate solution to the initial value problem with multiple Caputo fractional derivatives of generally non-commensurate orders. The multi-term differential equation is first transformed into a multi-order system and then into a system of recurrence relations for coefficients of formal fractional power series. The order of the fractional power series is discussed in relation to orders of derivatives appearing in the original equation. Application of the algorithm to an initial value problem results in a reliable and expected outcome.

Uncertain fractional differential equation (UFDE) is an excellent tool for describing complicated dynamic systems. This paper analyzes the valuation problems of currency option based on UFDE under optimistic value criterion. Firstly, an uncertain fractional currency model is formulated to describe the dynamics of foreign exchange rate. Then, the pricing formulae of European, American and Asian currency option are obtained under optimistic value criterion, respectively. Finally, numerical simulations are provided to illustrate the validity of our theoretical results.

In this paper, we consider with a class of system of differential equations whose argument transforms are involution. In this an initial value problem for a differential equation with involution is reduced to an initial value problem for a higher order ordinary differential equation. Than either two initial conditions are necessary for a solution, the equation is then reduced to a boundary value problem for a higher order ODE.

A mathematical expression to describe cation absorption of root is expressed with simulation results. The root cells selectively emit H+ ions with active transport consuming ATPs to establish electrical gradient. The gradient promotes external positive ions to flow into the roots, while carries negatively charged particles with symport. In this paper, a mathematical model whose independent variables are the concentrations of external and internal cation is presented. This differential equation is induced from Ohm’s law. The equation has terms for plant physiology, ion’s physical and electrical properties, growth of plant, and interaction between the root and the surroundings. Simulation showed that the physiology-related coefficient has important role on nutrition absorption.

In this paper we use the conformable fractional derivative to discuss some fractional linear differential equations with constant coefficients. By applying some similar arguments to the theory of ordinary differential equations, we establish a sufficient condition to guarantee the reliability of solving constant coefficient fractional differential equations by the conformable Laplace transform method. Finally, we analyze the analytical solution for a class of fractional models associated with Logistic model, Von Foerster model and Bertalanffy model is presented graphically for various fractional orders and solution of corresponding classical model is recovered as a particular case.

We consider a family of nonlinear singularly perturbed PDEs whose coefficients involve a logarithmic dependence in time with confluent Fuchsian singularities that unfold an irregular singularity at the origin and rely on a single perturbation parameter. We exhibit two distinguished finite sets of holomorphic solutions, so-called outer and inner solutions, by means of a Laplace transform with special kernel and Fourier integral. We analyze the asymptotic expansions of these solutions relatively to the perturbation parameter and show that they are (at most) of Gevrey order 1 for the first set of solutions and of some Gevrey order that hinges on the unfolding of the irregular singularity for the second.

The differential-equation of second order for space-time oscillation in flat tangential-spacetime of Minkowski-type can be written as an energy-equation. For this energy-equation the analogon of damping energy in two coupled Minkoswki-spacetimes with ftl can be calculated and formulated over Lamberts W-function or via Euler-equation.

The article outlines a creative approach by fusing ideas from fuzzy contractions, graph mappings, and Kannan mappings to create a brand-new notion known as "Kannan-graph-fuzzy contraction." This combination of concepts makes it possible to create a unique framework for resolving real-world issues controlled by nonlinear models. It is anticipated that using Kannan-graph-fuzzy contractions will offer a fresh approach to discussing solution existence in situations where these ideas overlap. This strategy might have a significant impact on a lot of real-world applications, proving how practical multidisciplinary techniques are for tackling issues in the real world originating from technology and engineering. The article includes several illustrated examples based on computer simulations, making the results more explicit.

This study focuses on the optimal deployment problem, and determines the optimal size of a military force to send to the battle field. The decision is optimized, based on an objective function, that considers the cost of deployment, the cost of the time it takes to win the battle, and the costs of killed and wounded soldiers with equipment. The cost of deployment is modeled as an explicit function of the number of deployed troops and the value of a victory with access to a free territory, is modeled as a function of the length of the time it takes to win the battle. The cost of lost troops and equipment, is a function of the size of the reduction of these lives and resources. An objective function, based on these values and costs, is optimized, under different parameter assumptions. The battle dynamics is modeled via the Lanchester differential equation system based on the principles of directed fire. First, the deterministic problem is solved analytically, via derivations and comparative statics analysis. General mathematical results are reported, including the directions of changes of the optimal deployment decisions, under the influence of alternative types of parameter changes. Then, the first order optimum condition from the analytical model, in combination with numerically specified parameter values, is used to determine optimal values of the levels of deployment in different situations. A concrete numerical case, based on documented facts from the Battle of Iwo Jima, during WW II, is analyzed, and the optimal US deployment decisions are determined under different assumptions. The known attrition coefficients of both armies, from USA and Japan, and the initial size of the Japanese force, are parameters. The analysis is also based on some parameters without empirical documentation, that are necessary to include to make optimization possible. These parameter values are motivated in the text. The optimal solutions are found via Newton- Raphson iteration. Finally, a stochastic version of the optimal deployment problem is defined. The attrition parameters are considered as stochastic, before the deployment decisions have been made. The attrition parameters of the two armies have the same expected values as in the deterministic analysis, are independent of each other, have correlation zero, and have relative standard deviations of 20%. All possible deployment decisions, with 5000 units intervals, from 0 to 150000 troops, are investigated, and the optimal decisions are selected. The analytical, and the two numerical, methods, all show that the optimal deployment level is a decreasing function of the marginal deployment cost, an increasing function of the marginal cost of the time to win the battle, an increasing function of the marginal cost of killed and wounded soldiers and lost equipment, an increasing function of the initial size of the opposing army, an increasing function of the efficiency of the soldiers in the opposing army and a decreasing function of the efficiency of the soldiers in the deployed army. With stochastic attrition parameters, the stochastic model also shows that the probability to win the battle is an increasing function of the size of the deployed army. When the optimal deployment level is selected, the probability of a victory is usually less than 100%, since it would be too expensive to guarantee a victory with 100% probability. Some of many results of relevance to the Battle of Iwo Jima, are the following: In the deterministic Case 0 analysis, the optimal US deployment level is 66200, the time to win the battle is 30 days and 14000 US soldiers are killed or wounded. If the marginal cost of the time it takes to win a victory doubles, the optimal deployment increases to 75400, the time to win a victory is reduced to 26 days, and less than 12000 soldiers are killed or wounded. In the stochastic Case 0 analysis, the optimal US deployment level is 65000, the expected time to win the battle is 46 days and almost 25000 US soldiers are expected to be killed or wounded. If the cost per killed or wounded soldier increases from 0 to 5 M $US, the optimal deployment level increases to 75000. Then, the victory is expected to appear after 35 days and 19900 US soldiers are expected to be killed or wounded.

We enhance Frobenius' method for solving linear ordinary differential equations about regular singular points. Key to Frobenius' approach is the exploration of the derivative with respect to a single parameter; this parameter is introduced through the powers of generalized power series. Extending this approach, we discover that tandem recurrence relations can be derived. These relations render coefficients for series occurring in logarithmic solutions. The method applies to the, practically important, exceptional cases in which the roots of the indicial equation are equal, or differ by a non-zero integer. We demonstrate the method on Bessel's equation and derive previously unknown tandem recurrence relations for coefficients of solutions of the second kind, for Bessel equations of all integer and half-integer order.

In the note, the authors find several simple and explicit forms for a family of inhomogeneous linear ordinary differential equations studied in "D. Lim, Differential equations for Daehee polynomials and their applications, J. Nonlinear Sci. Appl. 10 (2017), no. 4, 1303-1315; Available online at http://dx.doi.org/10.22436/jnsa.010.04.02".

Recently, there has been a growing interest in the production of electricity from renewable energy sources (RES). The RES investment is characterized by uncertainty, which is long-term, costly, depend on feed-in-tariff and support schemes. In this paper, we address the real option valuation (ROV) of a solar power plant investment. The real option framework is investigated. This framework considers the renewable certificate price, furthermore the cost of delay between establishing and operating the solar power plant. The optimal time of launching the project and assess the value of deferred option are discussed. The new three stage numerical methods are constructed, the Lobatto3C-Milstein (L3CM) methods. The numerical methods are integrated with concept of Black-Scholes option pricing theory, and applied in option valuation for solar energy investment with uncertainty. The numerical results of L3CM, finite difference and Monte Carlo methods are compared to show the efficiency of our methods. Our data set refers to the Arab Republic of Egypt.

Probability is an open question in the ontological interpretation of quantum mechanics. It has been discussed in some trajectory interpretations such as Bohmian mechanics and stochastic mechanics. New questions arise when the domain of probability extends to the complex space, including the generation of complex trajectories, the definition of the complex probability, the relation of the complex probability to the real quantum probability, and so on. The complex treatment proposed here applies the optimal quantum guidance law to derive the stochastic differential (SD) equation governing the particle’s random motions in the complex plane. The ensemble of the complex quantum random trajectories (CQRTs) solved from the complex SD equation is used to construct the probability distribution ρc(t,x,y) of the particle’s position over the complex plane z=x+iy. The correctness of the obtained complex probability is confirmed by the solution of the complex Fokker-Planck equation. The significant contribution of the complex probability is that it can be used to reconstruct both quantum probability and classical probability, and to clarify their relationship. Although quantum probability and classical probability are both defined on the real axis, they are obtained by projecting complex probability onto the real axis in different ways. This difference explains why the quantum probability cannot exactly converge to the classical probability when the quantum number is large.

In this paper, a class of nonlinear ordinary differential equations with impulses at variable times is considered. The existence and uniqueness of solution are given. At the same time, modifying the classical definitions of continuous dependence and Ga^teaux differentiability, some results on continuous dependence and Ga^teaux differentiable of solution relative to the initial value also are presented in new topology sense. For the autonomous impulsive system, the periodicity of solution is given. As an application, properties of solution for a type of controlled nonlinear ordinary differential equation with impulses at variable times is obtained. These results are foundation to study optimal control problems of systems governed by the differential equations with impulses at variable times.

This paper is concerned with variational methods for open quantum systems with Markovian dynamics governed by Hudson-Parthasarathy quantum stochastic differential equations. These QSDEs are driven by quantum Wiener processes of the external bosonic fields and are specified by the system Hamiltonian and system-field coupling operators. We consider the system response to perturbations in these operators and introduce a transverse Hamiltonian which encodes the propagation of the perturbations through the unitary system-field evolution. This approach provides an infinitesimal perturbation analysis tool which can be used for the development of optimality conditions in quantum control and filtering problems. Such settings employ, as performance criteria, quadratic (or more complicated) cost functionals of the system and field variables to be minimised over the energy and coupling parameters of system interconnections. We demonstrate an application of the transverse Hamiltonian variational approach to a mean square optimal coherent quantum filtering problem for a measurement-free field-mediated cascade connection of a quantum system with a quantum observer.

The paper's main goal is to introduce a novel type of interpolative convex contraction and build up some fresh findings for it using the interpolative convex contractions' progressive approach. We have established certain conclusions using orbitally S-complete and Suzuki type contractions in F-metric spaces. My research aims to examine the fixed point method's efficacy in solving fractional differential equation with boundary conditions.

The methods for constructing solutions to integro-differential equations of the Volterra type are considered. The equations are related to fractional conformable derivatives. Explicit solutions of homogeneous and inhomogeneous equations are constructed and a Cauchy-type problem is studied. It should be noted that the considered method is based on the construction of normalized systems of functions with respect to a differential operator of fractional order.

Here we are concerning with two problems of a coupled system of random and stochastic nonlinear differential equations with two coupled systems of nonlinear nonlocal random and stochastic integral conditions. The existence of solutions will be studied. The sufficient condition for the uniqueness of the solution will be given. The continuous dependence of the unique solution on the nonlocal conditions will be proved.

In this paper, we introduce some basic idea of Variational iteration method for short (VIM) to solve the Volterra’s integro-differential equations. The VIM is used to solve effectively, easily, and accurately a large class of non-linear problems with approximations which converge rapidly to accurate solutions. For linear problems, it’s exact solution can be obtained by only one iteration step due to the fact that the Lagrange multiplier can be exactly identified. It is to be noted that the Lagrange multiplier reduces the iteration on integral operator and also minimizes the computational time. The method requires no transformation or linearization of any forms. Two numerical examples are presented to show the effectiveness and efficiency of the method. Also, we compare the result with the result from Homotopy perturbation method (HPM). Finally, we investigate the absolute difference between variational iteration method and homotopy perturbation method and draw the graph of difference function by using Mathematica.

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the presence of both steady-state and transient shear and uniaxial elongational flows, supplemented by a variable link tension coefficient. Here, we provide the transient and stationary solutions of the tumbling-snake model under biaxial elongation both analytically, for small and large elongation rates, and via Brownian dynamics simulations, for the case of planar elongational flow over a wide range of rates, times, and the model parameters. We show that both the steady-state and transient first planar viscosity predictions are similar to their uniaxial counterparts, in accord with recent experimental data. The second planar viscosity seems to behave in all aspects similarly to the shear viscosity, if shear rate is replaced by elongation rate.

In this paper, we consider an epidemiological Susceptible-Infected-Susceptible (SIS) model in which the individuals cured from the infection return back into the group of people endangered by the infection. Delay differential equations were introduced to create more realistic models since many processes depend on past history. Many different mathematical methods, including the Susceptible-Infected-Recovered (SIR) model as well as derivatives thereof, can be used to prognose the spread of many various infectious diseases, including Covid-19.We use a nonlinear differential system with time delay as a mathematical formulation of this model. It is examined by the existence of a positive solution which is bounded by two functions $k_1, k_2\in C([t_0,\infty),(0,\infty))$ with constants $\lambda_1 > 0,\,\lambda_2 > 0$ where $0<\lambda_1\leq \lambda_2$. The conditions are put in place for existence of a positive solution that approaches zero in infinity, if $t\to\infty$. Infectious diseases accompany a person throughout their entire life, and in the history of mankind often appeared in frequent epidemics caused high mortality. The occurrence of infectious diseases in the population represents a serious health, social and economic problem. There is a lot of ambiguity concerning the time span and final spread of the epidemic thus making it very demanding for the governing bodies, healthcare institutions, and economic sector to precisely judge the future development.Several examples are shown at the end of this paper.

The paper proposes a new principle of finding and removing elements of mathematical model, redundant in terms of parametric identification of the model. It allows reducing computational and time complexity of the applications built on the model. Especially this is important for AI based systems, systems based on IoT solutions, distributed systems etc. Besides, the complexity reduction allows increasing an accuracy of mathematical models implemented. Despite the model order reduction methods are well known, they are extremely depended however on the problem area. Thus, proposed reduction principles can be used in different areas, what is demonstrated in this paper. The proposed method for the reduction of mathematical models of dynamic systems allows also the assessment of the requirements for the parameters of the simulator elements to ensure the specified accuracy of dynamic similarity. Efficiency of the principle is shown on the ordinary differential equations and on the neural network model. The given examples demonstrate efficient normalizing properties of the reduction principle for the mathematical models in the form of neural networks.

We consider a model of predator-prey interaction at fractional-order where the predation obeys the ratio-dependent functional response and the prey is linearly harvested. For the proposed model, we show the existence, uniqueness, non-negativity as well as the boundedness of the solutions. Conditions for the existence of all possible equilibrium points and their stability criteria, both locally and globally, are also investigated. The local stability conditions are derived using the Magtinon's theorem, while the global stability is proven by formulating an appropriate Lyapunov function. The occurance of Hopf bifurcation around the interior point is also shown analytically. At the end, we implement the Predictor-Corrector scheme to perform some numerical simulations.

In this paper, we define new kind of Daehee numbers, the degenerate Daehee numbers of the third kind, using degenerate log function as generating function. We obtain some identities for the degenerate Daehee numbers of the third kind associated with the Daehee, degenerate Daehee and degenerate Daehee numbers of the second kind. In addition, we derive a differential equation associated with the degenerate log function. We deduce some identities from the differential equation.

The Chebyshev collocation method (ChC) implemented as Chebfun is used in order to solve a class of second order one-dimensional singular and genuinely nonlinear boundary value problems. Efforts to solve these problems with conventional ChC have generally failed, and the outcomes obtained by finite differences or finite elements are seldom satisfactory. We try to fix this situation using the new Chebfun programming environment. However, for toughest problems we have to loosen the default Chebfun tolerance in Newton's solver as the ChC runs into trouble with ill-conditioning of the spectral differentiation matrices. Although in such cases the convergence is not quadratic the Newton updates decrease monotonically. This fact, along with the decreasing behaviour of Chebyshev coefficients of solutions, suggests that the outcomes are trustworthy, i.e., the collocation method has exponential (geometric) rate of convergence or at least an algebraic rate. We consider first a set of problems that have exact solutions or prime integrals and then another set of benchmark problems that do not possess these properties. Actually, for each test problem carried out we have determined how Chebfun solution converges, its length, the accuracy of the Newton method and especially how well the numerical results overlap with the analytical ones (existence and uniqueness).

The purpose of this paper is to find the time optimal control to a target set for semilinear stochastic functional differential equations involving time delays or memories under general conditions on a target set and nonlinear terms even though the equations contain unbounded principal operators. Our research approach is construct a fundamental solution for corresponding linear systems and establish variations of constant formula of solutions for given stochastic equations. The existence result of time optimal controls for one point target set governed by the given semilinear stochastic equation is also established.

In this article, we present a comparative study between the Adomian Decomposition Method (ADM) and Variational Iteration Method (VIM). The study outlines the significant features of the two methods, for solving nonlinear Volterra-Fredholm integro-differential equations. From the computational viewpoint, the VIM is more efficient, convenient and easy to use. Moreover, we proved the existence and uniqueness results and convergence of the solution. Finally, an example is included to demonstrate the validity and applicability of the proposed techniques.

This document focuses attention on the fundamental solution of an autonomous linear retarded functional differential equation (RFDE) along with its supporting cast of actors: kernel matrix, characteristic matrix, resolvent matrix, and the Laplace transform. The fundamental solution is presented in a form of the convolutional powers of the kernel matrix in the manner of a convolutional exponential matrix function. The fundamental solution combined with a solution representation gives an exact expression in explicit form for the solution of a RFDE. Algebraic graph theory is applied to the RFDE in the form of a weighted loop-digraph to illuminate the system structure and system dynamics and to identify the strong and weak components. Examples are provided in the document to elucidate the behavior of the fundamental solution. The paper introduces fundamental solutions of other functional differential equations.

To identify the exchange coefficients of a nonlinear polynomial tricompartmental general system of α+β , we follow the procedure; firstly, the recommended solution is to introduce an adequate time t*&gt;0 in a determinable manner. That is, after injecting quantity a into the main compartment, wait a moment for the exchange in the polynomial α+β order nonlinear general system to settle, and then compare this compartment to compartment 2 at t*. Secondly, applying the Taylor formula will linearize the system and identify the exchange coefficients. Finally, we will prove that the linearization method is stable.

We consider a Cauchy problem for differential equations in a Hilbert space $X$. The problem is stated in a time interval $I$, which can be finite or infinite. Weuse a fixed point argument for history-dependent operators to prove the unique solvability of the problem. Then, we state and prove convergence criteria for both a general fixed point problem and thecorresponding Cauchy problem. These criteria provide necessary and sufficient conditions on a sequence $\{u_n\}$ which guarantee its convergence to the solution of the corresponding problem, in the space of both continuous and continuously differentiable functions. We then specify our results in the study of a particular differential equation governed by two nonlinear operators. Finally, we provide an application in viscoelasticity and give a mechanical interpretation of the corresponding convergence result.

In a previous paper, we solved the partial differential equation of Mullins problem in the case of the evaporation-condensation in electronic devices and gave an exact solution relative to the geometric profile of the grain boundary grooving when materials are submitted to thermal and mechanical solicitation and fatigue effect. In this new research, new modelling of the grain groove profile was proposed and new analytical expressions of the groove profile, the derivative and the groove depth were obtained in the case of diffusion in thin polycrystalline films by the resolution of the fourth differential equation formulated by Mullins that supposed 〖y'〗^2≪1. The obtained analytical solution gave more accurate information on the geometric characteristics of the groove that were necessary to study the depth and the width of the groove. These new findings will open a new way to study with more accuracy the problem of the evaporation-condensation combined to the diffusion phenomenon on the material surfaces with the help of the analytical solutions.

In this paper, we compile the fractional power series method and the Laplace transform to design a new algorithm for solving the fractional Volterra integro-differential equation. For that, we assume the Laplace power series (LPS) solution in terms of power q=1m,m∈Z+, where the fractional derivative of order α=qγ for which γ∈Z+. This assumption will help us to write the integral, the kernel, and the nonhomogeneous terms as a LPS with the same power. The recurrence relations for finding the series coefficients can be constructed using this form. To demonstrate the algorithm's accuracy, the residual error is defined and calculated for several values of the fractional derivative. Two strongly nonlinear examples are discussed to provide the efficiency of the algorithm. The algorithm gains powerful results for this kind of problem. Under Caputo meaning the obtained results are illustrated numerically, and graphically. Geometrically, the behavior of the solution declares that the changing of the fractional derivative parameter values' in their domain alters the style of the attained solution in a symmetrical meaning and fully coinciding to the ordinary derivative value'. From these simulations, the results report that the recommended novel algorithm is a straightforward, accurate, and superb tool to generate analytic-approximate solutions for Integral, and integro-differental equations of fractional order.

The paper presents an adaptation of numerical solution of first order linear differential equation in fuzzy environment. The numerical method is re-established and studied with fuzzy concept to estimate its uncertain parameters whose values are not precisely known. Demonstrations of fuzzy solutions of the governing methods are carried out by the approaches, namely Modified Runge Kutta method and Runge Kutta Merson method. The results are compared with the exact solution which is found using generalized Hukuhara derivative (gH-derivative) concepts. Additionally, different illustrative examples and an application in industry of the methods are also undertaken with the useful table and graph to show the usefulness for attained to the proposed approaches.

In the paper, by virtue of techniques in combinatorial analysis, the authors simplify two families of nonlinear ordinary differential equations in terms of the Stirling numbers of the first kind.

The paper examines rational Darboux-Crum transforms (RDCTs) of the Jacobi-reference (JRef) potential on the line, with an emphasis on the transformations using only normalizable seed solutions. The Cooper-Ginocchio-Khare (CGK) potential constitutes the simplest example of this rich family of exactly solvable one-dimensional (1D) potentials. It was explicitly confirmed that that an arbitrary Darboux-Crum transformation (DCT) of the given Sturm-Liouville equation (SLE) can be represented as a chain of sequential Darboux deformations of the associated Liouville potentials. If the rational SLE (RSLE) under consideration is exactly solvable and the transformation function (TF) used at each step is nodeless then each RSLE constructed in such a way was proven being exactly solvable itself. The exact solvability of RDCTs of the JRef canonical SLE (CSLE) was also confirmed for the DCTs using pairs of juxtaposed eigenfunctions (an extension of the renowned Adler theorem to the RCSLEs of our interest). By re-expressing the reference polynomial fraction (RefPF) of the transformed CSLE in terms of the Krein determinant (KD), instead of Crum Wronskian (CW), it was shown that the eigenfunctions under consideration are expressible in terms of polynomial solutions of a generalized Heun equation with positions of the most poles dependent on the exponent parameters. The outlined technique was extended to two isospectral SUSY partners of the JRef potential (as well as to its third sibling with the inserted new lowest-energy level) which, by analogy with the CGK potential, are quantized by polynomial solutions of the 2-Heun equation with all four poles at some fixed positions on the real axis.

The susceptible-infected-recovered/removed-vaccinated (SIRV) epidemics model is an important generalization of the SIR epidemics model as it accounts quantitatively for the effects of vaccination campaigns on the temporal evolution of epidemics outbreaks. Additional to the time-dependent infection ($a(t)$) and recovery ($\mu (t)$) rates, regulating the transitions between the compartments $S\to I$ and $I\to R$, respectively, the time-dependent vaccination rate $v(t)$ accounts for the transition between the compartments $S\to V$ of susceptible to vaccinated fractions. An accurate analytical approximation is derived for arbitrary and different temporal dependencies of the rates which is valid for all times after the start of the epidemics for which the cumulative fraction of new infections $J(t)\ll 1$. As vaccination campaigns automatically reduce the rate of new infections by transferring susceptible to vaccinated persons, the limit $J(t)\ll 1$ is even better fulfilled than in the SIR-epidemics model. The comparison of the analytical approximation for the temporal dependence of the rate of new infections $\jt(t)=a(t)S(t)I(t)$, the corresponding cumulative fraction $J(t)$, and $V(t)$, respectively, with the exact numerical solution of the SIRV-equations for different illustrative examples proves the accuracy of our approach. The considered illustrative examples include the cases of stationary ratios with a delayed start of vaccinations, and an oscillating ratio of recovery to infection rate with a delayed vaccination at constant rate. The proposed analytical approximation is self-regulating as the final analytical expression for the cumulative fraction $\jinf $ after infinite time allows us to check the validity of the original assumption $J(t)\le \jinf \ll 1$.

This study develops a predictor-corrector algorithm for the numerical simulation of IVPs involving singular generalized fractional derivatives with Mittag-Leffler kernels. The proposed algorithm converts the considered IVP into a Volterra-type integral equation and then uses Trapezoidal rule to obtain approximate solutions. Numerical approximate solutions of some singular generalized fractional derivative with Mittag-Leffler kernels models have been presented to demonstrate the efficiency and accuracy of the proposed algorithm. The algorithm describes the influence of the fractional derivative parameters on the dynamics of the studied models. The suggested method is expected to be effectively employed in the field of simulating generalized fractional derivative models

In the present study, a mathematical model of non-steady partial differential equation from the process of oxygen mass transport in the human pulmonary circulation is proposed. Mathematical modelling of this kind of problems lead to a non-steady partial differential equation and for its numerical simulation, we have used finite differences. The aim of the process is the exact numerical analysis of the study, wherein consistency, stability and convergence is proposed. The necessity of doing the process is that, we would like to increase the order of numerical solution to a higher order scheme. An increment in the order of numerical solution makes the numerical simulation more accurate, also makes the numerical simulation being more complicated. In addition, the process of numerical analysis of the study in this order of solution needs more research work.

We introduce a novel stochastically correlated two-factor (i.e., bivariate) diffusion process under the square-root format, for which we obtained, analytically, the corresponding solutions for the conditional moment generating functions and conditional characteristic functions. Such solutions recover verbatim those of the uncorrelated case. We believe that the The aforementioned model encompasses a range of processes similar to that produced by a bivariate square-root process in which entries are correlated in the standard way, i.e., via a constant correlation coefficient. Note that, for the latter, closed-form solutions for the conditional characteristic functions are not available. We focus on the financial scenario of obtaining closed-form expressions for the exact price of a zero-coupon bond and Asian option prices via a Fourier cosine series method.

This paper delves into the dynamics of asset pricing within Bachelier’s market model (BMM), elucidating the representation of risky asset price dynamics and the definition of riskless assets. It highlights the fundamental differences between BMM and the Black-Scholes-Merton market model (BSMMM), including the extension of BMM to handle assets yielding a simple dividend. Our investigation further explores Bachelier’s term structure of interest rates (BTSIR), introducing a novel version of Bachelier’s Heath-Jarrow-Morton model and adapting the Hull-White interest rate model to fit BMM. The study concludes by examining the applicability of BMM in real-world scenarios, such as those involving environmental, social, and governance (ESG)-adjusted stock prices and commodity spreads.

This paper aims at computing feasible control strategies and the corresponding feasible state trajectories to drive an autonomous rear-axle bicycle robot from a given initial state to a final state such that the total running cost is minimized. Pontryagin’s Minimum Principle is applied and derives the optimality conditions from which the feasible control functions, expressed as functions of state and costate variables, are substituted into the combined state-costate system to obtain a new free-control state-costate nonlinear system of ordinary differential equations. A computer program was written in Scilab to solve the combined state-costate system and obtain the feasible state functions, the feasible costate functions and the feasible control functions. Associated Computational Simulations were provided to show the effectiveness and the reliability of the approach.

In the paper, the authors consider the generating functions of the Hermite polynomials and their squares, present explicit formulas for higher order derivatives of the generating functions of the Hermite polynomials and their squares, which can be viewed as ordinary differential equations or derivative polynomials, find differential equations that the generating functions of the Hermite polynomials and their squares satisfy, and derive explicit formulas and recurrence relations for the Hermite polynomials and their squares.

The classical Lotka-Volterra predator-prey model is globally stable and uniformly persistent. However, in real-life biosystems extinction of species due to stochastic effects is possible. In this paper, we consider the classical Lotka-Volterra predator-prey model under stochastic perturbations. For this model, using an analytical technique based on the direct Lyapunov method and a development of the ideas of R.Z. Khasminskii, we find the precise sufficient conditions for the stochastic extinction of one and both species and the precise necessary conditions for the stochastic system persistence. The stochastic extinction occurs via a process known as the stabilization by noise of the Khasminskii type, and, in order to establish the sufficient conditions for extinction, we found the conditions for the stabilization. The analytical results are illustrated by numerical simulations. {\bf Keywords:} stochastic perturbations, white noise, Ito's stochastic differential equation, the Lyapunov functions method, stability in probability, stabilization by noise, stochastic extinction, persistence.

Many problems in the fields of finance and actuarial science can be transformed into the problem of solving backward stochastic differential equations (BSDE) and partial differential equations (PDE) with jumps, which are often difficult to solve in high-dimensional cases. To solve this problem, this paper applies the deep learning algorithm to solve a class of high-dimensional nonlinear partial differential equations with jump terms and their corresponding backward stochastic differential equations (BSDE) with jump terms. Using the nonlinear Feynman-Kac formula, the problem of solving this kind of PDE is transformed into the problem of solving the corresponding backward stochastic differential equations with jump terms, and the numerical solution problem is turned into a stochastic control problem. At the same time, the gradient and jump process of the unknown solution are separately regarded as the strategy function and they are approximated respectively by using two multilayer neural networks as function approximators. Thus, deep learning-based method is used to overcome the “curse of dimensionality” caused by high-dimensional PDE with jump, and the numerical solution is obtained. In addition, this paper proposes a new optimization algorithm based on the existing neural network random optimization algorithm, and compares the results with the traditional optimization algorithm, and achieves good results. Finally, the proposed method is applied to three practical high-dimensional problems: Hamilton-Jacobi-Bellman equation, bond pricing under the jump Vasicek model and high-dimensional option pricing model with default risk. The proposed numerical method has obtained satisfactory accuracy and efficiency. The method has important application value and practical significance in investment decision-making, option pricing, insurance and other fields.

In the paper, the authors view some ordinary differential equations and their solutions from the angle of (the generalized) derivative polynomials and simplify some known identities for the Bernoulli numbers and polynomials, the Frobenius-Euler polynomials, the Euler numbers and polynomials, in terms of the Stirling numbers of the first and second kinds.

This paper develops feasible control strategies and associated system responses to bring an autonomous front-axle bicycle robot from specified initial conditions to final conditions such that a specific performance index is minimized. To solve the problem, the following approach is used: The feasible controls derived from the normal equations of optimality are substituted into the state and the costate systems and form a combined control-free state-costate system which is vectorized to enable and ease the application of a numerical method. A computer program written in Matlab computer programming language, codes a fourth-order Runge-Kutta numerical method and then solve the combined state-costate system of ordinary differential equations. The obtained results are the feasible bicycle robot trajectory, the feasible state functions, the feasible costate functions and the feasible control functions. Associated Computational Simulations are designed and provided to persuade on the effectiveness and the reliability of the approach.

The work is devoted to the study of a family of linear initial value problems of partial differential equations in the complex domain, dealing with two complex time variables. The use of a truncated Laplace-like transformation in the construction of the analytic solution allows to overcome a small divisor phenomenon arising from the geometry of the problem and represents an alternative approach to the one proposed in a recent work by the first two authors. The result leans on the application of a fixed point argument and the classical Ramis-Sibuya theorem.

This paper aims at computing feasible control strategies and the corresponding feasible state trajectories to drive an autonomous rear-axle bicycle robot from a given initial state to a final state such that the total running cost is minimized. Pontryagin’s Minimum Principle is applied and derives the optimality conditions from which the feasible control functions, expressed as functions of state and costate variables, are substituted into the combined state-costate system to obtain a new free-control state-costate nonlinear system of ordinary differential equations. A computer program was written in Scilab to solve the combined state-costate system and obtain the feasible state functions, the feasible costate functions and the feasible control functions. Associated Computational Simulations were provided to show the effectiveness and the reliability of the approach.

This paper is a slightly modified, abridged version of a previous work ``Parametric Gevrey asymptotics in two complex time variables through truncated Laplace transforms'' motivated by our contribution in the conference ``Formal and Analytic Solutions of Diff. (differential, partial differential, difference, q-difference, q-difference-differential) Equations on the Internet'' (FASnet20). It aims to clarify and give further detail at some crucial points concerning the asymptotic behavior of the solutions of the problems studied in that work.

This paper is a slightly modified, abridged version of a previous work "Parametric Gevrey asymptotics in two complex time variables through truncated Laplace transforms'' motivated by our contribution in the conference "Formal and Analytic Solutions of Diff. (differential, partial differential, difference, q-difference, q-difference-differential) Equations on the Internet'' (FASnet20). It aims to clarify and give further detail at some crucial points concerning the asymptotic behavior of the solutions of the problems studied in that work.

Hamiltonians describing the relativistic quantum dynamics of a particle with an arbitrary spin are shown to exhibit a supersymmetric structure when the even and odd elements of the Hamiltonian commute. For such supersymmetric Hamiltonians an exact Foldy-Wouthuysen transformation exits which brings it into a block-diagonal form separating the positive and negative energy subspaces. Here the supercharges transform between energy eigenstates of positive and negative energy. The relativistic dynamics of a charged particle in a magnetic field is considered for the case of a scalar (spin-zero) boson obeying the Klein-Gordan equation, a Dirac (spin one-half) fermion and a vector (spin-one) boson characterised by the Proca equation.

Recent advances in differential topology single out four-dimensions as being special, allowing for vast varieties of exotic smoothness (differential) structures, distinguished by their handlebody decompositions, even as the coarser algebraic topology is fixed. Should the spacetime we reside in takes up one of the more exotic choices, and there is no obvious reason why it shouldn't, apparent pathologies would inevitably plague calculus-based physical theories assuming the standard vanilla structure, due to the non-existence of a diffeomorphism and the consequent lack of a suitable portal through which to transfer the complete information regarding the exotic physical dynamics into the vanilla theories. An obvious plausible consequence of this deficiency would be the uncertainty permeating our attempted description of the microscopic world. We tentatively argue here, that a re-inspection of the key ingredients of the phenomenological particle models, from the perspective of exotica, could possibly yield interesting insights. Our short and rudimentary discussion is qualitative and speculative, because the necessary mathematical tools have only just began to be developed.

In this article, we consider the three dimensional $\alpha$-fractional nonlinear delay differential system of the form \begin{align*} D^{\alpha}\left(u(t)\right)&=p(t)g\left(v(\sigma(t))\right),\\D^{\alpha}\left(v(t)\right)&=-q(t)h\left(w(t))\right),\\D^{\alpha}\left(w(t)\right)&=r(t)f\left(u(\tau(t))\right),~ t \geq t_0, \end{align*} where $0 < \alpha \leq 1$, $D^{\alpha}$ denotes the Katugampola fractional derivative of order $\alpha$. We have established some new oscillation criteria of solutions of differential system by using generalized Riccati transformation and inequality technique. The obtained results are illustrated with suitable examples.

The comparison of group means in latent variable models plays a vital role in empirical research in the social sciences. The present article discusses extensions of invariance alignment and Haberman linking concerning the choice of linking functions for comparisons of many groups. Robust linking functions are proposed for invariance alignment and robust Haberman linking that are particularly suited to item response data under partial invariance. In a simulation study, it is shown that both linking approaches have comparable performance, and in some conditions, the newly proposed robust Haberman linking outperforms invariance alignment.

This study demonstrates a global, non-parametric, non-iterative optimization of time-mean value of a kind of index vibrated by time-varying forcing. It is based on the fact that the (steady) forced vibration of non-autonomous ordinary differential equation systems is well approximated by an analytical solution when the amplitude of forcing is sufficiently small and its base state without forcing is linearly stable and steady. It is applied to optimize a time-averaged heat-transfer rate on a two-dimensional thermal convection field in a square cavity with horizontal temperature difference, and the globally optimal way of vibrational forcing, i.e. the globally optimal, spatial distribution of vibrational heat and vorticity sources, is first obtained. The maximized vibrational thermal convection corresponds well to the state of internal gravity wave resonance. In contrast, the minimized thermal convection is weak, keeping the boundary layers on both sidewalls thick.

A new type of ordinary differential equation is introduced and discussed, namely, the time-dependent order ordinary differential equations. These equations can be solved via fractional calculus and are mapped into Volterra integral equations of second kind with singular integrable kernel. The solutions of the time-dependent order differential equations smoothly deforms solutions of the classical integer order ordinary differential equations into one-another, and can generate or remove singularities. An interesting symmetry of the solution in relation to the Riemann zeta function and Harmonic numbers was also proved.

The pollution caused by nitrogen dioxide is a major environmental problem in China. This study introduces a new type of atmospheric trace gas remote-sensing instrument, an airborne fiber imaging spectrometer. This spectrometer has a spectral range of 300–410 nm and works in push-broom mode with a 30° field of view on a flight path. Flight experiments were conducted on December 30, 2022 and January 5, 2023, covering heavily polluted areas of east of Beijing and Tangshan. This equipment obtained the density distribution of NO2 over the flight area. The results showed that pollution was mainly concentrated in the Caofeidian area and at the power station in the north, and the main source of pollution was anthropogenic. Satellite and airborne data near the pollution points were compared, and the two datasets showed a positive correlation with a correlation coefficient of 0.78 and 0.7, on the two days, respectively. This study demonstrates the capability of an airborne fiber imaging spectrometer for NO2 regional emission remote sensing and identifying the pollution points.

This research presents a solution to the problem of planning the optimum area for economic crops by developed mathematical models and developed an algorithm to solve the problem of planning the optimum area by considered economic value for the maximize profit of farmers. The data were collected from farmers in 8 provinces in the northeastern region of Thailand. The 3 economic crops studied were rice, cassava and sugarcane. The solving problem methods were 1) Created mathematical models and solved the problems with Lingo V.11. 2) Improved Differential Evolution algorithms (I-DE) to solve the problems, which had 3 local search methods included (Swap, Cyclic Move and K-variable moves). The results of this study showed that in the small and medium problems instances, Lingo V.11 and DE provided equal profit outcome but DE was faster but in the large size of test instances DE generated better solution than that of Lingo v.11 when Lingo simulation time is set to 250 hours and DE simulation time has set to maximum 21.82 minutes. 2) Comparing DE and I-DE , I-DE outperforms DE in finding the better solution for all size of test instances (small, medium and large).

Up to now, Schrödinger equation, Klein-Gordon equation (KGE) and Dirac equation are believed the fundamental equations of quantum mechanics. Schrödinger equation has a defect that there is no NKE solutions. Dirac equation has positive kinetic energy (PKE) and negative kinetic energy (NKE) branches. Both branches should have low momentum, or nonrelativistic, approximations: one is Schrödinger equation and the other is NKE Schrödinger equation. KGE has two problems: it is an equation of second time derivative, and calculated density is not definitely positive. To overcome the problems, it should be revised as PKE and NKE decoupled KGEs. The fundamental equations of quantum mechanics after the modification have at least two merits. They are of unitary in that everyone contains the first time derivative and are symmetric with respect to PKE and NKE. This reflects the symmetry of the PKE and NKE matters, as well as matter and dark matter, of our universe. The problems of one-dimensional step potentials are resolved by means of the modified fundamental equations for a nonrelativistic particle.

This paper presents a comprehensive investigation into the applicability and performance of two prominent growth models, namely the Verhulst model and the Montroll model, in the context of modeling tumor cell growth dynamics. Leveraging the power of Physics-Informed Neural Networks (PINNs), we aim to assess and compare the predictive capabilities of these models against experimental data obtained from the growth patterns of tumor cells. We employ a dataset comprising detailed measurements of tumor cell growth to train and evaluate the Verhulst and Montroll models. By integrating PINNs, we not only account for experimental noise but also embed physical insights into the learning process, enabling the models to capture the underlying mechanisms governing tumor cell growth. Our findings reveal the strengths and limitations of each growth model in accurately representing tumor cell proliferation dynamics. Furthermore, the study sheds light on the impact of incorporating physics-informed constraints on the model predictions. The insights gained from this comparative analysis contribute to advancing our understanding of growth models and their applications in predicting complex biological phenomena, particularly in the realm of tumor cell proliferation.

Summary: Studying biological systems generally relies on computational modeling and simulation, e.g., for model-driven discovery and hypothesis testing. Progress in standardization efforts led to the development of interrelated file formats to exchange and reuse models in systems biology, such as SBML, the Simulation Experiment Description Markup Language (SED-ML), or the Open Modeling EXchange format (OMEX). Conducting simulation experiments based on these formats requires efficient and reusable implementations to make them accessible to the broader scientific community and to ensure the reproducibility of the results. The Systems Biology Simulation Core Library (SBSCL) provides interpreters and solvers for these standards as a versatile open-source API in Java™. The library simulates even complex bio-models and supports deterministic Ordinary Differential Equations (ODEs); Stochastic Differential Equations (SDEs); constraint-based analyses; recent SBML and SED-ML versions; exchange of results, and visualization of in silico experiments; open modeling exchange formats (COMBINE archives); hierarchically structured models; and compatibility with standard testing systems, including the Systems Biology Test Suite and published models from the BioModels and BiGG databases. Availability and implementation: SBSCL is freely available at https://draeger-lab.github.io/SBSCL/.

In this paper, we investigate a new subclass of univalent functions defined by a generalized differential operator, and obtain some interesting properties of functions belonging to the class R^m λµ ( α, β, γ, υ ).

In this paper, we employ fuzzy Sumudu transform for solving system of linear fuzzy differential equations with fuzzy constant coefficients. The system with fuzzy constant coefficients is interpreted under strongly generalized differentiability. For this purpose, new procedures for solving the system are proposed. A numerical example is carried out for solving system adapted from fuzzy radioactive decay model. Conclusion is drawn in the last section and some potential research directions are given.

Electrically charged particles such as Electrons and Protons carry electric, E, and magnetic, B, fields. In addition to these fields, Quantum Mechanics (QM) endows these particles with an `arcane and spooky' field --- the wavefunction. This wavefunction of QM is not only assumed to be separate but distinct from the electromagnetic field. We herein upend this view by demonstrating otherwise. That is, we demonstrate that the four components of the Dirac wavefunction, can be shown to not only be an intimate, but, a direct function of the electromagnetic field carried by the particle in question. Insofar as unity, depth in our understanding and insight into both Dirac and Maxwell's equations as major pillars of Modern Physics, we believe that this work may very well inch us one-step-closer to the truth.

The paper is considering of the basic differential equations of hydrodynamics: the equation of continuity and motion. On the simplest example, it is shown that the equation of continuity in a system with the equation of motion leads to contradictions and erroneous results of modeling. A more correct form of the continuity equation is described. It is shown that the equations of motion can be written in the form of complete differentials. Three possible integral forms of the equations of motion are presented. As a conclusion, the existence and smoothness of the solution of the Navier-Stokes equations are considered.

We discuss a remote sensing system that is used to simultaneously detect range-resolved differential absorption LIDAR (light detection and ranging; DIAL) signals and integrated path DIAL signals (IP-DIAL) from aerosol targets for ranges up to 22 km. The DIAL/IP-DIAL frequency converter consists of an OPO pumped at 1064 nm to produce light at 1.6 μm and operates at 100 Hz pulse repetition frequency. The probe light is free space coupled to a movable platform that contains one transmitter and two receiver telescopes. Hybrid photon counting/current systems increase the dynamic range for detection by two orders of magnitude. Range resolved and column integrated dry-air CO2 and CH4 mixing ratios are obtained from line shape fits of CO2 and CH4 centered at 1602.2 nm and 1645.5 nm, respectively and measured at 10 different frequencies over ≈ 1.5 cm-1 bandwidth. The signal-to-noise ratios (SNR) of the IP-DIAL returns from cloud aerosols approach 1000:1 and the uncertainties in the mixing ratios weighted according to the integrated counts over the cloud segments range from 0.1 % to 1 %. The range averaged DIAL mixing ratios are in good agreement with the IP-DIAL mixing ratios at the 1 % to 2 % level for both CO2 and CH4. These results can serve as a validation method for future active and passive satellite observational systems.

In this paper, we established the existence of PC-mild solutions for non local fractional impulsive functional integro-differential equations with finite delay. The proofs are obtained by using the techniques of fixed point theorems, semi-group theory and generalized Bellman inequality. In this paper, we have used the distributed characteristic operators to define the mild solution of the system. Results obtained here improve and extend some known results.

We investigate the behavior of the solutions of the ray equation in isotropic media. We discuss local and global transversals and wavefronts and give examples of rays without wavefronts.

Shown is the derivation of Lorentz-Einstein k-factor in SRT as an amplitude-term of oscillation-differential equations of second order.This case is shown for classical Lorentz-factor as solution of an equation for undamped oscillation as well as the developed theorem as a second solution for advanced SRT of fourth order with an equation for damped oscillation-states.This advanced term allows a calculation for any velocities by real rest mass.Also accelerated coordinate -frames are discussed.

In the paper, by virtue of the Fa`a di Bruno formula, some properties of the Bell polynomials of the second kind, and the inversion formulas of binomial numbers and the Stirling numbers of the first and second kinds, the authors simplify meaningfully and significantly coefficients in two families of ordinary differential equations associated with higher order Frobenius–Euler numbers.

Physical processes are usually described using four-dimensional vector quantities - coordinate vector, momentum vector, current vector. But at the fundamental level they are characterized by spinors - coordinate spinors, impulse spinors, spinor wave functions. The propagation of fields and their interaction takes place at the spinor level, and since each spinor uniquely corresponds to a certain vector, the results of physical processes appear before us in vector form. For example, the relativistic Schrödinger equation and the Dirac equation are formulated by means of coordinate vectors, momentum vectors and quantum operators corresponding to them. In the Schrödinger equation the wave function is represented by a single complex quantity, in the Dirac equation a step forward is taken and the wave function is a spinor with complex components, but still coordinates and momentum are vectors. For a closed description of nature using only spinor quantities, it is necessary to have an equation similar to the Dirac equation in which momentum, coordinates and operators are spinors. It is such an equation that is presented in this paper. Using the example of the interaction between an electron and an electromagnetic field, we can see that the spinor equation contains more detailed information about the interaction than the vector equations. This is not new for quantum mechanics, since it describes interactions using complex wave functions, which cannot be observed directly, and only when measured goes to probabilities in the form of squares of the moduli of the wave functions. In the same way spinor quantities are not observable, but they completely determine observable vectors.In Section 2 of the paper, we analyze the quadratic form for an arbitrary four-component complex vector based on Pauli matrices. The form is invariant with respect to Lorentz transformations including any rotations and boosts. The invariance of the form allows us to construct on its basis an equation for a free particle combining the properties of the relativistic wave equation and the Dirac equation. For an electron in the presence of an electromagnetic potential it is shown that taking into account the commutation relations between the momentum and coordinate components allows us to obtain from this equation the known results describing the interactions of the electron spin with the electric and magnetic field. In section 3 of the paper this quadratic form is expressed through momentum spinors, which makes it possible to obtain an equation for the spinor wave function in spinor coordinate space by replacing the momentum spinor components by partial derivative operators on the corresponding coordinate spinor component. At the end of the paper the question on a possibility of second quantization of the electron field in the spinor coordinate space is touch upon. The statement that the electron and positron have the same, and exactly positive energy, and have opposite signs of charge and opposite signs of mass is also justified. Accordingly, in the process of annihilation their total energy, momentum, charge and mass do not change. It is shown that if we take as an axiom the conservation at interaction of the total mass of the system of particles taking into account its sign, then it is possible to explain the difference between bosons and fermions in the statistics to which they obey.

Calculus, as a starting point to everyone's study of mathematics, was developed along with classical mechanics. Though not rigorously formulated, the intuition for differential form and infinitesimal motion makes calculus not hard to understand. And hence it does not require great effort to grasp the principles of classical mechanics. However, the situation becomes quite different when it comes to the quantum world. Given a particle system, what of interest is always the physical quantity of this system, such as energy, position, momentum, entropy, etc. But it is no longer straightforward.

Molecular studies regarding regulatory elements such as small ncRNAs and their mechanisms are poorly understood in infectious diseases. Tuberculosis is one of the oldest infectious diseases of humanity, and it is still a challenge to prevent and treat it. The control of the infection as well as its diagnosis are still complex, and treatments used are linked to several side effects. This study aimed to investigate miRNA’s expression profile to identify possible biomarkers for tuberculosis. We applied NGS techniques to investigate miRNA’s global expression profile from blood samples of infected patients with tuberculosis, their respective healthy physicians, and external healthy individuals as controls. Samples from 22 individuals run through a differential expression, target genes, gene set enrichment, and miRNA-gene network analysis. We observed 153 altered miRNAs, among which, only three DEmiRNAs (hsa-let-7g-5p, hsa-miR-486-3p and hsa-miR-4732-5p) were found between the investigated patients and their respective physicians. These DEmiRNAs are suggested to play an important role in granuloma regulation and their immune physiopathology. Our results propose that miRNAs may be involved in immune modulation, regulating the repertoire of genes expressed in the immune system’s cells. Our findings encourage the application of miRNAs as potential biomarkers for tuberculosis.

As currently understood, the Dirac theory employs a 4 x1 type wavefunction. This 4x1 Dirac wavefunction is acted upon by a 4x4 Dirac Hamiltonian operator, in which process, four independent particle solutions result. Insofar as the real physical meaning and distinction of these four solutions, it is not clear what these solutions really mean. We demonstrate herein that these four independent particle solutions can be brought together under a single roof wherein the Dirac wavefunction takes a new form as a 4x4 wavefunction. In this new formation of the Dirac wavefunction, these four particle solutions precipitate into three distinct and mutuality dependent particles that are eternally bound in the same region of space. Given that Quarks are readily found in a mysterious threesome cohabitation-state eternally bound inside the Proton and Neutron, we make the suggestion that these Dirac particles might be Quarks. For the avoidance of speculation, we do not herein explore this idea further but merely present it as a very interesting idea worthy of further investigation. We however must say that, in the meantime, we are looking further into this very interesting idea, with the hope of making inroads in the immediate future.

We propose a new hierarchy of the vector derivative nonlinear Schr{\"o}dinger equations and consider the simplest multiphase solutions of this hierarchy. The study of the simplest solutions of these equations led to the following results. First, the three-leaf spectral curves $\Gamma=\{(\mu,\lambda)\}$ of the simplest multiphase solutions have a quite simple symmetry. They are invariant with respect to holomorphic involution $\tau$. The type of this involution depends on the genus of the spectral curve. Or the involution has the form $\tau:(\mu,\lambda)\to(\mu,-\lambda)$, or $\tau:(\mu,\lambda)\to(-\mu,-\lambda)$. The presence of symmetry leads to the fact that the dynamics of the solution is determined not by the entire spectral curve $\Gamma$, but by its factor $\Gamma/\tau$, which has a smaller genus. Secondly, it turned out that the dynamics of the two-component vector $\bp=(p_1,p_2)^t$ is determined, first of all, by the dynamics of its length $\abs{\bp}$. Independent equations determine the dependence of the direction of the vector $\bp$ from its length. In cases where the direction of the vector $\bp$ is fixed, the corresponding spectral curve splits into separate components. In conclusion, we note that, as in the case of the Manakov system, the equation of the spectral curve is invariant with respect to the orthogonal transformation of the vector solutions. I.e., the solution can be found from the spectral curve up to the orthogonal transformation. This fact indicates that the spectral curve does not depend on the individual components of the solution, but on their symmetric functions. Thus, the spectral data of multiphase solutions have two symmetries. These symmetries make it difficult to reconstruct signals from their spectral data. The work contains examples illustrating these statements.

In this study, a design of Morlet wavelet neural networks (MWNNs) is presented to solve the prediction differential model (PDM) using the global approximation capability of genetic algorithm (GA) and local quick interior-point algorithm scheme (IPAS), i.e., MWNN-GAIPAS. The famous PDM is known as a variant of functional differential system that works as an opposite of the historical delay differential models. A fitness function is optimized using the mean square error by applying the GA-IPAS for solving the PDM. Three PDM examples have been presented numerically to check the authenticity of the MWNN-GAIPAS. For the perfection of the designed MWNN-GAIPAS, the comparability of the obtained results and exact results is performed. Moreover, the neuron analysis is performed by taking the 3, 10 and 20 number of neurons. The statistical observations have been performed to authenticate the reliability of the MWNN-GAIPAS for solving the PDM.

We propose a notion of $(q,\sigma,\tau)$-differential graded algebra, which generalizes the notions of $(\sigma,\tau)$-differential graded algebra and $q$-differential graded algebra. We construct two examples of $(q,\sigma,\tau)$-differential graded algebra, where the first one is constructed by means of generalized Clifford algebra with two generators (reduced quantum plane), where we use a $(\sigma,\tau)$-twisted graded $q$-commutator. In order to construct the second example, we introduce a notion of $(\sigma,\tau)$-pre-cosimplicial algebra.

By expressing the Schrödinger wave function in the form $\psi =R{e}^{iS/\hslash }$ , where R and S are real functions, we have shown that the expectation value of S is conserved. The amplitude of the wave (R) is found to satisfy the Schrödinger equation while the phase (S) is related to the energy conservation. Besides the quantum potential that depends on R, viz., $V_Q=-\frac{{\hslash }^2}{2m}\frac{{\nabla }^{2}R}{R}$ , we have obtained a spin potential $V_S=-\frac{S{\nabla }^{2}S}{m}$ that depends on S which is attributed to the particle spin. The spin force is found to give rise to dissipative viscous force. The quantum potential may be attributed to the interaction between the two subfields S and R comprising the quantum particle. This results in splitting (creation/annihilation) of these subfields, each having a mass $m{c}^2$ with an internal frequency of $2m{c}^2/\hslash$ , satisfying the original wave equation and endowing the particle its quantum nature. The mass of one subfield reflects the interaction with the other subfield. If in Bohmian ansatz R satisfies the Klein-Gordon equation, then S must satisfies the wave equation. Conversely, if R satisfies the wave equation, then S yields the Einstein relativistic energy momentum equation.

The present work is concerned with the study of the generalized Langevin equation and its link to the physical theories of statistical mechanics and scale relativity. It is demonstrated that the form of the coefficients of Langevin equation depend critically on the assumption of continuity of the reconstructed trajectory. This in turn demands for the fluctuations of the diffusion term to be discontinuous in time. This paper further investigates the connection between the scale-relativistic and stochastic mechanics approaches, respectively, with the study of the Burgers equation, which in this case appears as a stochastic geodesic equation. By further demanding time reversibility of the drift the Langevin equation can also describe equivalent quantum-mechanical systems in a path-wise manner. The resulting statistical description obeys the Fokker-Plank equation of the probability density of the differential system, which can be readily estimated from Monte Carlo simulations of the random paths. Based on the Fokker-Plank formalism a new derivation of the transient probability densities is presented. Finally, stochastic simulations are compared to the theoretical results.

Antiangiogenic therapy (AAT) is an indirect oncological modality which is aimed at disruption of cancer cells nutrient supply. Invasive tumors have been shown to possess inherent resistant to this treatment, while compactly growing benign tumors react to it by shrinkage. It is generally accepted that AAT by itself is not curative. This study presents a mathematical model of non-invasive tumor growth with physiologically justified account of alteration of microvasculature and biomechanical aspects during tumor growth and AAT. In untreated setting the model reproduces tumor growth with saturation, where the maximum tumor volume depends on the level of angiogenesis. The outcomes of simulations of AAT depend on the tumor size at the moment of treatment initiation. If it is close to the stable size of avascular tumor grown in absence of angiogenesis, then the tumor is rapidly stabilized by AAT. Treatment of large tumors is accompanied by displacement of normal tissue due to the tumor shrinkage. During it, microvasculature undergoes distortion which degree depends on the displacement distance. As it affects tumor nutrient supply, the stable size of a tumor that undergoes AAT negatively correlates with its size at the beginning of treatment. For sufficiently large initial tumors, the long-term survival of tumor cells is compromised by competition with normal cells for severely limited inflow of nutrients, which makes AAT effectively curative.

Background: By using outcome prediction scores, it is possible to distinguish between good and poor performers with cochlear implants (CI) after CI implantation. The reasons for poor performance, despite good basic conditions, can be manifold. On the one hand, the postoperative fitting may be inadequate; on the other, neurophysiological disease processes may impair speech understanding with a CI. These disease processes are not yet fully understood. In acoustics, it is known that the auditory brainstem responses (ABR) and their latencies and amplitudes allow differential diagnosis based on reference values for normal-hearing individuals. The aim of this study was to provide reference values for electrically evoked brainstem responses (EABRs) in terms of rate-dependent latencies and amplitudes. Methods: 20 experienced adult CI recipients with a predicted and measured good postoperative word recognition score were recruited from the clinic’s patient pool. In the same stimulation mode and intensity we measured latencies and interpeak-latencies of EABRs and electrically evoked compound action potentials (ECAPs). With a defined supra-threshold stimulation intensity above the individual ECAP threshold, we applied stimulation at several rates between 11 and 91 stimuli per second. Results: We found rate dependences for EABR latency t3 and t5 in the order of 0.19 ms and 0.37 ms respectively, while ECAP was not affected by rate. Correspondingly, the interpeak intervalls’ rate dependences for t5-t1, t5-t3 and t3-t1 were of the order of 0.37 ms, 0.18 ms and 0.19 ms. Comparing the EABR amplitudes between the stimulation rates 11/s and 81/s, we found that at 81/s the amplitudes were significantly reduced down: to 73% for A3 and 81% for A5. These rate dependences of latency and amplitude in EABR have characteristics comparable to those of acoustic ABR. Conclusions: These data may serve to provide reference values for EABR and ECAP latencies, interpeak intervals and amplitudes with respect to stimulation rate. Altered response patterns of ECAPs and EABRs to normalised stimulation modes could be used in the future to describe and classify neuropathological processes in a better-differentiated way.

The core reasoning task for datalog engines is materialization, the evaluation of a datalog program over a database alongside its phys- ical incorporation into the database itself. The de-facto method of computing it, is through the recursive application of inference rules. Due to it being a costly operation, it is a must for datalog engines to provide incremental materialization, that is, to adjust the computation to new data, instead of restarting from scratch. One of the major caveats, is that deleting data is notoriously more involved than adding, since one has to take into account all possible data that has been entailed from what is being deleted. Differential Dataflow is a computational model that provides efficient incremen- tal maintenance, notoriously with equal performance between ad- ditions and deletions, and work distribution, of iterative dataflows. In this paper we investigate the performance of materialization with three reference datalog implementations, out of which one is built on top of a lightweight relational engine, and the two others are differential-dataflow and non-differential versions of the same rewrite algorithm, with the same optimizations.

This study aims to implement a SSDBI (single-stage differential boost inverter, SSDBI) applying in a single-stage BESS (battery energy storage system, BESS) topology, which can supply power from lower voltage battery module to AC load. Compared with the common two-stage topology, which has a two-stage converter and higher voltage battery module array, the single-stage topology can reduce the number of cells and components, and improve the power density. In addition, a modified SPWM control was proposed to reduce the control complexity of the SSDBI while improving the THD of the inverter. The modified SPWM control can reduce the duty ratio of the SSDBI and the stress on the components to improve the AC voltage output waveform and reduce the THD value about 2.654%.

Image segmentation and classification still represent an active area of research since no universal solution can be identified. Established segmentation algorithms like thresholding are problem specific, treat well the easy cases and mostly relied on single parameter i.e intensity. Machine learning approaches offer alternatives where predefined features are combined into different classifiers. On the other hand, the outcome of machine learning is only as good as the underlying feature space. Differential geometry can substantially improve the outcome of machine learning since it can enrich the underlying feature space with new geometrical objects, called invariants. In this way, the geometrical features form a high-dimensional feature space for each pixel, where original objects can be resolved. Alternatives based on the geometry of the image scale-invariant interest points have been exploited successfully in the field of computer vision. Here, we integrate geometrical feature extraction based on signal processing, machine learning, and input relying on domain knowledge. The approach is exemplified on the ISBI 2012 image segmentation challenge data set. As a second application, we demonstrate powerful image classification functionality based on the same principles, which was applied to the HeLa and HEp-2 data sets. Obtained results demonstrate that feature space enrichment properly balanced with feature selection functionality can achieve performance comparable to deep learning approaches.

Background: Traditional omic analysis relies on p-value and fold change as selection criteria. There is an ongoing debate on their effectiveness in delivering systemic and robust interpretation, due to their dependence on assumptions of conformity with various parametric distributions.Here, we propose a threshold-free selection method based on robust, non-parametric statistics, ensuring independence from the statistical distribution properties and broad applicability. Such methods could adapt to different initial data distributions, contrary to statistical techniques based on fixed thresholds. Methods: Our work extends the Rank Products methodology with a neutral selection method of high information-extraction capacity. We introduce the calculation of the RP distribution’s entropy to isolate the features of interest by their contribution to the distribution’s information content. The aim is a methodology performing threshold-free identification of the differentially expressed features, which are highly informative about the phenomenon under scrutiny. Conclusions: Applying the proposed method on microarray (transcriptomic and DNA methylation) and RNAseq count data of varying sizes and noise presence, we observe robust convergence for the different parameterisations to stable cutoff points. Functional analysis through BioInfoMiner and EnrichR was used to evaluate the information potency of the resulting feature lists. Overall, the derived functional terms provide a systemic description highly compatible with the results of traditional statistical hypothesis testing techniques. The methodology behaves consistently across different data types. The feature lists are compact and information-rich, indicating phenotypic aspects specific to the tissue and biological phenomenon i nvestigated. Selection by information content measures efficiently addresses problems, emerging from arbitrary thresholding, thus facilitating the full automation of the analysis.

An easier method for the calibration of differential drive robots is presented. Most of the calibration is done on-board and it is not necessary to expend too much time taking note of the robot's position. The calibration method does not need a big free space to perform the tests. The bigger space is just in a straight line, which is easy to find. Results with the proposed method are compared with those from UMB as a reference, and they show very little deviation while the proposed calibration is much simpler.

We extended the usage of the expansion method based on Sine-Gordon equation to the two dimensional Fisher equation. The relation between the trigonometric and hyperbolic functions are derived from the Sine-Gordon equation dened in two space dimension. The complex-valued traveling wave solutions to the two dimensional Fisher and Nagumo equations are set in forms a nite series of multiplications of powers of sech(.) and tanh(.) functions.