The asymptotic dynamics of the classes of rational difference equations (RDEs) of third order defined over the positive real-line as $$\displaystyle{x_{n+1}=\frac{x_{n}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-1}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-2}}{ax_n+bx_{n-1}+cx_{n-2}}}$$ and $$\displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-1}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-2}}}$$ is investigated computationally with theoretical discussions and examples. It is noted that all the parameters $a, b, c$ and the initial values $x_{-2}, x_{-1}$ and $x_0$ are all positive real numbers such that the denominator is always positive. Several periodic solutions with high periods of the RDEs as well as their inter-intra dynamical behaviours are studied.

Nanofluids are a novel class of heat transfer fluid that plays a vital role in industries. In mathematical investigations, these fluids are modeled in terms of traditional integer-order partial differential equations (PDEs). It is recognized that traditional PDEs cannot decode the complex behavior of physical flow parameters and memory effects. Therefore, this article intends to study the mixed convection heat transfer in nanofluid over an inclined vertical plate via fractional derivatives approach. The problem in hand is modeled in connection with Atangana-Baleanu fractional derivatives without singular and local kernel having strong memory. The human blood is considered as base fluid dispersing carbon nanotube (CNTs) (single-wall carbon nanotubes (SWCNTs) and multi-wall carbon nanotubes(MWCNTs )) into it to form blood-CNTs nanofluid. The nanofluids are considered to flow in a saturated porous medium under the influence of an applied magnetic field. The exact analytical expressions for velocity and temperature profiles are acquired using the Laplace transform technique and plotted in various graphs. The empirical results indicate that the memory effect decreases with increasing fractional parameters in the case of both temperature and velocity profiles. Moreover, the temperature profile is higher for blood-SWCNTs by reason of higher thermal conductivity whereas, this trend is opposite in case of velocity profile due densities difference.

The world is now undergoing through a global emergency due to COVID-19 which needs immediate remedies in order to strengthen the healthcare facility to save the nations. Looking towards to the remedies, research on different aspects including the genomic and proteomic level characterizations of the SARS-CoV2 are necessarily important. In this present study, the spatial representation/composition of twenty amino acids across the primary protein sequences of SARS-CoV2 have been looked into through different parameters viz. Shannon entropy, Hurst exponent in order to fetch the autocorrelation and amount of information over the spatial representations. Also frequency distribution of each of the amino acids over the protein sequences have been chalked out.

This paper analyzes the robust dissipativity of uncertain stochastic generalized neural networks (USGNNs) with Markovian jumping parameters and time-varying delays. In practical applications most of the systems refer to uncertainties, hence, the norm-bounded parameter uncertainties and stochastic disturbance are considered. Then, by constructing an appropriate Lyapunov-Krasovskii functional (LKF) and by employing integral inequalities LMI-based sufficient conditions of the considered systems are established. Numerical simulations are given to show the merit of the presented results.

The world is now undergoing through a global emergency due to COVID-19 which needs immediate remedies in order to strengthen the healthcare facility to save the nations. Looking towards to the remedies, research on different aspects including the genomic and proteomic level characterizations of the SARS-CoV2 are necessarily important. In this present study, the spatial representation/composition of twenty amino acids across the primary protein sequences of SARS-CoV2 have been looked into through different parameters viz. Shannon entropy, Hurst exponent in order to fetch the autocorrelation and amount of information over the spatial representations. Also frequency distribution of each of the amino acids over the protein sequences have been chalked out.

The purpose of this paper is to bring to light a method through which the global in time existence for arbitrary large in H¹ initial data of a strong solution to 3D periodic Navier-Stokes equations follows. The method consists of subdividing the time interval of existence into smaller sub-intervals carefully chosen. These sub-intervals are chosen based on the hypothesis that for any wavenumber m, one can find an interval of time on which the energy quantized in low-frequency components (up to m) of the solution u is lesser than the energy quantized in high-frequency components (down to m) or otherwise the opposite. We associate then a suitable number m to each one of the intervals and we prove that the norm ||u(t)||_H1 is bounded in both mentioned cases. The process can be continued until reaching the maximal time of existence T_max which yields the global in time existence of strong solution.

Comments are exploited by product vendors to measure satisfaction of consumers. With the advent of Natural Language Processing (NLP), comments on Google Play can be processed to extract knowledge on applications such as their reputation. Proposals in that direction are either informal or interested merely on functionality. Unlike, this work aims to determine reputation of Android applications in terms of confidentiality, integrity, availability and authentication (CIAA). This work proposes a model of assessing app reputation relying on sentiment analysis and text analysis of comments. While assuming that comments are reliable, we collect Google Play applications subject to comments which include security keywords. An in-depth analysis of keywords based on Naive Bayes classification is made to provide polarity of any comment. Based on comment polarity, reputation is evaluated for the whole application. Experiments made on real applications including dozens to billions of comments, reveal that developers lack to make efforts to guarantee CIAA services. A fine-grained analysis shows that not security reputed applications can be reputed in specific CIAA services. Results also show that applications with negative security polarities display in general positive functional polarities. This result suggests that security checking should include careful comment analysis to improve security of applications.

The purpose of this paper is to bring to light a method through which the global in time existence for arbitrary large in H¹ initial data of a strong solution to 3D periodic Navier-Stokes equations follows. The method consists of subdividing the time interval of existence into smaller sub-intervals carefully chosen. These sub-intervals are chosen based on the hypothesis that for any wavenumber m, one can find an interval of time on which the energy quantized in low-frequency components (up to m) of the solution u is lesser than the energy quantized in high-frequency components (down to m) or otherwise the opposite. We associate then a suitable number m to each one of the intervals and we prove that the norm ||u(t)||_H1 is bounded in both mentioned cases. The process can be continued until reaching the maximal time of existence T_max which yields the global in time existence of strong solution.

In this work, the problem is considered for eliminating a non-uniform rectilinear smearing of an image by mathematical and computer-based way, for example, a picture of several cars taken with a fixed camera and moving with different speeds. The problem is described by a set of 1-dimensional Fredholm integral equations (IEs) of the first kind of convolution type with a 1-dimensional point spread function (PSF) at uniform smearing and by a set of new 1-dimensional IEs of a general type (not convolution type) with a 2-dimensional PSF at non-uniform smearing. The problem is also described by one 2-dimensional IE of convolution type with a 2-dimensional PSF at uniform smearing and by a new 2-dimensional IE of a general type with a 4-dimensional PSF at non-uniform smearing. The problem for solving Fredholm IE of the first kind is ill-posed (unstable). Therefore, IEs of convolution type are solved by the Fourier transform (FT) method and Tikhonov's regularization (TR), and IEs of general type are solved by the quadrature/cubature and TR methods. Moreover, the magnitude of the image smear Δ is determined by the original “spectral method”, which increases the accuracy of image restoration. It is shown that the use of a set of 1-dimensional IEs is preferable to one 2-dimensional IE in the case of non-uniform smearing. In the inverse problem (image restoration), the Gibbs effect (the effect of false waves) in the image may occur. It can be edge or inner. The edge effect is well suppressed by the proposed technique “diffusing the edges". In the case of an inner effect, it is eliminated with difficulty, but the image smearing itself plays the role of diffusing and suppresses the inner Gibbs effect to a large extent. It is shown (in the presence of impulse noise in an image) that the well-known Tukey median filter can distort the image itself, and the Gonzalez adaptive filter also distorts the image (but to a lesser extent). We propose a modified adaptive filter. A software package was developed in MatLab and illustrative calculations are performed.

Periodic travelling waves are observed in various brain activities including visual, motor, language, sleep, and so on. There are several neural field models describing periodic waves assuming nonlocal interaction and, possibly, inhibition, time delay, or some other properties. In this work we study the influence of asymmetric connectivity functions and of time delay on the emergence of periodic waves and on their properties. Nonlinear wave dynamics is studied, including modulated and aperiodic waves. Multiplicity of waves for the same values of parameters is observed. External stimulation in order to restore wave propagation in a damaged tissue is discussed.

A regularization approach to model selection, within a generalized HJM framework, is introduced which learns the closest arbitrage-free model to a prespecified factor model. This optimization problem is represented as the limit of a one-parameter family of computationally tractable penalized model selection tasks. General theoretical results are derived and then specialized to affine term-structure models where new types of arbitrage-free machine learning models for the forward-rate curve are estimated numerically and compared to classical short-rate and the dynamic Nelson-Siegel factor models.

