Beyond the hypothesis of boundedness for the random coefficient of Airy, Hermite and Laguerre differential equations with uncertainties

[EN] In this work, we study the full randomized versions of Airy, Hermite and Laguerre differential equations, which depend on a random variable appearing as an equation coefficient as well as two random initial conditions. In previous contributions, the mean square stochastic solutions to the aforementioned random differential equations were constructed via the Frobenius method, under the assumption of exponential growth of the absolute moments of the equation coefficient, which is equivalent to its essential boundedness. In this paper we aim at relaxing the boundedness hypothesis to allow more general probability distributions for the equation coefficient. We prove that the equations are solvable in the mean square sense when the equation coefficient has finite moment-generating function in a neighborhood of the origin. A thorough discussion of the new hypotheses is included.This work has been supported by the Spanish Ministerio de Economia y Competitividad grant MTM2017-89664-P.Calatayud Gregori, J.; Cortés, J.; Jornet Sanz, M. (2020). Beyond the hypothesis of boundedness for the random coefficient of Airy, Hermite and Laguerre differential equations with uncertainties. Stochastic Analysis and Applications. 38(5):875-885. https://doi.org/10.1080/07362994.2020.1733017S875885385Neckel, T., & Rupp, F. (2013). Random Differential Equations in Scientific Computing. doi:10.2478/9788376560267Villafuerte, L., Braumann, C. A., Cortés, J.-C., & Jódar, L. (2010). Random differential operational calculus: Theory and applications. Computers & Mathematics with Applications, 59(1), 115-125. doi:10.1016/j.camwa.2009.08.061Cortés, J.-C., Jódar, L., Camacho, F., & Villafuerte, L. (2010). Random Airy type differential equations: Mean square exact and numerical solutions. Computers & Mathematics with Applications, 60(5), 1237-1244. doi:10.1016/j.camwa.2010.05.046Calbo, G., Cortés, J.-C., & Jódar, L. (2011). Random Hermite differential equations: Mean square power series solutions and statistical properties. Applied Mathematics and Computation, 218(7), 3654-3666. doi:10.1016/j.amc.2011.09.008Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation. Mediterranean Journal of Mathematics, 16(3). doi:10.1007/s00009-019-1338-6Calatayud, J., Cortés, J.-C., Jornet, M., & Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations, 2018(1). doi:10.1186/s13662-018-1848-8Gregori, J., López, J., & Sanz, M. (2018). Some Notes to Extend the Study on Random Non-Autonomous Second Order Linear Differential Equations Appearing in Mathematical Modeling. Mathematical and Computational Applications, 23(4), 76. doi:10.3390/mca23040076Calbo, G., Cortés, J.-C., & Jódar, L. (2010). Mean square power series solution of random linear differential equations. Computers & Mathematics with Applications, 59(1), 559-572. doi:10.1016/j.camwa.2009.06.007Calbo, G., Cortés, J.-C., Jódar, L., & Villafuerte, L. (2010). Analytic stochastic process solutions of second-order random differential equations. Applied Mathematics Letters, 23(12), 1421-1424. doi:10.1016/j.aml.2010.07.011CALBO SANJUÁN, G. (s. f.). Mean Square Analytic Solutions of Random Linear Models. doi:10.4995/thesis/10251/8721Jagadeesan, M. (2017). Simple analysis of sparse, sign-consistent JL. arXiv:1708.02966.Lin, G. D. (2017). Recent developments on the moment problem. Journal of Statistical Distributions and Applications, 4(1). doi:10.1186/s40488-017-0059-2Ernst, O. G., Mugler, A., Starkloff, H.-J., & Ullmann, E. (2011). On the convergence of generalized polynomial chaos expansions. ESAIM: Mathematical Modelling and Numerical Analysis, 46(2), 317-339. doi:10.1051/m2an/2011045Calbo, G., Cortés, J.-C., Jódar, L., & Villafuerte, L. (2011). Solving the random Legendre differential equation: Mean square power series solution and its statistical functions. Computers & Mathematics with Applications, 61(9), 2782-2792. doi:10.1016/j.camwa.2011.03.04

A mean square chain rule and its application in solving the random Chebyshev differential equation

