<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE ArticleSet PUBLIC "-//NLM//DTD PubMed 2.7//EN" "https://dtd.nlm.nih.gov/ncbi/pubmed/in/PubMed.dtd">
<ArticleSet>
<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>A gradient projection method for solving nonlinear optimal control problems with time-varying delays</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>363</FirstPage>
			<LastPage>378</LastPage>
			<ELocationID EIdType="pii">9002</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.30121.2691</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Seyed Mojtaba</FirstName>
					<LastName>Meshkani</LastName>
<Affiliation>Department of Mathematics, Payame Noor University, 19395-4697, Tehran, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Sohrab</FirstName>
					<LastName>Effati</LastName>
<Affiliation>Department of Applied Mathematics, Faculty of Mathematical Science, Ferdowsi University of Mashhad, Mashhad, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Aghileh</FirstName>
					<LastName>Heydari</LastName>
<Affiliation>Department of Mathematics, Payame Noor University, 19395-4697, Tehran, Iran</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>04</Month>
					<Day>05</Day>
				</PubDate>
			</History>
		<Abstract>An effective numerical method using gradient projection is proposed for solving an optimal control problems that involve time-varying delays in control and state variables. First, a variational inequality is established as necessary conditions. The main idea in variational inequality is to compute the gradient of the objective functional, taking into account time-dependent delays in control and state variables. Then, an iterative scheme utilizing a projection operator is presented, followed by a convergence analysis of the method for a coercive objective functional. At the end, several examples are provided to illustrate that the theoretical finding is efficient.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Nonlinear optimal control problems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Time delay systems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">variational inequality</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Time-varying delay</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Gradient projection method</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9002_88fdaf2f039648457e05a8b34416b31d.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>A robust fitted finite difference method for semi-linear two-parameter singularly perturbed PDEs</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>379</FirstPage>
			<LastPage>405</LastPage>
			<ELocationID EIdType="pii">9196</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.30916.2776</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Mekashaw Ali</FirstName>
					<LastName>Mohye</LastName>
<Affiliation>Department of Mathematics, College of Natural and Computational Science, Wolkite University, Wolkite, Ethiopia</Affiliation>

</Author>
<Author>
					<FirstName>Justin B.</FirstName>
					<LastName>Munyakazi</LastName>
<Affiliation>Department of Mathematics and Applied Mathematics University of  Western Cape</Affiliation>

</Author>
<Author>
					<FirstName>Tekle Gemechu</FirstName>
					<LastName>Dinka</LastName>
<Affiliation>Department of Applied Mathematics   Adama Science and Technology University</Affiliation>

</Author>
<Author>
					<FirstName>Yusuf Hussen</FirstName>
					<LastName>Haji</LastName>
<Affiliation>Department of Mathematics, Arsi University, Asela, Ethiopia</Affiliation>

</Author>
<Author>
					<FirstName>Abe NUra</FirstName>
					<LastName>Ware</LastName>
<Affiliation>Department of Mathematics, Arsi University, Asela, Ethiopia</Affiliation>

</Author>
<Author>
					<FirstName>Jemal Muhammed</FirstName>
					<LastName>Ahmed</LastName>
<Affiliation>Department of Mathematics, Oda Bultum University, Chiro, Ethiopia</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>06</Month>
					<Day>09</Day>
				</PubDate>
			</History>
		<Abstract>In this article, a new numerical approach is developed for nonlinear two-parameter singularly perturbed initial-boundary value problems. The implicit backward Euler discretization for the time derivative and the fitted operator technique in the spatial domain are employed. Newton&#039;s quasilinearization technique is applied to the nonlinear terms. An investigation of parameter-uniform error estimates shows that the developed approach is first-order accurate in both time and space. However, a temporal mesh refinement technique is introduced to improve the order of accuracy to two. Two examples are provided and implemented in Python to validate the applicability of the method, and the results are displayed in tables and graphs.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Nonlinear singularly perturbed problems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">quasilinearization</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">fitted operator finite difference method</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">uniform convergence</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9196_a1da3bcb9a2e740a13facab78cc434d7.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Modeling and analysis of deforestation and pollution dynamics induced by industrialization using the fractal-fractional Atangana-Baleanu derivative</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>407</FirstPage>
			<LastPage>430</LastPage>
			<ELocationID EIdType="pii">9198</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31414.2822</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Nobin</FirstName>
					<LastName>Daimary</LastName>
<Affiliation>Gauhati University, Assam, India</Affiliation>

</Author>
<Author>
					<FirstName>Ranu</FirstName>
					<LastName>Paul</LastName>