In this study, we first present a time-fractional L$\hat{e}$vy diffusion equation of the exponential option pricing models of European option pricing and the risk-neutral parameter. Then, we modify a particular L$\hat{e}$vy-time fractional diffusion equation of European-style options. Introduce a more general model from the models based on the L$\hat{e}$vy-time fractional diffusion equation and review some recent findings regarding of the Europe option pricing of risk-neutral free.

The ongoing outbreak of the novel coronavirus pneumonia (also known as COVID-19) has triggered a series of stringent control measures in China, such as city closure, traffic restrictions, contact tracing and household quarantine. These containment efforts often lead to changes in the contact pattern among individuals of the population. Many existing compartmental epidemic models fail to account for the effects of contact structure. In this paper, we devised a pairwise epidemic model to analyze the COVID-19 outbreak in China based on confirmed cases reported during the period February 3rd--17th, 2020. By explicitly incorporating the effects of family clusters and contact tracing followed by household quarantine and isolation, our model provides a good fit to the trajectory of COVID-19 infections and is useful to predict the epidemic trend. We obtained the average of the reproduction number $R=1.494$ ($95\%$ CI: $1.483-1.507$) for Hubei province and $R=1.178$ ($95\%$ CI: $1.145-1.158$) for China (except Hubei), suggesting that some existing studies may have overestimated the reproduction number by neglecting the dynamical correlations and clustering effects. We forecasted that the COVID-19 epidemic would peak on February 13th ($95\%$ CI: February $9-17$th) in Hubei and 6 days eariler in the regions outside Hubei. Moreover the epidemic was expected to last until the middle of March in China (except Hubei) and late April in Hubei. The sensitivity analysis shows that ongoing exposure for the susceptible and population clustering play an important role in the disease propagation. With the enforcement of household quarantine measures, the reproduction number $R$ effectively reduces and epidemic quantities decrease accordingly. Furthermore, we gave an answer to the public concern on how long the stringent containment strategies should maintain. Through numerical analysis, we suggested that the time for the resumption of work and production in China (except Hubei) and Hubei would be the middle of March and the end of April, 2020, respectively. These constructive suggestions may bring some immeasurable social-economic benefits in the long run.

In this study, by using Least Square Method, the dataset for the Gross Domestic Product per capita is modeled as a function satisfying the fractional differential equation. The function itself is assumed to be the finite summation of its previous values and the derivatives with unknown coefficients. Then, the prediction for the upcoming years is done by having an approach dividing the dataset into 4 regions corresponding to four different tasks. The mathematical model of the Gross Domestic Product (GDP) per capita of the countries (and union) which are Brazil, China, European Union (EU), India, Italy, Japan, UK, the USA, Spain, and Turkey is constructed with a new methodology called as the deep assessment method which comes from the expressing an arbitrary function modeling the dataset as the finite summation of its previous values and the derivatives with unknown coefficient. The method uses the fractional calculus properties combining with Least Square Method and is compared to Long short-term memory (LSTM) algorithm which is a special type of neural network used for time sequences in general.

A good recommender system can infer customers’ preferences based on their historical purchase records, and recommend products that the customers may be interested in, saving them a lot of time and energy. For enterprises, it is difficult to recommend accurately to each customer, and the bad recommendation may be counterproductive. Customer loyalty is an indicator that measures the preference relationship between customers and products in the field of marketing. A hypothesis is proposed in this study: if companies can divide customers into different groups based on customer loyalty, the recommendation effect on certain groups is better than that on overall customers. In this study, customer loyalty is measured by four features of the RFML model. All customers are viewed as points on a four-dimensional space, which are clustered by the k-means model. Two recommendation algorithms based on complex networks are tested: recommendation algorithm based on bipartite graph and PersonalRank (BGPR), and recommendation algorithm based on a single vertex set network and DeepWalk (SVDW). The experimental results show that customer loyalty has improved the effectiveness of the two algorithms over 14%, and the recommendation effect is the best on customer groups with a loyalty level of 4 (the highest level is 5). The recommendation algorithms with customer loyalty are better than using them alone.

In this paper, we study the existence of the random -pullback attractor of a non-autonomous local modified stochastic Swift-Hohenberg equation with multiplicative noise in stratonovich sense. It is shown that a random -pullback attractor exists in when its external force has exponential growth. Due to the stochastic term, the estimate are delicate, we overcome this difficulty by using the Ornstein-Uhlenbeck(O-U) transformation and its properties.

Rough volatility models are popularized by \cite{gatheral2018volatility}, where they have shown that the empirical volatility in the financial market is extremely consistent with rough volatility. Fractional Riccati equation as a part of computation for the characteristic function of rough Heston model is not known in explicit form as of now and therefore, we must rely on numerical methods to obtain a solution. In this paper, we give a short introduction to option pricing theory and an overview of the current advancements on the rough Heston model.

Primary energy consumption is one of the key drivers of global CO₂ emissions that, in turn, heavily depend on the efficiency of involved technologies. Either the improvement in technology efficiency or the expansion of non-fossil fuel consumption require large investments. The planning and financing of such investments, by policy makers or global energy firms, require, in turn, reliable measures of associated global spreads and their evolution in time. In this paper, our main contribution is the introduction of index measures for accessing global spreads (that is, measures of inequality or inhomogeneity in the statistical distribution of a related quantity of interest) of technology efficiency and CO₂ emission in primary energy consumption. These indexes are based on the Gini index, as used in economical sciences, and generalised entropy measures. Regarding primary energy sources, we consider petroleum, coal, natural gas and non-fossil fuels. Between our findings, we attest some stable relations in the evolution of global spreads of technology efficiency and CO₂ emission, and a positive relation between changes in global spreads of technology efficiency and use of non-fossil fuel.

In this paper, initially a mathematical model is formulated for transient frequency of power system considering time delays which occur while transmitting the control signals in open communication infrastructure. Time delay negligence in a power system leads to improper measurement of frequency variation in power system. The study of impact of time delays on the stability of power system is performed by estimating the decay rate of frequency wave form using Kalman Filter (KF). In power system, there is a possibility of multiple time delays. This paper also focusses on developing Interacting Multiple Model (IMM) Algorithm with multiple model space using Kalman Filter (KF) as state estimator tool. The multiple time delays in power system is considered as multiple model space. The result shows that KF provides better estimate of correct model for a particular input-set. The qualitative properties of Riccati difference equation (RDE) in terms of state error covariance of IMMKF are also analyzed and presented.

We introduce a new, and quite general variational model for opinion dynamics based on pairwise interaction potentials and a range of opinion evolution protocols ranging from random interactions to global synchronous flows in the opinion state space. The model supports the concept of topic “coupling”, allowing opinions held by individuals to be changed via indirect interaction with others on different subjects. Interaction topology is governed by a graph that determines interactions. Our model, which is really a family of variational models, has, as special cases, many of the previously established models for the opinion dynamics. After introducing the model, we study the dynamics of the special case in which the potential is either a tent function or a constructed bell-like curve. We find that even in these relatively simple potential function examples there emerges interesting behavior. We also present results of preliminary numerical explorations of the behavior of the model to motivate questions that can be explored analytically

In this age of Facebook, Instagram, and Twitter, there is rapidly growing interest in understanding network-enabled opinion dynamics in large groups of autonomous agents. The phenomena of opinion polarization, the spread of propaganda and fake news, and the manipulation of sentiment is of interest to large numbers of organizations and people. Whether it is the more nefarious players such as foreign governments that are attempting to sway elections or it is more open and above board, such as researchers who want to make large groups of people aware of helpful innovations, what is at stake is often significant. In this paper, we review opinion dynamics including the extensions of many classical models as well as some new models that deepen understanding. For example, we look at models that track the evolution of an individual’s power, that include noise, and that feature sequentially dependent topics, to name a few. While the first papers studying opinion dynamics appeared over 60 years ago, there is still a great deal of room for innovation and exploration. We believe that the political climate and the extraordinary (even unprecedented) events in the sphere of politics in the last few years will inspire new interest and new ideas. It is our aim to help those interested researchers understand what has already been explored in a significant portion of the field of opinion dynamics. We believe that in doing this, it will become clear that there is still much to be done.

