New setting is introduced to study the global powerful alliance. Global powerful alliance is about a set of vertices which are applied into the setting of neutrosophic graphs. Neighborhood has the key role to define this notion. Also, neighborhood is defined based on strong edges. Strong edge gets a framework as neighborhood and after that, too close vertices have key role to define global powerful alliance based on strong edges. The structure of set is studied and general results are obtained. Also, some classes of neutrosophic graphs excluding empty, path, star, and wheel and containing complete, cycle and r-regular-strong are investigated in the terms of set, minimal set, number, and neutrosophic number. Neutrosophic number is used in this way. It’s applied to use the type of neutrosophic number in the way that, three values of a vertex are used and they’ve same share to construct this number. It’s called “modified neutrosophic number”. Summation of three values of vertex makes one number and applying it to a set makes neutrosophic number of set. This approach facilitates identifying minimal set and optimal set which forms minimal-global-powerful-alliance number and minimal-global-powerful-alliance-neutrosophic number. Two different types of sets namely global-powerful alliance and minimal-global-powerful alliance are defined. Global-powerful alliance identifies the sets in general vision but minimal-global-powerful alliance takes focus on the sets which deleting a vertex is impossible. Minimal-global-powerful-alliance number is about minimum cardinality amid the cardinalities of all minimal-global-powerful alliances in a given neutrosophic graph. New notions are applied in the settings both individual and family. Family of neutrosophic graphs has an open avenue, in the way that, the family only contains same classes of neutrosophic graphs. The results are about minimal-global-powerful alliance, minimal-global-powerful-alliance number and its corresponded sets, minimal-global-powerful-alliance-neutrosophic number and its corresponded sets, and characterizing all minimal-global-powerful alliances, minimal-t-powerful alliance, minimal-t-powerful-alliance number and its corresponded sets, minimal-t-powerful-alliance-neutrosophic number and its corresponded sets, and characterizing all minimal-t-powerful alliances. The connections amid t-powerful-alliances are obtained. The number of connected components has some relations with this new concept and it gets some results. Some classes of neutrosophic graphs behave differently when the parity of vertices are different and in this case, cycle, and complete illustrate these behaviors. Two applications concerning complete model as individual and family, under the titles of time table and scheduling conclude the results and they give more clarifications and closing remarks. In this study, there’s an open way to extend these results into the family of these classes of neutrosophic graphs. The family of neutrosophic graphs aren’t study deeply and with more results but it seems that analogous results are determined. Slight progress is obtained in the family of these models but there are open avenues to study family of other models as same models and different models. There’s a question. How can be related to each other, two sets partitioning the vertex set of a graph? The ideas of neighborhood and neighbors based on strong edges illustrate open way to get results. A set is global powerful alliance when two sets partitioning vertex set have uniform structure. All members of set have more amount of neighbors in the set than out of set and reversely for non-members of set with less members in the way that the set is simultaneously t-offensive and(t-2)-defensive. A set is global if t=0. It leads us to the notion of global powerful alliance. Different edges make different neighborhoods but it’s used one style edge titled strong edge. These notions are applied into neutrosophic graphs as individuals and family of them. Independent set as an alliance is a special set which has no neighbor inside and it implies some drawbacks for these notions. Finding special sets which are well-known, is an open way to purse this study. Special set which its members have only one neighbor inside, characterize the connected components where the cardinality of its complement is the number of connected components. Some problems are proposed to pursue this study. Basic familiarities with graph theory and neutrosophic graph theory are proposed for this article.

Saif et al. (J. Math. Comput. Sci. 21 (2020) 127-135) considered modified Laplace transform and developed some of their certain properties and relations. Motivated by this work, in this paper, we define modified Sumudu transform and investigate many properties and relations including modified Sumudu transforms of the power function, sine, cosine, hyperbolic sine, hyperbolic cosine, exponential function, and function derivatives. Moreover, we attain two shifting properties and a scale preserving theorem for the modified Sumudu transform. We give modified inverse Sumudu transform and investigate some relations and examples. Furthermore, we show that the modified Sumudu transform is the theoretical dual transform to the modified Laplace transform.

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.

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.

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.

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.

The Thermodynamic Formalism provides a rigorous mathematical framework to study quantitative and qualitative aspects of dynamical systems. At its core there is a variational principle and corresponding, in its simplest form, to the Maximum Entropy principle, used as a statistical inference procedure to represent, by specific probability measures (Gibbs measures), the collective behaviour of complex systems. This framework has found applications in different domains of scienThe Thermodynamic Formalism provides a rigorous mathematical framework to study quantitative and qualitative aspects of dynamical systems. At its core there is a variational principle and corresponding, in its simplest form, to the Maximum Entropy principle, used as a statistical inference procedure to represent, by specific probability measures (Gibbs measures), the collective behaviour of complex systems. This framework has found applications in different domains of science, in particular, has been fruitful and influential in neurosciences. In this article, we review how the Thermodynamic Formalism can be exploited in the field of theoretical neuroscience, as a conceptual and operational tool, to link the dynamics of interacting neurons and the statistics of action potentials from either experimental data or mathematical models. We comment on perspectives and open problems in theoretical neuroscience that could be addressed within this formalism.ce, in particular, has been fruitful and influential in neurosciences. In this article, we review how the Thermodynamic Formalism can be exploited in the field of theoretical neuroscience, as a conceptual and operational tool, to link the dynamics of interacting neurons and the statistics of action potentials from either experimental data or mathematical models. We comment on perspectives and open problems in theoretical neuroscience that could be addressed within this formalism.