<Affiliation>Gauhati University, Assam, India</Affiliation>
<Identifier Source="ORCID">0009-0007-1614-3072</Identifier>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>15</Day>
				</PubDate>
			</History>
		<Abstract>This study presents a fractal-fractional model in the Atangana–Baleanu sense to investigate the dynamics of deforestation and pollution driven by industrialization. The model is analyzed for positivity and boundedness, and the existence and uniqueness of its solution are established using fixed-point theory. The system’s equilibrium points are identified, and the threshold parameter &lt;br /&gt;$\mathfrak{R_0}$ is determined, with local asymptotic stability confirmed for all equilibria. Sensitivity analysis highlights the key parameters influencing $\mathfrak{R_0}$, while Ulam–Hyers stability ensures robustness of the solution. Lagrangian polynomial interpolation is employed to approximate the solution, and phase portraits along with numerical simulations in Matlab illustrate the model’s dynamic behavior. The results demonstrate that the fractal-fractional approach provides a comprehensive framework for capturing complex environmental interactions, offering valuable insights into the effects of industrialization on deforestation and pollution.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Atangana-Baleanu fractal–fractional operator</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Stability</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Ulam-Hyres stability</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9198_9e76332a786da2f23b55998d929a734e.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Hunting cooperation in prey-predator models: spatiotemporal patterns and bifurcation analysis with holling type IV response</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>431</FirstPage>
			<LastPage>451</LastPage>
			<ELocationID EIdType="pii">9201</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31303.2805</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Atish Kumar</FirstName>
					<LastName>Sethy</LastName>
<Affiliation>Central University of Odisha, Koraput, Sunabeda N.A.D Post Office, 763004
Department of Mathematics</Affiliation>

</Author>
<Author>
					<FirstName>Jyotiska</FirstName>
					<LastName>Datta</LastName>
<Affiliation>Central University of Odisha, Koraput, Sunabeda N.A.D Post Office, 763004
Department of Mathematics</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>07</Month>
					<Day>30</Day>
				</PubDate>
			</History>
		<Abstract>This study examines a class of predator models that incorporate cooperative predation within&lt;br /&gt;specialized carnivore populations. The functional response is parameterized, and numerical simulations are employed to support the analytical investigation of pattern formation potential. The principal finding of this work is that stable Turing patterns, such as stripes, can emerge when predator distributions are more localized than those of their prey. Specialized predator groups that cooperate in hunting contribute to the formation of prey aggregation zones (roost patches), as cooperation enhances predation efficiency. The results demonstrate that although predators exhibit limited mobility, cooperative behavior during hunting promotes both successful predation and long-term coexistence with prey populations.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Predation cooperation</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">basin of attraction</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">limit cycle</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">turing pattern</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">global Stability</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9201_97e708fda3802908678f57fe0726a202.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Higher order numerical method for a class of singularly perturbed time dependent nonlinear reaction diffusion problems</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>453</FirstPage>
			<LastPage>470</LastPage>
			<ELocationID EIdType="pii">9215</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31159.2791</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Manikandan</FirstName>
					<LastName>Mariappan</LastName>
<Affiliation>Department of Mathematics, School of Engineering, Presidency University, Bengaluru - 560 064, Karnataka, India</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>07</Month>
					<Day>12</Day>
				</PubDate>
			</History>
		<Abstract>Nonlinear science plays an important role in modern technology. Because of the limitations over the linear theories and the chaotic nature of the problems in this technological era, investigation of nonlinear problems has become indispensable to analyse the dynamics of complicated and multi scale characteristics problems. This article aims at the analysis and implementation of a numerical method for a class of singularly perturbed time dependent nonlinear reaction diffusion problems. Together with classical finite difference operators, a piecewise uniform Shishkin mesh in the spatial direction and a uniform mesh in the temporal direction are used to formulate a new numerical method to solve the class of problems. The method is proved to be second order convergent in space and first order convergent in time uniformly with respect to the perturbation parameter. Numerical experiments are included which support the theoretical results.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Singularly perturbed time dependent nonlinear reaction diffusion problems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">boundary layers</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">finite difference scheme</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Shishkin mesh</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">parameter uniform convergence</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9215_64aae1dffcb8073d33edc7f4919f6ffc.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Convex optimization approach for the path-following problem of two collaborative robots</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>471</FirstPage>
			<LastPage>487</LastPage>
			<ELocationID EIdType="pii">9222</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31212.2798</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Maryam</FirstName>
					<LastName>Kia-Lashgami</LastName>
<Affiliation>Department of Applied Mathematics, Faculty of Mathematical Sciences, University of
Guilan, Rasht, Iran.</Affiliation>

</Author>
<Author>
					<FirstName>Mohammad</FirstName>
					<LastName>Keyanpour</LastName>