We develop a mathematical model to provide epidemic predictions for the 2019-nCov epidemic in Wuhan, China. We use reported case data up to January 31, 2020 from the Chinese Center for Disease Control and Prevention and the Wuhan Municipal Health Commission to parameterize the model. From the parameterized model we identify the number of unreported cases. We then use the model to project the epidemic forward with varying level of public health interventions. The model predictions emphasize the importance of major public health interventions in controlling 2019-nCov epidemics.

Deep coal beds have been suggested as possible usable underground geological locations for carbon dioxide storage. Furthermore, injecting carbon dioxide into coal beds can improve the methane recovery. Due to importance of this issue, a novel investigation has been done on adsorption of carbon dioxide on various types of coal seam. This study has proposed four types of Gaussian Process Regression (GPR) approaches with different kernel functions to estimate excess adsorption of carbon dioxide in terms of temperature, pressure and composition of coal seams. The comparison of GPR outputs and actual excess adsorption expresses that proposed models have interesting accuracy and also the Exponential GPR approach has better performance than other ones. For this structure, R2=1, MRE=0.01542, MSE=0, RMSE=0.00019 and STD=0.00014 have been determined. Additionally, the impacts of effective parameters on excess adsorption capacity have been studied for the first time in literature. According to these results, the present work has valuable and useful tools for petroleum and chemical engineers who dealing with enhancement of recovery and environment protection.

In this article, we investigate two types of double dispersion equations in two dierent dimensions. Double dispersion equation were derived to describe long nonlinear wave evolution in a thin hyperelastic rod. Conservation laws are obtained for these equations by the application of the multiplier method. Finally, travelling waves and line travelling waves are respectively considered for these two equations.

In the present study, a mathematical model of non-steady partial differential equation from the process of oxygen mass transport in the human pulmonary circulation is proposed. Mathematical modelling of this kind of problems lead to a non-steady partial differential equation and for its numerical simulation, we have used finite differences. The aim of the process is the exact numerical analysis of the study, wherein consistency, stability and convergence is proposed. The necessity of doing the process is that, we would like to increase the order of numerical solution to a higher order scheme. An increment in the order of numerical solution makes the numerical simulation more accurate, also makes the numerical simulation being more complicated. In addition, the process of numerical analysis of the study in this order of solution needs more research work.

Pressure fluctuations beneath hydraulic jumps downstream of Ogee spillways potentially damage stilling basin beds. This paper deals with the extreme pressures underneath free hydraulic jumps along a smooth stilling basin. The experiments were conducted in a laboratory flume. From the probability distribution of measured instantaneous pressures, the pressures with different non-exceedance probabilities (P*a%) could be determined. It was verified that the maximum pressure fluctuations, as well as the negative pressures, are located at the positions closest to the spillway toe. The minimum pressure fluctuations are located at the downstream of hydraulic jumps. It was possible to assess the cumulative curves of P*a% related to the characteristic points along the basin, and different Froude numbers. To benchmark, the results, the dimensionless forms of mean pressures, standard deviations, and pressures with different non-exceedance probabilities were assessed. It was found that an existing methodology can be used to interpret the present data, and pressure distribution in similar conditions, by using a new third-order polynomial relationship for the standard deviation (σ*X) with the determination coefficient (R2) equal to 0.717. It was verified that the new optimized adjustment gives more accurate results for the estimation of the maximum extreme pressures than the minimum extreme pressures.

A mathematical model given by a two - dimensional differential system is introduced in order to understand the transition process from the normal hematopoiesis to the chronic and accelerated acute stages in chronic myeloid leukemia. A previous model of Dingli and Michor is refined by introducing a new parameter in order to differentiate the bone marrow microenvironment sensitivities of normal and mutant stem cells. In the light of the new parameter, the system now has three distinct equilibria corresponding to the normal hematopoietic state, to the chronic state, and to the accelerated acute phase of the disease. A characterization of the three hematopoietic states is obtained based on the stability analysis. Numerical simulations are included to illustrate the theoretical results.

