University of GuilanJournal of Mathematical Modeling2345-394X11220230701A mixed algorithm for smooth global optimization207228620110.22124/jmm.2022.23133.2061ENRaoufZiadiLaboratory of Fundamental and Numerical Mathematics (LMFN), Department of Mathematics, University Ferhat Abbas Setif 1, 19000 Setif, Algeria0000-0002-6471-0445AbdelatifBencherif-MadaniLaboratory of Fundamental and Numerical Mathematics (LMFN), Department of Mathematics, University Ferhat Abbas Setif 1, 19000 Setif, AlgeriaJournal Article20221025This paper presents a covering algorithm for solving bound-constrained global minimization problems with a differentiable cost function. In the proposed algorithm, we suggest to explore the feasible domain using a one-dimensional global search algorithm through a number of parametric curves that are relatively spread and simultaneously scan the search space. To accelerate the corresponding algorithm, we incorporate a multivariate quasi-Newton local search algorithm to spot the lowest regions. The proposed algorithm converges in a finite number of iterations to an $\varepsilon$-approximation of the global minimum. The performance of the algorithm is demonstrated through numerical experiments on some typical test functions.https://jmm.guilan.ac.ir/article_6201_56e857ceeea5c3d0d27c810b087698e4.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Pricing American option under exponential Levy Jump-diffusion model using Random Forest instead of least square regression229244630010.22124/jmm.2022.21756.1909ENMohamedMaidoumiLAMAI, Cadi Ayyad University, Marrakech, Moroccomaidoumi-22@hotmail.frMehdiZahidLAMAI, Cadi Ayyad University, Marrakech, MoroccoBoubkerDaafiLAMAI, Cadi Ayyad University, Marrakech, MoroccoJournal Article20220211In this paper, we aim to propose a new hybrid version of the Longstaff and Schwartz algorithm under the exponential Levy Jump-diffusion model using Random Forest regression. For this purpose, we will build the evolution of the option price according to the number of paths. Further, we will show how this approach numerically depicts the convergence of the option price towards an equilibrium price when the number of simulated trajectories tends to a large number. In the second stage, we will compare this hybrid model with the classical model of the Longstaff and Schwartz algorithm (as a benchmark widely used by practitioners in pricing American options) in terms of computation time, numerical stability and accuracy. At the end of this paper, we will test both approaches on the Microsoft share “MSFT” as an example of a real market. https://jmm.guilan.ac.ir/article_6300_05e42c51116d4af34b2de83a55cadadc.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Investigation and solving of initial-boundary value problem including fourth order PDE by contour integral and asymptotic methods245255625110.22124/jmm.2023.23209.2069ENRezaDanaeiDepartment of Mathematics, Azarbaijan Shahid Madani University, Tabriz, Iran0000-0001-8485-2320MohammadJahanshahiDepartment of Mathematics, Azarbaijan Shahid Madani University, Tabriz, IranJournal Article20221107In this paper, we consider a fourth order mixed partial differential equation with some initial and boundary conditions which is unsolvable by classical methods such as Fourier, Fourier-Bircove and Laplace Transformation methods. For this problem we will apply the contour integral and asymptotic methods. The convergence of the appeared integrals, existence and uniqueness of solution, satisfying the solution and holding the given initial and boundary conditions are stablished by complex analysis theory and related contour integrals. Finally, the form of analytic and approximate solutions are given due to different cases of eigenvalues distributions in the complex plane.https://jmm.guilan.ac.ir/article_6251_f117c9029c5f3eac4a8c6d6dbd20115d.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701A new approach to solve weakly singular fractional-order delay integro-differential equations using operational matrices257275620710.22124/jmm.2023.23316.2080ENSaeedehRezabeykDepartment of Mathematics, Karaj Branch, Islamic Azad University, Karaj, IranSaeidAbbasbandyDepartment of Applied Mathematics, Faculty of Science, Imam Khomeini International University, Qazvin, 34149-16818, Iran0000-0003-3385-4152ElyasShivanianDepartment of Applied Mathematics, Faculty of Science, Imam Khomeini International University, Qazvin, 34149-16818, IranHesamDeriliDepartment of Mathematics, Karaj Branch, Islamic Azad University, Karaj, Iran0000-0001-5854-2249Journal Article20221125In this paper, we propose a new approach to solve weakly singular fractional delay integro-differential equations. In the proposed approach, we apply the operational matrices of fractional integration and delay function based on the shifted Chebyshev polynomials to approximate the solution of the considered equation. By approximating the fractional derivative of the unknown function as well as the unknown function in terms of the shifted Chebyshev polynomials and substituting these approximations into the original equation, we obtain a system of nonlinear algebraic equations. We present the convergence analysis of the proposed method. Finally, to show the accuracy and validity of the proposed method, we give some numerical examples.https://jmm.guilan.ac.ir/article_6207_f0eddb59636a94b8f6f3a8d3c7515f30.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Application of S-Boxes based on the chaotic Hindmarsh-Rose system for image encryption277300627810.22124/jmm.2023.23199.2068ENMeisamBavand SavadkouhiDepartment of Mathematics, Shahed University, Tehran, IranHamidHaj Seyyed JavadiDepartment of Mathematics and Computer Sciences, Shahed University, Tehran, IranMohammadAkbari TootkaboniDepartment of Pure Mathematics, Faculty of Mathematical Sciences, University of Guilan, IranJournal Article20221106The substitution box (S-Box) is a critical component in symmetric cipher algorithms. In this paper, we choose the Hindmarsh-Rose system to generate chaotic S-Boxes. We propose two S-Boxes based on the rotation algorithm relative to the rows (or columns) and the other based on the Zigzag transformation. The performance of the new S-Boxes is evaluated by bijective, nonlinearity, strict avalanche criterion (SAC), output bits independence criterion (BIC), differential approximation probability, linear approximation probability, and algebraic degree. The analysis results show that the presented S-Boxes have suitable cryptographic properties. Also, an image encryption algorithm based on two proposed S-Boxes, and a chaotic Hindmarsh-Rose system are presented. Experimental results show the recommended method has attained good security, and the suggested plan has potent resistance to different attacks.https://jmm.guilan.ac.ir/article_6278_6a82c8928f1ca101f71b85d237e07f83.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Symmetric-diagonal reductions as preprocessing for symmetric positive definite generalized eigenvalue solvers301322637410.22124/jmm.2023.23734.2120ENMoradAhmadnasabDepartment of Mathematics, Faculty of Science, University of Kurdistan, 66177-15175, Sanandaj, Iran0000-0002-3946-5599Journal Article20230201We discuss some potential advantages of the orthogonal symmetric-diagonal reduction in two main versions of the Schur-QR method for symmetric positive definite generalized eigenvalue problems. We also advise and use the appropriate reductions as preprocessing on the solvers, mainly the Cholesky-QR method, of the considered problems. We discuss numerical stability of the methods via providing upper bound for backward error of the computed eigenpairs and via investigating two kinds of scaled residual errors. We also propose and apply two kinds of symmetrizing which improve the stability and the performance of the methods. Numerical experiments show that the implemented versions of the Schur-QR method and the preprocessed versions of the Cholesky-QR method are usually more stable than the Cholesky-QR method. https://jmm.guilan.ac.ir/article_6374_f1b8d30f6e3adfda497477f8ae32625e.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701A new version of augmented self-scaling BFGS method323342653310.22124/jmm.2023.23425.2089ENMohamadJourakDepartment of Mathematics, Payame Noor University,
P.O. Box. 19395-3697, Tehran, IranSaeedNezhadhoseinDepartment of Mathematics, Payame Noor University,
P.O. Box. 19395-3697, Tehran, IranFarzadRahpeymaiiDepartment of Mathematics, Technical and Vocational University (TVU), Tehran, IranJournal Article20221210A new version of the augmented self-scaling memoryless BFGS quasi-Newton update, proposed in [Appl. Numer. Math. 167, 187--201, (2021)], is suggested for unconstrained optimization problems. To use the corresponding scaled parameter, the clustering of the eigenvalues of the approximate Hessian matrix about one point is applied with three approaches. The first and second approaches are based on the trace and the determinant of the matrix. The third approach is based on minimizing the measure function. The sufficient descent property is guaranteed for uniformly convex functions, and the global convergence of the proposed algorithm is proved both for the uniformly convex and general nonlinear objective functions, separately. Numerical experiments on a set of test functions of the CUTEr collection show that the proposed method is robust. In addition, the proposed algorithm is effectively applied to the salt and pepper noise elimination problem.https://jmm.guilan.ac.ir/article_6533_5e2bdae9bbc316cde5a82ba4fd3c6148.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Eigenvalue problem with fractional differential operator: Chebyshev cardinal spectral method343355668310.22124/jmm.2023.24239.2169ENAlirezaAfaridehDepartment of Mathematics, Tabriz Branch, Islamic Azad University, Tabriz, IranFarhadDastmalchi SaeiDepartment of Mathematics, Tabriz Branch, Islamic Azad University, Tabriz, Iran0000-0001-9829-2141BehzadNemati SarayDepartment of Mathematics, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan 45137-66731, Iran0000-0002-7655-8508Journal Article20230405In this paper, we intend to introduce the Sturm-Liouville fractional problem and solve it using the collocation method based on Chebyshev cardinal polynomials. To this end, we first provide an introduction to the Sturm-Liouville fractional equation. Then the Chebyshev cardinal functions are introduced along with some of their properties and the operational matrices of the derivative, fractional integral, and Caputo fractional derivative are obtained for it. Here, for the first time, we solve the equation using the operational matrix of the fractional derivative without converting it to the corresponding integral equation. In addition to efficiency and accuracy, the proposed method is simple and applicable. The convergence of the method is investigated, and an example is presented to show its accuracy and efficiency.https://jmm.guilan.ac.ir/article_6683_5a49be4f053702e71e18d4c4a590ade9.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Caputo fractional-time of a modified Cahn-Hilliard equation for the inpainting of binary images357373656910.22124/jmm.2023.22019.1938ENAnouarBen-loghfyryDepartment of mathematics, Faculty of Sciences and Technologies Mohammedia, University Hassan II, Casablanca, Morocco.0000-0003-3883-6957AbdelilahHakimLAMAI laboratory, university of Cadi Ayyad, Faculty of sciences and technology, Marrakesh, MoroccoJournal Article20220327In this work, we present a new version of the Cahn-Hilliard equation to deal with binary image inpainting. The proposed model is unique due to its memory effect ability implemented by the time fractional derivative. Also, this model has a new diffusion term that gives a topological reconnection and a well sharpness of edges and corners. We give an existence result with some numerical tests implemented by the convexity splitting to show the efficiency of the proposed model.https://jmm.guilan.ac.ir/article_6569_14fe58a8b8c38ca8a79b156afc8ca4ff.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Stability and bifurcation of stochastic chemostat model375394663710.22124/jmm.2023.24214.2165ENMehdiFatehi NiaDepartment of Mathematics, Yazd University, Yazd, IranNajmehKhajoeiDepartment of Mathematics, Yazd University, Yazd, IranJournal Article20230402The main purpose of this paper is to study dynamics of stochastic chemostat model. In this order, Taylor expansions, polar coordinate transformation and stochastic averaging method are our main tools. The stability and bifurcation of the stochastic chemostat model are considered. Some theorems provide sufficient conditions to investigate stochastic stability, $D$-bifurcation and $P$-bifurcation of the model. As a final point, to show the effects of the noise intensity and illustrate our theoretical results, some numerical simulations are presented.https://jmm.guilan.ac.ir/article_6637_e359d2e80aadd7c304c5e39ef76921fc.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701A fitted operator method of line scheme for solving two-parameter singularly perturbed parabolic convection-diffusion problems with time delay395410660110.22124/jmm.2023.23001.2039ENNaol TufaNegeroDepartment of Mathematics, Wollega University, Nekemte, Ethiopia0000-0003-1593-735XJournal Article20220928This paper presents a parameter-uniform numerical scheme for the solution of two-parameter singularly perturbed parabolic convection-diffusion problems with a delay in time. The continuous problem is semi-discretized using the Crank-Nicolson finite difference method in the temporal direction. The resulting differential equation is then discretized on a uniform mesh using the fitted operator finite difference method of line scheme. The method is shown to be accurate in $ O(\left(\Delta t \right)^{2} + N^{-2}) $, where $ N $ is the number of mesh points in spatial discretization and $ \Delta t $ is the mesh length in temporal discretization. The parameter-uniform convergence of the method is shown by establishing the theoretical error bounds. Finally, the numerical results of the test problems validate the theoretical error bounds.https://jmm.guilan.ac.ir/article_6601_796659fd522b120398d31fc34a7f5dd8.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X11220230701Robust computational technique for a class of singularly perturbed nonlinear differential equations with Robin boundary conditions411423664810.22124/jmm.2023.23515.2100ENManikandanMariappanDepartment of Mathematics, School of Engineering, Presidency University, Bengaluru - 560 064, Karnataka, India0000-0001-8268-7800IshwariyaRoseyDepartment of Science and Humanities, Amrita School of Engineering, Amrita Vishwa Vidyapeetham, Chennai-601103, Tamil Nadu, India0000-0001-5343-1187Journal Article20221228In this article, a class of singularly perturbed nonlinear differential equations with Robin boundary conditions is considered. A numerical method consists of the classical finite difference operator over a Shishkin mesh with two-mesh algorithm is constructed to solve the problems. The method is proved to be first order convergent uniformly with respect to the perturbation parameter. Experiments are carried out for two different types of Robin boundary conditions and Neumann boundary conditions as a special case of Robin boundary conditions.https://jmm.guilan.ac.ir/article_6648_6637163a609988f483a4aadd7287ee4e.pdf