Increasing rate of seepage velocity from several formation characteristics, such as permeability and porosity, in water aquifer environment greatly prompt pollution of water reservoirs within a short period of time. Considerably, migration rate of dissolved heavy metals from Solid Waste Dumpsites (SWD), such as municipal dumpsites and landfills, through heterogeneous aquifer environment, and finally into nearby water reservoirs are mainly influenced by variation of seepage velocity within the soil and water environment. This presents a dynamic system for water pollution that was studied using a formulated mathematical model to describe the transport process of dissolved heavy metals, mainly characterized by seepage velocities, within the water aquiferous environment. Permeability, porosity, fluid pressure and concentration of heavy metals in aquiferous environment were used as principal parameters that influence seepage velocity of the metals, in dissolved state, through the structural formation of water aquifers. The derived mathematical equations that constitute the model of this study were generated through Darcy’s law and the equation of continuity. The model was validated on structural river aquifer sediments, and it was solved using graphical method through matlab open-source software. The initial and boundary conditions were obtained by discretizing the geological setting of flow region so as to transform the gradient of the head, into the time domain.

An oil and vinegar scheme is a signature scheme based on multivariate quadratic polynomials over finite fields. The system of polynomials contains $n$ variables, divided into two groups: $v$ vinegar variables and $o$ oil variables. The scheme is called balanced (OV) or unbalanced (UOV), depending on whether $v = 0$ or not, respectively. These schemes are very fast and require modest computational resources, which make them ideal for low-cost devices such as smart cards. However, the OV scheme has been already proven to be insecure and the UOV scheme has been proven to be very vulnerable for many parameter choices. In this paper, we propose a new multivariate public key signature whose central map consists of a set of polynomials obtained from the multiplication of block matrices. Our construction is motivated by the design of the Simple Matrix Scheme for Encryption and the UOV scheme. We show that it is secure against the Separation Method, which can be used to attack the UOV scheme, and against the Rank Attack, which is one of the deadliest attacks against multivariate public-key cryptosystems. Some theoretical results on matrices with polynomial entries are also given, to support the construction of the scheme.

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.

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.

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.

New setting is introduced to study neutrosophic failed-independent number and failed independent neutrosophic-number arising neighborhood of different vertices. Neighbor is a key term to have these notions. Having all possible edges amid vertices in a set is a key type of approach to have these notions namely neutrosophic failed-independent number and failed independent neutrosophic-number. Two numbers are obtained but now both settings leads to approach is on demand which is finding biggest set which have all vertices which are neighbors. Let NTG : (V,E,σ,μ) be a neutrosophic graph. Then failed independent number I(NTG) for a neutrosophic graph NTG : (V,E,σ,μ) is maximum cardinality of a set S of vertices such that every two vertices of S are endpoints for an edge, simultaneously; failed independent neutrosophic-number In(NTG) for a neutrosophic graph NTG : (V,E,σ,μ) is maximum neutrosophic cardinality of a set S of vertices such that every two vertices of S are endpoints for an edge, simultaneously. As concluding results, there are some statements, remarks, examples and clarifications about some classes of neutrosophic graphs namely path-neutrosophic graphs, cycle-neutrosophic graphs, complete-neutrosophic graphs, star-neutrosophic graphs, complete-bipartite-neutrosophic graphs and complete-t-partite-neutrosophic graphs. The clarifications are also presented in both sections “Setting of Neutrosophic Failed-Independent Number,” and “Setting of Failed Independent Neutrosophic-Number,” for introduced results and used classes. Neutrosophic number is reused in this way. It’s applied to use the type of neutrosophic number in the way that, three values of a vertex are used and they’ve same share to construct this number to compare with other vertices. Summation of three values of vertex makes one number and applying it to a comparison. This approach facilitates identifying vertices which form neutrosophic failed-independent number and failed independent neutrosophic-number arising neighborhoods of vertices. In path-neutrosophic graphs, two neighbors, form maximal set but with slightly differences, in cycle-neutrosophic graphs, two neighbors forms maximal set. Other classes have same approaches. In complete-neutrosophic graphs, a set of all vertices leads us to neutrosophic failed-independent number and failed independent neutrosophic-number. In star-neutrosophic graphs, a set of vertices containing only center and one other vertex, makes maximal set. In complete-bipartite-neutrosophic graphs, a set of vertices including two vertices from different parts makes intended set but with slightly differences, in complete-t-partite-neutrosophic graphs, a set of t vertices from different parts makes intended set. In both settings, some classes of well-known neutrosophic graphs are studied. Some clarifications for each result and each definition are provided. Using basic set to extend this set to set of all vertices has key role to have these notions in the form of neutrosophic failed-independent number and failed independent neutrosophic-number arising neighborhood of vertices. The cardinality of a set has eligibility to neutrosophic failed-independent number but the neutrosophic cardinality of a set has eligibility to call failed independent neutrosophic-number. Some results get more frameworks and perspective about these definitions. The way in that, two vertices have connections amid each other, opens the way to do some approaches. A vertex could affect on other vertex but there’s no usage of edges. These notions are applied into neutrosophic graphs as individuals but not family of them as drawbacks for these notions. Finding special neutrosophic graphs which are well-known, is an open way to pursue this study. Some problems are proposed to pursue this study. Basic familiarities with graph theory and neutrosophic graph theory are proposed for this article.

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.

