The security of several full homomorphic encryption (FHE) schemes depends on the hardness of the approximate common divisor (ACD) problem. The analysis of attack and defense against the system is one of the frontiers of cryptography research. In this paper, the performance of existing algorithms, including orthogonal lattice, simultaneous diophantine approximation, multivariate polynomial and sample pre-processing are reviewed and analyzed for solving the ACD problem. Orthogonal lattice (OL) algorithms are divided into two categories (OL-$\land$ and OL-$\vee$) for the first time. And an improved algorithm of OL-$\vee$ is presented to solve the GACD problem. This new algorithm works well in polynomial time if the parameter satisfies certain conditions. Compared with Ding and Tao's OL algorithm, the lattice reduction algorithm is used only once, and when the error vector $\mathbf{r}$ is recovered in Ding et al.'s OL algorithm, the possible difference between the restored and the true value of $p$ is given. It is helpful to expand the scope of OL attacks.

The dynamics of the main photosynthetic structures is studied by mathematical modeling methods in this article. Chlorophyll portion variability in phytoplankton and the formation of energy-intensive substances in the process of photosynthesis underlie the models. These cellular components are considered in terms of their participation in the growth of specific biomass. Computational experiments are presented that simulate various degrees of environmental friendliness. The corresponding functions are built in accordance with seasonal fluctuations throughout the year in the Far East region of Russia. The stability of model solutions in long-term dynamics is also investigated. The models were tested for biological adequacy and their effectiveness was compared. For the model selected as a result of the comparison, the optimal control problem has been formulated and solved. It is shown that in this way it is possible to reduce the space of the initial components of the model system.

In this paper study, we develop a Bayesian Network model selection principle that address addresses the incommensurability of network features obtained from incongruous datasets and overcomes performance irregularities of the Minimum Description Length model selection principle. This is achieved (i) by approaching model evaluation as a classification problem, (ii) by estimating the effect that sampling error has on the satisfiability of conditional independence criterion, as reflected by Mutual Information, and (iii) by utilizing this error estimate to penalize uncertainty in the Minimum Uncertainty (MU) model selection principle. We validate our findings numerically and demonstrate the performance advantages of the MU criterion. Finally, we illustrate the advantages of the new model evaluation framework on a tRNA structural biology example.

This article develops duality principles and related convex dual formulations suitable for the local optimization of non-convex primal formulations for a large class of models in physics and engineering. The results are based on standard tools of functional analysis, calculus of variations and duality theory. In particular, we develop applications to a Ginzburg-Landau type equation.

Most studies modelling population mobility and the spread of infectious diseases, particularly using meta-population-multi-patched models, tend to focus on theoretical properties and numerical simulations of such models. There is relatively scanty literature published on fit, inference and uncertainty quantification on epidemic models with population mobility. In this research, we have used three estimation techniques to solve an inverse problem and quantify its uncertainty on a human mobility-based multi-patched epidemic model, using mobile phone sensing data and COVID-19 confirmed positive cases in Hermosillo, Mexico. First, we have utilized a Brownian bridge model using mobile phone GPS data to estimate residence and mobility parameters of the epidemic model. In the second step, we have estimated the optimal model epidemiological parameters by deterministically inverting the model using a Darwinian inspired evolutionary algorithm (EA) known as the genetic algorithm (GA). The third part of the analysis involves performing inference and uncertainty quantification on the epidemic model using two Bayesian Monte Carlo sampling methods: t-walk and Hamiltonian Monte Carlo (HMC). The results show that the estimated model parameters and incidence adequately fit the observed daily COVID-19 incidence in Hermosillo. Moreover, the estimated parameters from HMC result into large credible intervals, improving their coverage for the observed and predicted daily incidences. We also observe improved predictions when using multi-patch model with mobility against the single-patch model.

Given several nonnegative matrices with a single pattern of allocation among their zero/nonzero elements, the average matrix should have the same pattern, too. This is the first tenet of the pattern-multiplicative average (PMA) concept, while the second one suggests the multiplicative (or geometric) nature of averaging. The original concept of PMA was motivated by the practice of matrix population models as a tool to assess the population viability from long-term monitoring data. The task has reduced to searching for an approximate solution to an overdetermined system of polynomial equations for unknown elements of the average matrix (G), hence to a nonlinear constrained minimization problem for the matrix norm. Former practical solutions faced certain technical problems, which required sophisticated algorithms but returned acceptable estimates. Fresh idea to use the eigenvalue approximation ensues also from the basic equation of averaging, which determines the exact value of λ1(G), the dominant eigenvalue of matrix G, too, and the corre-sponding eigenvector. These are bound by the well-known linear equations, which reduce the task to a standard problem of linear programing (LP). The LP approach is realized for 13 fixed-pattern matrices gained in a case study of Androsace albana, an alpine short-lived perennial, monitored on permanent plots during 14 years. A standard software routine reveals the unique exact solution, rather than an approximate one, to the PMA problem, which turns the LP approach into ‘’the best of versatile optimization tools”. The exact solution turns out to be peculiar in reaching zero bounds for certain nonnegative entries of G, which deserves a modified problem formulation separating the lower bounds from zero.

Diabetes is a serious health problem throughout the world, including in Indonesia. The International Diabetes Federation (IDF) reports that the number of adults with diabetes is increasing every year. The Behavioral Risk Factor Surveillance System (BRFSS) is a survey conducted by the Centers for Disease Control and Prevention (CDC) in the United States. Classification methods in data mining techniques are used to classify diabetics and non-diabetics. The data mining process is carried out by preprocessing, feature selection, and dataset classification stages. In the preprocessing stage, data cleaning, data formatting, and data oversampling are carried out using the Synthetic Minority Over-sampling Technique (SMOTE). Next, the feature selection stage is carried out using the Particle Swarm Optimization (PSO) algorithm to find the best attributes. The dataset classification stage is carried out using the CART Model Decision Tree algorithm. The results of the performance evaluation of the CART algorithm are calculated using the confusion matrix and the MAE value, the results obtained for the CART algorithm without SMOTE and PSO obtained the best accuracy of 75.34% and the MAE value of 0.2466, while the CART algorithm using SMOTE and PSO can increase accuracy by 10 .94% to 86.28% and an MAE value of 0.1372.

In this paper, we extend the Merton model of investment in discrete time to the cases when there is a finite number of investors and the market is with frictions represented by convex penalty functions defined for each investor. In the main result of this paper, we proved the existence of optimal strategy of investment by using a new approach based on the formulation of an equivalent general equilibrium economy model.

The main objective of this study is to find out the influence of cooperation and intra-specific competition in escaping predation through refuge of the prey population and the effect of the two intra-specific interactions on the dynamics of prey-predator systems. For this purposes, two mathematical models with Holling type-II functional response function have been proposed and analyzed. The first model includes cooperation among prey populations whereas the second one incorporates intra-specific competition. Existence conditions and stability of different equilibrium points of both models have been carried out for the qualitative behaviour of the systems. Refuge through intra-specific competition has a stabilizing role, whereas cooperation has a destabilizing role on the system dynamics. Periodic oscillations are observed in both the systems through Hopf bifurcation. From the analytical and numerical findings, we conclude that intra-specific competition affects the prey population and continuously checks the refuge class under a critical value, and thus it never becomes too large to cause predator extinction due to food scarcity. Conversely, cooperation leads maximum individuals to escape predation through the refuge so that predators will suffer from low predation success.

This article develops duality principles and numerical results for a large class of non-convex variational models. The main results are based on fundamental tools of convex analysis, duality theory and calculus of variations. More specifically the approach is established for a class of non-convex functionals similar as those found in some models in phase transition. Finally, in the last section we present a concerning numerical example and the respective software.

The orthogonal matrix with cube be symmetric is a common class of matrices with important properties. We pay attention to this kind of matrices. By using the close relationship between the eigenvalues of a matrix and the trace of its power, we obtain the algorithm for its all possible different eigenvalues and multiplicities. The calculation formula is expressed only by the trace of a matrix and its power, avoiding solving the characteristic polynomial. The method is simple and practical. Furthermore, a new essential characterization for the sum of orthogonal matrix pairs being orthogonal is given as well.

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

Multicriteria decision aid requires a database as a decision matrix, in which two or more alternatives are evaluated according to two or more variables selected as decision criteria. Several problems of this nature use measures by Likert scales. Depending on the method, parameters from these data (e.g. means, modes or medians) are required for calculations. This parameterization of data in ordinal scales has fueled controversy for decades between authors who favor mathematical/statistical rigor and argue against the procedure, stating that ordinal scales should not be parameterized, and scientists from other areas who have shown gains from the process that compensate for this relaxation. The aim of this article is to demonstrate the advantages of the Composition of Probabilistic Preferences (CPP) method in multicriteria problems with data from Likert scales. The CPP is capable of allaying the protests raised and obtaining more accurate results than descriptive statistics or parametric models can bring. The proposed algorithm in R-code involves the use of CPP with empirical distributions and fitting histograms of data measured by Likert scales. Two case studies with simulated datasets having peculiar characteristics and a real case illustrate the advantages of the CPP.

The main goal of the paper is to use a generalized proportional Riemann-Liouville fractional derivative (GPRLFD) to model BAM neural networks and to study some stability properties of the equilibrium. Initially, several properties of the GPRLFD are proved such as the fractional derivative of a squared function. Also some comparison results for GPRLFD are provided. Two types of equilibrium of the BAM model with GPRLFD are defined. In connection with the applied fractional derivative and its singularity at the initial time the Mittag-Leffler exponential stability in time of the equilibrium is introduced and studied. An example is given illustrating the meaning of the equilibrium as well as its stability properties.

In this work, in spite of Milne’s recommendation using the three-point Newton–Cotes open formula (Milne’s rule) as a predictor rule and three-point Newton–Cotes closed formula (Simpson’s rule) as a corrector rule for 4-th differentiable functions with bounded derivatives. There is still a great need to introduce such formulas in other Lp spaces. Often, we need to approximate real integrals under the assumptions of the function involved. Because of that, the aim of this work is to introduce several Lp error estimates for the proposed perturbed Milne’s quadrature rule. Numerical experiments showing that our proposed quadrature rule is better than the classical Milne rule for certain types of functions are provided as well.

In this age of mass media and, in particular, social media-driven perception of reality, coupling disease and prophylactic opinion dynamics models can provide better insights into disease evolution than using a disease model alone. We develop in this work two disease-opinion dynamics models based on the epidemiology of the new coronavirus disease (COVID-19) and the availability or not of imperfect vaccines. We assume that susceptibility to infection decreases with the level of prophylactic attitude (personal hygiene, social distancing), and changes in prophylactic attitudes of susceptible individuals occur in response to perceived disease prevalence and vaccination coverage and efficacy in the population. We derive and discuss the disease-free equilibriums and reproduction numbers in the introduced models. We further assess the impacts of the distribution of opinions at disease introduction, the ability to detect presymptomatic, asymptomatic and symptomatic positive COVID-19 cases, the behavioural responses to the outbreak and the introduction of vaccination, and the effects of distortions of disease prevalence by public policy and mass media on disease dynamics. The insights highlighted from the proposed models are expected to make informative contributions to public policy in a context of opinion fluxes in response to perceived disease prevalence.

