Search

In order to solve the Sylvester equations more efficiently, a new four parameters positive and skew-Hermitian splitting (FPPSS) iterative method is proposed in this paper based on the previous research of the positive and skew-Hermitian splitting (PSS) iterative method. We prove that when coefficient matrix A and B satisfy certain conditions, the FPPSS iterative method is convergent in the parameter&rsquo, s value region. The numerical experiment results show that compared with previous iterative method, the FPPSS iterative method is more effective in terms of iteration number IT and runtime.

In this paper, the existence of the solution and its stability to the fractional boundary value problem (FBVP) were investigated for an implicit nonlinear fractional differential equation (VOFDE) of variable order. All existence criteria of the solutions in our establishments were derived via Krasnoselskii’s fixed point theorem and in the sequel, and its Ulam–Hyers–Rassias (U-H-R) stability is checked. An illustrative example is presented at the end of this paper to validate our findings.

In this paper, we construct and study a new family of multi-point Ehrlich-type iterative methods for approximating all the zeros of a uni-variate polynomial simultaneously. The first member of this family is the two-point Ehrlich-type iterative method introduced and studied by Trićković and Petković in 1999. The main purpose of the paper is to provide local and semilocal convergence analysis of the multi-point Ehrlich-type methods. Our local convergence theorem is obtained by an approach that was introduced by the authors in 2020. Two numerical examples are presented to show the applicability of our semilocal convergence theorem.

In this paper, we consider fractional ordinary differential equations and the fractional Euler–Poisson–Darboux equation with fractional derivatives in the form of a power of the Bessel differential operator. Using the technique of the Meijer integral transform and its modification, fundamental solutions to these equations are derived in terms of the Fox–Wright function, the Fox H-function, and their particular cases. We also provide some explicit formulas for the solutions to the corresponding initial-value problems in terms of the generalized convolutions introduced in this paper.

The paper presents k-version of the finite element method for boundary value problems (BVPs) and initial value problems (IVPs) in which global differentiability of approximations is always the result of the union of local approximations. The higher order global differentiability approximations (HGDA/DG) are always p-version hierarchical that permit use of any desired p-level without effecting global differentiability. HGDA/DG are true Ci, Cij, Cijk, hence the dofs at the nonhierarchical nodes of the elements are transformable between natural and physical coordinate spaces using calculus. This is not the case with tensor product higher order continuity elements discussed in this paper, thus confirming that the tensor product approximations are not true Ci, Cijk, Cijk approximations. It is shown that isogeometric analysis for a domain with more than one patch can only yield solutions of class C0. This method has no concept of finite elements and local approximations, just patches. It is shown that compariso of this method with k-version of the finite element method is meaningless. Model problem studies in R2 establish accuracy and superior convergence characteristics of true Cijp-version hierarchical local approximations presented in this paper over tensor product approximations. Convergence characteristics of p-convergence, k-convergence and pk-convergence are illustrated for self adjoint, non-self adjoint and non-linear differential operators in BVPs. It is demonstrated that h, p and k are three independent parameters in all finite element computations. Tensor product local approximations and other published works on k-version and their limitations are discussed in the paper and are compared with present work.

This paper provides several error estimations for total variation (TV) type regularization, which arises in a series of areas, for instance, signal and imaging processing, machine learning, etc. In this paper, some basic properties of the minimizer for the TV regularization problem such as stability, consistency and convergence rate are fully investigated. Both a priori and a posteriori rules are considered in this paper. Furthermore, an improved convergence rate is given based on the sparsity assumption. The problem under the condition of non-sparsity, which is common in practice, is also discussed, the results of the corresponding convergence rate are also presented under certain mild conditions.

Power electronic converters are mathematically represented by a system of ordinary differential equations discontinuous right-hand side that does not verify the conditions of the Cauchy-Lipschitz Theorem. More generally, for the properties that characterize their discontinuous behavior, they represent a particular class of systems on which little has been investigated over the years. The purpose of the paper is to prove the existence of at least one global solution in Filippov’s sense to the Cauchy problem related to the mathematical model of a power converter and also to calculate the error in norm between this solution and the integral of its averaged approximation. The main results are the proof of this theorem and the analytical formulation that provides to calculate the cited error. The demonstration starts by a proof of local existence provided by Filippov himself and already present in the literature for a particular class of systems and this demonstration is generalized to the class of electronic power converters, exploiting the non-chattering property of this class of systems. The obtained results are extremely useful for estimating the accuracy of the averaged model used for analysis or control of the effective system. In the paper, the goodness of the analytical proof is supported by experimental tests carried out on a converter prototype representing the class of power electronics converter.