New setting is introduced to study quasi-degree and quasi-co-degree arising from co-neighborhood. quasi-degree and quasi-co-degree is about a vertex which are applied into the setting of neutrosophic graphs. . The structure of set is studied and general results are obtained. Also, some classes of neutrosophic graphs namely path-neutrosophic graphs, cycle-neutrosophic graphs, complete-neutrosophic graphs and star-neutrosophic graphs, complete-bipartite-neutrosophic graphs and complete-multipartite-neutrosophic graphs are investigated in the terms of a vertex which is called either quasi-degree or quasi-co-degree. Neutrosophic number is reused in this way. It’s applied to use the type of neutrosophic number in the way that, three values of a vertex are used and they’ve same share to construct this number to compare with other vertices. Summation of three values of vertex makes one number and applying it to a comparison. This approach facilitates identifying vertices which form quasi-degree and quasi-co-degree. Quasi-degree is a value of a vertex which is maximum amid all values of vertices which are neighbors to a fixed vertex. Quasi-co-degree is a value of an edge which is maximum amid all values of edges which are neighbors to a fixed vertex but corresponded vertex is representative for this notion. Using different values which are related to a vertex inspire us to focus on edge and vertices which are corresponded to a fixed vertex. The notion of neighborhood is used to collect either vertices are titled neighbors or edges are incident to fixed vertex. In both settings, some classes of well-known neutrosophic graphs are studied. Some clarifications for each result and each definitions are provided. Using fixed vertex has key role to have these notions in the form of vertex or edge. The value of an edge has eligibility to call quasi-co-degree but the value of a vertex has eligibility to call quasi-degree. Some results get more frameworks and perspective about these definitions. The way in that, two vertices have connection together, open the way to define neighborhood and co-neighborhood. The maximum values in neighborhood and co-neighborhood introduces quasi-degree and quasi-co-degree, respectively. New name is chosen from degree. Since amid all vertices with different degrees, one vertex is chosen. In other words, one vertex is fixed and its degree turns out quasi-degree where two degrees could be assigned to a vertex. Degree of edges and degree of vertices. The number of edges which are incident to the vertex and the number of vertices which are neighbors to the vertex. Degree and co-degree are the notions which are transformed to use in quasi-style. Two neutrosophic values introduce two neutrosophic vertices separately in each settings. These notions are applied into neutrosophic graphs as individuals but not family of them as drawbacks for these notions. Finding special neutrosophic graphs which are well-known, is an open way to purse this study. Some problems are proposed to pursue this study. Basic familiarities with graph theory and neutrosophic graph theory are proposed for this article.

3D Zernike moments based on 3D Zernike polynomials have been successfully applied to the field of voxelized 3D shape retrieval and have attracted more attention in biomedical image processing. As the order of 3D Zernike moments increases, both computational efficiency and numerical accuracy decrease. Due to this phenomenon, a more efficient and stable method for computing high-order 3D Zernike moments was proposed in this study. The proposed recursive formula for computing 3D Zernike radial polynomials combines the recursive calculation of spherical harmonics to develop a voxel-based algorithm for the calculation of 3D Zernike moments. The algorithm was applied to the 3D shape Michelangelo's David with a size of 150×150×150 voxels. As compared to the method without additional acceleration, the proposed method uses a group action of order sixteen orthogonal group and saving unnecessary iterations, the factor of speed-up is 56.783±3.999 when the order of Zernike moments is between 10 and 450. The proposed method also obtained an accurate reconstructed shape with the error rate (normalized mean square error) of 0.00 (4.17×10^-3) when the reconstruction was computed for all moments up to order 450.

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.

In this study, we mainly focus on the behavior of an electron gas produced by argon plasma that is bounded by a moving rigid flat plate. The effects of an external magnetic field on the electrons collected with each other, with positive ions, and with neutral atoms in the plasma fluid are investigated. We used the BGK model of the Boltzmann kinetic equation, which is fundamental in the study of the gas dynamics from the continuum flow to the free-molecular regimes with Maxwellian velocity distribution functions of the various species. An exact analytical solution of the model equations for the unsteady flow is given using the moment and the traveling wave methods. The behavior of the mean velocity of electrons is illustrated, which compatible with the variation of the shear stress, viscosity coefficient, and the initial and boundary conditions. Besides, the thermodynamic prediction is investigated by applying the principles of irreversible thermodynamics and Gibbs formula. It is found that the thermodynamic variables verify the Boltzmann’s H–theorem, the principle of Le Chatelier, and the Thermodynamics Second Law for the non-equilibrium processes of the system. Qualitative agreements with previous related papers are introducing. 3-Dimensional graphics for the calculating variables are offering, and their behavior is deeply discussed.

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 present paper is concerned with the study of the quality of life index. Such an index has become an important index for measuring the well-being of individuals. However, the quality of life index is always a subjective, intangible, and often hard to quantify with precision due to the lack of quantitative models dealing with. The main goal of the present paper is thus to propose a mathematical, quantitative model for the measurement of a quality of life index. The main novelty is firstly the construction of a wavelet dynamic multiscale model to quantify and investigate the effect of time scale on the quality of life index measuring. The proposed procedure is acted empirically on a sample corresponding to Saudi Arabia as a case of study during the period from 2003 to 2020 as part of the 2030-vision plan. Saudi Arabia has implemented the so-called 2030-vision plan where the quality of life improvement is one of the main goals to be attempted. The findings show that wavelets are capable to localize the time-wise behavior of the index contrarily to classical studies which estimate a global view of the index. Moreover, the study shows the link between the quality of life behavior and many other indices.

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.

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 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.