This work - in two parts - tries to promote a new method in Mathematical Modeling, not ever used before, namely the use of the methods of Non-standard Analysis in Topoi. For the general mathematical theory of topoi (Topoi Theory) we refer to the book [10] . For non-standard analysis in the particular Boolean topos SET (the category of sets) we refer to the book [11]. This first part of the work is related, from an applicative point of view, mainly with Mathematical Neuroscience, more precisely a model of the human thinking, consciousness and subconsciousness. The connections between brain and mind will be also scketched, but a closer study also need some further topics to be studied in Part II: Artificial Intelligence, namely (Quantum) neural networks and (Quantum) Turing machine. However both parts (I. and II.) should be considered together, because they complement each other. The topoi model the intuitionistic logic (multi-valued) and have been used in Quantum Physics (see [9] and its references) while Non-Standard Analysis in SET (introduced by Abraham Robinson[31]) has been applied in Mathematical Economics (see [2] and its references); however, the combination was never used until now in Applied Mathematics. Even more, the Non-standard Analysis in Topoi is not too much studied from the point of view of Pure Mathematics either. One of the main objective of this work is to produce progress in this aria of Mathematics also, mainly for the SET-type topoi and exponential topoi, those topoi used in Quantum Physics and in which we need Non-stantard Analysis for our purposes. Based on the paper [13] we propose the logic of non-standard extensions in topoi as a model of the human thinking (based on infons), these theories representing top and very difficult results in Abstract Mathematics. We will analise the relation between brain and mind via the genetic-epigenetic interplay. More details will be provided in the Introduction, followed by a short presentation of the general theory, later connected with the (Quantum) Yang-Baxter equations in topoi. The connections between these equations and braid groups and knots, with intended applications in Neuroscience and Artificial Intelligence will be further analysed in Part II. Besides promoting a new idea of modeling, the present paper aims to promote the building of a Research Proposal for the European $''$Human Brain Project$''$ (https://www.humanbrainproject.eu/en/), and the NSF-Europe Program $''$Collaborative Research in Computational Neuroscience$''$ (https://www.nsf.gov/funding/pgm\_summ.jsp?pims\_id=5147), jointly with the new founded $''$National Centre for the Research of the Brain$''$ of the Romanian Academy of Science (https://acad.ro/centreAR/CNCC/CNCC.pdf, in Romanian). $''$Blue Brain Project$''$ (https://www.epfl.ch/research/domains/bluebrain/) could also be of big interest. We hope that the journals edited by MDPI (Axioms, in particular) and by the Simion Stoilow Institute of the Romanian Academy - Revue Roumaine de Mathematique Pures et Appliqu (http://imar.ro/journals/Revue\_Mathematique/home\_page.html) - among others - will host special issues related to the research subjects described in both parts of this work, starting with this one. For further details and some of the progress on this research subject, see www.ovidiufpasarescu.com.}

It is well known that the Continuous Galerkin Finite Element (CGFE) method is globally consistent with respect to the first law of thermodynamics. This means that, for any mesh, all obtained discrete solutions will conserve total energy. One might expect, that the method is, also, globally consistent with respect to the second law of thermodynamics. In this paper, we formally study if such conjecture is true. The heat conduction equation is used as the physical model for this analysis. In the present study it is proved that the conjecture is false: at least, for standard piecewise linear (1D and 2D) elements, the CGFE method is not always globally consistent with respect to the second law of thermodynamics. In other words, some obtained discrete solutions can violate the global postulate of the second law, which asserts that total entropy can never decrease.

In this paper, initially a mathematical model is formulated for transient frequency of power system considering time delays which occur while transmitting the control signals in open communication infrastructure. Time delay negligence in a power system leads to improper measurement of frequency variation in power system. The study of impact of time delays on the stability of power system is performed by estimating the decay rate of frequency wave form using Kalman Filter (KF). In power system, there is a possibility of multiple time delays. This paper also focuses on developing Interacting Multiple Model(IMM) Algorithm with multiple model space using Kalman Filter(KF) as state estimator tool. The multiple time delays in power system is considered as multiple model space The result shows that KF provides better estimate of correct model for a particular inputset. The qualitative properties of Riccati difference equation(RDE) in terms of state error covariance of IMMKF are also analyzed and presented.

This article aims to discuss, the stability and boundedness character of the solutions of the rational equation of the form \begin{equation}\label{eql21.1} y_{t+1}=\frac{\nu\epsilon^{-y_t}+\delta\epsilon^{-y_{t-1}}}{\mu+\nu y_t+\delta y_{t-1}},\quad t\in N(0). \end{equation} Here, $\epsilon>1, \nu,\delta,\mu\in (0,\infty)$ and $y_0, y_1$ are taken as arbitrary non-negative reals and $N(a)=\{a,a+1,a+2,\cdots \}$. Relevant examples are provided to validate our results. The exactness is tested using MATLAB.

In this study, we investigate spectral structure of conformable Sturm-Liouville problems and with this end, we obtain representation of solutions under different initial conditions and asymptotic formulas for eigenfunctions, eigenvalues, norming constants and normalized eigenfunctions. Consequently, we prove the existence of infinitely many eigenvalues. Also, we compare the solutions with graphics with different orders, different eigenvalues, different potentials and so, we observe the behaviors of eigenfunctions. We give an application to the α-orthogonality of eigenfunctions and reality of eigenvalues for conformable Sturm-Liouville problems defined by [15] in the last section.

The known three-component reaction-diffusion system modeling competition and co-existence of different language speakers is under study. A modification of this system is proposed, which is examined by Lie symmetry method; furthermore exact solutions in the form of traveling fronts are constructed and their properties are identified. Plots of the traveling fronts are presented and the relevant interpretation describing the language shift occurred in Ukraine during the Soviet times is suggested

This article presents a new method of segmenting images with gray levels. The method is based on determining several thresholds for separation of gray levels. The determination of these thresholds is done using the certainty of the neutrosophic information. The concept of this method can be stated simply: to choose the local maximums for the neutrosophic certainty.

The work is devoted to the investigation of virus quasispecies evolution and diversification due to mutations, competition for host cells, and cross-reactive immune responses. The model consists of a nonlocal reaction-diffusion equation for the virus density depending on the genotype considered as a continuous variable and on time. This equation contains two integral terms corresponding to the nonlocal effects of virus interaction with host cells and with immune cells. In the model, a virus strain is represented by a localized solution concentrated around some given genotype. Emergence of new strains corresponds to a periodic wave propagating in the space of genotypes. The conditions of appearance of such waves and their dynamics are described.

In this paper, we investigate the equilibrium points, stability of two equilibrium points, convergences of negative equilibrium point, periodic solutions, and existence of bounded or unbounded solutions of a system of nonlinear difference equations x_n+1 =x_n-1y_n - 1, y_n+1 = y_n-1x_n - 1 n = 0,1,..., where the initial values are real numbers. Additionally we present some numerical examples to verify our theoretical results.

We consider a model of predator-prey interaction at fractional-order where the predation obeys the ratio-dependent functional response and the prey is linearly harvested. For the proposed model, we show the existence, uniqueness, non-negativity as well as the boundedness of the solutions. Conditions for the existence of all possible equilibrium points and their stability criteria, both locally and globally, are also investigated. The local stability conditions are derived using the Magtinon's theorem, while the global stability is proven by formulating an appropriate Lyapunov function. The occurance of Hopf bifurcation around the interior point is also shown analytically. At the end, we implement the Predictor-Corrector scheme to perform some numerical simulations.

In this paper, the fixed-time stabilization problem for a class of uncertain chained system is addressed by using a novel nonsingular recursive terminal sliding mode control approach. A fixed-time controller and an adaptive law are designed to guarantee the uncertain chained form system both Lyapunov stable and fixed-time convergent within the settling time. The advantage of the controller based on the sliding mode is that the settling time does not depend on the system initial state. Furthermore, we use RBF neural network to estimate the uncertainty of the system. Finally, the simulation results demonstrate the performance of the control laws.

We address the dynamics of two-dimensional Navier-Stokes models with infinite delay and hereditary memory, whose kernels are a much larger class of functions than the one considered in the literature, on a bounded domain. We prove the existence and uniqueness of weak solutions by means of Faedo-Galerkin method. Moreover, we establish the existence of global attractor for the system with the existence of a bounded absorbing set and asymptotic compact property.

In light of Liu et al.'s original works, this research article revisits the solution of Burgers's nonlinear equation. The researcher found two exact and explicit solutions for groups $G_4$ and $G_6$, as well as a general solution. Their applications were conducted by using numerical calculations.

In light of Liu's original works, this paper revisits the solution of general Burgers's nonlinear equation. We obtain two exact and explicit solutions for group $G_4$ and $G_6$, and a most general solution as well. As applications, a numerical example is carried out.

This paper presents an extension of the resuts obtained in previous work by the author concerning the application of global optimization techniques to the design of finite strategic games with mixed strategies. In that publication the Fuzzy ASA global optimization method was applied to many examples of synthesis of strategic games with one previously specified Nash equilibrium, evidencing its ability in finding payoff functions whose respective games present those equilibria, possibly among others. That is to say, it was shown it is possible to establish in advance a Nash equilibrium for a generic finite state strategic game and to compute payoff functions that will make it feasible to reach the chosen equilibrium, allowing players to converge to the desired profile, considering that it is an equilibrium of the game as well. Going beyond this state of affairs, the present article shows that it is possible to "impose" multiple Nash equilibria to finite strategic games by following the same reasoning as before, but with a slight change: using the same fundamental theorem of Richard D. McKelvey, modifying the original prescribed objective function and globally minimizing it. The proposed method, in principle, is able to find payoff functions that result in games featuring an arbitrary number of Nash equiibria, paving the way to a substantial number of potential applications.

This paper investigates an inverse problem of option pricing in the extended Black--Scholes model. We identify the model coefficients from the measured data and attempt to find arbitrage opportunities in financial markets using a Bayesian inference approach. The posterior probability density function of the parameters is computed from the measured data. The statistics of the unknown parameters are estimated by Markov Chain Monte Carlo (MCMC), which explores the posterior state space. The efficient sampling strategy of MCMC enables us to solve inverse problems by the Bayesian inference technique. Our numerical results indicate that the Bayesian inference approach can simultaneously estimate the unknown drift and volatility coefficients from the measured data.

In this study, based on our previous study, we examined the mathematical properties, especially the stability of the equilibrium for our proposed mathematical model. By means of the results of the stability in this study, we also used actual data representing transient booms and resurgent booms, and conducted parameter estimation for our proposed model using Bayesian inference. In addition, we conducted a model fitting to five actual data. By this study, we reconfirmed that we can express the resurgences or minute vibrations of actual data by means of our proposed model.

This paper investigates the inverse option problems (IOP) in the extended Black--Scholes model arising in financial markets. We identify the volatility and the drift coefficient from the measured data in financial markets using a Bayesian inference approach, which is presented as an IOP solution. The posterior probability density function of the parameters is computed from the measured data. The statistics of the unknown parameters are estimated by a Markov Chain Monte Carlo (MCMC) algorithm, which exploits the posterior state space. The efficient sampling strategy of the MCMC algorithm enables us to solve inverse problems by the Bayesian inference technique. Our numerical results indicate that the Bayesian inference approach can simultaneously estimate the unknown trend and volatility coefficients from the measured data.

In this paper, the author obtains an analytic exact form of Heaviside function, which is also known as Unit Step function and constitutes a fundamental concept of the Operational Calculus.In particulat, this function is explicitly expressed in a very simple manner by the aid of purely algebraic representations. The novelty of this work is that the proposed explicit formula is not performed in terms of non – elementary special functions, e.g. Dirac delta function or Error function and also is neither the limit of a function, nor the limit of a sequence of functions with point wise or uniform convergence. Hence, it may be much more appropriate and useful in the computational procedures which are inserted into Operational Calculus techniques and other engineering practices.

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In the first part of this two-paper series, we proposed and evaluated the theory of the 2D discrete Fourier Transform (DFT) in polar coordinates. The theory of the actual manipulated quantities was shown, including the standard set of shift, modulation, multiplication, and convolution rules. In this second part of the series, we address the computational aspects of the 2D DFT in polar coordinates. Specifically, we demonstrate how the decomposition of the 2D DFT as a DFT, Discrete Hankel Transform (DHT) and inverse DFT sequence can be exploited for efficient code. We also demonstrate how the proposed 2D DFT can be used to approximate the continuous forward and inverse Fourier transform in polar coordinates in the same manner that the 1D DFT can be used to approximate its continuous counterpart.

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In this paper, we propose and evaluate the theory of the 2D discrete Fourier Transform (DFT) in polar coordinates. This discrete theory is shown to arise from discretization schemes that have been previously employed with the 1D DFT and the discrete Hankel Transform (DHT). The proposed transform possesses orthogonality properties, which leads to invertibility of the transform. In the first part of this two-part paper, the theory of the actual manipulated quantities is shown, including the standard set of shift, modulation, multiplication, and convolution rules. Parseval and modified Parseval relationships are shown, depending on which choice of kernel is used. Similar to its continuous counterpart, the 2D DFT in polar coordinates is shown to consist of a 1D DFT, DHT and 1D inverse DFT.

The key issue in the post-genomic era is how to systematically describe the association between small molecule transcripts or translations inside cells. With the rapid development of high-throughput “omics” technologies, the achieved ability to detect and characterize molecules with other molecule targets opens up the possibility of investigating the relationships between different molecules from a global perspective. In this article, a Molecular Associations Network(MAN) is constructed and comprehensively analyzed by integrating the associations among miRNA, lncRNA, protein, drug, and disease, in which any kind of potential associations can be predicted. More specifically, each node in MAN can be represented as a vector by combining two kinds of information including the attributes of the node itself (e.g. sequences of ncRNAs and proteins, semantics of diseases and molecular fingerprints of drugs) and the manner of the node in the complex network (associations with other nodes). Random Forest classifier is trained to classify and predict new interactions or associations between biomolecules. In the experiment, the proposed method achieves a superb performance with 0.9735 AUC in 5-fold cross-validation, which show that the proposed method can provide new insight for exploration of the molecular mechanisms of disease and valuable clues for disease treatment.

The feedback collection and analysis has remained an important subject matter since long. The traditional techniques for student feedback analysis are based on questionnaire-based data collection and analysis. However, the student expresses their feedback opinions on online social media sites, which need to be analyzed. This study aims at the development of fuzzy-based sentiment analysis system for analyzing student feedback and satisfaction by assigning proper sentiment score to opinion words and polarity shifters present in the input reviews. Our technique computes the sentiment score of student feedback reviews and then applies fuzzy-logic module to analyze and quantify student’s satisfaction at the fine-grained level. The experimental results reveal that the proposed work has outperformed the baseline studies as well as state-of-the-art machine learning classifiers.

In this study, we investigate the form of solutions, stability character and asymptotic behavior of the following four rational difference equations x_{n+1} = (1/(x_{n}(x_{n-1}±1)±1)), x_{n+1} = ((-1)/(x_{n}(x_{n-1}±1)∓1)), such that their solutions are associated with Tribonacci numbers.

This paper deals with the form, the stability character, the periodicity and the global behavior of solutions of the following four rational difference equations x_{n+1} = ((±1)/(x_{n}(x_{n-1}±1)-1)) x_{n+1} = ((±1)/(x_{n}(x_{n-1}∓1)+1)).

Nowadays industrial dryers are used instead of traditional methods for drying. In designing dryers suitable for controlling the process of drying and reaching a high quality product, it is necessary to predict the instantaneous moisture loss during drying. For this purpose, ten mathematical-experimental models with a neural network model based on the kinetic data of pistachio drying are studied. The data obtained from the cabinet dryer will be evaluated at four temperatures of inlet air and different air velocities. The pistachio seeds will be placed in a thin layer on an aluminum sheet on a drying tray and weighed by a scale attached to the computer at different times. In the neural network, data are divided into three parts: educational (60%), validation (20%) and test (20%). Finally, the best mathematical-experimental model using genetic algorithm and the best neural network structure for predicting instantaneous moisture are selected based on the least squared error and the highest correlation coefficient.

One challenging problem is the representation of three-dimensional datasets that vary with time. These datasets can be though as a cloud of points that gradually deforms. But point-wise variations lack of information about the overall deformation pattern, and more importantly, about the extreme deformation locations inside the cloud. The present article applies a technique in computational mechanics to derive the strain-rate state of a time-dependent and three-dimensional data distribution, by which one can characterize its main trends of shift. Indeed, the tensorial analysis methodology is able to determine the global deformation rates in the entire dataset. With the use of this technique, one can characterize the significant fluctuations in a reduced multivariate description of an urban system and identify the possible causes of those changes: calculating the strain-rate state of a PCA-based multivariate description of an urban system, we are able to describe the clustering and divergence patterns between the districts of the city and to characterize the temporal rate in which those variations happen.

Smartphone devices, particularly android devices used by billions of people everywhere in the world. Similarly, this increasing rate attracts mobile botnet attackers that is a network of interconnected nodes operated by command and control (C&C) method to expand malicious activities. At present, mobile botnet attacks carried Distributed denial of services (DDoS) that causes to steal sensitive data, remote access, spam generation, etc. Consequently, various approaches are defined in the literature to detect mobile botnet using static or dynamic analysis. In this paper, we have proposed a novel hybrid model, which is a combination of static and dynamic method that relies on machine learning to identify android botnet applications having C&C capability. The results evaluated through machine learning classifiers in which Random forest classifier outperform other ML techniques, i.e. Naïve Bayes, Support Vector Machine, and Simple logistics. Our proposed framework can achieve 97.48% accuracy in detecting such harmful applications. Furthermore, we highlight some research directions and possible solutions regarding botnet attacks for the entire community.

In this paper, we propose a new mathematical model based on the association between susceptible and recovered individual, where the association between susceptible and recovered individual is disturbed by white noise. This model is based on demographic changes and is used for long term behavior. We study the stability of equilibria of the deterministic model and prove the conditions for the extinction of diseases. Then, we investigate and obtain the critical condition of the stochastic epidemic model for the extinction and the permanence in mean of the disease with the white noise. To verify our results, we present some numerical simulations for real data related to disease.

Causal relationships lie at the very core of our scientific description of biophysical phenomena.  Nevertheless, observable facts such as chaotic dynamics, symmetry breaks, long-range collisionless neural interactions, zero-value energy singularities, particle/wave duality, ecological approaches to visual perception elude cause effect inductive explanations.  We illustrate how simple topological claims, seemingly useless and far away from scientific inquiry  (e.g., “given at least some wind on Earth, there must at all times be a cyclone or anticyclone somewhere”; “if one stirs to dissolve a lump of sugar in a cup of coffee, it appears there is always a point without motion”; “at any moment, there is always a pair of antipodal points on the Earth’s surface with equal temperatures and barometric pressures”) describe a real world, made of projections and mappings unsupplied with causality.  Unavoidable, necessary, passive, spontaneous changes in systems’ shapes, dimensions and symmetries incarnate biophysical quantifiable counterparts of abstract mathematical concepts.  Real features arising from topological mechanisms do not require causal powers to be produced, because are independent of exerted physical forces, energetic fields’ actions and thermodynamic constraints.  Providing examples that shed novel interpretative light on phenomena such as double slit experiments, cellular mechanisms and brain function, we conclude as follows: countless scientific matters can be accessed through methodological weapons which not just break causality, replace temporal correlations and probabilistic a priori knowledge of previous cases, but also permit augments in available information devoid of entropy increases and energy losses.

The aim of this expository article is to present recent developments in the centuries-old discussion on the interrelations between several types of domination in graphs. However, the novelty even more prominent in the newly discovered simplified presentations of several older results. Domination can be seen as arising from real-world application and extracting classical results as first described by this article.The main part of this article, concerning a new domination and older one, is presented in a narrative that answers two classical questions: (i) To what extend must closing set be dominating? (ii) How strong is the assumption of domination of a closing set? In a addition, we give an overview of the results concerning domination. The problem asks how small can a subset of vertices be and contain no edges or, more generally how can small a subset of vertices be and contain other ones. Our work was as elegant as it was unexpected being a departure from the tried and true methods of this theory that had dominated the field for one fifth a century. This expository article covers all previous definitions. The inability of previous definitions in solving even one case of real-world problems due to the lack of simultaneous attentions to the worthy both of vertices and edges causing us to make the new one. The concept of domination in a variety of graphs models such as crisp, weighted and fuzzy, has been in a spotlight. We turn our attention to sets of vertices in a fuzzy graph that are so close to all vertices, in a variety of ways, and study minimum such sets and their cardinality. A natural way to introduce and motivate our subject is to view it as a real-world problem. In its most elementary form, we consider the problem of reducing waste of time in transport planning. Our goal here is to first describe the previous definitions and the results, and then to provide an overview of the flows ideas in their articles. The final outcome of this article is twofold: (i) Solving the problem of reducing waste of time in transport planning at static state; (ii) Solving and having a gentle discussions on problem of reducing waste of time in transport planning at dynamic state. Finally, we discuss the results concerning holding domination that are independent of fuzzy graphs. We close with a list of currently open problems related to this subject. Most of our exposition assumes only familiarity with basic linear algebra, polynomials, fuzzy graph theory and graph theory.

Present study offers mathematical calculations of the roll-coating procedure lubricated with upper Convected Maxwell Fluid. An incompressible isothermal viscoelastic fluid is considered with roll and the porous web having uniform velocities. By employing Lubrication Approximation Theory, the desired equations of motion for the fluid concerned over porous web are modelled and analyzed. The suction rate on the web and injection rate at the roll surface are anticipated proportionate. Results for velocity profile and pressure gradient are obtained analytically. Fluid parameters of industrial significance i.e. detachment point, pressure, sheet /roll separating force, power contribution and coating thickness are also calculated numerically. Substantial and monotonic increase is witnessed in these quantities with the increase of flow parameters.

We consider a control problem for a diffusive PDE model of heterogeneous population of thermostatically controlled loads (TCLs) aiming to balance the aggregate power consumption within a given amount of time. Using the Green’s function approach, the problem is formulated as an approximate controllability problem for a residue depending on control parameters nonlinearly. A sufficient condition for approximate controllability is derived in terms of initial temperature distribution, operation time of TCLs and threshold value of the aggregate power consumption. Case studies allow to reveal the advantages of the proposed solution from numerical calculations point of view.

The aim of this expository article is to present recent developments in the centuries-old discussion on the interrelations between several types of domination in graphs. However, the novelty even more prominent in the newly discovered simplified presentations of several older results. The main part of this article, concerning a new domination and older one, is presented in a narrative that answers two classical questions: (i) To what extend must closing set be dominating? (ii) How strong is the assumption of domination of a closing set? In a addition, we give an overview of the results concerning domination. The problem asks how small can a subset of vertices be and contain no edges or, more generally how can small a subset of vertices be and contain other ones. Our work was as elegant as it was unexpected being a departure from the tried and true methods of this theory that had dominated the field for one fifth a century. This expository article covers all previous definitions. The inability of previous definitions in solving even one case of real-world problems due to the lack of simultaneous attentions to the worthy both of vertices and edges causing us to make the new one. The concept of domination in a variety of graphs models such as crisp, weighted and fuzzy, has been in a spotlight. We turn our attention to sets of vertices in a fuzzy graph that are so close to all vertices, in a variety of ways, and study minimum such sets and their cardinality. A natural way to introduce and motivate our subject is to view it as a real-world problem. In its most elementary form, we consider the problem of reducing waste of time in transport planning. Our goal here is to first describe the previous definitions and the results, and then to provide an overview of the flows ideas in their articles. The final outcome of this article is twofold: (i) Solving the problem of reducing waste of time in transport planning at static state; (ii) Solving and having a gentle discussions on problem of reducing waste of time in transport planning at dynamic state. Finally, we discuss the results concerning holding domination that are independent of fuzzy graphs. We close with a list of currently open problems related to this subject. Most of our exposition assumes only familiarity with basic linear algebra, polynomials, fuzzy graph theory and graph theory.

A non-Newtonian fluid model is used to investigate the two-dimensional pulsatile blood flow through a tapered artery with a stenosis. The mixed convection effects of heat and mass transfer are also taken into account. By applying non-dimensionalization and radial coordinate transformation, we simplify the system in a tube. Under the finite difference scheme, numerical solutions are calculated for velocity, temperature concentration, resistance, impedance, wall shear stress and shearing stress at the stenosis throat. Finally, Quantitative analysis is carried out.

In the article, by applied the concept of strongly convex function and one known identity, we establish several Ostrowski type inequalities involving conformable fractional integrals. As applications, some new error estimations for the midpoint formula are provided as well.

In this work, we establish some new Hermite-Hadamard type inequalities for convex functions via conformable fractional integral. Moreover, we show that through the conformable fractional integral we can find some new Hermite-Hadamard type inequalities for convex functions via the classical integrals.

Using an algebraic method for solving the wave equation in quantum mechanics, we encountered a new class of orthogonal polynomials on the real line. One of these is a four-parameter polynomial with a discrete spectrum. Another that appeared while solving a Heun-type equation has a mix of continuous and discrete spectra. Based on these results and on our recent study of the solution space of an ordinary differential equation of the second kind with four singular points, we introduce a modification of the hypergeometric polynomials in the Askey scheme. Up to now, all of these polynomials are defined only by their three-term recursion relations and initial values. However, their other properties like the weight function, generating function, orthogonality, Rodrigues-type formula, etc. are yet to be derived analytically. This is an open problem in orthogonal polynomials.

In this paper, we study the well-posedness and asymptotics of a one-dimensional thermoelastic laminated beam system either with or without structural damping, where the heat conduction is given by Fourier's law effective in the rotation angle displacements. We show that the system is well-posed by using Lumer-Philips theorem, and prove that the system is exponentially stable if and only if the wave speeds are equal, by using the perturbed energy method and Gearhart-Herbst-Prüss-Huang theorem. Furthermore, we show that the system with structural damping is polynomially stable provided that the wave speeds are not equal, by using the second-order energy method.

In this paper, we work with a diffusion-perturbed risk model comprising a surplus generating process and an investment return process. The investment return process is of standard Black-Scholes type, that is, it comprises a single risk-free asset that earns interest at a constant rate and a single risky asset whose price process is modelled by a geometric Brownian motion. Additionally, the company is allowed to purchase noncheap proportional reinsurance priced via the expected value principle. Using the Hamilton-Jacobi-Bellman approach, we derive a second-order Volterra integrodifferential equation which we transform into a linear Volterra integral equation of the second kind. We proceed to solve this integral equation numerically using the block-by-block method for the optimal reinsurance retention level that minimizes the ultimate ruin probability. The numerical results based on light- and heavy-tailed distributions show that proportional reinsurance and investments play a vital role in enhancing the survival of insurance companies. But the ruin probability exhibits sensitivity to the volatility of the stock price.

This article deals with the conformable double Laplace transforms and their some properties with examples and also the existence Condition for the conformable double Laplace transform is studied. Finally, in order to obtain the solution of nonlinear fractional problems, we present a modified conformable double Laplace that we call conformable double Laplace decomposition methods (CDLDM). Then, we apply it to solve, Regular and singular conformable fractional coupled burgers equation illustrate the effectiveness of our method some examples are given.

In this paper, we work with a diffusion-perturbed risk model comprising a surplus generating process and an investment return process. The investment return process is of standard Black-Scholes type, that is, it comprises a single risk-free asset that earns interest at a constant rate and a single risky asset whose price process is modelled by a geometric Brownian motion. Additionally, the company is allowed to purchase noncheap proportional reinsurance priced via the expected value principle. Using the Hamilton-Jacobi-Bellman approach, we derive a second-order Volterra integrodifferential equation which we transform into a linear Volterra integral equation of the second kind. We proceed to solve this integral equation numerically using the block-by-block method for the optimal reinsurance retention level that minimizes the ultimate ruin probability. The numerical results based on light- and heavy-tailed distributions show that proportional reinsurance and investments play a vital role in enhancing the survival of insurance companies. But the ruin probability exhibits sensitivity to the volatility of the stock price.

An SIRS (Susceptible–Infected–Removed-Susceptible) mathematical model for the transmission dynamics of the Transfusion–Transmitted Malaria (TTM) model with optimal control pair u₁(t) and u₂(t) was developed and studied in this research work. The model Transfusion–Transmitted Malaria disease–free equilibrium and endemic equilibriums points were determined. The model exhibited two equilibriums; disease-free and endemic equilibrium. It is shown that the disease–free equilibrium was locally asymptotically stable if the associated basic reproduction numbers R₀ is less than unity while the disease persists if R₀ is greater than unity. The global stability of the Transfusion–Transmitted Malaria model at the disease-free equilibrium was established using the comparison method. The optimality system was derived and an optimal control model of blood screening and drug treatment for the Transfusion–Transmitted Malaria model was investigated. Conditions for the optimal control were considered using Pontryagin’s Maximum Principle and solved numerically using the Forward and Backward Finite Difference Method (FBDM). Numerical results obtained are in perfect agreement with our analytical results.

We investigate the exact travelling wave solutions of the Tzitzéica-Dodd–Bullough–Mikhailov (TDBM) and Tzitzéica-type (TT) equations, appear in nonlinear optics, using simple ansatz approach. We use the Painlevé and traveling wave transformations to convert the mentioned nonlinear equations into nonlinear ordinary differential equations. Then, the simple ansatz approach has been employed to procure wave solutions, in terms of bright, singular soliton and trigonometric functions, for the Tzitzéica-Dodd–Bullough–Mikhailov (TDBM) and Tzitzéica-type (TT) equations. The study shows that the simple ansatz approach is straightforward and easy to implement for solving nonlinear problems.

Many various using of this new-born fuzzy model for solving real-world problems and urgent requirements involve introducing new concept for analyzing the situations which leads to solve them by proper, quick and ecient method based on statistical data. This gap between the model and its solution cause that we introduce nikfar domination in neutrosophic graphs as creative and eective tool for studying a few selective vertices of this model instead of all ones by using special edges. Being special selection of these edges aect to achieve quick and proper solution to these problems. Domination hasn't ever been introduced. So we don't have any comparison with another denitions. The most used graphs which have properties of being complete, empty, bipartite, tree and like stu and they also achieve the names for themselves, are studied as fuzzy models for getting nikfar dominating set or at least becoming so close to it. We also get the relations between this special edge which plays main role in doing dominating with other special types of edges of graph like bridges. Finally, the relation between this number with other special numbers and characteristic of graph like order are discussed.

We introduce a new variation on the domination theme which we call nikfar domination as reducing waste of time in transportation planning. We determine the nikfar domination number  v for several classes of fuzzy graphs. The bounds is obtained for it. We prove both of the Vizing's conjecture and the Grarier-Khelladi's conjecture are monotone decreasing fuzzy graph property for nikfar domination. We obtain Nordhaus-Gaddum (NG) type results for these parameters. Finally, we discuss about nikfar dominating set of a fuzzy tree by using the bridges and -strong edges equivalence.

A method of solution to solve the compressible unsteady 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity equation in cylindrical coordinates is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the boundary layer of the tube. A reduction to a single partial differential equation is possible and integral calculus methods are applied for the case of a body force directed to the centre of the tube to obtain an integral form of the Hunter-Saxton equation. Also an extension for a more general body force is shown where in addition there is a rotational force applied.

The manuscript surveys the special functions of the Fox-Wright type. These functions are generalizations of the hypergeometric functions. Notable representatives of the type are the Mittag-Leffler functions and the Wright function. The integral representations of such functions are given and the conditions under which these function can be represented by simpler functions are demonstrated. The connection with generalized fractional differential and integral operators is demonstrated and discussed.

Exact and soliton type solutions have great importance in propagation of surface waves, fluid dynamics, optics, and many other elds of nonlinear sciences. In this study, the explicit and exact soliton type solutions for two space-time fractional Equal- Width (FEW) equations with conformable derivative are procured via the hyperbolic function approach. The wave type solutions are represented in some hyperbolic and trigonometric functions.

DGJ (Daftardar-Gejii-Jafaris) method is used to obtain numerical solution of the third order fractional differential equation. Providing the DGJ method converges, the approximate solution is a good and effective numerical result which is close to the exact solution or the exact solution. For this,the examples of the explaning the method are presented. The proposed method is implemented for the approximation solution of the third order nonlinear fractional partial differential equations. The method was shown to be unsuitable and inconsistent for an example of a nonlinear fractional partial differential equation depend on initial-boundary value conditions. The fact that these numerical results are not consistent can be explained by the fact that the method is not convergent.

A method of solution to solve the compressible 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity equation in cylindrical coordinates is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the boundary layer of the tube. It is concluded that the total divergence of the flow can be expressed as the integral with respect to time of the line integral of the dot product of inertial and azimuthal velocity. The line integral is evaluated on a contour that is annular and traces the boundary layer as time increases in the flow. Finally it's explicit dependence on the gradient of a density function and Laplacian of pressure is shown with complete dependence on the vorticity and it's rate of change for a flow associated with the Navier Stokes equations.

The paper is devoted to the investigation of regional observability with constraints on the state of semilinear parabolic systems. The purpose is to reconstruct the initial state between two prescribed functions only on an internal subregion $\omega$ of the system evolution domain $\Omega$. The proofs use two approaches, the subdifferential and HUM approach. Finally, a numerical example is provided to verify the effectiveness of our theory results.

In this paper, Lie symmetry and Painlevé Techniques are applied to the SIRD (Susceptible, Infected, Recovered and Dead) model. A demonstration of the integrability of the model is provided to present an exact solution. The study revealed that nonlinear system passes Painlevé test and does not possesses complex chaotic behaviour. However, the system fails to pass the Painlevé test while constraints reach values equivalent to the corresponding complex chaotic behaviour. The two-dimensional Lie symmetry algebra and the commutator table of the infinitesimal generators are obtained. Lie symmetry analysis serves to linearize the nonlinear system and find the corresponding invariant solution.

We apply utilized the extended form of the auxiliary equation method to obtain extensively reliable exact travelling wave solutions of perturbed Gerdjikov–Ivanov equation (GIE) that is widely used as a model in the field theory of quanta and non-linear optics. The method is based on a simple first order second degree ODE. The new form of the approach gives more solutions to the governing equation efficiently.

By using the new recently proposed extended rational methods, we set some solitary and multi wave solutions of the coupled Higgs system in the present study. Some new families of hyperbolic and trigonometric solutions of the coupled Higgs equations are successfully obtained. Some solutions are real as some are complex valued functions. Particular forms of some solutions derived by choice of some parameters are demonstrated in three dimensional spaces. The complex valued solutions are also depicted by plotting both real and complex component.

A new mathematical model is presented to study the effects of macrophages on the bone fracture healing process. The model consists of a system of nonlinear ordinary differential equations that represents the interactions among classically and alternatively activated macrophages, mesenchymal stem cells, osteoblasts, and pro- and anti-inflammatory cytokines. A qualitative analysis of the model is performed to determine the equilibria and their corresponding stability properties. Numerical simulations are also presented to support the theoretical results and to monitor the evolution of a broken bone for different types of fractures under various medical interventions. The model can be used to guide clinical experiments and to explore possible medical treatments that accelerate the bone fracture healing process either by surgical interventions or drug administrations.

In this paper, by extending some results of informational genomics, we present a new randomness test based on the empirical entropy of strings and some properties of repeatability and unrepeatability of substrings of certain lengths. We give the theoretical motivations of our method and some experimental results of its application to a wide class of strings: decimal representations of real numbers, roulette outcomes, logistic maps, linear congruential generators, quantum measurements, natural language texts, and genomes. It will be evident that the evaluation of randomness resulting from our tests does not distinguish among the different sources of randomness (natural, or pseudo-casual).

Information-related measures are useful tools for multi-variable data analysis, as measures of dependence among variables, and as descriptions of order and disorder in biological and physical systems.  Measures, like marginal entropies, mutual / interaction / multi -information, have long been used in a number of fields including descriptions of systems complexity and biological data analysis.  The mathematical relationships among these measures are therefore of significant inherent interest.  Relations between common information measures include the duality relations based on Möbius inversion on lattices.  These are the direct consequence of the symmetries of the lattices of the sets of variables (subsets ordered by inclusion).  While these relationships are of significant interest there has been, to our knowledge, no systematic examination of the full range of relationships of this diverse range of functions into a unifying formalism as we do here.  In this paper we define operators on functions on these lattices based on the Möbius inversions that map functions into one another (Möbius operators).  We show that these operators form a simple group isomorphic to the symmetric group S₃.  Relations among the set of functions on the lattice are transparently expressed in terms of the operator algebra, and, applied to the information measures, can be used to derive a wide range of relationships among diverse information measures.  The Möbius operator algebra is naturally generalized which yields extensive new relationships.  This formalism now provides a fundamental unification of information-related measures, and the isomorphism of all distributive lattices with the subset lattice implies an even broader application of these results.

Pollution of sub-surface water reservoirs mainly rivers and streams through contaminated water point sources (CWPS) was studied. The objective was to formulate a discrete time delay mathematical model which describes the dynamics of reservoir pollution using mixing-problem processes that involve single species contaminants such as nitrates, phosphorous and detergents. The concentration  of pollutants was expressed as a function of the inflow and outflow rates using the principle for the conservation of mass. Systems of ODEs generated from principles of mixing problems were refined into a system of DDEs so that the concentration of pollutant leaving the reservoir at time would be determined at some earlier instant, for the delay. The formulated model is a mathematical discrete time delay model which would be used to describe the dynamics of sub-surface water reservoir pollution. The results from the validation of the model were analyzed   to determine how time delays in the mixing processes affect the rate of particle movement in water reservoirs.

A new type of ordinary differential equation is introduced and discussed, namely, the time-dependent order ordinary differential equations. These equations can be solved via fractional calculus and are mapped into Volterra integral equations of second kind with singular integrable kernel. The solutions of the time-dependent order differential equations smoothly deforms solutions of the classical integer order ordinary differential equations into one-another, and can generate or remove singularities. An interesting symmetry of the solution in relation to the Riemann zeta function and Harmonic numbers was also proved.

In this paper, we established the existence of PC-mild solutions for non local fractional impulsive functional integro-differential equations with finite delay. The proofs are obtained by using the techniques of fixed point theorems, semi-group theory and generalized Bellman inequality. In this paper, we have used the distributed characteristic operators to define the mild solution of the system. Results obtained here improve and extend some known results.

The Colebrook equation is a popular model for estimating friction loss coefficients in water and gas pipes. The model is implicit in the unknown flow friction factor . To date, the captured flow friction factor  can be extracted from the logarithmic form analytically only in the term of the Lambert -function. The purpose of this study is to find an accurate and computationally efficient solution based on the shifted Lambert -function also known as the Wright -function. The Wright -function is more suitable because it overcomes the problem with the overflow error by switching the fast growing term  of the Lambert -function to the series expansions that further can be easily evaluated in computers without causing overflow run-time errors. Although the Colebrook equation transformed through the Lambert -function is identical to the original expression in term of accuracy, a further evaluation of the Lambert -function can be only approximate. Very accurate explicit approximations of the Colebrook equation that contains only one or two logarithms are shown. The final result is an accurate explicit approximation of the Colebrook equation with the relative error of no more than 0.0096%. The presented approximations are in the form suitable for everyday engineering use, they are both accurate and computationally efficient.

In this article, we consider the three dimensional $\alpha$-fractional nonlinear delay differential system of the form \begin{align*} D^{\alpha}\left(u(t)\right)&=p(t)g\left(v(\sigma(t))\right),\\D^{\alpha}\left(v(t)\right)&=-q(t)h\left(w(t))\right),\\D^{\alpha}\left(w(t)\right)&=r(t)f\left(u(\tau(t))\right),~ t \geq t_0, \end{align*} where $0 < \alpha \leq 1$, $D^{\alpha}$ denotes the Katugampola fractional derivative of order $\alpha$. We have established some new oscillation criteria of solutions of differential system by using generalized Riccati transformation and inequality technique. The obtained results are illustrated with suitable examples.