<Affiliation>Department of Applied Mathematics, Faculty of Mathematical Sciences,
University of Guilan, P.O. Box 1914, Rasht, Iran.</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>07</Month>
					<Day>19</Day>
				</PubDate>
			</History>
		<Abstract>This paper investigates the time-optimal path-tracking problem for a collaborative robotic system, considering some limitations and dynamic characteristics. This problem is formulated for a robotic system consisting of two-link planar manipulators with and without bar cases along a predetermined geometric path in minimum time. The main challenges are to satisfy both the co-position and co-time conditions of the end-effector movement, as well as the physical limitations of the applied torque to the joints. Through discretization and convexification, we convert the problem into a convex cone optimization problem. The numerical example confirms the effectiveness of the method.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Convex problem</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Mathematical modeling</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Collaborative robots</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Optimization</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Control</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9222_b482a06a3e5d0e335f0284dbeb87add2.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Mittag-Leffler wavelet-based numerical method for fractional pantograph delay differential equations</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>489</FirstPage>
			<LastPage>508</LastPage>
			<ELocationID EIdType="pii">9233</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31321.2807</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Arezoo</FirstName>
					<LastName>Ghasempour</LastName>
<Affiliation>Department of Mathematics, Faculty of Mathematical Sciences, Alzahra University, Tehran, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Yadollah</FirstName>
					<LastName>Ordokhani</LastName>
<Affiliation>Department of Mathematics, Faculty of Mathematical Sciences, Alzahra University, Tehran, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Mohsen</FirstName>
					<LastName>Razzaghi</LastName>
<Affiliation>Department of Mathematics and Statistics, Mississippi State University, MS, USA</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>04</Day>
				</PubDate>
			</History>
		<Abstract>This paper proposes a robust numerical framework for solving fractional pantograph delay differential equations. The approach leverages the Riemann–Liouville fractional integral operator, represented through Mittag-Leffler wavelet functions within a collocation-based scheme. To facilitate computation, an operational matrix is constructed, enabling the transformation of the fractional differential system into a system of algebraic equations. The proposed method’s accuracy, stability, and convergence are rigorously validated through comprehensive numerical experiments.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Fractional pantograph differential equations</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Mittag-Leffler wavelets</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">operational matrix</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9233_c9128f6e51cdf96bde74edba9cb69641.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>An accelerated method for solving constrained multi-objective optimization</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>509</FirstPage>
			<LastPage>532</LastPage>
			<ELocationID EIdType="pii">9240</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.30721.2753</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Niloofar</FirstName>
					<LastName>Salehi Mokari</LastName>
<Affiliation>Department of Mathematics, Faculty of Mathematical Sciences and Computer, Shahid Chamran University of Ahvaz, Ahvaz, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Hadi</FirstName>
					<LastName>Basirzadeh</LastName>
<Affiliation>Department of Mathematics, Faculty of Mathematical Sciences and Computer, Shahid Chamran University of Ahvaz, Ahvaz, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Vahid</FirstName>
					<LastName>Morovati</LastName>
<Affiliation>Department of Mathematics, University of Hormozgan, Bandarabbas, Iran</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>05</Month>
					<Day>20</Day>
				</PubDate>
			</History>
		<Abstract>A novel non-parametric algorithm is introduced for solving constrained multi-objective optimization problems. At each iteration, a convex sub-problem is solved to determine the search direction, while a non-monotone line search technique is used to determine the step size. An adaptive acceleration term, computed from changes in the search directions, is incorporated to scale the step and dynamically enhance convergence performance. The algorithm’s effectiveness relies on a diverse set of initial feasible solutions to accurately approximate the non-dominated boundary. Benchmark tests validate the approach, with Pareto fronts compared to those obtained using the Zoutendijk method. Numerical evaluations demonstrate superior performance in terms of convergence rate and solution quality. The algorithm is also applied to a real-world engineering design problem involving speed reduction, highlighting its computational efficiency and robustness in practical applications.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Constrained Multi-objective Optimization Problems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Feasible Direction Methods</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Line-Search Techniques</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Pareto Critical Point</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9240_6325aaf7c4dccebd2cd9ddacf487bb28.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>The method based on quintic B-spline functions for addressing time-fractional advection-dispersion equations</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>533</FirstPage>
			<LastPage>556</LastPage>
			<ELocationID EIdType="pii">9243</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31501.2833</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Mohamed</FirstName>
					<LastName>Adel</LastName>
<Affiliation>Department of Mathematics, Faculty of Science, Islamic University of Madinah,
 Medina, KSA.</Affiliation>

</Author>
<Author>
					<FirstName>Somayieh</FirstName>
					<LastName>Abdi-mazraeh</LastName>
<Affiliation>Department of Applied Mathematics, Faculty of Mathematics, Statistics and
 Computer Sciences,
 University of Tabriz, Tabriz, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Golamreza</FirstName>
					<LastName>Zaki</LastName>
<Affiliation>Department of Applied Mathematics, Faculty of Mathematics, Statistics and
 Computer Sciences,
 University of Tabriz, Tabriz, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Safar</FirstName>
					<LastName>Irandust-Pakchin</LastName>