We investigate, whereas the SU(4)/SU(2) model, is viable for describing strong interactions. The existing problems of confinement and asymptotic freedom, two phenomena that can not be described by the SU(3) model, might indicate that we need something "larger" than the SU(3) model. The considered coset SU(4)/SU(2) or the Stiefel manifold V^2(C_4) contains an su(3) algebra, plus additional degrees of freedom that resembes the Feddev-Poppov concept. The richer structure of this coset, give us enough room, to seek for new phenomena, as its dimensionality is 12. The consideration of the SU(4)/SU(2) model flavors firstly the unification of nuclear fields (strong and weak) in an SU(4) model.

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.

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.

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.

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.

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.

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.

The main idea of this paper is to investigate the exact solutions and dynamic properties in optical nanofibers, which is modeled by space-time fractional perturbed nonlinear schr\"odinger equation involving Kerr law nonlinearity with conformal fractional derivative. Firstly, by the complex fractional traveling wave transformation, the traveling wave system of the original equation is obtained, then a conserved quantity, namely the Hamiltonian is constructed, and the qualitative analysis of this system is conducted via this quantity by classifying the equilibrium points. Moreover, the prior estimate of the existence of the soliton and periodic solution is established via the bifurcation method. Furthermore, all exact traveling wave solutions are constructed to illustrate our results explicitly by the complete discrimination system for polynomial method.

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 the present paper, the effect of the non-linear thermal radiation on the neutral gas mixture in the unsteady state is investigated for the first time. The unsteady BGK technique of the Boltzmann kinetic equations for a neutral non-homogenous gas is solved. The solution of the unsteady case makes the problem more general significance than the stationary one. For this purpose, the moments' method, together with the traveling wave method, is applied. The temperature and concentration are calculated for each gas component and mixture for the first time.Furthermore, the study is held for aboard range of temperatures ratio parameter and a wide range of the molar fraction. The distribution functions are calculated for each gas component and the gas mixture. The significant non-equilibrium irreversible thermodynamic characteristics the entire system is acquired analytically. That technic allows us to investigate the consistency of Boltzmann's H-theorem, Le Chatelier principle, and thermodynamics laws. Moreover, the ratios among the different participation of the internal energy alteration are evaluated via the Gibbs formula of total energy. The final results are utilized to the argon-helium non-homogenous gas at different magnitudes of radiation force strength and molar fraction parameters. 3D-graphics are presented to predict the behavior of the calculated variables, and the obtained results are theoretically discussed.

In this paper it is shown that mathematical description of strain, constitutive law and dynamics obtained by direct replacement of integer order derivatives with Caputo or Riemann-Liouville fractional order partial derivatives, having integral representation on finite interval, in case of a guitar string, is nonobjective. The basic idea is that different observers, using this type of descriptions, obtain different results which cannot be reconciled, i.e. transformed into each other using only formulas that link the coordinates of the same point in two fixed orthogonal reference frames and formulas that link the numbers representing the same moment of time in two different choices of the origin of time measuring. This is not an academic curiosity! It is rather a problem: which one of the obtained results is correct?

Energy efficiency of Data Center (DC) operations heavily relies on IT and cooling systems performance. A reliable and efficient cooling system is necessary to produce a persistent flow of cold air to cool servers that are subjected to constantly increasing computational load due to the advent of IoT- enabled smart systems. Consequently, increased demand for computing power will bring about increased waste heat dissipation in data centers. In order to bring about a DC energy efficiency, it is imperative to explore the thermal characteristics analysis of an IT room (due to waste heat). This work encompasses the employment of an unsupervised machine learning modelling technique for uncovering weaknesses of the DC cooling system based on real DC monitoring thermal data. The findings of the analysis result in the identification of areas for energy efficiency improvement that will feed into DC recommendations. The methodology employed for this research includes statistical analysis of IT room thermal characteristics, and the identification of individual servers that frequently occur in the hotspot zones. A critical analysis has been conducted on available big dataset of ambient air temperature in the hot aisle of ENEA Portici CRESCO6 computing cluster. Clustering techniques have been used for hotspots localization as well as categorization of nodes based on surrounding air temperature ranges. The principles and approaches covered in this work are replicable for energy efficiency evaluation of any DC and thus, foster transferability. This work showcases applicability of best practices and guidelines in the context of a real commercial DC that transcends the set of existing metrics for DC energy efficiency assessment.

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.

New setting is introduced to study co-neighborhood, neutrosophic t-neighborhood, neutrosophic quasi-vertex set, neutrosophic quasi-order, neutrosophic neighborhood, neutrosophic co-t-neighborhood, neutrosophic quasi-edge set, neutrosophic quasi-size, Neutrosophic number, neutrosophic co-neighborhood, co-neutrosophic number, quasi-number and quasi-co-number. Some classes of neutrosophic graphs are investigated.

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.

A data-driven analysis method known as dynamic mode decomposition (DMD) approximates the linear Koopman operator on projected space. In the spirit of Johnson-Lindenstrauss Lemma, we will use random projection to estimate the DMD modes in reduced dimensional space. In practical applications, snapshots are in high dimensional observable space and the DMD operator matrix is massive. Hence, computing DMD with the full spectrum is infeasible, so our main computational goal is estimating the eigenvalue and eigenvectors of the DMD operator in a projected domain. We will generalize the current algorithm to estimate a projected DMD operator. We focus on a powerful and simple random projection algorithm that will reduce the computational and storage cost. While clearly, a random projection simplifies the algorithmic complexity of a detailed optimal projection, as we will show, generally the results can be excellent nonetheless, and quality understood through a well-developed theory of random projections. We will demonstrate that modes can be calculated for a low cost by the projected data with sufficient dimension.