In the paper we show that the general relativity in recent Einstein-Palatini formulation is equivalent to a gauge field. We begin with a bit of information of the Einstein-Palatini formulation and derive Einstein field equations from it. In the next section, we consider general relativity with a positive cosmological constant in terms of the corrected curvature. We show that in terms of the corrected curvature general relativity takes the form typical for a gauge field. Finally, we give a geometrical interpretation of the corrected curvature.

It is shown that if a non-constant meromorphic function f(z) is of finite order and shares certain values with its shifts/difference operators then f(z) coincides with that particular shift/difference operator.

Let G[F,Vk, Huv] be the graph with k pockets, where F is a simple graph of order n ≥ 1,Vk= {v1,v2,··· ,vk} is a subset of the vertex set of F and Hvis a simple graph of order m ≥ 2,v is a specified vertex of Hv. Also let G[F,Ek, Huv] be the graph with k edge pockets, where F is a simple graph of order n ≥ 2, Ek= {e1,e2,···ek} is a subset of the edge set of F and Huvis a simple graph of order m ≥ 3, uv is a specified edge of Huvsuch that Huv− u is isomorphic to Huv− v. In this paper, we derive closed-form formulas for resistance distance and Kirchhoff index of G[F,Vk, Hv] and G[F,Ek, Huv] in terms of the resistance distance and Kirchhoff index F, Hv and F, Huv, respectively.