<Affiliation>Department of Mathematics</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>25</Day>
				</PubDate>
			</History>
		<Abstract>This paper introduces a numerical method designed to address the fractional time advection-dispersion equation. Initially, the time dimension is discretized by employing the L1 method. Subsequently, quintic B-spline functions are utilized for the discretization of the spatial dimension. This &lt;br /&gt;approach yields a system of algebraic equations that can be efficiently solved. The proposed method is proven to be unconditionally stable. Numerical experiments provide compelling evidence of the method’s eﬃciency and eﬀectiveness</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Advection-dispersion equation</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Caputo fractional derivative</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">quintic B-spline functions</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9243_9250861472f7a03bbcd991807a6d0570.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Numerical pricing of American options under a nonlinear Black-Scholes framework with mixed fractional Brownian motion</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>557</FirstPage>
			<LastPage>577</LastPage>
			<ELocationID EIdType="pii">9245</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.30891.2773</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Afshin</FirstName>
					<LastName>Babaei</LastName>
<Affiliation>Faculty of MAthematical sciences, University of Mazandaran, Babolsar, Iran.</Affiliation>

</Author>
<Author>
					<FirstName>Maryam</FirstName>
					<LastName>Rezaei</LastName>
<Affiliation>Faculty of Mathematical Sciences, University of Mazandaran, Babolsar, Iran.</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>06</Month>
					<Day>05</Day>
				</PubDate>
			</History>
		<Abstract>Transaction costs significantly impact option pricing and trading strategies in financial markets‎. ‎This study investigates the valuation of American options under transaction costs‎, ‎modeled as a linear function of the underlying asset price‎. ‎To capture long-range dependence in asset returns‎, ‎the underlying dynamics are described by a mixed fractional Brownian motion (fBm)‎. ‎The model incorporates dividend-paying stocks‎, ‎along with time-varying interest and dividend rates‎. ‎A compact finite difference scheme is developed to solve the resulting nonlinear Black-Scholes equation‎, ‎ensuring numerical stability and accuracy‎. ‎The proposed framework offers an efficient approach for pricing American options in realistic market conditions.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Mixed fractional Brownian motion‎</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">‎Transaction costs‎</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">‎American options‎</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">‎Compact difference scheme</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9245_9a406f97ece14eddc4d6d2bf17bca9b5.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Blood flow through a stenosed artery under MHD two-phase conditions with Hall effect, radiation and diffusion</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>579</FirstPage>
			<LastPage>594</LastPage>
			<ELocationID EIdType="pii">9247</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31551.2840</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Surendra Kumar</FirstName>
					<LastName>Agarwal</LastName>
<Affiliation>IIS Deemed to be University, Jaipur
MANSAROVAR</Affiliation>

</Author>
<Author>
					<FirstName>Gauri</FirstName>
					<LastName>Sethi</LastName>
<Affiliation>Department of Mathematics IIS Deemed to be University Jaipur Rajasthan</Affiliation>

</Author>
<Author>
					<FirstName>Astha</FirstName>
					<LastName>Pareek</LastName>
<Affiliation>IIS Deemed to be University, Jaipur
MANSAROVAR</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>09</Month>
					<Day>01</Day>
				</PubDate>
			</History>
		<Abstract>This study investigates the combined effects of heat and mass transfer on two-phase blood flow through a stenosed artery, considering the influence of the Hall current. Blood is modeled as a Newtonian fluid in both the central core and plasma regions. A system of differential equations governing momentum, temperature, and concentration is developed separately for each region. The analysis incorporates magnetic field, thermal radiation, and Hall current effects. Key flow characteristics, including flow resistance, total volumetric flow rate, and wall shear stress, are evaluated for varying magnetic field strengths, radiation parameters, and Hall current intensities. The findings indicate that stronger magnetic fields and radiation levels lead to a reduction in blood flow velocity and temperature. The inclusion of the Hall current introduces a cross-flow component due to the induced electric field, further modifying the velocity distribution, particularly in the plasma region. Moreover, an increase in the Schmidt number enhances the concentration profiles in both the core and plasma regions. Overall, the Hall effect significantly alters the magnetic interaction with the flow, impacting the distribution of mechanical and thermal quantities throughout the arterial segment. The results have potential biomedical applications in magnetic resonance imaging (MRI), targeted drug delivery, and electromagnetic therapy, where controlled magnetic and electric fields influence blood flow and heat transfer in vascular systems.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Blood flow</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Lorentz force</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Artery</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Stenosis</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Heat-Mass transfer</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9247_d446f8e6bed93c2cb19b1f055ee8aa10.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>A fast and cheap approach for strengthening Lagrangian bound for the generalized Celis-Dennis-Tapia subproblem</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>595</FirstPage>
			<LastPage>608</LastPage>
			<ELocationID EIdType="pii">9264</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31942.2884</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Temadher Alassiry</FirstName>
					<LastName>Almaadeed</LastName>
<Affiliation>Qatar University-College of Arts and Sciences-Dept Mathematics and Statistics, P.O. Box 2713
Qatar University, Doha- Qatar</Affiliation>

</Author>
<Author>
					<FirstName>Abdelouahed</FirstName>
					<LastName>Hamdi</LastName>