A problem involving the government's inspection of coins in Lower Slobbovia is discussed in the American Mathematical Monthly. It is desirable to classify and correct counterfeit coins by utilizing the fewest possible measures. However, previous article only proves the existence of the solution, which does not give a practical solution to the problem. In this article, we provide an effective method to achieve the identification of coin heaps using Dirac delta functions, Singularity functions (and their selectivity), unit-order construction of continuum mappings and construction of discretization of unique basis vectors without attempted expansions in S-base space.Finally, possible applications of this mathematical approach in gene mutation detection are also presented.

The study considers building a dynamic model of profit maximization for the car-sharing system and its verification on the example of implementing such a system in Astana (Republic of Kazakhstan). The region, bounded by the administrative boundaries of Astana, was divided into subregions by covering the region with regular hexagons placed side by side. A dataset was built with information on 1168 trips to Astana from January to March 2023. The Kepler visualization service constructed maps of the beginning and end of trips in the region and a map of trips binding to the hexagonal grid cells. Each cell of the grid corresponds to a specific sub-region, for which the quantitative parameters necessary for solving the profit maximization problem in the car-sharing system are calculated. Stations with cars will be placed in the cells of the grid, which will be available to car-sharing service customers. Based on the collected data, a dynamic (four periods per day) and static profit maximization model in the car-sharing system was built. Modeling was carried out based on the built models in the example of Astana. It was established that using a dynamic model of profit maximization in the car-sharing system increases the profit by 3.7%.

This paper is concerned with the optimal linear quadratic Gaussian (LQG) control problem for discrete time-varying system with multiple input delays and multiplicative noises. The main contributions are two-fold. Firstly, when the state variables can be observed exactly, we obtain a necessary and suffcient condition for the multiple-delays system in terms of the non-homogeneous relationship between the state and costate, which is the solution to the coupled forward and backward stochastic difference equations. Secondly, when the state variables are partially observed, we derive a suboptimal linear output feedback controller for the discrete-time system based on the obtained results of the optimal LQG control. Numerical examples are shown to illustrate the proposed algorithm.

In recent years, intensive work has been carried out to achieve maximum energy production in photovoltaic power plants. Today, this challenge has to be solved for arrays of solar panels that move with the object on which they are installed (vehicles, trucks, ships, etc.). The present work tries to find the best solution for both constant voltage systems and for cases of maximum intensity. Once found, these solutions will be applied in real time for the dynamic recharging of battery packs, trying to achieve energetically autonomous vehicles.To solve these situations, a new computational method of voltage and amperage calculation has been developed in this paper, based on Dijkstra’s minimum path finding algorithm on graph theory, adapted to electrical circuits. For this purpose, panel performance analysis sensors, developed at the University of Malaga, have been combined with different electronic solutions (wifi relay devices using esp8266 or similar chips) achieving an accurate and sufficiently fast solution at very low cost. With this, series or parallel configurations can be carried out depending on the type of energy generation sought. The theoretical solutions using Minkowski paths, achieved in the past, have been simulated and subsequently constructed in this article, indicating the diagrams necessary for their realisation.

Modern computational biology makes widespread use of mathematical models of biological systems, in particular systems of ordinary differential equations, as well as models of dynamic systems described in other formalisms, such as agent-based models. Parameters are numerical values of quantities reflecting certain properties of a modeled system and affecting model solu-tions. At the same time, depending on parameter values, different dynamic regimes - stationary, oscillatory, established as a result of transient modes of various types - can be implemented in the system. Predicting changes in the solution dynamics type depending on changes of model parameters is an important scientific task. Nevertheless, this problem does not have an analytical solution for all formalisms in a general case. The routinely used method of performing a series of computa-tional experiments, i.e. solving a series of direct problems with various sets of parameters fol-lowed by expert analysis of solution plots is labor-intensive with a large number of parameters and a decreasing step of the parametric grid. In this regard, the development of methods allow-ing to obtain and analyze information on a set of computational experiments in an aggregate form is relevant. This work is devoted to developing a method for visualization and classification of various dy-namic regimes of a model using a composition of the dynamic time warping (DTW-algorithm) and principal coordinates analysis (PCoA) methods. This method allows to obtain a qualitative visualisation of the results of the set of solutions of a mathematical model and to carry out the correspondence between the values of the model parameters and the type of dynamic regimes of its solutions. This method has been tested on the Lotka-Volterra model and artificial sets of var-ious dynamics

Fluctuating stock prices make it difficult for investors to see investment opportunities. One tool that can help investors overcome this is forecasting techniques. Long Short-Term Memory (LSTM) is one of deep learning methods used in forecasting time series. The training and success of deep learning is strongly influenced by the selection of hyperparameters. This research uses a hybrid method between the Genetic Algorithm (GA) and LSTM to find a suitable model for predicting stock prices. GA is used in optimizing the architecture such as the number of epochs, window size, and the number of LSTM units in the hidden layer. Tuning optimizer is also carried out using several optimizers to achieve the best value. From method that has been applied, it shows that the method has a good level of accuracy with MAPE values below 10% in every optimizer used. The error rate generated is quite low, in case-1 with a minimum RMSE value of 93.03 and 94.40, & in case-2 with an RMSE value of 104.99 and 150.06 during training and testing. A fairly stable and small value is generated by setting it using the Adam Optimizer.

In this article, the multitouch option, also called the n  touch option (or the “baseball” option when n  3 ) is analyzed and valued in closed form. This is a kind of barrier option that sets a gradual knockout / knock-in mechanism based on the number of times the underlying asset has crossed a predefined barrier in various time intervals before expiry. The higher the number of predefined time intervals during which the barrier has been touched, the lower the value of a knock-out contract at expiry, and conversely for a knock-in one. Multitouch options can be viewed as an extension of step barrier options, preserving the ability of the latter to adjust the exposure to risk over time, while eliminating the notorious danger of “sudden death” that holders of step barrier options are faced with. Unlike occupation time derivatives, the payoff at expiry does not depend on the amount of time spent outside the authorized range, but on the number of passages beyond the authorized range.

More than fifty million people over the age of fifty have osteoporosis issues due to low bone mass, disease contribution, and not exercising continuously. Certain most important risk factors and causes of osteoporosis contain age, gender, hormones, visits to inexperienced doctors, and certain medical conditions. In this manuscript, we diagnosed certain frank operational laws under the availability of complex intuitionistic fuzzy (CIF) information. Furthermore, using the evaluated operational laws, we pioneered the theory of CIF frank weighted averaging (CIFFWA), CIF frank ordered weighted averaging (CIFFOWA), CIF frank hybrid averaging (CIFFHA), CIF frank weighted geometric (CIFFWG), CIF frank ordered weighted geometric (CIFFOWG), CIF frank hybrid geometric (CIFFHG) operators, and described their beneficial and valuable results with certain useful properties. Moreover, we aim to evaluate the most dangerous type of osteoporosis based on their symptoms, causes, and risk factors using diagnosed approaches. Finally, to enhance the worth of the evaluated operators, we compare the diagnosed result with certain prevailing results and illustrated their geometrical representation with the help of MATLAB.

In this work, an analytical and numerical analysis of the transition to chaos in five nonlinear systems of ordinary and partial differential equations, which are models of the autocatalytic chemical processes and the numbers of interacting populations, is carried out. It is analytically and numerically shown that in all considered systems of equations, further complication of the dynamics of solutions and the transition to chemical and biological turbulence is carried out in full accordance with the universal Feigenbaum-Sharkovsky-Magnitskii bifurcation theory through subharmonic and homoclinic cascades of bifurcations of stable limit cycles. In this case, irregular (chaotic) attractors in all cases are exclusively singular attractors in the sense of the FShM theory.

Grid refinement is used to save computing costs while preserving high precision of fluid simulation. In the Lattice Boltzmann Method, grid refinement often uses interpolated values. We found a method, where interpolation in space and time is not required. For this purpose, we used the moment matching condition and rescaling the non-equilibrium part of populations, and developed the recalibration procedure that allows to transfer information between different LBM stencils in the simulation domain. Then, we build non-uniform lattice which uses stencils with different shapes on the transition. The resulting procedure is verified on the 2D Poisselle flow and the advected vortex benchmark. In prospect, grids with adaptive geometry can be built with the use of the proposed method.

The option pricing problem is equivalent with the hedging problem of the option, i.e. what the writer should do in order to hedge the risk that she undertakes selling a contract and moreover what is the probability of profit selling at a specific price. In other words, the writer first designs a way to hedge the risk and then prices the contract according to the amount needed to hedge that risk. The probability of profit is also a useful information for the buyer. The hedging strategy should be practically possible for the writer otherwise has no meaning. In this note we will discuss the option pricing problem and in particular the effect of the volatility on the binomial model which is a way to hedge practically a specific option in contrast to every pricing model that assumes rebuilding of the replicating portfolio continuously in time. In order to use the binomial model we have to modified it accordingly as we have seen in a previous paper. We also point out four open problems regarding the hedging problem.

This paper presents a concise mathematical framework for investigating both feed-forward and backward process, during the training to learn model weights, of an artificial neural network (ANN). Inspired from the idea of the two-step rule for backpropagation, we define a notion of F_adjoint which is aimed at a better description of the backpropagation algorithm. In particular, by introducing the notions of F-propagation and F-adjoint through a deep neural network architecture, the backpropagation associated to a cost/loss function is proven to be completely characterized by the F-adjoint of the corresponding F-propagation relatively to the partial derivative, with respect to the inputs, of the cost function.

The main goal of this work is to develop a novel boundary element method (BEM) model for analyzing ultrasonic wave propaga-tion in three-temperature anisotropic viscoelastic porous media. Because of strong nonlinearity of ultrasonic wave propagation in three-temperature porous media problems, the analytical or numerical solutions to the problems under consideration are always challenging, necessitating the development of new computational techniques. As a result, we use a new BEM model to solve such problems. A time stepping procedure based on the linear multistep method is obtained after solving the discretized boundary in-tegral equation with the quadrature rule. The calculation of a double integral is required to obtain fundamental solutions, but this increases the total BEM computation time. Our proposed BEM technique is used to solve the current problem and improve the formulation efficiency. The numerical results are graphed to demonstrate the effects of viscosity and anisotropy on the nonlinear ultrasonic stress waves in three-temperature porous media. The validity, accuracy, and efficiency in the proposed methodology were demonstrated by comparing obtained results to the corresponding solution obtained from finite difference method (FDM).