In this paper, we investigate a new $q$-analogue of the higher order degenerate Changhee polynomials and numbers, which are called the Witt-type formula for the $q$-analogue of degenerate Changhee polynomials of order $r$. We can derive some new interesting identities related to the degenerate $(h,q)$-Changhee polynomials and numbers.

We investigate the value distribution of dierence polynomials of entire and meromorphic functions, which can be viewed as the Hayman's conjecture. And we also study the uniqueness of dierence polynomials sharing a common value.

This paper attempts to propose a Riemann spring model via an analogy between the Riemann Zeta function of a complex number and the elastic spring in series.

The generalized Kuramoto–Sivashinsky equation is investigated using the modified Kudrayshov method for the exact analytical solution. The modified Kudrayshov method converts the nonlinear partial differential equation to algebraic equations, as a result of various steps, which on solving the so obtained equation systems yields the analytical solution. By this way various exact solutions including complex structures are found and drawn their behaviour in complex plane by Maple to compare the uniqueness of solutions.

In this paper, we investigate an analytical solution of a class of nonlinear fractional integro-differential equation base on a generalized fractional power series expansion. The fractional derivatives are described in the conformable's type. The new approach is a modified form of the well-known Taylor series expansion. The illustrative examples are presented to demonstrate the accuracy and effectiveness of the proposed method