Numerical methods, including machine learning methods, are now actively used in the applications related to guided wave propagation. The method proposed in this study for material properties characterization is based on the algorithm of the clustering of multivariate data series obtained as a result of the application of the matrix pencil method to the experimental data. In the proposed technique, multi-objective optimization is employed to improve the accuracy of particular parameter identification. At the first stage, the computationally efficient method based on the calculation of the Fourier transform of Green's matrix is employed iteratively and the obtained solution is used for the filter construction with decreasing bandwidth, which allows us to obtain nearly noise-free classified data (with mode separation). The filter provides data separation between all guided waves in a natural way, which is needed at the second stage, where the slower method based on the minimization of the slowness residuals is applied to the data. The method might be applied for material properties identification in plates with thin coatings/interlayers, multi-layered anisotropic laminates etc.

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.

This work probes the combined effects of magnetic field and viscous dissipation on heat field and examine the second law analysis (entropy generation) in an electrically conducting fluid under the effect of wall mass transfer over continuous stretched non-isothermal surface with variable viscosity. The viscosity of the fluid is assumed to be an inverse linear function of temperature. The governing equation for the problem are changed to dimensionless ordinary differential equations by using similarity transformation and solved numerically by using Rung Kutta and Shooting technique. Velocity, concentration and temperature distribution are obtained and used to compute the entropy generation and the Bejan number in the flow field. The effect of variable viscosity, Schmidt number, Hartman and Reynolds number on the velocity, concentration, temperature, entropy generation and Bejan number are studied and discussed.

We compute explicitly trace of one-dimensional diffusion process which can be regarded as a Dirichlet form on graphs and we study its conservativeness property. Finally, we give an example of a discretization of one dimensional diffusion of Bessel's process with order $\nu$.

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.

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.

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.

In this paper, we investigate the dynamics of following higher order difference equation x_{n+1}=A+B((x_{n})/(x_{n-m}²)) with A,B and initial conditions are positive numbers, and m∈{2,3,⋯}. Especially we study the boundedness, periodicity, semi-cycles, global asymptotically stability and rate of convergence of solutions of related higher order difference equations.

By partly generalizing the Lipschitz condition of existing results to the generalized Lipschitz one, the author utilizes a fixed point theorem, variational method and Lyapunov function method to derive the unique existence of globally asymptotical input-to-state stability of positive stationary solution for Gilpin-Ayala competition model with diffusion and delayed feedback under Dirichlet zero boundary value. Remarkably, it is the first paper to derive the unique existence of the stationary solution of reaction-diffusion Gilpin-Ayala competition model, which is globally asymptotical input-to-state stability. And numerical examples illuminate the effectiveness and feasibility of the proposed methods.

In this paper we provide an ``expository overview" of classic epidemic projection models. Starting with the simple case of an epidemic that grows exponentially we then investigate ``compartmental" models. These assume that the growth of an infected population is limited endogenously by the size of the underlying pool of susceptibles. We then describe a new family of so-called "Exo-r" statistical models, which hinge on an exogenously driven growth rate of the infected population. This family, which can be used to model both infections and deaths, captures parsimoniously both the depletion of susceptibles and the effect of interventions such as lockdowns and ``social distancing". The model is used to fit numbers of Covid-19 infections in China. It is also used to model and project deaths in the United States. Results are used to inform a discussion on i) the challenges at hand and ii) the extent to which epidemic projection models may be useful despite being wrong.

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.

In a one price economy, the Fundamental Theorem of Asset Pricing (FTAP) establishes that no-arbitrage is equivalent to the existence of an equivalent martingale measure. Such an equivalent measure can be derived as the normal unit vector of the hyperplane that separates the attainable gain subspace and the convex cone representing arbitrage opportunities. However, in two-price financial models (where there is a bid-ask price spread), the set of attainable gains is not a subspace anymore. We use convex optimization, and the conic property of this region to characterize the “No-Arbitrage” principle in financial models with the bid-ask price spread present. This characterization will lead us to the generation of a set of price factor random variables. Under such a set, we can find the lower and upper bounds (supper-hedging and sub-hedging bounds) for the price of any future cash flow. We will show that for any given cash flow, for which the price is outside the above range, we can build a trading strategy that provides one with an arbitrage opportunity. We will generalize this structure to any two-price finite-period financial model.

We review a modern differential geometric description of the fluid isotropic motion and featuring it the diffeomorphism group structure, modelling the related dynamics, as well as its compatibility with the quasi-stationary thermodynamical constraints. There is analyzed the adiabatic liquid dynamics, within which, following the general approach, there is explained in detail, the nature of the related Poissonian structure on the fluid motion phase space, as a semidirect Banach groups product, and a natural reduction of the canonical symplectic structure on its cotangent space to the classical Lie-Poisson bracket on the adjoint space to the corresponding semidirect Lie algebras product. We also present a modification of the Hamiltonian analysis in case of the isotermal liquid dynamics. We study the differential-geometric structure of the adiabatic magneto-hydrodynamic superfluid phase space and its related motion within the Hamiltonian analysis and invariant theory. In particular, we construct an infinite hierarchies of different kinds of integral magneto-hydrodynamic invariants, generalizing those before constructed in the literature, and analyze their differential-geometric origins. A charged liquid dynamics on the phase space invariant with respect to an abelian gauge group transformation is also investigated, some generalization of the canonical Lie-Poisson type bracket is obtained.

In this manuscript a new approach into analyzing the local stability of equilibrium points of non-linear Caputo fractional planar systems is shown. It is shown that the equilibrium points of such systems can be a stable focus or unstable focus. In addition, it is proposed that previous results regarding the stability of equilibrium points have been incorrect, the results here attempt to correct such results. Lastly, it is proposed that a Caputo fractional planar system cannot undergo a Hopf bifurcation, contrary to previous results prior. Though, it is shown that such systems can undergo a Hopf bifurcation (topologically).

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.