The influence of neutrosophy in the previous period is constantly growing in many areas of science and technology. Moreover, various applications of the neutrosophic approach have become more common in recent years. Our goal in this research is to utilize the neutrosophy to improve the performance of the Dai-Liao conjugate gradient (CG) method. Specifically, in this research, we propose and investigate a new neutrosophic logic system to calculate the key parameter t involved in the Dai–Liao CG iterations. Theoretical analysis and numerical experience indicate that the efficiency and robustness of the new rule for determining t. Combining the neutrosophy and the Dai-Liao conjugate gradient method, we propose and explore a new Dai-Liao CG iterations for solving large-scale unconstrained optimization models. The global convergence is established under common assumptions and the backtracking line search. Finally, by conducting numerical experiments, computational evidence demonstrates that the new fuzzy neutrosophic Dai-Liao conjugate gradient method is computationally effective and robust.

Positive invariant set is an important concept of dynamic systems. The purpose of this paper is to study the sufficient and necessary conditions that the set of ellipsoids and the Lorenz cone are positive invariant sets of discrete time systems. By means of nonlinear programming and induced norm, the problem of positive invariance is formulated as an optimization problem, and the equivalent dual form optimization is also presented using the dual optimization method. Our results provide more alternative methods for determining the positive invariance of quadratic form convex sets from the point of view of optimization and dual optimization. The effectiveness of this method is demonstrated by numerical examples.

Theoretical approaches were applied to study the effect of magnetic field and heat transfer on the flow of blood plasma through an asymmetric arterial segment. The plasma was considered to be unsteady, laminar and an incompressible fluid through non-uniform arterial segment in a two-dimensional flow. Axial velocity, pressure gradient, stream function and pressure shift per wavelength were evaluated in blood plasma flow in arteries by use of coupled linear partial differential equations, solved with the help of Finite difference method. The corresponding initial and boundary conditions were obtained by discretizing the flow channel. The effect of magnetic field on blood plasma in arteries was determined by creating a magnetic field gradient through the application of a varying strength of magnetic field on the flow. Heat transfer characteristics due to applied magnetic field and viscosity of blood was obtained by use of linear partial differential equations to determine varying temperature conditions in heat transfer characteristics. Numerical results for velocity profiles, magnetic profiles and temperature profiles were obtained to characterize blood plasma viscosity by use of various non-dimensional parameters. The results were graphically presented using MATLAB software. The results obtained helped in analyzing theoretically the effects of magnetic field and heat transfer in arterial plasma flow.

Optimization of Magnetohydrodynamics parameters on the velocity profile and temperature distribution of incompressible fluid flow on a porous channel was evaluated. The fluid flow was considered to be unsteady, incompressible and flowing in 2-D in porous channel. The effects of magnetic parameter, Darcy’s number and fluid pressure on velocity profiles, and the effect of variable viscosity and Eckert number on temperature distribution in incompressible fluid flow were determined. The flow was considered to be in a channel running along the x-axis on which the magnetic field exist and finally along the y-axis on the porous part of the channel. The resulting model equations were solved by Finite Difference Method (FDM) in MATLAB software. Analysis of results indicated that increasing Darcy’s number and fluid pressure leads to increase in fluid velocity profile, while increasing magnetic parameter decreases fluid velocity profile. Also, it was observed that increase in both Eckert number and fluid viscosity lead to increase in temperature distribution. Optimization in temperature was achieved by increasing the magnetic field while viscosity was optimized by increasing the length of the porous part of the channel. This study will helpful to contribute alternative equations and methodology to engineering and in factories where getting the MHD parameters optimally is the main objective, particularly on temperature, velocity and pressure.

Hessenberg differential algebraic equations (Hessenberg-DAEs) with high index play a critical role in the modeling of mechanical systems and multibody dynamics. Motivated by the widely used Lie Group Differential Algebraic Equation (LGDAE) method which only handles index 2 systems, we propose a Modified Extended Lie Group Differential Algebraic Equation (MELGDAE) method for solving index 3 Hessenberg-DAEs, and provide theoretical analysis to deepen the foundation of the MELGDAE method.The performance of the MELGDAE method is compared with the standard methods RADAU and MEBDF on index 2 and 3 DAE systems, and it is demonstrated that the MELGDAE integrator exhibits competitive performance in terms of high accuracy and the preservation of algebraic constraints. In particular, all differential variables in index 3 Hessenberg DAEs achieve second-order convergence using the MELGDAE method, which suggests potential for extension to Hessenberg-DAEs with an index of 4 or higher.

Remanufacturing is a way for product recovery initiative and maintaining sufficient production flow to satisfy customer demand through circular economy. Recently, manufacturers have been progressively incorporating remanufacturing processes, making their supply chain vulnerable to disruptions. Among the disruptions that usually occur is the shortage in spare parts supply, which causes unexpected delay in production and possible loss of sales. In the face of potential disruption, it is essential for remanufacturing facilities to manage their supply chains in an optimum manner as to decrease the adverse impact of disruptions to their business. In this study, a two-stage production-inventory system is analysed by developing a cost-minimisation model focusing on the recovery schedule after a disruption occurrence in sourcing spare parts for a remanufacturer’s production cycle. The developed model was solved using the branch and bound algorithm. Through numerical experiments, the results indicate that the optimal recovery schedule and the number of recovery cycles are significantly dependant on disruption time, lost sales and backorder costs. Moreover, a sensitivity analysis shows that the lost sales option seems to be more effective than the backorder option in optimising the system’s overall cost. The proposed model could assist managers in deciding the optimal production strategy while providing interesting managerial insights into vital spare parts recovery issues when disruption strikes.

Cancer diseases lead to the second-highest death rate all over the world. The dynamics of invasion of cancer cells into the human body tissues and metastasis are the main causes of death in patients with cancer. This study deals with theoretical investigation of the dynamics of invasion of cancer cells for tumour growths in human body tissues using discretized Cahn-Hilliard, concentration and reaction-diffusion equations which were solved by Finite Difference Method with the aid of MATLAB computer software. A Crank-Nicolson numerical scheme was developed for the discretized model equations. The numerical result obtained was used to describe the dynamics of cancer invasion of tissues with respect to cancer cells density on tumour growth, turbulence and mobility and equilibrium between charge and discharge of cancer cells. The results of the study provide new insights into combating cancer disease by providing mitigating and intervention measures to this major health problem.

A new low-pass filter (LPF) circuit with independent adjustment of various, incl. digitally controlled, pole frequency resistors, pole quality factor and transmission coefficient. The proposed low-pass filter is based on the use of the properties of a multi-differential operational amplifier that performs the functions of a signal adder. The peculiarity of the filter is that it has three inputs, with respect to which different transmission coefficients are implemented, incl. inverting (-1) and non-inverting (+1). To check these properties of a low-pass filter in the Micro-Cap environment, computer simulation of a specific circuit was performed on a multi-differential operational amplifier AD 8130. Mathematical expressions are given for the main parameters of the proposed low-pass filter, which allow parametric synthesis of elements of a specific circuit under given restrictions on the used element base.

This article develops applications of the generalized method of lines to numerical solutions of the time-independent, incompressible Navier-Stokes system in fluid mechanics. We recall that for such a method, the domain of the partial differential equation in question is discretized in lines (or more generally in curves), and the concerning solutions are written on these lines as functions of the boundary conditions and the domain boundary shape.

Background: Articles published in scholarly journals form part of the scientific evidence base. It is the responsibility of the scientific community to maintain its integrity. In 2011 the BMJ commissioned a feature article to draw attention to an article that had appeared in another journal- The Lancet 13 years previously. The Lancet had already retracted the article. These actions exemplify the best traditions of scientific record-keeping. Objective: This submission examines whether the main claims summary made in the BMJ were factual. Method: We examine what was published in the Lancet against what was published in the BMJ and verify against the findings in the GMC hearings transcripts and verdict of the UK High Court. Results: The 6 points highlighted in BMJ had errors and need to be corrected. Conclusions: There are significant differences between what was reported in the Lancet paper and what was alleged to be there by the BMJ. This article aims only to point to errors in the BMJ article, to set the record straight. It does not show there was a causal association between MMR vaccination and autism.

In recent years, there has been a growing interest in the preoperative modeling of fluid-structure interaction in the treatment of cerebral aneurysms. In this study, we investigate two cases involving laminar incompressible fluid flow interacting with hyperelastic materials (a vertical flap and aneurysm walls) under the effect of fluid flow. We present a finite element method (FEM) for these prototypical two-dimensional ($2D$) configurations, taking into account the complex flows and deformation of the variant's structure models. We use an Arbitrary Lagrangian-Eulerian (ALE) formulation in a continuum, fully monolithic coupled way, and discretize the fluid and solid domains using the quadratic LBB stable $P_2P_1$ finite element pair to approximate the displacement, velocity, and pressure spaces independently. The resulting discretized form of the nonlinear algebraic system is linearized using a variant of Newton's procedure, and the Jacobian matrices are approximated via the divided difference method. The resulting linear systems are solved using a direct solver MUMPS. We evaluate hydrodynamic forces such as drag and lift coefficients for nonlinear elastic material models, including Saint-Venant Kirchoff, Neo-Hookean, and Mooney-Rivlin (2 parameters) separately. To gain more physical insight into the problem, we verify the computed results by comparing the velocity, viscosity, and pressure fields. Our aim is to qualitatively analyze the changes in the mechanics of the wall behavior of the elasticity of the flap vs. the complexity of deformation, interface, details for prototypical flow situations, corresponding displacement, and deformation. We also provide a physical interpretation of blood flow magnitude and basic hemodynamic phenomena of wall shear stresses (WSS) in aneurysma and the range of deformation. Overall, this study presents a comprehensive approach to modeling fluid-structure interaction in cerebral aneurysms, which could help in the development of more effective treatment strategies.

The information data accuracy affects the working reliability of production management system directly. Mobile equipment is an important part of integrated manufacturing system, in the course of its works, due to the complex working environment and the dynamic change of its position, it lead to the important problem of information detection distortion. Because the application of fractional differential operator has the dual functions of enhancing signal strength and reducing data error, we puts forward a data processing method of obtaining production information data of various accuracy by changing the parameters of fractional differential operator. Establish a mobile equipment detection data fusion processing model based on fractional order partial differential equations and apply it in the processing experiment of mobile equipment testing information detection data in an integrated manufacturing system, the feasibility and effectiveness of the method are verified. We conclude that the fractional order partial differential equations used in the processing of mobile equipment monitoring information data of an integrated manufacturing system has the function of obtaining production information data of various accuracy, which are of great significance to improve the working reliability in integrated manufacturing system. which is of great significance to improve the working reliability and scientific decision of integrated manufacturing system.

The buildup of plaque in the arteries characterizes atherosclerosis, which causes the walls of the arteries to thicken, the lumen to narrow, and the wall to thin in certain areas. These changes can lead to alterations in blood flow, potentially resulting in aneurysms and heart attacks if left untreated. This paper presents a phenomenological model to explain the mechanics of plaque rupture in stenosed bifurcated elastic arteries. The model considers the interaction between the plaque and artery wall, blood flow, mechanical properties of the artery wall and plaque, and hemodynamic forces in the system. Using the Navier-Stokes equations to describe blood flow and elastic properties of artery walls, our study shows that blood flow can become turbulent, leading to backflow, vortices, and possible stagnation. Certain regions can become highly vulnerable and result in elevated heat transfer between blood and arterial walls, which can lead to the rupture of the plaque cap. The study focuses on blood flow features such as velocity profiles and wall displacement on fluid-structure interaction, which are consistent with the literature. Finally, we calculate the wall shear stress (WSS) for minimum and maximum times while considering elastic walls. Our findings may provide valuable insights into the mechanisms of plaque rupture and inform the development of improved diagnostic and therapeutic approaches.