[EN] In this paper a new version of the chain rule for calculat- ing the mean square derivative of a second-order stochastic process is proven. This random operational calculus rule is applied to construct a rigorous mean square solution of the random Chebyshev differential equation (r.C.d.e.) assuming mild moment hypotheses on the random variables that appear as coefficients and initial conditions of the cor- responding initial value problem. Such solution is represented through a mean square random power series. Moreover, reliable approximations for the mean and standard deviation functions to the solution stochastic process of the r.C.d.e. are given. Several examples, that illustrate the theoretical results, are included.This work was completed with the support of our TEX-pert.Cortés, J.; Villafuerte, L.; Burgos-Simon, C. (2017). A mean square chain rule and its application in solving the random Chebyshev differential equation. Mediterranean Journal of Mathematics. 14(1):14-35. https://doi.org/10.1007/s00009-017-0853-6S1435141Calbo, G., Cortés, J.C., Jódar, L., Villafuerte, L.: Analytic stochastic process solutions of second-order random differential equations. Appl. Math. Lett. 23(12), 1421–1424 (2010). doi: 10.1016/j.aml.2010.07.011El-Tawil, M.A., El-Sohaly, M.: Mean square numerical methods for initial value random differential equations. Open J. Discret. Math. 1(1), 164–171 (2011). doi: 10.4236/ojdm.2011.12009Khodabin, M., Maleknejad, K., Rostami, K., Nouri, M.: Numerical solution of stochastic differential equations by second order Runge Kutta methods. Math. Comp. Model. 59(9–10), 1910–1920 (2010). doi: 10.1016/j.mcm.2011.01.018Santos, L.T., Dorini, F.A., Cunha, M.C.C.: The probability density function to the random linear transport equation. Appl. Math. Comput. 216(5), 1524–1530 (2010). doi: 10.1016/j.amc.2010.03.001González Parra, G., Chen-Charpentier, B.M., Arenas, A.J.: Polynomial Chaos for random fractional order differential equations. Appl. Math. Comput. 226(1), 123–130 (2014). doi: 10.1016/j.amc.2013.10.51El-Beltagy, M.A., El-Tawil, M.A.: Toward a solution of a class of non-linear stochastic perturbed PDEs using automated WHEP algorithm. Appl. Math. Model. 37(12–13), 7174–7192 (2013). doi: 10.1016/j.apm.2013.01.038Nouri, K., Ranjbar, H.: Mean square convergence of the numerical solution of random differential equations. Mediterran. J. Math. 12(3), 1123–1140 (2015). doi: 10.1007/s00009-014-0452-8Villafuerte, L., Braumann, C.A., Cortés, J.C., Jódar, L.: Random differential operational calculus: theory and applications. Comp. Math. Appl. 59(1), 115–125 (2010). doi: 10.1016/j.camwa.2009.08.061Øksendal, B.: Stochastic differential equations: an introduction with applications, 6th edn. Springer, Berlin (2007)Soong, T.T.: Random differential equations in science and engineering. Academic Press, New York (1973)Wong, B., Hajek, B.: Stochastic processes in engineering systems. Springer Verlag, New York (1985)Arnold, L.: Stochastic differential equations. Theory and applications. John Wiley, New York (1974)Cortés, J.C., Jódar, L., Camacho, J., Villafuerte, L.: Random Airy type differential equations: mean square exact and numerical solutions. Comput. Math. Appl. 60(5), 1237–1244 (2010). doi: 10.1016/j.camwa.2010.05.046Calbo, G., Cortés, J.C., Jódar, L.: Random Hermite differential equations: mean square power series solutions and statistical properties. Appl. Math. Comp. 218(7), 3654–3666 (2011). doi: 10.1016/j.amc.2011.09.008Calbo, G., Cortés, J.C., Jódar, L., Villafuerte, L.: Solving the random Legendre differential equation: Mean square power series solution and its statistical functions. Comp. Math. Appl. 61(9), 2782–2792 (2010). doi: 10.1016/j.camwa.2011.03.045Cortés, J.C., Jódar, L., Company, R., Villafuerte, L.: Laguerre random polynomials: definition, differential and statistical properties. Utilit. Math. 98, 283–293 (2015)Cortés, J.C., Jódar, L., Villafuerte, L.: Mean square solution of Bessel differential equation with uncertainties. J. Comp. Appl. Math. 309, 383–395 (2017). doi: 10.1016/j.cam.2016.01.034Golmankhaneh, A.K., Porghoveh, N.A., Baleanu, D.: Mean square solutions of second-order random differential equations by using homotopy analysis method. Romanian Reports Physics 65(2), 1237–1244 (2013)Khalaf, S.L.: Mean square solutions of second-order random differential equations by using homotopy perturbation method. Int. Math. Forum 6(48), 2361–2370 (2011)Khudair, A.R., Ameen, A.A., Khalaf, S.L.: Mean square solutions of second-order random differential equations by using Adomian decomposition method. Appl. Math. Sci. 5(49), 2521–2535 (2011)Agarwal, R.P., O’Regan, D.: Ordinary and partial differential equations. Springer, New York (2009

Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties

[EN] In this paper we study random non-autonomous second order linear differential equations by taking advantage of the powerful theory of random difference equations. The coefficients are assumed to be stochastic processes, and the initial conditions are random variables both defined in a common underlying complete probability space. Under appropriate assumptions established on the data stochastic processes and on the random initial conditions, and using key results on difference equations, we prove the existence of an analytic stochastic process solution in the random mean square sense. Truncating the random series that defines the solution process, we are able to approximate the main statistical properties of the solution, such as the expectation and the variance. We also obtain error a priori bounds to construct reliable approximations of both statistical moments. We include a set of numerical examples to illustrate the main theoretical results established throughout the paper. We finish with an example where our findings are combined with Monte Carlo simulations to model uncertainty using real data.This work has been supported by the Spanish Ministerio de Economia y Competitividad grant MTM2017-89664-P. Marc Jornet acknowledges the doctorate scholarship granted by Programa de Ayudas de Investigacion y Desarrollo (PAID), Universitat Politecnica de Valencia.Calatayud-Gregori, J.; Cortés, J.; Jornet-Sanz, M.; Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations. (3):1-29. https://doi.org/10.1186/s13662-018-1848-8S1293Apostol, T.M.: Mathematical Analysis, 2nd edn. Pearson, New York (1976)Boyce, W.E.: Probabilistic Methods in Applied Mathematics I. Academic Press, New York (1968)Calbo, G., Cortés, J.C., Jódar, L.: Random Hermite differential equations: mean square power series solutions and statistical properties. Appl. Math. Comput. 218(7), 3654–3666 (2011)Calbo, G., Cortés, J.C., Jódar, L., Villafuerte, L.: Solving the random Legendre differential equation: mean square power series solution and its statistical functions. Comput. Math. Appl. 61(9), 2782–2792 (2011)Casabán, M.C., Cortés, J.C., Navarro-Quiles, A., Romero, J.V., Roselló, M.D., Villanueva, R.J.: Computing probabilistic solutions of the Bernoulli random differential equation. J. Comput. Appl. Math. 309, 396–407 (2017)Casabán, M.C., Cortés, J.C., Romero, J.V., Roselló, M.D.: Solving random homogeneous linear second-order differential equations: a full probabilistic description. Mediterr. J. Math. 13(6), 3817–3836 (2016)Cortés, J.C., Jódar, L., Camacho, J., Villafuerte, L.: Random Airy type differential equations: mean square exact and numerical solutions. Comput. Math. Appl. 60(5), 1237–1244 (2010)Cortés, J.C., Jódar, L., Company, R., Villafuerte, L.: Laguerre random polynomials: definition, differential and statistical properties. Util. Math. 98, 283–295 (2015)Cortés, J.C., Jódar, L., Villafuerte, L.: Random linear-quadratic mathematical models: computing explicit solutions and applications. Math. Comput. Simul. 79(7), 2076–2090 (2009)Cortés, J.C., Jódar, L., Villafuerte, L.: Mean square solution of Bessel differential equation with uncertainties. J. Comput. Appl. Math. 309(1), 383–395 (2017)Cortés, J.C., Sevilla-Peris, P., Jódar, L.: Analytic-numerical approximating processes of diffusion equation with data uncertainty. Comput. Math. Appl. 49(7–8), 1255–1266 (2005)Díaz-Infante, S., Jerez, S.: Convergence and asymptotic stability of the explicit Steklov method for stochastic differential equations. J. Comput. Appl. Math. 291(1), 36–47 (2016)Dorini, F., Cunha, M.: Statistical moments of the random linear transport equation. J. Comput. Phys. 227(19), 8541–8550 (2008)Dorini, F.A., Cecconello, M.S., Dorini, M.B.: On the logistic equation subject to uncertainties in the environmental carrying capacity and initial population density. Commun. Nonlinear Sci. Numer. Simul. 33, 160–173 (2016)Golmankhaneh, A.K., Porghoveh, N.A., Baleanu, D.: Mean square solutions of second-order random differential equations by using homotopy analysis method. Rom. Rep. Phys. 65(2), 350–362 (2013)Grimmett, G.R., Stirzaker, D.R.: Probability and Random Processes. Clarendon Press, Oxford (2000)Henderson, D., Plaschko, P.: Stochastic Differential Equations in Science and Engineering. Cambridge Texts in Applied Mathematics. World Scientific, Singapore (2006)Hussein, A., Selim, M.M.: A developed solution of the stochastic Milne problem using probabilistic transformations. Appl. Math. Comput. 216(10), 2910–2919 (2009)Hussein, A., Selim, M.M.: Solution of the stochastic transport equation of neutral particles with anisotropic scattering using RVT technique. Appl. Math. Comput. 213(1), 250–261 (2009)Hussein, A., Selim, M.M.: Solution of the stochastic radiative transfer equation with Rayleigh scattering using RVT technique. Appl. Math. Comput. 218(13), 7193–7203 (2012)Khodabin, M., Maleknejad, K., Rostami, M., Nouri, M.: Numerical solution of stochastic differential equations by second order Runge–Kutta methods. Math. Comput. Model. 53(9–10), 1910–1920 (2011)Khodabin, M., Rostami, M.: Mean square numerical solution of stochastic differential equations by fourth order Runge–Kutta method and its application in the electric circuits with noise. Adv. Differ. Equ. 2015, 62 (2015)Khudair, A.K., Ameen, A.A., Khalaf, S.L.: Mean square solutions of second-order random differential equations by using Adomian decomposition method. Appl. Math. Sci. 51(5), 2521–2535 (2011)Khudair, A.K., Haddad, S.A.M., Khalaf, S.L.: Mean square solutions of second-order random differential equations by using the differential transformation method. Open J. Appl. Sci. 6, 287–297 (2016)Lesaffre, E., Lawson, A.B.: Bayesian Biostatistics. Statistics in Practice. Wiley, New York (2012)Li, X., Fu, X.: Stability analysis of stochastic functional differential equations with infinite delay and its application to recurrent neural networks. J. Comput. Appl. Math. 234(2), 407–417 (2010)Licea, J.A., Villafuerte, L., Chen-Charpentier, B.M.: Analytic and numerical solutions of a Riccati differential equation with random coefficients. J. Comput. Appl. Math. 309(1), 208–219 (2013)Liu, S., Debbouche, A., Wang, J.: On the iterative learning control for stochastic impulsive differential equations with randomly varying trial lengths. J. Comput. Appl. Math. 312, 47–57 (2017)Loève, M.: Probability Theory. Vol. I. Springer, Mineola (1977)Lord, G.J., Powell, C.E., Shardlow, T.: An Introduction to Computational Stochastic PDEs. Cambridge Texts in Applied Mathematics. Dover, New York (2014)Nouri, K., Ranjbar, H.: Mean square convergence of the numerical solution of random differential equations. Mediterr. J. Math. 12(3), 1123–1140 (2015)Rencher, A.C., Schaalje, G.B.: Linear Models in Statistics, 2nd edn. Wiley, New York (2008)Santos, L.T., Dorini, F.A., Cunha, M.C.C.: The probability density function to the random linear transport equation. Appl. Math. Comput. 216(5), 1524–1530 (2010)Seber, G.A.F., Wild, C.J.: Nonlinear Regression. Cambridge Texts in Applied Mathematics. Wiley, New York (2003)Smith, R.C.: Uncertainty Quantification: Theory, Implementation, and Applications. SIAM, Philadelphia (2014)Soheili, A.R., Toutounian, F., Soleymani, F.: A fast convergent numerical method for matrix sign function with application in SDEs (Stochastic Differential Equations). J. Comput. Appl. Math. 282, 167–178 (2015)Soong, T.T.: Random Differential Equations in Science and Engineering. Academic Press, New York (1973)Villafuerte, L., Braumann, C.A., Cortés, J.C., Jódar, L.: Random differential operational calculus: theory and applications. Comput. Math. Appl. 59(1), 115–125 (2010)Xiu, D.: Numerical Methods for Stochastic Computations. A Spectral Method Approach. Cambridge Texts in Applied Mathematics. Princeton University Press, New York (2010

