Open Journal of Mathematical Sciences (OMS)

In this paper, we investigate the dynamical behavior of a four-dimensional fractional-order chaotic system using a Caputo fractional derivative. The examination of local stability reveals that the system undergoes Hopf bifurcations, leading to oscillatory states. Numerical findings based on Lyapunov exponents, bifurcation diagrams, and phase portraits validate the existence of chaos and hyperchaos over extensive parameter ranges. Furthermore, the system exhibits diverse chaotic phenomena, including self-excited attractors produced from unstable equilibria and coexisting attractors that emerge under different initial conditions. For numerical simulations, the Adams-Bashforth-Moulton predictor-corrector method is employed. Also, analytical calculations are carried out using the Maple program, and the MATLAB software program is used to illustrate the results. The results demonstrate that the proposed fractional-order system accurately represents multistability, coexisting chaotic attractors, and complex dynamics depending on parameters and derivative order.

This study presents a new extension of the length-biased truncated Lomax distribution by incorporating the sine Topp-Leone family, resulting in a flexible model called the sine Topp-Leone length-biased truncated Lomax distribution. The proposed distribution demonstrates remarkable adaptability in modeling data with increasing hazard rates. Key statistical and reliability properties of the new model are thoroughly examined, including the survival function, hazard rate, reversed hazard rate, quantile function, moments, incomplete moments, and Renyi and Tsallis entropy measures. Parameter estimation is conducted using both classical and Bayesian methods, considering symmetric and asymmetric loss functions. Due to the computational challenges of Bayesian estimation, Markov Chain Monte Carlo techniques with independent gamma priors are employed. Simulation studies confirm the consistency of the proposed estimators, showing improved accuracy with increasing sample size, with Bayesan estimates provide lower absolute biases and mean squared errors. Finally, applications to two real datasets demonstrate the superior flexibility and effectiveness of the proposed distribution compared to existing alternatives.

This work investigates a porous elastic system under the Green–Naghdi heat conduction theory. We provide an explicit characterization of the decay rate, which may be exponential or polynomial, depending on the relationship between the wave propagation speeds. To verify the asymptotic behavior of the solution, we present an algorithm for obtaining a numerical solution using finite element methods and finite differences in time.

Dynamic boundary conditions provide a stringent verification setting for physics-informed neural networks (PINNs), because the approximation must satisfy an interior parabolic equation together with an evolution law posed on the boundary. This paper develops a transparent benchmark study for the one-dimensional heat equation with Wentzell-type dynamic boundary conditions. Two manufactured solutions are employed. Benchmark A is a cosine solution for which homogeneous boundary forcing is admissible only under the compatibility condition \(k=\alpha\pi^2\); it is intentionally retained as a reduced Wentzell test because the endpoint derivatives vanish and the boundary-gradient term is inactive. Benchmark B is a fully active test with nonzero endpoint derivatives at both boundaries, so that the \(\beta u_x\) contribution is explicitly exercised. The numerical evidence is based on directly computed norms, residuals, five independent random-feature seeds, an implemented finite-difference reference solver, and a targeted ablation study. The reported neural solver is a least-squares physics-informed random-feature network with hard initial-condition enforcement, \(250\) hidden features, \(3980\) interior collocation points, and \(401\) boundary times per endpoint. Over five seeds, the mean relative errors are \(E_{L^2}=(8.65\pm1.27)\times10^{-6}\) and \(E_{L^{\infty}}=(1.12\pm0.20)\times10^{-5}\) for Benchmark A, and \(E_{L^2}=(9.45\pm2.42)\times10^{-8}\) and \(E_{L^{\infty}}=(2.77\pm0.99)\times10^{-7}\) for Benchmark B. The finite-difference baseline displays approximately second-order convergence under the parabolic scaling \(\Delta t=2h^2\). These results show that the two-case manufactured suite provides a useful verification framework for smooth one-dimensional Wentzell-type dynamics, while the scope of the claims remains restricted to this benchmark setting rather than general neural-solver superiority.

In this paper, we present new sharp forms of Iyengar-type integral inequalities for differentiable and twice-differentiable functions. Using exact integral identities and the Cauchy-Schwarz inequality, we derive optimal \(L^2\)-bounds for the deviation between the midpoint value and the integral mean of a function. These results refine and extend several classical Hermite-Hadamard and Iyengar inequalities.