We consider the problem of scheduling tasks on a single machine subject to an uncertain breakdown to minimize the flowtime. Assume the machine is unavailable during the breakdown; the starting time of the breakdown is a random variable s with distribution function D(s) and the terminating time of the breakdown has no any other information; jobs are non-resumable. Under the assumptions, starting from the perspective of statistical optimization, we first define the expected competitive ratio of an algorithm to find the optimized solution with the considered problem. Further, we propose and prove a certain result on the expected competitive ratio of the SPT rule. In particular, we show that under a quite loose condition the expected competitive ratio of the SPT rule is no more than 5/4. Meanwhile, we also make some discussions about our studies. What we have done will rich and improve the studying results on the area of scheduling to minimize flowtime and advance the development of online optimization and statistical optimization.

The quality of underground space information has become a major problem that endangers the safety of underground spaces. Currently, the main methods for the high-precision and long-distance transmission of detection information are radar and optical methods. However, in practical applications, we found that the radar method has the shortcomings of large energy loss and poor anti-jamming ability, which limit the accuracy of information data transmission and distance. The optical method has the shortcomings that the weather has a great impact on its accuracy and can only be applied to static objects above ground; therefore, it has the limitation of application objects and use environment. More importantly, the current high-precision information remote detection methods are limited to the detection of overground space objects and are not applicable to the detection of various information data in underground space. In this study, we analyze the spectral properties of the fractional differential operator and find that it is suitable for studying non-linear, non-causal, and non-stationary signals. The theory of fractional calculus is applied to the field of data processing, and a mathematical model of remote transmission and high-precision detection of information based on fractional difference is established, which realizes the functions of high-precision and remote detection of information. By fusing the information data to detect the mathematical model over a long distance and with high accuracy, a mathematical model for stratum data processing used to provide long-distance and high-accuracy data was established. Through application in engineering practice, the effectiveness of this method for underground space information data detection was verified.

In this paper, an analytical form of the Unit Step Function (or Heaviside Step function) is presented. This important function constitutes a fundamental concept of Operational Calculus and is also involved in many other fields of applied and engineering mathematics. In particular, Heaviside Step Function is performed in a very simple manner by the use of a finite number of standard operations, as the summation of six inverse tangent functions. The novelty of this work when compared with other analytical expressions, is that the proposed exact formula contains two arbitrary single - valued continuous functions which satisfy only one restriction. In addition, the proposed representation is not exhibited in terms of miscellaneous special functions, e.g. Bessel functions, Beta function, Logistic function. Besides, it is neither the limit of a function, nor the limit of a sequence of functions with point – wise or uniform convergence. Hence, this formula may be much more practical, flexible and useful in the computational procedures which are inserted into Operational Calculus techniques and other engineering practices.

Malaria is a critical fevered illness caused by Plasmodium parasites transmitted among people through the bites of infected female Anopheles mosquitoes. Public awareness about the disease is important for the control of disease. This article proposes a mathematical model to study the dynamics of malaria disease transmission with the influence of awareness-based control interventions. We found two equilibria of the model, namely the disease-free and endemic equilibrium. Disease-free equilibrium is stable globally if basic reproduction number (R0) is less than unity (R0<1). Some basic mathematical properties of the proposed model, such as nonnegativity and boundedness of solutions, the feasibility of the equilibrium points and their stability properties, have been studied. Finally, we adopted optimal control to minimize the cost of disease control and solve the problem by formulating Hamiltonian functional. The optimal use of insecticides for controlling the mosquito population is determined. Numerical simulations have been provided for the confirmation of the analytical results. We conclude that media awareness with optimal control approach is best for cost-effective malaria disease management.

Animal models of diseases in particular, the mice, are considered to be the cornerstone for translational research in immunology. The aim of the present study is to model the geometry and analyze the network structure of murine lymphatic system (LS). The algorithm for building the graph model of the LS makes use of anatomical data. To identify the edge directions of the graph model, a mass balance approach to lymph dynamics based on Hagen-Poiseuille equation is applied. It is the first study in which a geometric model of the murine LS has been developed and characterized in terms of its structural organization and the lymph transfer function. Our study meets the demand for quantitative mechanistic approaches in the growing field of immunoengineering to utilize or exploit the lymphatic system for immunotherapy.

In engineering practice, various types of information data are affected by many factors during the collection process. For example, information data measurement errors are caused by equipment performance and the working environment. During the transmission of detection information, the signal distortion caused by energy loss and signal interference causes unpredictable detection errors in the collected information data. Through the study of fractional calculus theory, it was found that it is suitable for studying nonlinear, non-causal, and non-stationary signals, and has the dual functions of improving detection information and enhancing signal strength. Therefore, under the influence of many factors, we applied the fractional difference algorithm to the field of information data processing. Multi-sensor detection data fusion algorithm based on the fractional partial differential equation was adopted to establish its online detection data. A multi-sensor detection data fusion algorithm based on a fractional partial differential equation is established, which effectively fuses the information data detection errors caused by various influencing factors and greatly improves the detection accuracy of information data. The effectiveness of this method is proved through its application in an experiment.