Improving the approximation of the first and second order statistics of the response stochastic process to the random Legendre differential equation

[EN] In this paper, we deal with uncertainty quantification for the random Legendre differential equation, with input coefficient A and initial conditions X-0 and X-1. In a previous study (Calbo et al. in Comput Math Appl 61(9):2782-2792, 2011), a mean square convergent power series solution on (-1/e, 1/e) was constructed, under the assumptions of mean fourth integrability of X-0 and X-1, independence, and at most exponential growth of the absolute moments of A. In this paper, we relax these conditions to construct an L-p solution (1 <= p <= infinity) to the random Legendre differential equation on the whole domain (-1, 1), as in its deterministic counterpart. Our hypotheses assume no independence and less integrability of X-0 and X-1. Moreover, the growth condition on the moments of A is characterized by the boundedness of A, which simplifies the proofs significantly. We also provide approximations of the expectation and variance of the response process. The numerical experiments show the wide applicability of our findings. A comparison with Monte Carlo simulations and gPC expansions is performed.This work has been supported by the Spanish Ministerio de Economia y Competitividad Grant MTM2017-89664-P. Marc Jornet acknowledges the doctorate scholarship granted by Programa de Ayudas de Investigacion y Desarrollo (PAID), Universitat Politecnica de Valencia.Calatayud-Gregori, J.; Cortés, J.; Jornet-Sanz, M. (2019). Improving the approximation of the first and second order statistics of the response stochastic process to the random Legendre differential equation. Mediterranean Journal of Mathematics. 16(3):1-14. https://doi.org/10.1007/s00009-019-1338-6S114163Soong, T.T.: Random Differential Equations in Science and Engineering. Academic Press, New York (1973)Strand, J.L.: Random ordinary differential equations. J. Differ. Equ. 7(3), 538–553 (1970)Smith, R.C.: Uncertainty quantification. Theory, implementation, and application. SIAM Comput. Sci. Eng. New York (2013) ISBN 9781611973211Fishman, G.: Monte Carlo: Concepts, Algorithms, and Applications. Springer, Berlin (2013)Cortés, J.-C., Romero, J.-V., Roselló, M.-D., Santonja, F.-J., Villanueva, R.-J.: Solving continuous models with dependent uncertainty: a computational approach. Abstr. Appl. Anal. 2013, 983839 (2013). https://doi.org/10.1155/2013/983839Xiu, D.: Numerical Methods for Stochastic Computations. A Spectral Method Approach. Cambridge Texts in Applied Mathematics. Princeton University Press, New York (2010)El-Tawil, M.A.: The approximate solutions of some stochastic differential equations using transformations. Appl. Math. Comput. 164(1), 167–178 (2005)Cortés, J.-C., Sevilla-Peris, P., Jódar, L.: Constructing approximate diffusion processes with uncertain data. Math. Comput. Simul. 73(1–4), 125–132 (2006)Cortés, J.-C., Jódar, L., Villafuerte, L., Villanueva, R.-J.: Computing mean square approximations of random diffusion models with source term. Math. Comput. Simul. 76(1–3), 44–48 (2007)Khodabin, M., Maleknejad, K., Rostami, M., Nouri, M.: Numerical solution of stochastic differential equations by second order Runge–Kutta methods. Math. Comput. Model. 53(9–10), 1910–1920 (2011)Nouri, K., Ranjbar, H.: Mean square convergence of the numerical solution of random differential equations. Mediterran. J. Math. 12(3), 1123–1140 (2015)Nouri, N.: Study on stochastic differential equations via modified Adomian decomposition method. U.P.B. Sci. Bull. Ser. A 78(1), 81–90 (2016)Khodabin, M., Rostami, M.: Mean square numerical solution of stochastic differential equations by fourth order Runge–Kutta method and its application in the electric circuits with noise. Adv. Differ. Equ. 623, 1–19 (2015)Díaz-Infante, S., Jerez, S.: Convergence and asymptotic stability of the explicit Steklov method for stochastic differential equations. J. Comput. Appl. Math. 291(1), 36–47 (2016)Soheili, Ali R, Toutounian, F., Soleymani, F.: A fast convergent numerical method for matrix sign function with application in SDEs (Stochastic Differential Equations). J. Comput. Appl. Math. 282, 167–178 (2015)Øksendal, B.: Stochastic Differential Equations. Springer, Berlin (2003)Villafuerte, L., Braumann, C.A., Cortés, J.-C., Jódar, L.: Random differential operational calculus: theory and applications. Comput. Math. Appl. 59(1), 115–125 (2010)Licea, J., Villafuerte, L., Chen-Charpentier, B.M.: Analytic and numerical solutions of a Riccati differential equation with random coefficients. J. Comput. Appl. Math. 309(1), 208–219 (2013)Cortés, J.-C., Jódar, L., Camacho, J., Villafuerte, L.: Random Airy type differential equations: mean square exact and numerical solutions. Comput. Math. Appl. 60(5), 1237–1244 (2010)Calbo, G., Cortés, J.-C., Jódar, L.: Random Hermite differential equations: mean square power series solutions and statistical properties. Appl. Math. Comput. 218(7), 3654–3666 (2011)Calbo, G., Cortés, J.-C., Jódar, L., Villafuerte, L.: Solving the random Legendre differential equation: mean square power series solution and its statistical functions. Comput. Math. Appl. 61(9), 2782–2792 (2011)Cortés, J.C., Jódar, L., Villafuerte, L.: Mean square solution of Bessel differential equation with uncertainties. J. Comput. Appl. Math. 309(1), 383–395 (2017)Golmankhaneh, A.K., Porghoveh, N.A., Baleanu, D.: Mean square solutions of second-order random differential equations by using homotopy analysis method. Roman. Rep. Phys. 65(2), 350–362 (2013)Khudair, A.K., Ameen, A.A., Khalaf, S.L.: Mean square solutions of second-order random differential equations by using Adomian decomposition method. Appl. Math. Sci. 51(5), 2521–2535 (2011)Khudair, A.K., Haddad, S.A.M., Khalaf, S.L.: Mean square solutions of second-order random differential equations by using the differential transformation method. Open J. Appl. Sci. 6, 287–297 (2016)Norman, L., Kotz, S., Balakrishnan, N.: Continuous Univariate Distributions, vol. 1. Wiley, Oxford (1994)Ernst, O.G., Mugler, A., Starkloff, H.-J., Ullmann, E.: On the convergence of generalized polynomial chaos expansions. ESAIM Math. Modell. Num. Anal. 46(2), 317–339 (2012)Shi, W., Zhang, C.: Error analysis of generalized polynomial chaos for nonlinear random ordinary differential equations. Appl. Num. Math. 62(12), 1954–1964 (2012)Calatayud, J., Cortés, J.-C., Jornet, M.: On the convergence of adaptive gPC for non-linear random difference equations: theoretical analysis and some practical recommendations. J. Nonlinear Sci. Appl. 11(9), 1077–1084 (2018