An SIR (Susceptible – Infected – Removed) mathematical model for the transmission dynamics of the Transfusion –Transmitted Malaria (TTM) model with optimal control pair and was developed and studied in this research work. The model Transfusion –Transmitted Malaria disease – free equilibrium and endemic equilibriums points were determined. The model exhibited two equilibriums; disease-free and endemic equilibrium. It was shown that the disease – free equilibrium was locally asymptotically stable if the associated basic reproduction numbers  is less than unity while the disease persists if  is greater than unity. The global stability of the Transfusion –Transmitted Malaria model at the disease – free equilibrium was established using the comparison method. The optimality system was derived and an optimal control model of blood screening and drug treatment for the Transfusion –Transmitted Malaria model was investigated. Conditions for the optimal control were considered using Pontryagin’s Maximum Principle and solved numerically using the Forward and Backward Finite Difference Method (FBDM). Numerical results obtained are in perfect agreement with our analytical results.

Malaria and Pneumonia are leading causes of serious illness in children and adults worldwide with their death rate and prevalence on the rise. Such alarming statistics may retard the milestones so far achieved in meeting the Millennium Development Goals 4 and 6 whose targets are to improve child survival and reverse the high prevalence of diseases such as pneumonia and malaria respectively. Two sub-models of malaria-pneumonia co-infection namely malaria model and pneumonia model were considered first and then followed by the full malaria-pneumonia co-infection model. The malaria model, pneumonia model and co-infection model basic reproduction numbers denoted by R_m, R_p and R_mp respectively was obtained using the Next Generation Matrix method. The model disease free equilibrium’s local and global stability was analysed using Descartes’ Rule of signs and Comparison method. The bifurcation analysis for the malaria, pneumonia and co-infection models was studied using the Centre Manifold Theory. The sensitivity indices of the model basic reproduction numbers R_m, R_p and R_mp to the parameters in the models were calculated. Optimal control theory was applied using the Pontryagins’ Maximum Principle to investigate optimal strategies for controlling the spread of malaria, pneumonia and co-infection models using insecticide treated bed nets (u₁(t)) spraying of mosquitoes insecticides (u₂(t)), sanitation (u₃(t)), vaccination (u₄(t)), anti-malaria drugs (u₅(t)), anti-pneumonia drugs (u₆(t)), both anti-malaria drugs and anti-pneumonia drugs (u₇(t)) as the system time control variables.  Numerical simulations using a set of parameter values were provided to validate the analytical results.