The operation status of manufacturing equipment is directly related to the reliability of the operation of manufacturing equipment and the continuity of operation of the production system. Based on the analysis of the operation status of manufacturing equipment and its characteristics, it is proposed that the concept of assessing the operation status of manufacturing equipment can be realized by applying real-time acquisition of accurate inspection data of important parts of weak motion units and comparing them with their motion status evaluation criteria. Through the study of the application characteristics of fractional-order calculus theory, a differential data fusion model based on the fractional-order differential operator is established. The advantages of Internet of Things (IoT）technology and fractional order differential fusion algorithm are integrated to obtain real-time high-precision data of the operating parameters of manufacturing equipment, and the research objective of operating condition assessment of manufacturing equipment is realized. The feasibility and effectiveness of the method are verified by the application of the method in machining center operation status assessment.

Symbolic pole/zero analysis is an important step when designing an analog operational amplifier. Generally, a simplified symbolic analysis of analog circuits suffers from NP-hardness, i.e., an exponential growth of the number of symbolic terms of the transfer function with the circuit size. In this study, we present a mathematical model combined with a heuristic-metaheuristic solution method for the symbolic pole/zero simplification in operational transconductance amplifiers (OTA). At first, the circuit is symbolically solved and an improved root splitting method is applied to extract symbolic poles/zeroes from the exact expanded transfer function. Then, a hybrid algorithm based on heuristic information and a metaheuristic technique using simulated annealing is performed for the simplification of the derived symbolic pole/zero expressions. The developed method has been tested on three analog OTAs. The obtained results show the effectiveness of the proposed method to achieve accurate simplified symbolic pole/zero expressions with the least complexity.

In this work we derive a generalized series expansion of the acrtangent function by using the enhanced midpoint integration (EMI). Algorithmic implementation of the generalized series expansion utilizes a simple two-step iteration. This approach significantly improves the convergence and requires no surd numbers in computation of the arctangent function.

This article develops dual variational formulations for a large class of models in variational optimization. The results are established through basic tools of functional analysis, convex analysis and duality theory. The main duality principle is developed as an application to a Ginzburg-Landau type system in superconductivity in the absence of a magnetic field. In the first sections, we develop new general dual convex variational formulations, more specifically, dual formulations with a large region of convexity around the critical points which are suitable for the non-convex optimization for a large class of models in physics and engineering. Finally, in the last section we present some numerical results concerning the generalized method of lines applied to a Ginzburg-Landau type equation.

This paper presents PrognosEase; a software that provides an easier way to produce different types of run-to-failure data mimicking real-world conditions to simplify prognosis studies in terms of data collection and improvement in ML degradation modelling process. Different types of degradation types made available to meet different types of applications. Besides, some primary ML tests were performed to ensure that complexity patterns of real systems could be observed in the training/testing predictions attitude. This paper also presents the impacts, limitations and potential improvements of the data generator.

In this paper we give approximations to the Mittag-Leffler functions in terms of elementary functions using different methods. This allowed us to establish a practical method we called integerization principle. This principle states that many fractional nonlinear oscillators may be solved by means of the solution to some integer-order Duffing oscillator equation. The accuracy of the obtained results is illsutrated in concrete examples. Formulas for estimating the errors in the approximations are also provided.

The adaptive identification method developed to evaluate the parameters of the Bouc-Wen hysteresis (BWH). The adaptive approach based on the use of adaptive observers. We synthesize adaptive identification algorithms using the second Lyapunov method. The boundedness of adaptive system processes shown in coordinate and parametric spaces. We prove exponential dissipativity of processes in an adaptive system. Estimating method proposed for signaling uncertainty in the system.

This article develops a duality principle and a related convex dual formulation suitable for a large class of models in physics and engineering. The results are based on standard tools of functional analysis, calculus of variations and duality theory. In particular, we develop applications to a model in non-linear elasticity.

The system identification problem with multiple nonlinearities is relevant. Its decision depends on many factors. These include: feedbacks, the method of connecting nonlinear links, signal properties. They affect the identifiability of the system parameters. We introduced a condition for the excitation constancy for state variables, which considers the S-identifiability of the system. We propose system decomposition by measuring input to identify parameters. Each subsystem has an implicit identification representation. It guarantees obtaining estimates of subsystem parameters based on experimental data. The trajectories boundedness of adaptive system proved in parametric and coordinate spaces. Conditions guaranteeing exponential stability of the system obtained. Systems of self-oscillation generation and nonlinear correction of a nonlinear system consider. Conditions for the trajectories boundedness of the adaptive system obtained for these cases. The influence of nonlinearity and feedback on the system performance estimated.

(1) Background: As Patient-reported outcome measures face challenges with low response rate on surveys, different incentive mechanism have been proposed to achieve a higher response rate. However, it seems that monetary incentives are the only mechanism with proven effect. Nevertheless, it is less likely to be used than other mechanisms due to the monetary cost. (2) Methods: In this research work, a cryptographic scheme for rewarding patients with cryptographic tokens is developed and implemented on the Ethereum test network. (3) Results: The model is able to distribute decentralised tokens to patients who complete the PROMs in a fair, private, decentralised approach. The token can be further used by patients to exchange more healthcare services, encouraging more patients to participate in PROMs. At the same time, an IER detection method is built to improve the quality of PROMs, avoiding the healthcare provider paying for meaningless PROMs from patients who barely participate for incentives. (4) Conclusions: This work provides an privacy-preserving incentive model to increase the response rate for PROMs surveys. Our model prevents patients providing invalid responses to gain rewards.

In this paper we study a conjecture proposed by P.Duxbury , C.laror , L.Leduino de Salles Neto in 2021\cite{conjecture} on the Golumb Ruler Problem which is a classical optimization model in discrete case . In \cite{conjecture} the authors constucted a continuous model for the Golumb Ruler Problem associated to the discrete case and conjectured that the optimal value of both cases are equal . We deal with this conjecture via algebraic geometrical methods .

In this research, new setting is introduced for new SuperHyperNotions, namely, SuperHyperDefensive SuperHyperAlliances and Neutrosophic SuperHyperDefensive SuperHyperAlliances. Two different types of SuperHyperDefinitions are debut for them but the research goes further and the SuperHyperNotion, SuperHyperUniform, and SuperHyperClass based on that are well-defined and well-reviewed. The literature review is implemented in the whole of this research. For shining the elegancy and the significancy of this research, the comparison between this SuperHyperNotion with other SuperHyperNotions and fundamental SuperHyperNumbers are featured. The definitions are followed by the examples and the instances thus the clarifications are driven with different tools. The applications are figured out to make sense about the theoretical aspect of this ongoing research. The “Cancer’s Recognitions” are the under research to figure out the challenges make sense about ongoing and upcoming research. The special case is up. The cells are viewed in the deemed ways. There are different types of them. Some of them are individuals and some of them are well-modeled by the group of cells. These types are all officially called “SuperHyperVertex” but the relations amid them all officially called “SuperHyperEdge”. The frameworks “SuperHyperGraph” and “neutrosophic SuperHyperGraph” are chosen and elected to research about “Cancer’s Recognitions”. Thus these complex and dense SuperHyperModels open up some avenues to research on theoretical segments and “Cancer’s Recognitions”. Some avenues are posed to pursue this research. It’s also officially collected in the form of some questions and some problems. Assume a SuperHyperGraph. An “SuperHyperAlliance” is a minimal SuperHyperSet of SuperHyperVertices with minimum cardinality such that either of the following expressions hold for the cardinalities of SuperHyperNeighbors of s∈S:,|S∩N(s)|>|S∩(V \N(s))|,and|S∩N(s)|<|S∩(V \N(s))|.Thefirst Expression, holds if S is SuperHyperOffensive. And the second Expression, holds if S is “SuperHyperDefensive”. It’s useful to define “neutrosophic” version of SuperHyperDefensive SuperHyperAlliances. Since there’s more ways to get type-results to make SuperHyperDefensive SuperHyperAlliances more understandable. For the sake of having neutrosophic SuperHyperDefensive SuperHyperAlliances, there’s a need to “redefine” the notion of “SuperHyperDefensive SuperHyperAlliances”. The SuperHyperVertices and the SuperHyperEdges are assigned by the labels from the letters of the alphabets. In this procedure, there’s the usage of the position of labels to assign to the values. Assume a SuperHyperAlliance. It’s redefined neutrosophic SuperHyperAlliance if the mentioned Table holds, concerning, “The Values of Vertices, SuperVertices, Edges, HyperEdges, and SuperHyperEdges Belong to The Neutrosophic SuperHyperGraph” with the key points, “The Values of The Vertices & The Number of Position in Alphabet”, “The Values of The SuperVertices&The Minimum Values of Its Vertices”, “The Values of The Edges&The Minimum Values of Its Vertices”, “The Values of The HyperEdges&The Minimum Values of Its Vertices”, “The Values of The SuperHyperEdges&The Minimum Values of Its Endpoints”. To get structural examples and instances, I’m going to introduce the next SuperHyperClass of SuperHyperGraph based on SuperHyperDefensive SuperHyperAlliances. It’s the main. It’ll be disciplinary to have the foundation of previous definition in the kind of SuperHyperClass. If there’s a need to have all SuperHyperConnectivities until the SuperHyperDefensive SuperHyperAlliances, then it’s officially called “SuperHyperDefensive SuperHyperAlliances” but otherwise, it isn’t SuperHyperDefensive SuperHyperAlliances. There are some instances about the clarifications for the main definition titled “SuperHyperDefensive SuperHyperAlliances”. These two examples get more scrutiny and discernment since there are characterized in the disciplinary ways of the SuperHyperClass based on SuperHyperDefensive SuperHyperAlliances. For the sake of having neutrosophic SuperHyperDefensive SuperHyperAlliances, there’s a need to “redefine” the notion of “neutrosophic SuperHyperDefensive SuperHyperAlliances” and “neutrosophic SuperHyperDefensive SuperHyperAlliances”. The SuperHyperVertices and the SuperHyperEdges are assigned by the labels from the letters of the alphabets. In this procedure, there’s the usage of the position of labels to assign to the values. Assume a neutrosophic SuperHyperGraph. It’s redefined “neutrosophic SuperHyperGraph” if the intended Table holds. And SuperHyperDefensive SuperHyperAlliances are redefined “neutrosophic SuperHyperDefensive SuperHyperAlliances” if the intended Table holds. It’s useful to define “neutrosophic” version of SuperHyperClasses. Since there’s more ways to get neutrosophic type-results to make neutrosophic SuperHyperDefensive SuperHyperAlliances more understandable. Assume a neutrosophic SuperHyperGraph. There are some neutrosophic SuperHyperClasses if the intended Table holds. Thus SuperHyperPath, SuperHyperCycle, SuperHyperStar, SuperHyperBipartite, SuperHyperMultiPartite, and SuperHyperWheel, are “neutrosophic SuperHyperPath”, “neutrosophic SuperHyperCycle”, “neutrosophic SuperHyperStar”, “neutrosophic SuperHyperBipartite”, “neutrosophic SuperHyperMultiPartite”, and “neutrosophic SuperHyperWheel” if the intended Table holds. A SuperHyperGraph has “neutrosophic SuperHyperDefensive SuperHyperAlliances” where it’s the strongest [the maximum neutrosophic value from all SuperHyperDefensive SuperHyperAlliances amid the maximum value amid all SuperHyperVertices from a SuperHyperDefensive SuperHyperAlliances.] SuperHyperDefensive SuperHyperAlliances. A graph is SuperHyperUniform if it’s SuperHyperGraph and the number of elements of SuperHyperEdges are the same. Assume a neutrosophic SuperHyperGraph. There are some SuperHyperClasses as follows. It’s SuperHyperPath if it’s only one SuperVertex as intersection amid two given SuperHyperEdges with two exceptions; it’s SuperHyperCycle if it’s only one SuperVertex as intersection amid two given SuperHyperEdges; it’s SuperHyperStar it’s only one SuperVertex as intersection amid all SuperHyperEdges; it’s SuperHyperBipartite it’s only one SuperVertex as intersection amid two given SuperHyperEdges and these SuperVertices, forming two separate sets, has no SuperHyperEdge in common; it’s SuperHyperMultiPartite it’s only one SuperVertex as intersection amid two given SuperHyperEdges and these SuperVertices, forming multi separate sets, has no SuperHyperEdge in common; it’s SuperHyperWheel if it’s only one SuperVertex as intersection amid two given SuperHyperEdges and one SuperVertex has one SuperHyperEdge with any common SuperVertex. The SuperHyperModel proposes the specific designs and the specific architectures. The SuperHyperModel is officially called “SuperHyperGraph” and “Neutrosophic SuperHyperGraph”. In this SuperHyperModel, The “specific” cells and “specific group” of cells are SuperHyperModeled as “SuperHyperVertices” and the common and intended properties between “specific” cells and “specific group” of cells are SuperHyperModeled as “SuperHyperEdges”. Sometimes, it’s useful to have some degrees of determinacy, indeterminacy, and neutrality to have more precise SuperHyperModel which in this case the SuperHyperModel is called “neutrosophic”. In the future research, the foundation will be based on the “Cancer’s Recognitions” and the results and the definitions will be introduced in redeemed ways. The recognition of the cancer in the long-term function. The specific region has been assigned by the model [it’s called SuperHyperGraph] and the long cycle of the move from the cancer is identified by this research. Sometimes the move of the cancer hasn’t be easily identified since there are some determinacy, indeterminacy and neutrality about the moves and the effects of the cancer on that region; this event leads us to choose another model [it’s said to be neutrosophic SuperHyperGraph] to have convenient perception on what’s happened and what’s done. There are some specific models, which are well-known and they’ve got the names, and some SuperHyperGeneral SuperHyperModels. The moves and the traces of the cancer on the complex tracks and between complicated groups of cells could be fantasized by a neutrosophic SuperHyperPath(-/SuperHyperCycle, SuperHyperStar, SuperHyperBipartite, SuperHyperMultipartite, SuperHyperWheel). The aim is to find either the longest SuperHyperDefensive SuperHyperAlliances or the strongest SuperHyperDefensive SuperHyperAlliances in those neutrosophic SuperHyperModels. For the longest SuperHyperDefensive SuperHyperAlliances, called SuperHyperDefensive SuperHyperAlliances, and the strongest SuperHyperCycle, called neutrosophic SuperHyperDefensive SuperHyperAlliances, some general results are introduced. Beyond that in SuperHyperStar, all possible SuperHyperPaths have only two SuperHyperEdges but it’s not enough since it’s essential to have at least three SuperHyperEdges to form any style of a SuperHyperCycle. There isn’t any formation of any SuperHyperCycle but literarily, it’s the deformation of any SuperHyperCycle. It, literarily, deforms and it doesn’t form. A basic familiarity with SuperHyperGraph theory and neutrosophic SuperHyperGraph theory are proposed.

Magnetic Nanoparticles (MNPs) is a promising technique to cure brain diseases. On the one hand, by serving as drug carriers, they can bypass the blood-brain barrier and deliver drug molecules to the brain parenchyma; on the other hand, their transport trajectory can be manipulated by applying an external magnetic field. However, due to the complex microstructure of brain tissues, e.g. the anisotropy of white matter (WM), how to achieve desired drug distribution patterns, e.g. uniform distribution, by tuning the drug delivery system is largely unknown. Here, in this study, by adopting a mathematical model capable of capturing the diffusion trajectories of MNPs in the microstructures, we systematically investigated the effects of key parameters in the MNPs delivery system on the equivalent diffusion coefficient of MNPs in the microenvironment of brain WM. The results show that uniform distribution of MNPs in anisotropic tissues can be achieved by adjusting the particle size and magnetic field. We have not only obtained a deeper understanding on how to optimise the MNPs delivery system, it can also be anticipated that an improved mathematical model could even help to achieve complex drug distribution patterns in the complicated brain environment by designing an appropriate combination of the key parameters.

: In this article, the authors propose an algorithmic approach to building a model of the dynamics of economic and, in particular, innovation processes. The approach under consideration is based on a complex algorithm that includes (1) decomposition of the time series into components using singular spectrum analysis; (2) recognition of the optimal component model based on fuzzy rules, and (3) creation of statistical models of individual components with their combination. It is shown that this approach corresponds to the high uncertainty characteristic of the tasks of the dynamics of innovation processes. The proposed algorithm makes it possible to create effective models that can be used both for analysis and for predicting the future states of the processes under study. The advantage of this algorithm is the possibility to expand the base of rules and components used for modeling. This is an important condition for improving the algorithm and its applicability for solving a wide range of problems.

Human behaviour was tipped as the mainstay in the control of further SARS-CoV-2 (COVID-19) spread especially after the lifting of restrictions by many countries. Countries in which restrictions were lifted soon after the first wave, had subsequent waves of the COVID-19 infections and it remains to be seen whether there will be more waves. In this paper, we formulate a deterministic model for COVID-19 incorporating dynamic non-pharmaceutical interventions, dubbed social dynamics. The model steady states are determined and their stability analysed. Numerical simulations are carried out to determine the pack of various parameters that influence the social dynamics. In South Africa, the first wave was the only wave in which the only interventions rested solely on human behavior. The model is thus fitted to COVID-19 data on the first wave in South Africa. The results presented in this paper have implications on the trajectory of the pandemic in the presence of changing social processes.

A proxy war, between a coalition of countries, BLUE and a country, RED, is considered. RED wants to increase the size of the RED territory. BLUE wants to involve more regions in trade and other types of cooperation. GREEN is a small and independent nation that wants to become a member of BLUE. RED attacks GREEN and tries to invade. BLUE decides to give optimal arms support to GREEN. This support can help GREEN in the war against RED and simultaneously reduce the military power of RED, which is valuable to BLUE, also outside this proxy war, since RED may confront BLUE also in other regions. The optimal control problem of dynamic arms support, from the BLUE perspective, is defined in general form. The objective function is a weighted sum of the present value of the free GREEN territory and the present value to BLUE of the net loss of military resources in the RED army. The net loss of RED at a particular point in time is a function of the location of the front line and the size of the mobile GREEN forces behind the RED line. First, it is assumed that the expected net loss is proportional to the force ratio behind the red front line. It is proved that the net loss function is a strictly concave quadratic function of x, the location of the front line. It is also proved that the unique maximum of the net loss function occurs at the same front location, also if the net loss function is proportional to the strength ratio behind the RED lines, raised to any strictly positive exponent. From the BLUE perspective, there is an optimal position of the front. This position is a function of the weights in the objective function and all other parameters. Optimal control theory is used to determine the optimal dynamic BLUE strategy, conditional on a RED strategy, which is observed by BLUE military intelligence. The optimal arms support strategy for BLUE is to initially send a large volume of arms support to GREEN, to rapidly move the front to the optimal position. Then, the support should be almost constant during most of the war, keeping the war front location stationary. In the final part of the conflict, when RED will have almost no military resources left and tries to retire from the GREEN territory, BLUE should strongly increase the arms support and make sure that GREEN rapidly can regain the complete territory and end the war

: Seasonal influenza (a.k.a flu) is responsible for considerable morbidity and mortality across the globe. The three recognized pathogens that cause epidemics during the winter season are influenza A, B and C. The influenza virus is particularly dangerous due to its mutability. Vaccines are an effective tool in preventing seasonal influenza, and their formulas are updated yearly according to WHO recommendations. However, in order to facilitate decision-making in the planning of the intervention, policymakers need information on the projected costs and quantities related to introducing the influenza vaccine, in order to help governments obtain an optimal allocation of the vaccine each year. In this paper, an approach based on a Controlled Elitism Non-Dominated Sorting Genetic Algorithm (CENSGA) model is introduced to optimize the allocation of influenza vaccination. A bi-objective model is formulated to control the infection volume, and reduce the unit cost of the vaccination campaign. An SIR (Susceptible–Infected–Recovered) model is employed for representing a potential epidemic. The model constraints are based on the epidemiological model, time management, and vaccine quantity. A two-phase optimization process is proposed: guardian control followed by contingent controls. The proposed approach is an evolutionary metaheuristic multi-objective optimization algorithm with a local search procedure based on a hash table. Moreover, in order to optimize the scheduling of a set of policies over a predetermined time to form a complete campaign, an extended CENSGA is introduced with a variable-length chromosome (VLC) along with mutation and crossover operations. To validate the applicability of the proposed CENSGA, it is compared with the classical Non-Dominated Sorting Genetic Algorithm (NSGA-II). The results are analyzed using graphical and statistical comparisons in terms of cardinality, convergence, distribution and spread quality metrics, illustrating that the proposed CENSGA is effective and useful for determining the optimal vaccination allocation campaigns.

In this paper mathematical analysis of the HIV/AIDS deterministic model exposed in [Espitia, C. et. al. Mathematical Model of HIV/AIDS Considering Sexual Preferences Under Antiretroviral Therapy, a Case Study in San Juan de Pasto, Colombia, Journal of Computational Biology 29 (2022) 483–493] is made. The objective is to gain insight into the qualitative dynamics of the model determining the conditions for the persistence or effective control of the disease in the community through the study of basic properties such as positiveness and boundedness, calculus of basic reproduction number, stationary points such as disease free equilibrium (DFE), boundary equilibrium (BE) and endemic equilibrium (EE) are calculated, local stability (LAS) of disease free equilibrium. It research allow to conclude that the best way to reduce contagion and consequently to reach a DFE is thought to be the reduction of homosexual partners rate as they are the most affected population by the virus, and are therefore the most likely to become infected and to spread the infection. Increasing the departure rate of infected individuals, leads to a decrease in untreated infected heterosexual men and untreated infected women.

The conceptual and theoretical backbones of symmetry/asymmetry detections or similarity/dissimilarity, and identity/unidentity study are automorphism or isomorphism respectively. However, the development of equations and methods for symmetry/asymmetry detections, similarity/dissimilarity, and identity/unidentity measures deviates from these backbones. In this article, an equation was proposed for symmetry/asymmetry detections, similarity/dissimilarity, and identity/unidentity measures, and proved that its isoreflective pairs-points are functionally bijective and inverse. The proposal, called Kabirian-based optinalysis, is based on the conceptual and theoretical frameworks of automorphism and isomorphism. The Kabirian-based optinalysis is also proven and characterized as invariant (robust) under translation (i.e., scaling and location shift), and rotation or reflection. Computing codes were written in python language for Kabirian-based optinalysis to serve as working codes for application and verification.

After the great success of Mobile wallet, the Internet of Things (IoT) leaves the door wide open for consumers to use their connected devices to access their bank accounts and perform routine banking activities from anywhere, anytime and with any device. However, consumers need to feel safe when interacting with IoT-based payment systems, and their personal information should be protected as much as possible. Unlike as usually done in the literature, in this paper, we introduce two lightweight and secure IoT-based payment protocols based on an identity-based signature scheme. We adopt a server-aided verification technique to construct the first scheme. This technique allows to outsource the heavy computation overhead on the sensor node to a cloud server while maintaining the user's privacy. The second scheme is built upon a pairing-free ECC-based security protocol to avoid the heavy computational complexity of bilinear pairing operations. The security reduction results of both schemes are held in the Random Oracle Model (ROM) under the discrete logarithm and computational Diffie-Hellman assumptions. Finally, we experimentally compare the proposed schemes against each other and against the original scheme on the most commonly used IoT devices: a smartphone, a smartwatch and the embedded device Raspberry Pi. Compared with existing schemes, our proposed schemes achieve significant efficiency in the term of communication and computational overheads

This research investigates a model for the COVID-19 outbreak in Indonesia by analyzing comorbid diseases, self-quarantine, government-provided quarantine, and immunization regimens. Some model parameters are derived from academic journals, others are extrapolated from data collected by the Indonesian Ministry of Health in 2020, and others are assumed. Evaluation of the model reveals non-endemic equilibrium points, endemic equilibrium points, and the basic reproduction number (BRN). The contact rate and infected probability factors impact the increase in the number of COVID-19-infected individuals with comorbid disease within a group. According to the sensitivity analysis of the BRN, the key parameters impacting the spread of COVID-19 are the rate of susceptible recruitment, the rate of contact, the COVID-19 infection death rate, and the likelihood of infection. In addition, sensitivity analyses are provided to examine the effect of parameter changes on subpopulations. We discovered that the natural death rate is the most sensitive parameter based on the sensitivity index after reaching equilibrium. Symmetry aspects appear in some of the visualization of the model’s solution and the BRN and parameters’ sensitivity.

Effective prediction of wastewater treatment is beneficial for precise control of wastewater treatment processes. The nonlinearity of pollutant indicators such as COD and TP makes the model difficult to fit and has low prediction accuracy. The classical deep learning methods have been shown to perform nonlinear modeling. However, there are enormous numerical differences between multi-dimensional data in the prediction problem of wastewater treatment, such as COD above 3000 mg/L and TP around 30 mg/L. It will make current normalization methods challenging to handle effectively, leading to the training failing to converge and the gradient disappears or exploding. This paper proposes a multi-factor prediction model based on deep learning. The model consists of a combined normalization layer and a codec. The combined normalization layer combines the advantages of three normalization calculation methods: z-score, Interval, and Max, which can realize the adaptive processing of multi-factor data, fully retain the characteristics of the data, and finally cooperate with the codec to learn the data characteristics and output the prediction results. Experiments show that the proposed model can overcome data differences and complex nonlinearity in predicting industrial wastewater pollutant indicators and achieve better prediction accuracy than classical models.

Chimeric Antigen Receptor (CAR)-T cell immunotherapy revolutionized cancer treatment and consists of the genetic modification of T lymphocytes with a CAR gene, aiming to increase their ability to recognize and kill antigen-specific tumor cells. The dynamics of CAR-T cell responses in patients presents a multiphasic kinetics with distribution, expansion, contraction, and persistence phases. The characteristics and duration of each phase depend on the tumor type, the infused product, and on patient-specific characteristics. We present a mathematical model which describes the multiphasic CAR-T cell dynamics resulting from the interplay between CAR-T and tumor cells, considering patient and product heterogeneities. The CAR-T cell population is divided into functional (distributed and effector), memory, and exhausted CAR-T cell phenotypes. The model is able to describe the diversity of CAR-T cell dynamic behaviors in different patients and hematological cancers as well as their therapy outcomes. Our results indicate that the joint assessment of the area under the concentration-time curve in the first 28 days and the corresponding fraction of non-exhausted CAR-T cells may be considered as potential markers to classify therapy responses. Overall, the analysis of different CAR-T cell phenotypes can be a key aspect for a better understanding of the whole CAR-T cell dynamics.

We present a differentially private extension of the block coordinate descent based on objective perturbation. The algorithm iteratively performs linear regression in a federated setting on vertically partitioned data. In addition to a privacy guarantee, the algorithm also offers a utility guarantee; a tolerance parameter indicates how much the differentially private regression may deviate from an analysis without differential privacy. The algorithm’s performance is compared with the standard block coordinate descent algorithm and the trade-off between utility and privacy is studied. The performance is studied using both artificial test data and the forest fires data set. We find that the algorithm is fast and able to generate practical predictions with single-digit privacy budgets, albeit with some accuracy loss.

In this work, we explore the extreme searching of multidimensional functions by natural gradient descent based on Dirichlet and generalized Dirichlet distributions. The natural gradient is based on describing multidimensional surface with probability distributions, which allows us to reduce changing the accuracy of gradient and step-size. In this article, we propose an algorithm of natural gradient descent based on Dirichlet and generalized Dirichlet distributions. We demonstrate that the natural gradient descent with step-size adaptation with Dirichlet and generalized Dirichlet distributions has higher accuracy and does not take a large number of iterations for minimizing test functions than gradient descent and Adam.

In this note we describe a new approach to the option pricing problem by introducing the notion of the safe (and acceptable) for the writer price of an option, in contrast to the fair price used in the Black-Scholes model. Our starting point is that the option pricing problem is closely related with the hedging problem by practical techniques. Recalling that the Black - Scholes model does not give us the price of the option but the initial value of a replicating portfolio we observe easily that has a serious disadvantage because assumes the building of this replicating portfolio continuously in time and this is a disadvantage of any model that assumes such a construction. Here we study the problem from the practical point of view concerning mainly the over the counter market. This approach is not affected by the number of the underlying assets and is particularly useful for incomplete markets. In the usual Black-Scholes or binomial or some other approaches one assumes that one can invest or borrow at the same risk free rate $r>0$ which is not true in general. Even if this is the case one can immediately observes that this risk free rate is not a universal constant but is different among different people or institutions. So, the fair price of an option is not so much fair! Moreover, the two sides are not, in general, equivalent against the risk therefore the notion of a fair price has no meaning at all. We also define a variant of the usual binomial model, by estimating safe upward and downward rates $u,d$, trying to give a cheaper safe or acceptable price for the option.

A novel kind of power spectrum is constructed, the inter-spike spectrum, which transforms any signal into its spike-frequency domain. This method clearly shows the apparent cycles in the data and overcomes the problems for spike-train-like signals when using the obvious idea of Fourier-transforming it. We invent this instructive approach with the idea of transforming the τ-recurrence rate of a recurrence plot (RP), which often has a spiky appearance. The τ-recurrence rate is the density of recurrence points along diagonals of the RP, which are parallel to the main diagonal with a distance of τ. In this context the inter-spike spectrum can be interpreted as a nonlinear power spectrum of a potentially high dimensional system which constitutes the RP. The proposed measure is robust to noise and is able to detect and analyze bifurcations.

Brain tumors are abnormal cells growth in the brain tissues that can be cancerous, or not. In any case, they could be a very aggressive disease that should be detected as early as possible. Usually, magnetic resonance imaging (MRI) is the main tool commonly adopted by neurologists and radiologists to identify and classify any possible anomalies present in the brain anatomy. In the present work, an automatic unsupervised method called Z2-$\gamma$, based on the use of adaptive finite-elements and suitable pre-processing and post-processing techniques, is introduced. In particular, the proposed methodology is able to automatically classify whether a given MR image represents a healthy or a diseased brain and in this latter case, is able to locate the tumor area.

The fractional Hilbert transform, a generalization of the Hilbert transform, has been extensively studied in the literature because of its extensive use in optics, engineering, and signal processing. In the present work, we aimed to expand the fractional Hilbert transform to a space of generalized functions known as Boehmians. We introduce a new fractional convolution operator for the fractional Hilbert transform to prove a convolution theorem similar to the classical Hilbert transform and also to extend the fractional Hilbert transform to Boehmians. We also construct a suitable Boehmian space on which the fractional Hilbert transform exists. Further, we investigate convergence of the fractional Hilbert transform for the class of Boehmians and discuss the continuity of the extended fractional Hilbert transform.

Pancreatic cancer represents one of the difficult problems of contemporary medicine. The development illness evolves very slowly, takes place in special sake (stroma) and manifests clinically close to a final stage. Another feature of this pathology is a coexistence (symbiotic) effect between cancer cells and normal cells inside stroma. All these aspects make it difficult to understand the pathogenesis of pancreatic cancer and develop a proper therapy. The emergence of pancreatic pre-cancer and cancer cells represents a branching stochastic process engaging populations of 64 cells differing in the number of acquired mutations. In this study we formulate and calibrate the mathematical model of pancreatic cancer using the quasispecies framework. The mathematical model incorporates the mutation matrix, fineness landscape matrix and the death rates. Each element of the mutation matrix presents the probability of appearing a specific mutation in the branching sequence of cells representing the accumulation of mutations. The model incorporates the cancer cell elimination by effect CD8 T cells (CTL). The down-regulation of the effector function of CTLs and exhaustion are parameterized. The symbiotic effect of coexistence of normal and cancer cells is considered. The computational predictions obtained with the model are consistent with empirical data. The modelling approach can be used to investigate other types of cancers and examine various treatment procedures.

In this work we develop a new algorithm for efficient computation of the Voigt/complex error function. In particular, in this approach we propose a two-domain scheme where number of the grid-points is dependent on the input parameter $y$. The error analysis we performed shows that the MATLAB implementation meets the requirements for radiative transfer applications involving the HITRAN spectroscopic database. The run-time test shows that this MATLAB implementation provides rapid computation especially at smaller range of the parameter x.

The paper investigates new robust code constructions based on bent functions and spline–wavelet transformation. Implementation of bent functions in codes construction increases the probability of error detection in the data channel and cryptographic devices. Meanwhile, use of spline–wavelets theory for constructing the codes gives the possibility to increase system security from the actions of an attacker. Presented constructions combines spline–wavelets functions and bent functions. Developed robust codes, compared to existing ones, have a higher parameter of maximum error masking probability. Illustrated codes are ensuring the security of transmitted information. Some of the granted constructions were implemented on FPGA.

This short paper demonstrates that the galactic spiral equations developed from ROTASE model can be used to morphologically simulate the phyllotactic spiral-ring patterns of plants. It is also noticed during this project that although the numbers in the Fibonacci sequence dominate the phyllotactic numbers, there are still significant phyllotactic numbers are not in the Fibonacci sequence, however, most of those numbers are well fit into a Fibonacci sibling sequence starting with numbers 0 and 2. Many plants show the phyllotactic number pairs of 10-16, 16-26 and 26-42 from the Fibonacci sibling sequence, which is the second dominated sequence in the phyllotaxis.

Working memory refers to the capability of the nervous system to selectively retain short-term memories in an active state. The long-standing viewpoint is that neurons play an indispensable role and working memory is encoded by synaptic plasticity. Furthermore, some recent studies have shown that calcium signaling assists the memory processes and the working memory might be affected by the astrocyte density. Over the last few decades, growing evidence has also revealed that astrocytes exhibit diverse coverage of synapses which are considered to participate in neuronal activities. However, very little effort has yet been made to attempt to shed light on the potential correlations between these observations. Hence, in this article we will leverage a computational neuron-astrocyte model to study the short-term memory performance subject to various astrocytic coverage and we will demonstrate that the short-term memory is susceptible to this factor. Our model may also provide plausible hypotheses for the various sizes of calcium events as they are reckoned to be correlated with the astrocytic coverage.

In this paper, we will consider an autonomous two-dimensional ODE Kolmogorov type 1 system with three parameters, which is a particular system of the general predator–prey systems with 2 a Holling type II. By reformulating this system as a set of two second order differential equations, we 3 will investigate the nonlinear dynamics of the system from the Jacobi stability point of view, using 4 the Kosambi-Cartan-Chern (KCC) geometric theory. We will determine the nonlinear connection, the 5 Berwald connection and the five KCC-invariants which express the intrinsic geometric properties 6 of the system, including the deviation curvature tensor. Furthermore, we will obtain necessary and 7 sufficient conditions on the parameters of the system in order to have the Jacobi stability near the 8 equilibrium points and we will point out these on a few examples.

Nowadays the practice of developing algorithms to maintain the confidentiality of data shows that there is a lack of some features, such as velocity, predictability, etc. Generating pseudorandom numbers is one such problem that lies in the basement of many algorithms, even in hardware microprograms. An unreliable generator can cause cyberattacks on it, despite the security in the upper layers. At the same time, the algorithm should be fast enough to provide uninterrupted circuit work for the entire system. The paper presents a new algorithm generating pseudorandom numbers on cellular automata, which is not only fast and easy-repeating, but unpredictable enough and can be used in cryptographic systems. Using the NIST statistical test suite for random and pseudorandom number generators (PRNG), it is shown that the presented algorithm is more than three times superior to the state-of-the-art methods and algorithms in terms of ? − ?????. A high level of the presented algorithm’s parallelization allows for implementation it effectively on calculators with parallel structure. CPU-based architecture, FPGA-based architecture, CUDA- based architecture of PRNG and different PRNG implementations are presented to confirm high performance of the proposed solution.

A calibrated mathematical model of antiviral immune response to SARS-CoV-2 infection is developed. The model considers the innate and antigen-specific responses to SARS-CoV-2 infection. Recently published data sets from human challenge studies with SARS-CoV-2 were used for parameter estimation. Understanding the regulation of multiple intertwined reaction components of the immune system is necessary for linking the clinical phenotypes of COVID-19 with the kinetics of immune responses. Consideration of multiple immune reaction components in a single calibrated mathematical model allowed us to address some fundamental issues related to pathogenesis of COVID-19, i.e. sensitivity of the peak viral load to parameters characterizing the specific response components, the kinetic coordination of the individual responses, and the factors favoring a prolonged viral persistence. The model provides a tool for predicting the infectivity of patients, i.e. the amount of virus which is transmitted via droplets from the person infected with SARS-CoV-2, depending on the time of infection. The thresholds in the relative unbalance between innate and adaptive response parameters which lead to a prolonged persistence of SARS-CoV-2 due to the loss of a kinetic response synchrony/coordination were identified.

In this paper it will be considered three deterministic models for the study of the interaction between the human immune system and a virus: the logistic model, the Gompertz model, and the generalized logistic model (or Richards model). A qualitative analysis of this models based on dynamical systems theory will be performed, by studying the local behavior of the equilibria and obtaining local dynamics properties from the linear stability point of view. Also, we will compare this models in order to understanding which is more appropriate to models the interaction between the human immune system and a virus. It will be obtained some natural medical interpretations which are available for all three models and can be useful to the medical community.

In the present study we introduce a deterministic mathematical model in order to study the interaction between the human immune system and a virus, like COVID 19. The mathematical analysis is based on the tools of dynamical systems theory, by modeling the interactions between the immune system and the virus, using a predator-prey method and following the ideas of G. Moza, from \cite{TGV}. It will be obtained some conclusions with medical relevance and the main three of them are the followings: 1) A deficiency of a single type of immunity in the early stages of virus proliferation, may lead to a large virus multiplication and the human body can loses the fight against this virus; 2) If the level of all components of the immunity system are at the normal threshold from the first moment of the infection and the immune system kill the virus at higher rate than the rate of virus reproduction, then the human body has the ability to stop the multiplication of the virus and liquidate it, that means the virus will be destroyed; 3) If the level of at least one type of immunity can be increased beyond the normal threshold, by medical interventions (like vaccination) in the early stages of virus infection, then the immune system has a better chance to win the fight with the virus.