Mean square convergent non-standard numerical schemes for linear random differential equations with delay

[EN] In this paper, we are concerned with the construction of numerical schemes for linear random differential equations with discrete delay. For the linear deterministic differential equation with discrete delay, a recent contribution proposed a family of non-standard finite difference (NSFD) methods from an exact numerical scheme on the whole domain. The family of NSFD schemes had increasing order of accuracy, was dynamically consistent, and possessed simple computational properties compared to the exact scheme. In the random setting, when the two equation coefficients are bounded random variables and the initial condition is a regular stochastic process, we prove that the randomized NSFD schemes converge in the mean square (m.s.) sense. M.s. convergence allows for approximating the expectation and the variance of the solution stochastic process. In practice, the NSFD scheme is applied with symbolic inputs, and afterward the statistics are explicitly computed by using the linearity of the expectation. This procedure permits retaining the increasing order of accuracy of the deterministic counterpart. Some numerical examples illustrate the approach. The theoretical m.s. convergence rate is supported numerically, even when the two equation coefficients are unbounded random variables. M.s. dynamic consistency is assessed numerically. A comparison with Euler's method is performed. Finally, an example dealing with the time evolution of a photosynthetic bacterial population is presented.This work has been supported by the Spanish Ministerio de Economia, Industria y Competitividad (MINECO), the Agencia Estatal de Investigacion (AEI) and Fondo Europeo de Desarrollo Regional (FEDER UE) grant MTM2017-89664-P.Calatayud, J.; Cortés, J.; Jornet, M.; Rodríguez, F. (2020). Mean square convergent non-standard numerical schemes for linear random differential equations with delay. Mathematics. 8(9):1-17. https://doi.org/10.3390/math8091417S11789Bocharov, G. A., & Rihan, F. A. (2000). Numerical modelling in biosciences using delay differential equations. Journal of Computational and Applied Mathematics, 125(1-2), 183-199. doi:10.1016/s0377-0427(00)00468-4Jackson, M., & Chen-Charpentier, B. M. (2017). Modeling plant virus propagation with delays. Journal of Computational and Applied Mathematics, 309, 611-621. doi:10.1016/j.cam.2016.04.024Chen-Charpentier, B. M., & Diakite, I. (2016). A mathematical model of bone remodeling with delays. Journal of Computational and Applied Mathematics, 291, 76-84. doi:10.1016/j.cam.2014.11.025Kyrychko, Y. N., & Hogan, S. J. (2010). On the Use of Delay Equations in Engineering Applications. Journal of Vibration and Control, 16(7-8), 943-960. doi:10.1177/1077546309341100Harding, L., & Neamţu, M. (2016). A Dynamic Model of Unemployment with Migration and Delayed Policy Intervention. Computational Economics, 51(3), 427-462. doi:10.1007/s10614-016-9610-3Mickens, R. E. (2005). Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations. Journal of Difference Equations and Applications, 11(7), 645-653. doi:10.1080/10236190412331334527Patidar, K. C. (2016). Nonstandard finite difference methods: recent trends and further developments. Journal of Difference Equations and Applications, 22(6), 817-849. doi:10.1080/10236198.2016.1144748García, M. A., Castro, M. A., Martín, J. A., & Rodríguez, F. (2018). Exact and nonstandard numerical schemes for linear delay differential models. Applied Mathematics and Computation, 338, 337-345. doi:10.1016/j.amc.2018.06.029Castro, M. Á., García, M. A., Martín, J. A., & Rodríguez, F. (2019). Exact and Nonstandard Finite Difference Schemes for Coupled Linear Delay Differential Systems. Mathematics, 7(11), 1038. doi:10.3390/math7111038Villafuerte, L., Braumann, C. A., Cortés, J.-C., & Jódar, L. (2010). Random differential operational calculus: Theory and applications. Computers & Mathematics with Applications, 59(1), 115-125. doi:10.1016/j.camwa.2009.08.061Cortés, J.-C., Jódar, L., Roselló, M.-D., & Villafuerte, L. (2012). Solving initial and two-point boundary value linear random differential equations: A mean square approach. Applied Mathematics and Computation, 219(4), 2204-2211. doi:10.1016/j.amc.2012.08.066Calatayud, J., Cortés, J.-C., Jornet, M., & Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations, 2018(1). doi:10.1186/s13662-018-1848-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation. Mediterranean Journal of Mathematics, 16(3). doi:10.1007/s00009-019-1338-6Licea, J. A., Villafuerte, L., & Chen-Charpentier, B. M. (2013). Analytic and numerical solutions of a Riccati differential equation with random coefficients. Journal of Computational and Applied Mathematics, 239, 208-219. doi:10.1016/j.cam.2012.09.040Burgos, C., Calatayud, J., Cortés, J.-C., & Villafuerte, L. (2018). Solving a class of random non-autonomous linear fractional differential equations by means of a generalized mean square convergent power series. Applied Mathematics Letters, 78, 95-104. doi:10.1016/j.aml.2017.11.009Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Random differential equations with discrete delay. Stochastic Analysis and Applications, 37(5), 699-707. doi:10.1080/07362994.2019.1608833Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). $$\mathrm {L}^p$$-calculus Approach to the Random Autonomous Linear Differential Equation with Discrete Delay. Mediterranean Journal of Mathematics, 16(4). doi:10.1007/s00009-019-1370-6Cortés, J. C., & Jornet, M. (2020). Lp-Solution to the Random Linear Delay Differential Equation with a Stochastic Forcing Term. Mathematics, 8(6), 1013. doi:10.3390/math8061013Caraballo, T., Cortés, J.-C., & Navarro-Quiles, A. (2019). Applying the random variable transformation method to solve a class of random linear differential equation with discrete delay. Applied Mathematics and Computation, 356, 198-218. doi:10.1016/j.amc.2019.03.048Cortés, J. C., Jódar, L., & Villafuerte, L. (2007). Numerical solution of random differential equations: A mean square approach. Mathematical and Computer Modelling, 45(7-8), 757-765. doi:10.1016/j.mcm.2006.07.017Cortés, J. C., Jódar, L., & Villafuerte, L. (2007). Mean square numerical solution of random differential equations: Facts and possibilities. Computers & Mathematics with Applications, 53(7), 1098-1106. doi:10.1016/j.camwa.2006.05.030El-Tawil, M. A. (2005). The approximate solutions of some stochastic differential equations using transformations. Applied Mathematics and Computation, 164(1), 167-178. doi:10.1016/j.amc.2004.04.062Calatayud, J., Cortés, J.-C., Díaz, J. A., & Jornet, M. (2019). Density function of random differential equations via finite difference schemes: a theoretical analysis of a random diffusion-reaction Poisson-type problem. Stochastics, 92(4), 627-641. doi:10.1080/17442508.2019.1645849Calatayud, J., Cortés, J.-C., Díaz, J. A., & Jornet, M. (2020). Constructing reliable approximations of the probability density function to the random heat PDE via a finite difference scheme. Applied Numerical Mathematics, 151, 413-424. doi:10.1016/j.apnum.2020.01.012Burgos, C., Cortés, J.-C., Villafuerte, L., & Villanueva, R.-J. (2020). Mean square convergent numerical solutions of random fractional differential equations: Approximations of moments and density. Journal of Computational and Applied Mathematics, 378, 112925. doi:10.1016/j.cam.2020.112925Strand, J. . (1970). Random ordinary differential equations. Journal of Differential Equations, 7(3), 538-553. doi:10.1016/0022-0396(70)90100-2Buckwar, E. (2000). Introduction to the numerical analysis of stochastic delay differential equations. Journal of Computational and Applied Mathematics, 125(1-2), 297-307. doi:10.1016/s0377-0427(00)00475-1Antonio Dorini, F., & Sampaio, R. (2012). Some Results on the Random Wear Coefficient of the Archard Model. Journal of Applied Mechanics, 79(5). doi:10.1115/1.400645