Previous works were dedicated to the functional k-Nearest Neighbors (kNN) and the local linearity method estimations of a regression operator. In this paper, a sequence pair of (Xi,Yi)i=1,…,n of functional mixing observations are considered. We treat the local linear estimation of the cumulative function of Yi given functional input variable Xi. Precisely, we combine the kNN method with the local linear algorithm to construct a new and fast efficiency estimator of the conditional distribution function. The main purpose of this paper is to prove the strong convergence of the constructed estimator under mixing conditions. An application to the functional times series prediction is used to compare our proposed estimator with the existing competitive estimators, and show its efficiency and superiority.

In previous papers, a technique has been suggested to avoid order reduction when integrating initial boundary value problems with several kinds of exponential methods. The technique implies in principle to calculate additional terms at each step from those already necessary without avoiding order reduction. The aim of the present paper is to explain the surprising result that, many times, in spite of having to calculate more terms at each step, the computational cost of doing it through Krylov methods decreases instead of increases. This is very interesting since, in that way, the methods improve not only in terms of accuracy, but also in terms of computational cost.

Here, we review some results of fractional volatility models, where the volatility is driven by fractional Brownian motion (fBm). In these models, the future average volatility is not a process adapted to the underlying filtration, and fBm is not a semimartingale in general. So, we cannot use the classical Itô’s calculus to explain how the memory properties of fBm allow us to describe some empirical findings of the implied volatility surface through Hull and White type formulas. Thus, Malliavin calculus provides a natural approach to deal with the implied volatility without assuming any particular structure of the volatility. The aim of this paper is to provides the basic tools of Malliavin calculus for the study of fractional volatility models. That is, we explain how the long and short memory of fBm improves the description of the implied volatility. In particular, we consider in detail a model that combines the long and short memory properties of fBm as an example of the approach introduced in this paper. The theoretical results are tested with numerical experiments.