This paper reports a new global Mittag-Leffler synchronization criterion with regard to fractional-order hyper-chaotic financial systems by designing the suitable impulsive control and the state feedback controller. The significance of this impulsive synchronization lies in the fact that the backward economic system can synchronize asymptotically with the advanced economic system under the effective impulse macroeconomic management means. Matlab LMI-toolbox is utilized to deduce the feasible solution in numerical example, which shows the effectiveness of the proposed methods. It is worth mentioning that the LMI-based criterion usually requires the activation function of the system to be Lipschitz, but the activation function in this paper is fixed and truly nonlinear, which cannot be assumed to be Lipschitz continuous. This is another mathematical difficulty overcome in this paper.

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

The development of synthetic biology has enabled us to make massive progress on biotechnology and to approach research questions from a brand new perspective. In particular, the design and study of gene regulatory networks in vitro, in vivo and in silico, have played an increasingly indispensable role in understanding and controlling biological phenomena. Among them, it is of great interest to understand how associative learning is formed at the molecular circuit level. Noticeably, mathematical models have been increasingly used to predict the behaviors of molecular circuits. The Fernando’s model, which is thought to be one of the first works in this line of research using the Hill equation, attempted to design a synthetic circuit that mimics Hebbian learning in the neural network architecture. In this article, we carry out in-depth computational analysis of the model and demonstrate that the reinforcement effect can be achieved by choosing the proper parameter values. We also construct a novel circuit that can demonstrate forced dissociation, which was not observed in the Fernando’s model. Our work can be readily used as reference for synthetic biologists who consider implementing the circuits of this kind in biological systems.