<Affiliation>Qatar University-College of Arts and Sciences-Dept Mathematics and Statistics, P.O. Box 2713,  Qatar University, Doha- Qatar</Affiliation>

</Author>
<Author>
					<FirstName>Akram</FirstName>
					<LastName>Taati</LastName>
<Affiliation>Faculty of Mathematical Sciences, University of Guilan, Rasht, Iran.</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>10</Month>
					<Day>11</Day>
				</PubDate>
			</History>
		<Abstract>In this paper, we consider the generalized Celis-Dennis-Tapia problem&lt;br /&gt;which is the problem of minimizing a nonconvex quadratic function subject to&lt;br /&gt;two quadratic inequality constraints, one of which being convex. When there is&lt;br /&gt;a positive duality gap, by exploiting an equivalent form of the dual Lagrangian&lt;br /&gt;problem, we propose to improve the dual bound by adding one or two linear cuts&lt;br /&gt;to the Lagrangian relaxation. The present work is motivated by and generalizes&lt;br /&gt;the results of [14] for the problem with two strictly convex quadratic constraints.&lt;br /&gt;Our main contribution is to show that one can include the feasible region in a con-&lt;br /&gt;vex set and then follow the approach in [14] to construct the linear cuts based on&lt;br /&gt;supporting hyperplanes of the convex set. Numerical experiments are conducted&lt;br /&gt;to assess the quality of the proposed bounds.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">quadratically constrained quadratic programming</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Celis-Dennis-Tapia problem</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">dual Lagrangian bound</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Supporting hyperplane</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9264_3bb0c7a429b19885a3efa780276fd8df.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Mathematical model and hybrid meta-heuristic solution approaches for hub location problem with hybrid drone-airplane delivery mode</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>609</FirstPage>
			<LastPage>646</LastPage>
			<ELocationID EIdType="pii">9269</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31720.2857</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Mehrnaz</FirstName>
					<LastName>Mohebbi</LastName>
<Affiliation>Shiraz University of Technology</Affiliation>

</Author>
<Author>
					<FirstName>Hamid Reza</FirstName>
					<LastName>Maleki</LastName>
<Affiliation>Faculty of Mathematics, Shiraz University of Technology, Shiraz, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Sadegh</FirstName>
					<LastName>Niroomand</LastName>
<Affiliation>Shiraz University of Technology</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>09</Month>
					<Day>19</Day>
				</PubDate>
			</History>
		<Abstract>This study addresses the integrated hub location and drone delivery problem, an area with few5&lt;br /&gt;prior investigations. We propose a bi-objective integer linear programming model to minimize total cost and total drone route time. A novel three-zone structure allows drone transfers between zones via cargo planes, enhancing realism and complexity. Drone capacities are categorized as light and heavy, improving allocation flexibility. Due to the model’s complexity, several metaheuristic algorithms including Genetic Algorithm, Differential Evolution, Simulated Annealing, and their hybrid versions (SA-GA and SA-DE) are developed and compared. Parameter tuning is performed using the Taguchi method. Computational experiments on various instances show that hybrid algorithms outperform classical methods and scale effectively for larger problems, providing a practical and integrated framework for hub location and drone delivery planning.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Hub location problem</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Drone delivery problem</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Meta-heuristic algorithm</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Hybrid meta-heuristic algorithm</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Taguchi Method</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9269_5940ad57fb95609dea042874deb754db.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Two-stage multi-objective technology portfolio planning under resource constraints (case study: Iranian technology development fund)</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>645</FirstPage>
			<LastPage>671</LastPage>
			<ELocationID EIdType="pii">9289</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31547.2839</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Marzieh</FirstName>
					<LastName>Shaverdi</LastName>
<Affiliation>Department of STI Financing and Economics,
National Research Institution for Science Policy (NRISP), Tehran, Iran.</Affiliation>

</Author>
<Author>
					<FirstName>Khatereh</FirstName>
					<LastName>Ghorbani-Moghadam</LastName>