This paper describes a number of simple but quite effective methods for constructing exact solutions of nonlinear partial differential equations that involve a relatively small amount of intermediate calculations. The methods employ two main ideas: (i) simple exact solutions can serve to construct more complex solutions of the equations under consideration and (ii) exact solutions of some equations can serve to construct solutions of other, more complex equations. In particular, we propose a method for constructing complex solutions from simple solutions using translation and scaling. We show that in some cases, rather complex solutions can be obtained by adding one or more terms to simpler solutions. There are situations where nonlinear superposition allows us to construct a complex composite solution using similar simple solutions. We also propose a few methods for constructing complex exact solutions to linear and nonlinear PDEs by introducing complex-valued parameters into simpler solutions. The effectiveness of the methods is illustrated by a large number of specific examples (over 30 in total). These include nonlinear heat equations, reaction–diffusion equations, wave type equations, Klein–Gordon type equations, equations of motion through porous media, hydrodynamic boundary layer equations, equations of motion of a liquid film, equations of gas dynamics, Navier–Stokes equations, and some other PDEs. Apart from exact solutions to `ordinary’ partial differential equations, we also describe some exact solutions to more complex nonlinear delay PDEs. Along with the unknown function at the current time, u=u(x,t), these equations contain the same function at a past time, w=u(x,t−τ), where τ&gt, 0 is the delay time. Furthermore, we look at nonlinear partial functional-differential equations of the pantograph type, which, in addition to the unknown u=u(x,t), also contain the same functions with dilated or contracted arguments, w=u(px,qt), where p and q are scaling parameters. We propose an efficient approach to construct exact solutions to such functional-differential equations. Some new exact solutions of nonlinear pantograph-type PDEs are presented. The methods and examples in this paper are presented according to the principle “from simple to complex”.

In this paper, we present a numerical solution of a two-phase fractional Stefan problem with time derivative described in the Caputo sense. In the proposed algorithm, we use a special case of front-fixing method supplemented by the iterative procedure, which allows us to determine the position of the moving boundary. The presented method is an extension of a front-fixing method for the one-phase problem to the two-phase case. The novelty of the method is a new discretization of the partial differential equation dedicated to the second phase, which is carried out by introducing a new spatial variable immobilizing the moving boundary. Then, the partial differential equation is transformed to an equivalent integro-differential equation, which is discretized on a homogeneous mesh of nodes with a constant spatial and time step. A new convergence criterion is also proposed in the iterative algorithm determining the location of the moving boundary. The motivation for the development of the method is that the analytical solution of the considered problem is impossible to calculate in some cases, as can be seen in the figures in the paper. Moreover, the change of the boundary conditions makes obtaining a closed analytical solution very problematic. Therefore, creating new numerical methods is very valuable. In the final part, we also present some examples illustrating the comparison of the analytical solution with the results received by the proposed numerical method.

This paper is concerned with the orthogonal polynomials. Upper and lower bounds of Legendre polynomials are obtained. Furthermore, entropies associated with discrete probability distributions is a topic considered in this paper. Bounds of the entropies which improve some previously known results are obtained in terms of inequalities. In order to illustrate the results obtained in this paper and to compare them with other results from the literature some graphs are provided.

In this paper a new method called the least squares differential quadrature method (LSDQM) is introduced as a straightforward and efficient method to compute analytical approximate polynomial solutions for nonlinear partial differential equations with fractional time derivatives. LSDQM is a combination of the differential quadrature method and the least squares method and in this paper it is employed to find approximate solutions for a very general class of nonlinear partial differential equations, wherein the fractional derivatives are described in the Caputo sense. The paper contains a clear, step-by-step presentation of the method and a convergence theorem. In order to emphasize the accuracy of LSDQM we included two test problems previously solved by means of other, well-known methods, and observed that our solutions present not only a smaller error but also a much simpler expression. We also included a problem with no known exact solution and the solutions computed by LSDQM are in good agreement with previous ones.

In this paper, we study the oscillatory behavior of solutions for a type of generalized proportional fractional differential equations with forcing and damping terms. Several oscillation criteria are established for the proposed equations in terms of Riemann-Liouville and Caputo settings. The results of this paper generalize some existing theorems in the literature. Indeed, it is shown that for particular choices of parameters, the obtained conditions in this paper reduce our theorems to some known results. Numerical examples are constructed to demonstrate the effectiveness of the our main theorems. Furthermore, we present and illustrate an example which does not satisfy the assumptions of our theorem and whose solution demonstrates nonoscillatory behavior.

In this survey paper, we present both some basic properties of the four-parameters Wright function and its applications in Fractional Calculus. For applications in Fractional Calculus, the four-parameters Wright function of the second kind is especially important. In the paper, three case studies illustrating a wide spectrum of its applications are presented. The first case study deals with the scale-invariant solutions to a one-dimensional time-fractional diffusion-wave equation that can be represented in terms of the Wright function of the second kind and the four-parameters Wright function of the second kind. In the second case study, we consider a subordination formula for the solutions to a multi-dimensional space-time-fractional diffusion equation with different orders of the fractional derivatives. The kernel of the subordination integral is a special case of the four-parameters Wright function of the second kind. Finally, in the third case study, we shortly present an application of an operational calculus for a composed Erdélyi-Kober fractional operator for solving some initial-value problems for the fractional differential equations with the left- and right-hand sided Erdélyi-Kober fractional derivatives. In particular, we present an example with an explicit solution in terms of the four-parameters Wright function of the second kind.

The purpose of this paper is to introduce the split combination of variational inclusion problem which combines the concept of the modified variational inclusion problem introduced by Khuangsatung and Kangtunyakarn and the split variational inclusion problem introduced by Moudafi. Using a modified Halpern iterative method, we prove the strong convergence theorem for finding a common solution for the hierarchical fixed point problem and the split combination of variational inclusion problem. The result presented in this paper demonstrates the corresponding result for the split zero point problem and the split combination of variation inequality problem. Moreover, we discuss a numerical example for supporting our result and the numerical example shows that our result is not true if some conditions fail.

Let Z(s)=(Z1(s),…,Zp(s))t be an isotropic second-order stationary multivariate spatial process. We measure the statistical association between the p random components of Z with the correlation coefficients and measure the spatial dependence with variograms. If two of the Z components are correlated, the spatial information provided by one of them can improve the information of the other. To capture this association, both within components of Z(s) and across s, we use a cross-variogram. Only two robust cross-variogram estimators have been proposed in the literature, both by Lark, and their sample distributions were not obtained. In this paper, we propose new robust cross-variogram estimators, following the location estimation method instead of the scale estimation one considered by Lark, thus extending the results obtained by García-Pérez to the multivariate case. We also obtain accurate approximations for their sample distributions using saddlepoint techniques and assuming a multivariate-scale contaminated normal model. The question of the independence of the transformed variables to avoid the usual dependence of spatial observations is also considered in the paper, linking it with the acceptance of linear variograms and cross-variograms.

The present paper concretizes the models proposed by S. Ri and N. Secelean. S. Ri proposed the construction of the fractal interpolation function (FIF) considering finite systems consisting of Rakotch contractions, but produced no concretization of the model. N. Secelean considered countable systems of Banach contractions to produce the fractal interpolation function. Based on the abovementioned results, in this paper, we propose two different algorithms to produce the fractal interpolation functions both in the affine and non-affine cases. The theoretical context we were working in suppose a countable set of starting points and a countable system of Rakotch contractions. Due to the computational restrictions, the algorithms constructed in the applications have the weakness that they use a finite set of starting points and a finite system of Rakotch contractions. In this respect, the attractor obtained is a two-step approximation. The large number of points used in the computations and the graphical results lead us to the conclusion that the attractor obtained is a good approximation of the fractal interpolation function in both cases, affine and non-affine FIFs. In this way, we also provide a concretization of the scheme presented by C.M. Păcurar.

In this paper, we introduce a new family of efficient and optimal iterative methods for finding multiple roots of nonlinear equations with known multiplicity ( m &ge, 1 ) . We use the weight function approach involving one and two parameters to develop the new family. A comprehensive convergence analysis is studied to demonstrate the optimal eighth-order convergence of the suggested scheme. Finally, numerical and dynamical tests are presented, which validates the theoretical results formulated in this paper and illustrates that the suggested family is efficient among the domain of multiple root finding methods.

This paper is devoted to the study of the Cubic B-splines to find the numerical solution of linear and non-linear 8th order BVPs that arises in the study of astrophysics, magnetic fields, astronomy, beam theory, cylindrical shells, hydrodynamics and hydro-magnetic stability, engineering, applied physics, fluid dynamics, and applied mathematics. The recommended method transforms the boundary problem to a system of linear equations. The algorithm we are going to develop in this paper is not only simply the approximation solution of the 8th order BVPs using Cubic-B spline but it also describes the estimated derivatives of 1st order to 8th order of the analytic solution. The strategy is effectively applied to numerical examples and the outcomes are compared with the existing results. The method proposed in this paper provides better approximations to the exact solution.

This paper constructs the generalized viscosity implicit midpoint rule for nonexpansive mappings in Banach space. It obtains strong convergence conclusions for the proposed algorithm and promotes the related results in this field. Moreover, this paper gives some applications. Finally, the paper gives six numerical examples to support the main results.

In the fields of continuous casting and the roll stepped cooling, the heat transfer coefficient is piecewise linear. However, few papers discuss the solution of the backward heat conduction problem in this situation. Therefore, the aim of this paper is to solve the backward heat conduction problem, which has the piecewise linear heat transfer coefficient. Firstly, the ill-posed of this problem is discussed and the truncated regularized optimization scheme is introduced to solve this problem. Secondly, because the regularization parameter is the key factor for the regularization method, this paper presents an improved method for choosing the regularization parameter to reduce the iterative number and proves the fourth-order convergence of this method. Furthermore, the numerical simulation experiments show that, compared with other methods, the improved method of fourth-order convergence effectively reduces the iterative number. Finally, the truncated regularized optimization scheme is used to estimate the initial temperature, and the results of numerical simulation experiments illustrate that the inverse values match the exact values very well.

This paper is concerned with a class of the discrete Mackey–Glass model that describes the process of the production of blood cells. Prior to proceeding to the main results, we prove the boundedness and extinction of its solutions. By means of the contraction mapping principle and under appropriate assumptions, we prove the existence of almost periodic positive solutions. Furthermore and by the implementation of the discrete Lyapunov functional, sufficient conditions are established for the exponential convergence of the almost periodic positive solution. Examples, as well as numerical simulations are illustrated to demonstrate the effectiveness of the theoretical findings of the paper. Our results are new and generalize some previously-reported results in the literature.

This paper uses the viscosity implicit midpoint rule to find common points of the fixed point set of a nonexpansive mapping and the zero point set of an accretive operator in Banach space. Under certain conditions, this paper obtains the strong convergence results of the proposed algorithm and improves the relevant results of researchers in this field. In the end, this paper gives numerical examples to support the main results.

The fractional calculus is useful in describing the natural phenomena with memory effect. This paper addresses the fractional form of Ambartsumian equation with a delay parameter. It may be a challenge to obtain accurate approximate solution of such kinds of fractional delay equations. In the literature, several attempts have been conducted to analyze the fractional Ambartsumian equation. However, the previous approaches in the literature led to approximate power series solutions which converge in subdomains. Such difficulties are solved in this paper via the Homotopy Perturbation Method (HPM). The present approximations are expressed in terms of the Mittag-Leffler functions which converge in the whole domain of the studied model. The convergence issue is also addressed. Several comparisons with the previous published results are discussed. In particular, while the computed solution in the literature is physical in short domains, with our approach it is physical in the whole domain. The results reveal that the HPM is an effective tool to analyzing the fractional Ambartsumian equation.

This paper further improves the Lie group method with Magnus expansion proposed in a previous paper by the authors, to solve some types of direct singular Sturm&ndash, Liouville problems. Next, a concrete implementation to the inverse Sturm&ndash, Liouville problem algorithm proposed by Barcilon (1974) is provided. Furthermore, computational feasibility and applicability of this algorithm to solve inverse Sturm&ndash, Liouville problems of higher order (for n=2,4) are verified successfully. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied. In conclusion, this work provides a method that can be adapted successfully for solving a direct (regular/singular) or inverse Sturm&ndash, Liouville problem (SLP) of an arbitrary order with arbitrary boundary conditions.

In the paper, a survey of the main results concerning univariate and multivariate exponential power (EP) distributions is given, with main attention paid to mixture representations of these laws. The properties of mixing distributions are considered and some asymptotic results based on mixture representations for EP and related distributions are proved. Unlike the conventional analytical approach, here the presentation follows the lines of a kind of arithmetical approach in the space of random variables or vectors. Here the operation of scale mixing in the space of distributions is replaced with the operation of multiplication in the space of random vectors/variables under the assumption that the multipliers are independent. By doing so, the reasoning becomes much simpler, the proofs become shorter and some general features of the distributions under consideration become more vivid. The first part of the paper concerns the univariate case. Some known results are discussed and simple alternative proofs for some of them are presented as well as several new results concerning both EP distributions and some related topics including an extension of Gleser’s theorem on representability of the gamma distribution as a mixture of exponential laws and limit theorems on convergence of the distributions of maximum and minimum random sums to one-sided EP distributions and convergence of the distributions of extreme order statistics in samples with random sizes to the one-sided EP and gamma distributions. The results obtained here open the way to deal with natural multivariate analogs of EP distributions. In the second part of the paper, we discuss the conventionally defined multivariate EP distributions and introduce the notion of projective EP (PEP) distributions. The properties of multivariate EP and PEP distributions are considered as well as limit theorems establishing the conditions for the convergence of multivariate statistics constructed from samples with random sizes (including random sums of random vectors) to multivariate elliptically contoured EP and projective EP laws. The results obtained here give additional theoretical grounds for the applicability of EP and PEP distributions as asymptotic approximations for the statistical regularities observed in data in many fields.

In the present paper, we consider an important problem from the point of view of application in sciences and engineering, namely, a new class of nonlinear Love-equation with infinite memory in the presence of source term that takes general nonlinearity form. New minimal conditions on the relaxation function and the relationship between the weights of source term are used to show a very general decay rate for solution by certain properties of convex functions combined with some estimates. Investigations on the propagation of surface waves of Love-type have been made by many authors in different models and many attempts to solve Love&rsquo, s equation have been performed, in view of its wide applicability. To our knowledge, there are no decay results for damped equations of Love waves or Love type waves. However, the existence of solution or blow up results, with different boundary conditions, have been extensively studied by many authors. Our interest in this paper arose in the first place in consequence of a query for a new decay rate, which is related to those for infinite memory &piv, in the third section. We found that the system energy decreased according to a very general rate that includes all previous results.

The aim of this paper is to study the Lagrange interpolation on the unit circle taking only into account the separation properties of the nodal points. The novelty of this paper is that we do not consider nodal systems connected with orthogonal or paraorthogonal polynomials, which is an interesting approach because in practical applications this connection may not exist. A detailed study of the properties satisfied by the nodal system and the corresponding nodal polynomial is presented. We obtain the relevant results of the convergence related to the process for continuous smooth functions as well as the rate of convergence. Analogous results for interpolation on the bounded interval are deduced and finally some numerical examples are presented.

In this paper, first, we introduce a path for a convex combination of a pseudocontractive type of mappings with a perturbed mapping and prove strong convergence of the proposed path in a real reflexive Banach space having a weakly continuous duality mapping. Second, we propose two modified implicit iterative methods with a perturbed mapping for a continuous pseudocontractive mapping in the same Banach space. Strong convergence theorems for the proposed iterative methods are established. The results in this paper substantially develop and complement the previous well-known results in this area.

The convergence using the fuzzy semi-metric and dual fuzzy semi-metric is studied in this paper. The infimum type of dual fuzzy semi-metric and the supremum type of dual fuzzy semi-metric are proposed in this paper. Based on these two types of dual fuzzy semi-metrics, the different types of triangle inequalities can be obtained. We also study the convergence of these two types of dual fuzzy semi-metrics.

In this paper, we present the application of the homotopy analysis method for solving integrodifferential equations. In this method, a series is created, the successive elements of which are determined by calculating the appropriate integral of the previous element. In this elaboration, we prove that, if this series is convergent, then its sum is the solution of the objective equation. We formulate and prove the sufficient condition of this convergence, and we give also the estimation of error of an approximate solution obtained by taking the partial sum of the considered series. Moreover, we present in this paper the example of using the investigated method for determining the vibrations of the freely supported reinforced concrete beam as well as for solving the equation of movement of the electromagnet jumper mechanical system.

In this paper, we consider a resolvent operator which depends on the composition of two mappings with &oplus, operation. We prove some of the properties of the resolvent operator, that is, that it is single-valued as well as Lipschitz-type-continuous. An existence and convergence result is proven for a generalized implicit set-valued variational inclusion problem with &oplus, operation. Some special cases of a generalized implicit set-valued variational inclusion problem with &oplus, operation are discussed. An example is constructed to illustrate some of the concepts used in this paper.

In this paper, we investigate a class of nonlinear Langevin equations involving two fractional orders with nonlocal integral and three-point boundary conditions. Using the Banach contraction principle, Krasnoselskii’s and the nonlinear alternative Leray Schauder theorems, the existence and uniqueness results of solutions are proven. The paper was appended examples which illustrate the applicability of the results.

Probabilistic evolution theory (PREVTH) forms a framework for the solution of explicit ODEs. The purpose of the paper is two-fold: (1) conversion of multinomial right-hand sides of the ODEs to purely second degree multinomial right-hand sides by space extension, (2) decrease the computational burden of probabilistic evolution theory by using the condensed Kronecker product. A first order ODE set with multinomial right-hand side functions may be converted to a first order ODE set with purely second degree multinomial right-hand side functions at the expense of an increase in the number of equations and unknowns. Obtaining purely second degree multinomial right-hand side functions is important because the solution of such equation set may be approximated by probabilistic evolution theory. A recent article by the authors states that the ODE set with the smallest number of unknowns can be found by searching. This paper gives the details of a way to search for the optimal space extension. As for the second purpose of the paper, the computational burden can be reduced by considering the properties of the Kronecker product of vectors and how the Kronecker product appears within the recursion of PREVTH: as a Cauchy product structure.

The Colebrook-White equation is often used for calculation of the friction factor in turbulent regimes; it has succeeded in attracting a great deal of attention from researchers. The Colebrook–White equation is a complex equation where the computation of the friction factor is not direct, and there is a need for trial-error methods or graphical solutions; on the other hand, several researchers have attempted to alter the Colebrook-White equation by explicit formulas with the hope of achieving zero-percent (0%) maximum deviation, among them Dejan Brkić and Pavel Praks. The goal of this paper is to discuss the results proposed by the authors in their paper:” Accurate and Efficient Explicit Approximations of the Colebrook Flow Friction Equation Based on the Wright ω-Function” and to propose more accurate formulas.

This paper focuses on the problem of finite-time stabilization of homogeneous, non-Lipschitz systems with dilations. A key contribution of this paper is the design of a virtual recursive Hölder, non-Lipschitz state feedback which renders the non-Lipschitz systems in the special case dominated by a lower-triangular nonlinear system, finite-time stable. The proof is based on a recursive design algorithm developed recently, to construct the virtual Hölder continuous, finite-time stabilizer as well as a C1 positive definite and proper Lyapunov function that guarantees finite-time stability of the non-Lipschitz non-linear systems.

This is a survey paper illuminating the distinguished role of the Mittag-Leffler function and its generalizations in fractional analysis and fractional modeling. The content of the paper is connected to the recently published monograph by Rudolf Gorenflo, Anatoly Kilbas, Francesco Mainardi and Sergei Rogosin.

The rate of change of any function versus its independent variables was defined as a derivative. The fundamentals of the derivative concept were constructed by Newton and l’Hôpital. The followers of Newton and l’Hôpital defined fractional order derivative concepts. We express the derivative defined by Newton and l’Hôpital as an ordinary derivative, and there are also fractional order derivatives. So, the derivative concept was handled in this paper, and a new definition for derivative based on indefinite limit and l’Hôpital’s rule was expressed. This new approach illustrated that a derivative operator may be non-linear. Based on this idea, the asymptotic behaviors of functions were analyzed and it was observed that the rates of changes of any function attain maximum value at inflection points in the positive direction and minimum value (negative) at inflection points in the negative direction. This case brought out the fact that the derivative operator does not have to be linear; it may be non-linear. Another important result of this paper is the relationships between complex numbers and derivative concepts, since both concepts have directions and magnitudes.

Using some of the extended fixed point results for Geraghty contractions in b-metric spaces given by Faraji, Savić and Radenović and their idea to apply these results to nonlinear integral equations, in this paper we present some existence and uniqueness conditions for the solution of a nonlinear Fredholm&ndash, Volterra integral equation with a modified argument.

In Mombeni et al. (2019), Birnbaum-Saunders and Weibull kernel estimators were introduced for the estimation of cumulative distribution functions (c.d.f.s) supported on the half-line $[0,\infty)$. They were the first authors to use asymmetric kernels in the context of c.d.f. estimation. Their estimators were shown to perform better numerically than traditional methods such as the basic kernel method and the boundary modified version from Tenreiro (2013). In the present paper, we complement their study by introducing five new asymmetric kernel c.d.f. estimators, namely the Gamma, inverse Gamma, lognormal, inverse Gaussian and reciprocal inverse Gaussian kernel c.d.f. estimators. For these five new estimators, we prove the asymptotic normality and we find asymptotic expressions for the following quantities: bias, variance, mean squared error and mean integrated squared error. A numerical study then compares the performance of the five new c.d.f. estimators against traditional methods and the Birnbaum-Saunders and Weibull kernel c.d.f. estimators from Mombeni et al. (2019). By using the same experimental design, we show that the lognormal and Birnbaum-Saunders kernel c.d.f. estimators perform the best overall, while the other asymmetric kernel estimators are sometimes better but always at least competitive against the boundary kernel method., 38 pages, 2 tables, 9 figures

In this paper, we study the dynamics of a fractional-order epidemic model with general nonlinear incidence rate functionals and time-delay. We investigate the local and global stability of the steady-states. We deduce the basic reproductive threshold parameter, so that if R0&lt, 1, the disease-free steady-state is locally and globally asymptotically stable. However, for R0&gt, 1, there exists a positive (endemic) steady-state which is locally and globally asymptotically stable. A Holling type III response function is considered in the numerical simulations to illustrate the effectiveness of the theoretical results.

In this paper, we introduce a new class of multi-dimensional robust optimization problems (named (P)) with mixed constraints implying second-order partial differential equations (PDEs) and inequations (PDIs). Moreover, we define an auxiliary (modified) class of robust control problems (named (P)(b¯,c¯)), which is much easier to study, and provide some characterization results of (P) and (P)(b¯,c¯) by using the notions of normal weak robust optimal solution and robust saddle-point associated with a Lagrange functional corresponding to (P)(b¯,c¯). For this aim, we consider path-independent curvilinear integral cost functionals and the notion of convexity associated with a curvilinear integral functional generated by a controlled closed (complete integrable) Lagrange 1-form.

Compared with the linear Black–Scholes model, nonlinear models are constructed through taking account of more practical factors, such as transaction cost, and so it is difficult to find an exact analytical solution. Combining the Haar wavelet integration method, which can transform the partial differential equation into the system of algebraic equations, the homotopy perturbation method, which can linearize the nonlinear problems, and the variational iteration method, which can solve the large system of algebraic equations efficiently, a novel numerical method for the nonlinear Black–Scholes model is proposed in this paper. Compared with the traditional methods, it has higher efficiency and calculation precision.

The accurate description of the flow of water in porous media is of the greatest importance due to its numerous applications in several areas (groundwater, soil mechanics, etc.). The nonlinear Richards equation is often used as the governing equation that describes this phenomenon and a large number of research studies aimed to solve it numerically. However, due to the nonlinearity of the constitutive expressions for permeability, it remains a challenging modeling problem. In this paper, the stationary form of Richards’ equation used in saturated soils is solved by two numerical methods: generalized finite differences, an emerging method that has been successfully applied to the transient case, and a finite element method, for benchmarking. The nonlinearity of the solution in both cases is handled using a Newtonian iteration. The comparative results show that a generalized finite difference iteration yields satisfactory results in a standard test problem with a singularity at the boundary.

The paper deals with a two-person zero-sum differential game for a dynamical system described by differential equations with the Caputo fractional derivatives of an order α∈(0,1) and a Bolza-type cost functional. A relationship between the differential game and the Cauchy problem for the corresponding Hamilton–Jacobi–Bellman–Isaacs equation with fractional coinvariant derivatives of the order α and the natural boundary condition is established. An emphasis is given to construction of optimal positional (feedback) strategies of the players. First, a smooth case is studied when the considered Cauchy problem is assumed to have a sufficiently smooth solution. After that, to cope with a general non-smooth case, a generalized minimax solution of this problem is involved.

This paper focuses on a method to construct wavelet Riesz bases with homogeneous boundary condition and use them to a kind of second-order elliptic equation. First, we construct the splines on the interval [0,1] and consider their approximation properties. Then we define the wavelet bases and illustrate the condition numbers of stiffness matrices are small and bounded. Finally, several numerical examples show that our approach performs efficiently.

For the aggregation problem of attributes with a correlation relationship, it is often necessary to take the correlation factor into account in order to make the decision results more objective and reasonable. The Heronian mean is an aggregation operator which reflects the interaction between attributes. It is of great theoretical and practical significance to study and popularize the multiple attribute decision-making methods based on the Heronian mean operator. In this paper, we first give a new three-parameter generalized weighted Heronian mean (TPGWHM), which has a series of excellent properties such as idempotency, monotonicity and boundedness. At the same time, the relationship between the TPGWHM and the existing aggregation operators is given. Then, we propose the intuitionistic fuzzy three-parameter generalized weighted Heronian mean (IFTPGWHM) and give its idempotency, monotonicity, boundedness and limit properties. On this basis, a multiple attribute decision-making method based on the TPGWHM and a multiple attribute decision-making method based on the IFTPGWHM are given, and corresponding examples are given and analyzed.

Linear independence of the blending functions is a necessary requirement for T-spline in isogeometric analysis. The main work in this paper focuses on the analysis about T-splines of degree one, we demonstrate that all the blending functions of such T-spline of degree one are linearly independent. The advantage owned by one degree T-spline is that it can avoid the problem of judging whether the model is analysis-suitable or not, especially for occasions that need a quick response from the analysis results. This may provide a new way of using T-spline for a CAD and CAE integrating scenario, since one degree T-spline still guarantees the topology flexibility and is compatible with the spline-based modeling system. In addition, we compare the numerical approximations of isogeometric analysis and finite element analysis, and the experiment indicates that isogeometric analysis using T-spline of degree one can reach a comparable result with classical method.