The importance of statistics in biological research is inevitable. The appliance of statistics is the most powerful tool to support the scientific hypothesis and to give credibility to biological research methodology, to interpret convoluted explanations and conclusions based on the research findings. The evidences from literatures point the persistent statistical errors, selection of non-proper statistical test with the consequence of misinterpretation of the scientific results which are published in international journals. One of the most extensively used statistical tests in the field of biology (preclinical and clinical) is analysis of variance (ANOVA). ANOVA test is in place of multiple T-tests to compare the means of more than two groups at a time; and there are some important points which biologists should be aware of them to avoid possible misinterpretation of results with the consequence of wrong conclusion. Accordingly, the aim of the current review is to help biologists to understand the basic concepts of ANOVA test and reach to a more valid interpretation of achieved results.

Citizen participation has always been of vital importance in decision-making processes and, likewise, joint work in the formulation and execution of projects that respond to territorial needs. This in order to ensure the best impact of public or private investment, optimizing the result and community articulation. For this reason, this investigation proposes a methodology of evaluation and prioritization of projects based on an approach integrated between the modelling in dynamics systems (SDM) and the Analytical process network (ANP), in which the citizen factors are used as qualitative and quantitative variables inside the posed balance sheet model. This is how the model developed by the Stella Architect software, it offers to the decision makers approximate information about the reality of the community, especially those who were affected by the Colombian conflict and where the expectations of a comprehensive repair could disrupt the socioeconomic system and reveal more pressing interests than those of the central government. The methodology is described by the information gotten through the participatory workshops done in the California county that belongs to the Magdalena region, this is a very affected region by the conflict in Colombia. The investigation results confirm the useful of, first of all, the use of the modelling as an information-generating tool by allowing the simulation of relevant variables in the evaluation and prioritization of projects at different moments of time and, in second place, use the knowledge of the experts that expose the ANP method for establishing the importance of each variable compared with other variables and, likewise of each project compared to the rest of the alternatives in each of the different instants of time in which the evaluation is carried out, given that the results could change in any time, the official ranking is showed after an optimization process with partial results.