<Affiliation>Mosaheb Institute of Mathematics</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>09</Month>
					<Day>01</Day>
				</PubDate>
			</History>
		<Abstract>This study proposes a novel two-stage multi-objective framework for optimal technology portfolio planning and resource allocation under constraints, specifically for Technology Development Funds (TDFs). The integrated methodology combines the Analytic Network Process (ANP) for prioritizing strategic technology fields with a multi-period Mixed-Integer Linear Programming (MILP) model, solved using a Revised Multi-Choice Goal Programming (RMCGP) approach. The model’s objectives are to maximize technological export potential, maximize international technological cooperation, and minimize financial risk, while incorporating critical real-world mechanisms such as staged financing contingent on Technology Readiness Level (TRL) progress, loan moratorium, and repayment periods.&lt;br /&gt;The framework was validated through a real-world case study of an Iranian Technology Development Fund (ITDF), involving eight technology fields and up to 30 projects per field. Key quantitative results demonstrate the model’s efficacy: by reducing the risk objective’s weight from 0.3 to 0.1, the number of approved projects increased over fivefold (from 12 to 65), and the total allocated resources surged nearly tenfold (from $22.2 million to $217.5 million). Sensitivity analysis revealed that fields with high export potential and collaboration capacity (e.g., Advanced Machinery) received the highest funding, while the staged financing mechanism successfully identified and terminated 25% of projects for insufficient technical progress after the first stage. The proposed model provides a robust decision-support tool for policymakers to enhance the strategic impact and financial efficiency of national technology investments.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Technology prioritization</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Technology portfolio selection</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Staged financing</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Moratorium period</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Reinvestment strategy</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9289_c5c9c33fed7df383c3c4c61b0fe3db6b.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Full and partial controllability of the Kermack-Mckendrick system with time- varying incidence rates</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>673</FirstPage>
			<LastPage>697</LastPage>
			<ELocationID EIdType="pii">9295</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31390.2818</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Hamza</FirstName>
					<LastName>El Mahjour</LastName>
<Affiliation>MASI research Team
Department of Information Systems and Communication
ENSAT, Abdelmalek Esaadi University, Morocco</Affiliation>

</Author>
<Author>
					<FirstName>Aadil</FirstName>
					<LastName>Lahrouz</LastName>
<Affiliation>LAM research laboratory
Department of Mathematics
FSTT, Abdelmalek Essaadi University</Affiliation>

</Author>
<Author>
					<FirstName>Omar</FirstName>
					<LastName>Zakary</LastName>
<Affiliation>Statistics and Modelling research team
Department of Mathematics
FSBM, University of Hassan II</Affiliation>

</Author>
<Author>
					<FirstName>Mariam</FirstName>
					<LastName>Redouane</LastName>
<Affiliation>LAM research laboratory
Department of Mathematics
FSTT, Abdelmalek Essaadi University</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>12</Day>
				</PubDate>
			</History>
		<Abstract>This study contributes to epidemic control literature by introducing a time-varying inci-&lt;br /&gt;dence rate and establishing global controllability of the nonlinear SIR system, offering a&lt;br /&gt;practical framework for adaptive control strategies. We derive explicit solutions for partial&lt;br /&gt;controllability, demonstrating the feasibility of controlling the infected population, pro-&lt;br /&gt;viding guidance for outbreak management. Numerical methods exploiting an algorithmic&lt;br /&gt;approach achieve full control, targeting a desired state (Sd, Id). A novel hybrid method&lt;br /&gt;integrates analytical solutions with algorithmic optimization, leveraging explicit expres-&lt;br /&gt;sions for I(t) and S(t) to enhance precision and efficiency of epidemic control strategies,&lt;br /&gt;advancing adaptive management approaches</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Epidemic model</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Varying infection rate</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Full Control</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Partial Control</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Hybrid Method</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9295_fd6f0670afa9641d22cff7f55adf61dd.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Dual analysis of myocardial infarction using fractional mathematical modeling and machine learning</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>699</FirstPage>
			<LastPage>722</LastPage>
			<ELocationID EIdType="pii">9299</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31534.2837</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Tharmalingam</FirstName>
					<LastName>Gunasekar</LastName>
<Affiliation>Department of Mathematics, Vel Tech Rangarajan Dr. Sagunthala R\&amp;D Institute of Science and Technology, Chennai - 600062, Tamil Nadu, India</Affiliation>

</Author>
<Author>
					<FirstName>Sumaiya</FirstName>
					<LastName>Banu</LastName>
<Affiliation>Department of Mathematics, Vel Tech Rangarajan Dr. Sagunthala R\&amp;D Institute of Science and Technology, Chennai - 600062, Tamil Nadu, India</Affiliation>

</Author>
<Author>
					<FirstName>Shyam Sundar</FirstName>
					<LastName>Santra</LastName>
<Affiliation>JIS College of Engineering, Kalyani, Nadia, West Bengal</Affiliation>

</Author>
<Author>
					<FirstName>Dumitru</FirstName>
					<LastName>Baleanu</LastName>