Financial system is essentially chaotic and unstable if there is not any external inputs. By means of Lyapunov function method, design of switching law, novel fuzzy assumption, $L^p$ estimation technique and Laplace semigroup theory, the author presents the boundedness and LMI-based (globally) asymptotical input-to-state stability criteria of financial systems. Particularly, the globally asymptotical stability in the meaning of switching implies that when the time $t$ is big enough, the dynamic of any subsystem must approach its unique equilibrium point. Besides, the global financial crisis often erupts periodically, which illuminates that the global stability in the classical sense is actually meaningless. So the stability in the meaning of switching proposed in this paper is suitable and appropriate. Numerical examples illuminate the effectiveness of the obtained results.

This paper aims to investigate the global asymptotic stability of following system of higher order difference equations with quadratic terms: x_{n+1}=A+B((y_{n})/(y_{n-m}²)),y_{n+1}=A+B((x_{n})/(x_{n-m}²)) where A and B are positive numbers and the initial values are positive numbers. We also study the rate of convergence and oscillation behaviour of the solutions of related system.

In the upper charged layers of the atmosphere, the plasma is very rarefied. The collisions between its molecules are almost non-existent, and the driving forces behind them are the Lorentz forces resulting from the electric and magnetic fields. For this reason, we are interested in studying the behavior of non-collision plasmas because of its essential applications, such as the movement of satellites in the charged atmosphere. In this paper, the flow problem of collisionless gaseous plasma is examined. For that propose, we solve the unsteady Vlasov-Maxwell system of non-linear partial differential equations analytically. Methods of moments and traveling wave parameters are used to acquire an exact solution. Specific macroscopic properties of collisionless gaseous plasma are calculated along with electrical and magnetic fields. Further, thermodynamic estimation, such as entropy and entropy production, is presented. Those calculations allow us to measure the consistency with the laws of non-equilibrium thermodynamics. Relations between internal energy modification participations are predicted using Gibbs' equation for collisionless plasma. The modification effect of internal energies due to electro-magnetic fields is found to be small compared with the internal energy change due to the effect of entropy. That is because these fields are self-induced by plasma particles due to the sudden movement of the rigid plane plate. The results are accomplished according to the typical argon gaseous plasma model. Three-dimensional diagrams showing the measured variables are drawn to investigate and discuss their behavior. The problem has many commercial applications for the movement of objects in the charged atmosphere.

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.

The purpose of this paper is threefold: (i) to introduce and study the Euler-Lagrange prolongations of flatness PDEs solutions (best approximation of flatness) via associated least squares Lagrangian densities and integral functionals on Riemannian manifolds; (ii) to analyze some decomposable multivariate dynamics represented by Euler-Lagrange PDEs of least squares Lagrangians generated by flatness PDEs and Riemannian metrics; (iii) to give examples of explicit flat extremals and non-flat approximations.

Applying a procedure similar to that of E.S. Bring, by using a 4^th degree Tschirnhaus transformation, it was possible to transform the Bring-Jerrard normal quintic (BJQ) equation into a De Moivre form (DMQ), so that it could be solved by radicals. The general solution by radicals of the De Moivre equations of any degree is presented. By the same procedure the BJSx (normal sextic) equation was taken to another one without the 2^nd, 4^th and 6^th terms which was transformed into a cubic (solvable) equation. By applying a 6^th degree Tschirnhaus transformation to the BJSp (normal septic) equation its Bring-Jerrard binormal (without the 2^nd, 3^rd, 4^th and 5^th terms) form was obtained.

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.

Autonomous differential equations of fractional order and non-singular kernel are solved. While solutions can be obtained through numerical, graphical, or analytical solutions, we seek an implicit analytical solution.

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.

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.

Radiation Fireballs are singular phenomena which involve severe thermal radiation and consequently, they need to be duly assessed and prevented. Although the radiative heat transfer produced by a sphere is relatively well known, the shadowing measures implemented to control the fireball’s devastating effects have frequently posed a difficult analytical instance, mainly due to its specific configuration. In this article, since the usual solving equations for the said cases are impractical, the authors propose a novel graphic-algorithm method that sorts the problem efficiently for different kinds of obstructions and relative positions of the fireball and the defenses. Adequate application of this method may improve the safety of a significant number of facilities exposed to such risks.

In the current era of ‘big data’, scientists are able to quickly amass enormous amount of data in a limited number of experiments. The investigators then try to hypothesize about the root cause based on the observed trends for the predictors and the response variable. This involves identifying the discriminatory predictors that are most responsible for explaining variation in the response variable. In the current work, we investigated three related multivariate techniques: Principal Component Regression (PCR), Partial Least Squares or Projections to Latent Structures (PLS), and Orthogonal Partial Least Squares (OPLS). To perform a comparative analysis, we used a publicly available dataset for Parkinson’ disease patien ts. We first performed the analysis using a cross-validated number of principal components for the aforementioned techniques. Our results demonstrated that PLS and OPLS were better suited than PCR for identifying the discriminatory predictors. Since the X data did not exhibit a strong correlation, we also performed Multiple Linear Regression (MLR) on the dataset. A comparison of the top five discriminatory predictors identified by the four techniques showed a substantial overlap between the results obtained by PLS, OPLS, and MLR, and the three techniques exhibited a significant divergence from the variables identified by PCR. A further investigation of the data revealed that PCR could be used to identify the discriminatory variables successfully if the number of principal components in the regression model were increased. In summary, we recommend using PLS or OPLS for hypothesis generation and systemizing the selection process for principal components when using PCR.rewordexplain later why MLR can be used on a dataset with no correlation

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