In this paper, the functional variable method is used to establish solitary wave solutions and periodic wave solutions of the loaded quadratic non-linear Klein-Gordon equation, the loaded cubic non-linear Klein-Gordon equation and the loaded coupled non-linear Klein-Gordon equation. All solutions of these equations have been examined and three dimensional graphics of the obtained solutions have been drawn by using the Matlab software. The main advantage of the proposed functional variable method over other methods is that it provides more new exact traveling wave solutions along with additional free parameters. The graphical representations of the soliton solutions and the periodic wave solutions by using distinct values of random parameter are demonstrated to better understand their physical features. The exact solutions have its great importance to reveal the internal mechanism of the physical phenomena.

Digital twins are the foundation of autonomous off-road vehicles. Dynamic reduction methods are one of several ways to develop digital twins for off-road vehicles. The article commences with a comprehensive overview of the most widely used dynamic reduction methods and then introduces performance metrics for assessing their efficacies in the context of digital twins. The paper additionally includes a detailed mathematical derivation of the state-space representation for reduced-order finite element models. The state-space representation of the reduced finite element models facilitates their export to problem-solving environments for dynamic analysis. The state-space models are eventually solved utilizing the built-in libraries of numerical solvers in textual and graphical programming platforms. In addition, the article identifies the set of solvers that best suit the simulation of virtual models for off-road vehicles. This article also includes an evaluation of the simulation results of digital models with modes ranging from 0 to 30 Hz. In addition, the article demonstrates the lower bound of the frequency range necessary and sufficient to be retained in off-road vehicle virtual models. Finally, the paper presents the simulation outcomes for digital models of commercial off-road vehicles with custom-built control, electrical, and control systems.

The application of transformations of the state equations of continuous linear and bilinear systems to the canonical form of controllability allows one to simplify the computation of Gramians of these systems. In the paper we developed the method and obtain algorithms for computation of the controllability and observability Gramians of continuous linear and bilinear stationary systems with many inputs and one output, based on the method of spectral expansion of the Gramians and the iterative method for computing the bilinear systems Gramians. An important feature of the concept is the idea of separability of the Gramians expansion: sep333arate computation of the scalar and matrix parts of the spectral Gramian expansion reduces the sub-Gramian matrices computation to calculation of numerical sequences of their elements. For the continuous linear systems with one output the method and the algorithm of the spectral decomposition of the controllability Gramian are developed in the form of Xiao matrices. Analytical expressions for the diagonal elements of the Gramian matrices are obtained, making use of which the rest elements can be calculated. For continuous linear systems with many outputs the spectral decompositions of the Gramians in the form of generalized Xiao matrices are obtained, which allows to significantly reduce the number of calculations. The obtained results are generalized for continuous bilinear systems with one output. Iterative spectral algorithms for computation of elements of Gramians of these systems are proposed. Examples are given that illustrate theoretical results.

This paper presents a cooperative parallel tabu search meta-heuristic for the vehicle routing problem with time windows. It is based on the scheme in which several search threads cooperate by asynchronously exchanging information on the best solutions identified. The exchanges are performed through a mechanism called adaptive memory which holds and manages a pool of solutions. This enforces the asynchronous strategy of information exchanges and ensures the independence of the individual search threads. Each of these independent threads implements a tabu search meta-heuristic. Comprehensive computational experiments and comparisons to best known solutions show that the proposed cooperative parallel tabu search algorithm is able to achieve 48 new best solutions on VRPTW benchmark instances.

Power series expansions are useful in approximation theory and mathematical physics. The manuscript presents several types of fractional Taylor expansions of sufficiently smooth functions. This is achieved by employing an incremental regularization procedure to the computation of the derivative. The series are constructed in an algorithmic way, which allows for its implementation in computer algebra systems.

A correct system design can be systematically obtained from a specification model of a real-time system that integrates hybrid measurements in a realistic industrial environment, this has been carried out through complete Matlab / Simulink / Stateflow models. However, there is a widespread interest in carrying out that modeling by resorting to Machine Learning models, which can be understood as Automated Machine Learning for Real-time systems that present some degree of hybridization. An induction motor controller which must be able to maintain a constant air flow through a filter is one of these systems and it is discussed in the paper as a study case of closed-loop control system. The article discusses a practical application of ML methods that demonstrates how to replace such closed loop in industrial control systems with a Simulink block generated from neural networks to show how the proposed procedure can be applied to derive complete hybrid system designs with artificial neural networks (ANN). In the proposed ANN-based method to design a real-time hybrid system with continuous and discrete components, we use a typical design of a neural network, in which we define the usual phases: training, validation, and testing. The generated output of the model is made up of reference variables values of the cyber-physical system, which represent the functional and dynamic aspects of model. They are used to feed Simulink/Stateflow blocks in the real target system.

The process of revitalising quiescent cells in order for them to proliferate plays a pivotal role in the repair of warn out tissues as well as for tissue homeostasis. This process is also crucial in the growth, development and well-being of higher multi-cellular organisms such as mammals. Deregulation of quiescent-proliferation transition is related to many diseases, in particular cancer. Recent studies have revealed that this process is regulated tightly by the Rb−E2F bistable switch mechanism. Based on experimental observations, in this study, we formulate a mathematical model to examine the effect of the growth factor concentration on the proliferation-quiescence transition in human cells. Working with a non-dimensionalised model, we prove positivity, boundedness and uniqueness of solutions. To understand model solution behaviour close to bifurcation points, we carry out bifurcation analysis, which is further verified by use of numerical bifurcation analysis, sensitivity analysis and numerical simulations. Indeed, bifurcation and numerical analysis of the model predicted a transition between stable, bistable and stable states, which are dependent on the growth factor concentration parameter (GF). The derived predictions confirm experimental observations.

The capacitated team orienteering problem with time windows (CTOPTW) is a NP-hard combinatorial optimization problem. In the CTOPTW, a set of customers is given each with a profit, a demand, a service time and a time window. A homogeneous fleet of vehicles is available for serving customers and collecting their associated profits. Each vehicle is constrained by a maximum tour duration and a limited capacity. The CTOPTW is concerned with the determination of a preset number of vehicle tours that begin and end at a depot, visit each customer no more than once while satisfying the time duration, time window and vehicle capacity constraints on each tour. The objective is to maximize the total profit collected. In this study we propose an iterated local search (ILS) algorithm to deal with the CTOPTW. ILS is a single solution based meta-heuristic that successively invokes a local search procedure to explore the solution space. A perturbation operator is used to modify the current local optimum solution in order to provide a starting solution for the local search procedure. As different problems and instances have different characteristics, the success of the ILS is highly dependent on the local search procedure, the perturbation operator(s) and the perturbation strength. The basic ILS uses a single perturbation operator and the perturbation strength remains the same during the optimization process. To address these issues, we use three different perturbation operators and a varying perturbation strength which changes as the algorithm progresses. The idea is to assign a larger perturbation strength in the early stages of the search in order to focus on exploring the search space. The perturbation strength is gradually decreased so that we focus more on exploitation. The computational results show that the proposed ILS algorithm is able to generate high quality solutions on the CTOPTW benchmark instances taken from the scientific literature, demonstrating its efficiency in terms of both the solution quality and computational time. Moreover, the proposed ILS produces 21 best known results and 5 new best solutions.