<Affiliation>Department of Computer Science and Mathematics, Lebanese American University, Beirut - 11022801, Lebanon</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>30</Day>
				</PubDate>
			</History>
		<Abstract>This paper presents a novel fractional-order mathematical model of myocardial infarction in women who are users of combined oral contraceptive pill and who also develop comorbidity due to various reasons. The system of equations incorporate Caputo fractional derivative to capture memory effects of the model. Existence and uniqueness of solution of the mathematical model is derived. Numerical simulations were rigorously conducted on the math model with varying fractional order namely, $0.3$, $0.5$ and $0.8$ using Euler&#039;s method. The numerical results thus obtained are simulated by Adam&#039;s method for 200 days period. The output from these simulations form the dataset of the Bayesian regularization neural network (BRNN) with dataset split for training, testing and validatating the computational model. Bayesian regularization is incorporated to handle overfitting efficiently. Root Mean Square Error (RMSE) are computed for all three fractional orders respectively. Regression analysis is conducted which yielded perfect correlation \((R=1)\) accross the all datasets. The combined mathematical and computational analysis form a strong layout in myocardial infarction risk prediction, diagnosis and treatment in young women.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Myocardial Infarction</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Fractional Mathematical Modeling</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Caputo derivative</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Combined Oral Contraceptive Pill</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Bayesian Regularization Neural Network</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9299_6593aeb65997ccb0968e730a257ed0b3.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Classification of flow behavior near generalized equilibrium points in piecewise smooth systems</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>723</FirstPage>
			<LastPage>739</LastPage>
			<ELocationID EIdType="pii">9301</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.30519.2738</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Zohreh</FirstName>
					<LastName>Zarei</LastName>
<Affiliation>Department of Math., University of Hormozgan, Bandar Abbas, Iran</Affiliation>

</Author>
<Author>
					<FirstName>Majid</FirstName>
					<LastName>Karimi Amaleh</LastName>
<Affiliation>Department of Math, University of Hormozgan, Bandar Abbas, Iran</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>06</Month>
					<Day>01</Day>
				</PubDate>
			</History>
		<Abstract>The aim of this paper is to classify the various states of flow behavior for piecewise smooth systems near generalized equilibrium points. Seven categories are introduced based on the sign of the vector field across the discontinuity boundary, each encompassing distinct dynamical configurations. We investigate how a small perturbation parameter $\mu$ influences the existence, type, and stability of generalized singular points in planar piecewise linear systems. Starting with a one-dimensional example to illustrate core mechanisms, we extend the analysis to two dimensions, providing a detailed classification grounded in the signs of the system’s components. Our results yield a comprehensive framework for understanding how generalized singular points govern local dynamics, including bifurcations induced by parameter variation. This work contributes to the theoretical foundation for analyzing discontinuity-induced phenomena such as sliding modes and non-smooth bifurcations.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">generalized equilibrium point</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">piecewise smooth systems</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">bifurcation analysis</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">perturbation parameter</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Stability</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9301_fe3017a88aba3c657d183b055fb140ae.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Dynamical analysis of novel minimal SEIR model incorporating asymptomatic transmission</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>741</FirstPage>
			<LastPage>760</LastPage>
			<ELocationID EIdType="pii">9302</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31474.2832</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Viralkumar D.</FirstName>
					<LastName>Patel</LastName>
<Affiliation>Department of Mathematics, Sardar Patel University, Vallabh  Vidyanagar, Anand, Gujarat, India.

Department of Mathematics, Shri Alpesh N. Patel Post Graduate Institute of Science &amp;  Research, Anand.</Affiliation>

</Author>
<Author>
					<FirstName>Abdulvahid H.</FirstName>
					<LastName>Hasmani</LastName>
<Affiliation>Department of Mathematics, Sardar Patel University, Vallabh Vidyanagar, Gujarat, India.</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>25</Day>
				</PubDate>
			</History>
		<Abstract>We propose a modified SEIR model that includes asymptomatic transmission directly in the infection term, avoiding the need for a separate asymptomatic compartment while keeping the model realistic. The basic reproduction number ($\mathcal{R}_0$) is calculated to measure the potential for disease spread. Local stability analysis shows that the disease-free equilibrium is stable when $\mathcal{R}_0 &lt; 1$ and unstable when $\mathcal{R}_0 &gt; 1$, while the endemic equilibrium is locally stable in the latter case. A forward bifurcation at $\mathcal{R}_0 = 1$ is identified, indicating a smooth transition from the disease-free equilibrium to a unique endemic equilibrium without coexistence of the two equilibria. Global stability results show that the disease-free state is globally asymptotically stable for $\mathcal{R}_0 \leq 1$, and the endemic state is globally asymptotically stable for $\mathcal{R}_0 &gt; 1$. Simulations using early COVID‑19 data support these findings, showing that higher asymptomatic transmission prolongs outbreaks, increases peaks, and delays elimination. Evaluation of control strategies reveals that isolation is more effective than testing alone, and their combination produces the greatest overall reduction in disease spread under appropriate assumptions.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Modified $SEIR$ Model</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Asymptomatic Infection</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Stability analysis</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">forward bifurcation</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">COVID-19 simulation</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9302_30a58393c2cfbe438486b4983c339a78.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Novel Legendre-Jaiswal functions for solving time-space fractional partial differential equations</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>761</FirstPage>
			<LastPage>786</LastPage>
			<ELocationID EIdType="pii">9305</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2025.31438.2827</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Santoshi</FirstName>
					<LastName>Tarei</LastName>
<Affiliation>Department of Mathematics
National Institute of Technology Rourkela</Affiliation>

</Author>
<Author>
					<FirstName>Ankur</FirstName>
					<LastName>Kanaujiya</LastName>