So called Higuchi’s method of fractal dimension estimation is widely used and the term Higuchi’s fractal dimension even occurs in many publications. This paper deals with this method from mathematical point of view. Terms distance and dimension and its basic properties are explained and Higuchi’s dimension according the original source is defined. Definition of Higuchi’s dimension was comparated with mathematical definition of the distance and dimension. It is showed, that the definition of the Higuchi’s dimension does not satisfy axioms of distance and dimension. So called Higuchi’s method and Higuchi’s dimension are mathematically incorrect. Therefore, all results achieved by this method are scientifically unreliable.

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.

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 real financial market, the delayed market feedback and the delayed effect of government macro-control are inevitable. And both the delay of market feedback and the delay of macro-control effect bring about the mathematical difficulties in studying stabilization and synchronization of the hyper-chaotic financial system. However, employing Lyapunov function method, differential mean value theorem, suitable bounded hypotheses and pulse control technology results in the globally asymptotical stabilization and synchronization criteria. It is the first paper to drive the stabilization and synchronization criteria under the assumptions of the double delays. Finally, numerical examples illuminate the effectiveness of the proposed methods.

Instead of the conventional 0 and 1 values, bipolar reasoning uses -1, 0, +1 to describe double-sided judgements in which neutral elements are halfway between positive and negative evaluations (e.g., “uncertain” lies between “impossible” and “totally sure”). We discuss the state-of-the-art in bipolar logics and recall two medieval forerunners, i.e., William of Ockham and Nicholas of Autrecourt, who embodied a bipolar mode of thought that is eminently modern. Starting from the trivial observation that “once a wheat sheaf is sealed and tied up, the packed down straws display the same orientation”, we work up a new theory of the bipolar nature of networks, suggesting that orthodromic (i.e., feedforward, bottom-up) projections might be functionally coupled with antidromic (i.e., feedback, top-down) projections via the mathematical apparatus of presheaves/globular sets. When an entrained oscillation such as a neuronal spike propagates from A to B, changes in B might lead to changes in A, providing unexpected antidromic effects. Our account points towards the methodological feasibility of novel neural networks in which message feedback is guaranteed by backpropagation mechanisms endowed in the same feedforward circuits. Bottom-up/top-down transmission at various coarse-grained network levels provides fresh insights in far-flung scientific fields such as object persistence, memory reinforcement, visual recognition, Bayesian inferential circuits and multidimensional activity of the brain. Implying that axonal stimulation by external sources might backpropagate and modify neuronal electric oscillations, our theory also suggests testable previsions concerning the optimal location of transcranial magnetic stimulation’s coils in patients affected by drug-resistant epilepsy.

An actual problem of inertial navigation related to the creation of a multifunctional device that performs the functions of a meter for both angular velocities and linear accelerations is considered. The possibility of constructing such a device based on a hybrid-type modulation micromechanical gyroscope has shown. The ability to measure linear accelerations has provided by the presence of a small symmetrical distance between the axes of the elastic suspension relative to the center of mass of the sensor element. A mathematical model of the device ("heavy" gyroscope) as a high-quality three-dimensional oscillatory system is constructed. It has shown numerically that the reaction of the system to the movement of an object has, along with precession, the observed nutation. This allows you to develop an algorithm for separating information about the angular and linear motion of the base and implement a two-component angular velocity meter and a two-component linear acceleration meter in one device.

Four discrete models using the exact spectral derivative discretization finite difference (ESDDFD) method are proposed for a chaotic five-dimensional, conformable fractional derivative financial system incorporating ethics and market confidence. Since the system considered was recently studied using the conformable Euler finite difference (CEFD) method and found to be hyperchaotic, and the CEFD method was recently shown to be valid only at fractional index , the source of the hyperchaos is in question. Through numerical experiments, illustration is presented that the hyperchaos previously detected is in part an artifact of the CEFD method as it is absent from the ESDDFD models.

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.

This article presents a method of segmenting images with gray levels that uses Onicescu's information energy calculated in the context of the neutrosophic theory. Starting from the information energy calculation for complete neutrosophic information, it is shown how to extend its calculation for incomplete and inconsistent neutrosophic information. The segmentation method is based on calculation of thresholds for separating the gray levels using the local maximum points of the Onicescu information energy.

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.

New ideas on the framework of Neutrosophic SuperHyperGraph for different styles of Neutrosophic SuperHyper-Bipartite and Neutrosophic SuperHyper-Path are introduced. More instances and more clarifications alongside sufficient references.

In today's world it is very much important to maintain the security of information and its risks. The biometric-based techniques are very much useful in these problems. Among the several kinds of biometric-based technique, face detection is much complex and much more important. Due to the age and several other problems, a human face structure changes over time, again a human has lots of expressions. Sometimes due to the lighting condition or the variation of the angle of an input device, the pattern of a human face structure also changed. As a result, the face cannot be detected properly. In this paper, a method is proposed that can detect the human faces both automatically and manually very efficiently. In manual mode, a user can select the input faces referred by the system according to their choice. In automated mode, the system detected all possible face areas using the Kohonen Self-Organizing Feature Map technique. This method reduced the complex color image into a vector quantized image with desired colors. Then a color segmentation technique is used to detect the possible face skin areas from the vector quantized image. Then the Histogram Oriented Gradient technique used to detect the feature from the faces and K-Nearest Neighbour Classifier is used to compare both face images detected by the two modes. The automated method prosed better accuracy than the manual method.