This study introduces the Secant–Weibull Autoregressive Conditional Duration (SW–ACD) model and its extension with exogenous calendar effects (SW–ACD–X) . The primary innovation is the integration of the Secant-Weibull distribution as the innovation law, which allows the framework to capture non-monotonic intensity shapes, such as unimodal and bathtub patterns, that are typically inaccessible to standard monotonic models. A significant methodological contribution is the SW–ACD–X model, which endogenously incorporates intraday seasonality into the conditional mean equation. This joint estimation strategy provides an integrated alternative to traditional two-step pre-filtering by simultaneously capturing the interaction between deterministic diurnal patterns and stochastic duration clustering. The numerical properties of the proposed models are assessed through Monte Carlo simulations, which demonstrate asymptotic consistency while highlighting inherent identification challenges in small-sample regimes. Model estimation is implemented using a dual approach: Frequentist Maximum Likelihood and Bayesian Hamiltonian Monte Carlo (HMC) via the No-U-Turn Sampler (NUTS) in RStan. Empirical application to high-frequency IBM transaction data shows that the SW–ACD–X exhibits promising fit advantages over established benchmarks, including the W–ACD–X, LW–ACD–X, G-ACD-X , and Lomax–ACD–X models. Comprehensive model selection based on AIC, BIC, WAIC, and LOOIC confirms that the proposed model is a robust tool for analyzing market microstructure, liquidity dynamics, and the complex patterns of high-frequency durations.

In this paper, we introduce and study the class of modified \((p, h)\)-convex stochastic processes, which unifies and extends several existing notions of convexity in the stochastic setting. We establish fundamental arithmetic properties of this class and derive Hermite–Hadamard-type inequalities using classical and fractional Katugampola-type operators. We also investigated Ostrowski-type and Jensen-type inequalities in a simple and unified framework. Our results generalize and unify many existing results in the literature, providing a comprehensive framework for the study of convex stochastic processes.

For all positive non-square integer multipliers, there are infinitely many triangular numbers that are multiples of other triangular numbers. With a simple change of variables, one obtains a Pell equation, whose odd solutions provide the indices of the many infinitely triangular numbers multiple of other triangular numbers. General algebraic expressions of fundamental solutions of Pell equations are found for the multiplier expressed in function of the closest integer square. Finally, recurrent relations yielding the triangular numbers and their multiples and indices are calculated for non-square multipliers.

We introduce a class of mixed general bivariational inequalities in a real Hilbert space and show that several known models, including mixed variational inequalities, bivariational inequalities, general variational inequalities, and complementarity problems, arise as special cases. An auxiliary-principle framework is then used to derive predictor–corrector type iterative schemes. A basic descent estimate is established under a g-partially relaxed monotonicity assumption, and a convergence theorem is obtained under natural continuity and uniqueness hypotheses. The presentation has been streamlined to make the algorithmic steps explicit, and a scalar example is included to illustrate the resolvent formulation.

Let \(\eta\) be a fixed positive integer. Let \(S\) be a subset of \(\mathbb{Z}\), \(\star:S\times S\to \mathbb{Z}\) be a binary function, and \(\zeta_{\eta}:\{\xi\in \mathbb{Z}:\gcd(\xi,\eta)=1\}\to \{0,1\}\) be a function. For a simple connected graph \(G\) of order \(n\), a bijective function \(f:V(G)\to S\) (where \(|S|=n\)) is called an arithmetic cordial labeling modulo \(\eta\) under the arithmetic structure \(\langle S,\zeta_\eta,\star\rangle\) if the induced function \(f_\eta^*:E(G)\to \{0,1\}\), defined by \(f_\eta^*(ab)=1\) whenever \(\gcd(f(a)\star f(b),\eta)= 1\) and \(\zeta_\eta(f(a)\star f(b))=1\); otherwise, \(f_\eta^*(ab)=0\), satisfies the condition \(|e_{f_\eta^*}(0)-e_{f_\eta^*}(1)|\leq 1\), where \(e_{f_\eta^*}(i)\) is the number of edges with label \(i\) (\(i=0,1\)). In this paper, we explore the arithmetic cordial labeling of some graphs under conditions imposed on the function \(\zeta_\eta\). The graphs included are star graphs, ladder graphs, alternate cycle snake graphs, join graphs, corona graphs, and tensor product graphs.