<Affiliation>Department of Mathematics, NIT Rourkela, Odisha, India</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>08</Month>
					<Day>18</Day>
				</PubDate>
			</History>
		<Abstract>This paper examines a new fractional function based on Legendre and Jaiswal polynomials to solve linear and nonlinear time-space fractional partial differential equations of linear and nonlinear class. The Caputo sense is applied while using the fractional derivative. These problems can be solved using the collocation method, operational, and pseudo-operational matrices of integer and fractional-order integration. Using operational matrices, pseudo-operational matrices, and the collocation method, the problem is transformed into a system of algebraic equations. An upper bound on the error of the fractional-order integral operational matrix is computed. Furthermore, a detailed stability and convergence analysis of the collocation scheme presented to validate the robustness of the numerical approach. The applicability and effectiveness of the approach are demonstrated through several benchmark examples, including linear and non-linear fractional convection-diffusion, convection-diffusion-reaction, and nonlinear Fisher&#039;s equation. The numerical results confirm that the proposed method is stable, rapidly convergent, and highly accurate, outperforming several existing techniques in both efficiency and precision.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Legendre polynomial</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Jaiswal polynomial</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Fractional partial differential equation</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">error analysis</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9305_65082bf065e0cbe54d0cc6155545d721.pdf</ArchiveCopySource>
</Article>

<Article>
<Journal>
				<PublisherName>University of Guilan</PublisherName>
				<JournalTitle>Journal of Mathematical Modeling</JournalTitle>
				<Issn>2345-394X</Issn>
				<Volume>14</Volume>
				<Issue>2</Issue>
				<PubDate PubStatus="epublish">
					<Year>2026</Year>
					<Month>05</Month>
					<Day>01</Day>
				</PubDate>
			</Journal>
<ArticleTitle>Modeling log-volatility with zero returns: empirical evidence for asymmetric SV and log-GARCH models</ArticleTitle>
<VernacularTitle></VernacularTitle>
			<FirstPage>787</FirstPage>
			<LastPage>801</LastPage>
			<ELocationID EIdType="pii">9400</ELocationID>
			
<ELocationID EIdType="doi">10.22124/jmm.2026.31979.2889</ELocationID>
			
			<Language>EN</Language>
<AuthorList>
<Author>
					<FirstName>Abdeljalil</FirstName>
					<LastName>Settar</LastName>
<Affiliation>Laboratory of mathematics, computer science and applications (LMCSA), FST Mohammedia,
Hassan II university, 20650 Casablanca, Morocco</Affiliation>

</Author>
<Author>
					<FirstName>Sara</FirstName>
					<LastName>Chegdal</LastName>
<Affiliation>Laboratory of mathematics, computer science and applications (LMCSA), FST Mohammedia,
Hassan II university, 20650 Casablanca, Morocco</Affiliation>

</Author>
<Author>
					<FirstName>Mustapha</FirstName>
					<LastName>Kabil</LastName>
<Affiliation>Laboratory of mathematics, computer science and applications (LMCSA), FST Mohammedia,
Hassan II university, 20650 Casablanca, Morocco</Affiliation>

</Author>
<Author>
					<FirstName>Ghassane</FirstName>
					<LastName>Benrhmach</LastName>
<Affiliation>Department of Mathematics and Statistics, College of Engineering, Abu Dhabi University, 59911 Abu
Dhabi, UAE</Affiliation>

</Author>
</AuthorList>
				<PublicationType>Journal Article</PublicationType>
			<History>
				<PubDate PubStatus="received">
					<Year>2025</Year>
					<Month>10</Month>
					<Day>15</Day>
				</PubDate>
			</History>
		<Abstract>In this work, we address the challenges posed by zero returns in both stochastic volatility (SV) and log-GARCH models in their asymmetric form. Building upon EM imputation for handling zero returns, we propose a unified approach that enhances parameter estimation robustness for both model classes. Specifically, we employ the Quasi-Maximum Likelihood (QML) estimation, incorporating the Kalman filter for both asymmetric SV and asymmetric log-GARCH models, to ensure robust parameter estimation even in the presence of zero returns. By comparing the performance of these models under our proposed framework, we provide new insights into their relative strengths in capturing the asymmetric volatility dynamics in the presence of zero returns. This contribution extends the existing literature by proposing a computational framework applicable to such models, based on a logarithmic specification of volatility.</Abstract>
		<ObjectList>
			<Object Type="keyword">
			<Param Name="value">Stochastic volatility</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">log-GARCH</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">zero returns</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">quasi-maximum likelihood</Param>
			</Object>
			<Object Type="keyword">
			<Param Name="value">Kalman filter</Param>
			</Object>
		</ObjectList>
<ArchiveCopySource DocType="pdf">https://jmm.guilan.ac.ir/article_9400_bd9ba5bd36544b3dee3440a32fed3017.pdf</ArchiveCopySource>
</Article>
</ArticleSet>