Lp-solution to the random linear delay differential equation with stochastic forcing term

[EN] This paper aims at extending a previous contribution dealing with the random autonomous-homogeneous linear differential equation with discrete delay tau > 0, by adding a random forcing term f(t) that varies with time: x'(t) = ax(t) + bx(t-tau) + f(t), t >= 0, with initial condition x(t) = g(t), -tau <= t <= 0. The coefficients a and b are assumed to be random variables, while the forcing term f(t) and the initial condition g(t) are stochastic processes on their respective time domains. The equation is regarded in the Lebesgue space L-p of random variables with finite p-th moment. The deterministic solution constructed with the method of steps and the method of variation of constants, which involves the delayed exponential function, is proved to be an L-p-solution, under certain assumptions on the random data. This proof requires the extension of the deterministic Leibniz's integral rule for differentiation to the random scenario. Finally, we also prove that, when the delay tau tends to 0, the random delay equation tends in L-p to a random equation with no delay. Numerical experiments illustrate how our methodology permits determining the main statistics of the solution process, thereby allowing for uncertainty quantification.This work has been supported by the Spanish Ministerio de Economia, Industria y Competitividad (MINECO), the Agencia Estatal de Investigacion (AEI) and Fondo Europeo de Desarrollo Regional (FEDER UE) grant MTM2017-89664-P.Cortés, J.; Jornet, M. (2020). Lp-solution to the random linear delay differential equation with stochastic forcing term. Mathematics. 8(6):1-16. https://doi.org/10.3390/math8061013S11686Xiu, D., & Karniadakis, G. E. (2004). Supersensitivity due to uncertain boundary conditions. International Journal for Numerical Methods in Engineering, 61(12), 2114-2138. doi:10.1002/nme.1152Casabán, M.-C., Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 32, 199-210. doi:10.1016/j.cnsns.2015.08.009Strand, J. . (1970). Random ordinary differential equations. Journal of Differential Equations, 7(3), 538-553. doi:10.1016/0022-0396(70)90100-2Villafuerte, L., Braumann, C. A., Cortés, J.-C., & Jódar, L. (2010). Random differential operational calculus: Theory and applications. Computers & Mathematics with Applications, 59(1), 115-125. doi:10.1016/j.camwa.2009.08.061Cortés, J.-C., Jódar, L., Roselló, M.-D., & Villafuerte, L. (2012). Solving initial and two-point boundary value linear random differential equations: A mean square approach. Applied Mathematics and Computation, 219(4), 2204-2211. doi:10.1016/j.amc.2012.08.066Calatayud, J., Cortés, J.-C., Jornet, M., & Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations, 2018(1). doi:10.1186/s13662-018-1848-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation. Mediterranean Journal of Mathematics, 16(3). doi:10.1007/s00009-019-1338-6Licea, J. A., Villafuerte, L., & Chen-Charpentier, B. M. (2013). Analytic and numerical solutions of a Riccati differential equation with random coefficients. Journal of Computational and Applied Mathematics, 239, 208-219. doi:10.1016/j.cam.2012.09.040Burgos, C., Calatayud, J., Cortés, J.-C., & Villafuerte, L. (2018). Solving a class of random non-autonomous linear fractional differential equations by means of a generalized mean square convergent power series. Applied Mathematics Letters, 78, 95-104. doi:10.1016/j.aml.2017.11.009Nouri, K., & Ranjbar, H. (2014). Mean Square Convergence of the Numerical Solution of Random Differential Equations. Mediterranean Journal of Mathematics, 12(3), 1123-1140. doi:10.1007/s00009-014-0452-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Random differential equations with discrete delay. Stochastic Analysis and Applications, 37(5), 699-707. doi:10.1080/07362994.2019.1608833Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). $$\mathrm {L}^p$$-calculus Approach to the Random Autonomous Linear Differential Equation with Discrete Delay. Mediterranean Journal of Mathematics, 16(4). doi:10.1007/s00009-019-1370-6Caraballo, T., Cortés, J.-C., & Navarro-Quiles, A. (2019). Applying the random variable transformation method to solve a class of random linear differential equation with discrete delay. Applied Mathematics and Computation, 356, 198-218. doi:10.1016/j.amc.2019.03.048Zhou, T. (2014). A Stochastic Collocation Method for Delay Differential Equations with Random Input. Advances in Applied Mathematics and Mechanics, 6(4), 403-418. doi:10.4208/aamm.2012.m38Shi, W., & Zhang, C. (2017). Generalized polynomial chaos for nonlinear random delay differential equations. Applied Numerical Mathematics, 115, 16-31. doi:10.1016/j.apnum.2016.12.004Khusainov, D. Y., Ivanov, A. F., & Kovarzh, I. V. (2009). Solution of one heat equation with delay. Nonlinear Oscillations, 12(2), 260-282. doi:10.1007/s11072-009-0075-3Shaikhet, L. (2016). Stability of equilibrium states of a nonlinear delay differential equation with stochastic perturbations. International Journal of Robust and Nonlinear Control, 27(6), 915-924. doi:10.1002/rnc.3605Benhadri, M., & Zeghdoudi, H. (2018). Mean square asymptotic stability in nonlinear stochastic neutral Volterra-Levin equations with Poisson jumps and variable delays. Functiones et Approximatio Commentarii Mathematici, 58(2). doi:10.7169/facm/1657Santonja, F.-J., & Shaikhet, L. (2012). Analysing Social Epidemics by Delayed Stochastic Models. Discrete Dynamics in Nature and Society, 2012, 1-13. doi:10.1155/2012/530472Liu, L., & Caraballo, T. (2018). Analysis of a Stochastic 2D-Navier–Stokes Model with Infinite Delay. Journal of Dynamics and Differential Equations, 31(4), 2249-2274. doi:10.1007/s10884-018-9703-xLupulescu, V., & Abbas, U. (2011). Fuzzy delay differential equations. Fuzzy Optimization and Decision Making, 11(1), 99-111. doi:10.1007/s10700-011-9112-7Krapivsky, P. L., Luck, J. M., & Mallick, K. (2011). On stochastic differential equations with random delay. Journal of Statistical Mechanics: Theory and Experiment, 2011(10), P10008. doi:10.1088/1742-5468/2011/10/p10008GARRIDO-ATIENZA, M. J., OGROWSKY, A., & SCHMALFUSS, B. (2011). RANDOM DIFFERENTIAL EQUATIONS WITH RANDOM DELAYS. Stochastics and Dynamics, 11(02n03), 369-388. doi:10.1142/s0219493711003358Cortés, J.-C., Villafuerte, L., & Burgos, C. (2017). A Mean Square Chain Rule and its Application in Solving the Random Chebyshev Differential Equation. Mediterranean Journal of Mathematics, 14(1). doi:10.1007/s00009-017-0853-6Cortés, J. C., Jódar, L., & Villafuerte, L. (2007). Numerical solution of random differential equations: A mean square approach. Mathematical and Computer Modelling, 45(7-8), 757-765. doi:10.1016/j.mcm.2006.07.017Braumann, C. A., Cortés, J.-C., Jódar, L., & Villafuerte, L. (2018). On the random gamma function: Theory and computing. Journal of Computational and Applied Mathematics, 335, 142-155. doi:10.1016/j.cam.2017.11.045Khusainov, D. Y., & Pokojovy, M. (2015). Solving the Linear 1D Thermoelasticity Equations with Pure Delay. International Journal of Mathematics and Mathematical Sciences, 2015, 1-11. doi:10.1155/2015/47926