When presenting special relativity, it is customary to single out the so-called paradoxes. One of these paradoxes is the formal occurrence of speeds exceeding the speed of light. An essential part of such paradoxes arises from the incompleteness of the relativistic calculus of velocities. In special relativity, the additive group is used for velocities. However, the use of only group operations imposes artificial restrictions on possible computations. Naive expansion to vector space is usually done by using non-relativistic operations. We propose to consider arithmetic operations in the special theory of relativity in the framework of the Cayley–Klein model for projective spaces. We show that such paradoxes do not arise in the framework of the proposed relativistic extension of algebraic operations.

We use a system biology approach to translate the interaction of Bacillus Calmette-Gurin (BCG) + interleukin 2 (IL-2) for the treatment of bladder cancer into a mathematical model. The model is presented as a system of differential equations with the following variables: number of tumor cells, bacterial cells, immune cells, and cytokines involved in the tumor-immune response. This work investigates the delay effect induced by the proliferation of tumor antigen-specific effector cells after the immune system destroys BCG-infected urothelium cells following BCG and IL-2 immunotherapy in the treatment of bladder cancer. For the proposed model, three equilibrium states are found analytically. The stability of all equilibria is analyzed using the method of Lyapunov functionals construction and the method of linear matrix inequalities (LMIs).

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.

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.

In this paper, the existence of two nontrivial stationary solutions for the nonlinear Gilpin Ayala two species competition model is given by using the mountain pass lemma, and the local stability criterion of the trivial solution is given by using Lyapunov function method. Based on the local stability criterion, we give some suggestions on how to avoid the population extinction. This is, when the population is on the verge of extinction, we should try our best to avoid the diffusion behavior and reduce the diffusion coefficient, otherwise the species are easy to go extinct. Numerical example shows the effectiveness of the proposed method.

In this paper, we consider a make-to-order supply chain which satisfies demand that is dependent on both price and quoted lead -time. The manufacturer chooses the lead -time and the order quantity, and the retailer sets the revenue shares. The interactions between the manufacturer and the retailer are modelled as a Nash Game, and the existence and uniqueness of pure strategy equilibrium are demonstrated. A mechanism that enables the supply chain to coordinate the decisions of the members is developed. Lastly, we also analyze how the supply chain system parameters impact the optimal supply chain decisions and the supply chain performance.

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.

The aim of this work is to introduce a mathematical model representing the evolution of the temperature in a vegetation cover and the ground underneath it. Vegetation, and its interaction with soil, plays a very important role in the protection of soil surface from the action of sun and precipitations. A reduction in the vegetated mass increase the risk of desertification, soil erosion or surface runoffs which which can give rise to soil loss and sediment retention. These processes can favour climate change and global warming, which are major concerns nowadays. The mathematical model presented takes into account the main processes involved in vegetation cover and the interaction with the soil, among which, we can mention the Leaf Area Index, which is a dimensionless quantity defined as the one-sided green leaf area per unit ground surface area, or albedo and co-albedo which are clearly influenced by the vegetation. It is also considered a nonlinear heat capacity in the soil which incorporates the latent heat of fusion, when the phase change takes place. The numerical technique used to solve the mathematical model is based on a finite volume scheme with Weighted Essentially Non Oscillatory technique for spatial reconstruction and the third order Runge-Kutta Total Variation Diminishing numerical scheme is used for time integration. Some numerical examples are solved to obtain the distribution of temperature both in the vegetation cover and the soil.

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.

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.

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.

Increased data acquisition by uncalibrated, heterogeneous digital sensor systems such as smartphones present new challenges. Binary metrics are proposed for the quantification of cyber-physical signal characteristics and features, and a standardized constant-Q variation of the Gabor atom is developed for use with wavelet transforms. Two different continuous wavelet transform (CWT) reconstruction formulas are presented and tested under different signal to noise ratio (SNR) conditions. A sparse superposition of Nth order Gabor atoms worked well against a synthetic blast transient using the wavelet entropy and an entropy-like parametrization of the SNR as the CWT coefficient-weighting functions. The proposed methods should be well suited for sparse feature extraction and dictionary-based machine learning across multiple sensor modalities.

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.

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 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.

In this paper, the authors employ Mountain Pass Lemma, the method of weak solution regularization and Lyapunov function method to derive the unique existence of globally exponential stable positive stationary solution of a single-species model with diffusion and delayed feedback. The obtained stability criterion illuminates that under some suitable conditions, a certain internal competition is conducive to the overall stability of the population, and a certain amount of family planning is conducive to the overall stability of the population. A numerical example and three tables show the effectiveness of the proposed methods.

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).

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.

With the revolution of computing and biology technology, data sets containing information could be huge and complex that sometimes are difficult to handle. Dynamic computing is an efficient approach to solve some of the problems. Since neighborhood multigranulation rough sets(NMGRS) were proposed, few papers focused on how to calculate approximations in NMGRS and how to update them dynamically. Here we propose approaches for computing approximations in NMGRS and updating them dynamically. First, static approaches for computing approximations in NMGRS are proposed. Second, search region in data set for updating approximations in NMGRS is shrunk. Third, matrix-based approaches for updating approximations in NMGRS while decreasing or increasing neighborhood classes are proposed. Fourth, incremental algorithms for updating approximations in NMGRS while decreasing or increasing neighborhood classes are designed. Finally, the efficiency and validity of the designed algorithms are verified by experiments.

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.

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&lt;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.