Some notes to extend the study on random non-autonomous second order linear differential equations appearing in Mathematical Modeling

The objective of this paper is to complete certain issues from our recent contribution [J. Calatayud, J.-C. Cort\'es, M. Jornet, L. Villafuerte, Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties, Advances in Difference Equations, 2018:392, 1--29 (2018)]. We restate the main theorem therein that deals with the homogeneous case, so that the hypotheses are clearer and also easier to check in applications. Another novelty is that we tackle the non-homogeneous equation with a theorem of existence of mean square analytic solution and a numerical example. We also prove the uniqueness of mean square solution via an habitual Lipschitz condition that extends the classical Picard Theorem to mean square calculus. In this manner, the study on general random non-autonomous second order linear differential equations with analytic data processes is completely resolved. Finally, we relate our exposition based on random power series with polynomial chaos expansions and the random differential transform method, being the latter a reformulation of our random Fr\"obenius method.Comment: 15 pages, 0 figures, 2 table

Maximal LpL^p-regularity for stochastic evolution equations

We prove maximal $L^p$-regularity for the stochastic evolution equation \{{aligned} dU(t) + A U(t)\, dt& = F(t,U(t))\,dt + B(t,U(t))\,dW_H(t), \qquad t\in [0,T], U(0) & = u_0, {aligned}. under the assumption that $A$ is a sectorial operator with a bounded $H^\infty$-calculus of angle less than $\frac12\pi$ on a space $L^q(\mathcal{O},\mu)$. The driving process $W_H$ is a cylindrical Brownian motion in an abstract Hilbert space $H$. For $p\in (2,\infty)$ and $q\in [2,\infty)$ and initial conditions $u_0$ in the real interpolation space \XAp we prove existence of unique strong solution with trajectories in L^p(0,T;\Dom(A))\cap C([0,T];\XAp), provided the non-linearities F:[0,T]\times \Dom(A)\to L^q(\mathcal{O},\mu) and B:[0,T]\times \Dom(A) \to \g(H,\Dom(A^{\frac12})) are of linear growth and Lipschitz continuous in their second variables with small enough Lipschitz constants. Extensions to the case where $A$ is an adapted operator-valued process are considered as well. Various applications to stochastic partial differential equations are worked out in detail. These include higher-order and time-dependent parabolic equations and the Navier-Stokes equation on a smooth bounded domain \OO\subseteq \R^d with $d\ge 2$. For the latter, the existence of a unique strong local solution with values in (H^{1,q}(\OO))^d is shown.Comment: Accepted for publication in SIAM Journal on Mathematical Analysi

A comparative linear mean-square stability analysis of Maruyama- and Milstein-type methods

In this article we compare the mean-square stability properties of the Theta-Maruyama and Theta-Milstein method that are used to solve stochastic differential equations. For the linear stability analysis, we propose an extension of the standard geometric Brownian motion as a test equation and consider a scalar linear test equation with several multiplicative noise terms. This test equation allows to begin investigating the influence of multi-dimensional noise on the stability behaviour of the methods while the analysis is still tractable. Our findings include: (i) the stability condition for the Theta-Milstein method and thus, for some choices of Theta, the conditions on the step-size, are much more restrictive than those for the Theta-Maruyama method; (ii) the precise stability region of the Theta-Milstein method explicitly depends on the noise terms. Further, we investigate the effect of introducing partially implicitness in the diffusion approximation terms of Milstein-type methods, thus obtaining the possibility to control the stability properties of these methods with a further method parameter Sigma. Numerical examples illustrate the results and provide a comparison of the stability behaviour of the different methods.Comment: 19 pages, 10 figure