Search

We congratulate the authors Bertsimas, Pauphilet and van Parys (hereafter BPvP) and Hastie, Tibshirani and Tibshirani (hereafter HTT) for providing fresh and insightful views on the problem of variable selection and prediction in linear models. Their contributions at the fundamental level provide guidance for more complex models and procedures.

The aim of this paper is to discuss the main result in the paper by D.Y. Gao and X. Lu [On the extrema of a nonconvex functional with double-well potential in 1D, Z. Angew. Math. Phys. (2016) 67:62]. More precisely we provide a detailed study of the problem considered in that paper, pointing out the importance of the norm on the space $C^{1}[a,b]$; because no norm (topology) is mentioned on $C^{1}[a,b]$ we look at it as being a subspace of $W^{1,p}(a,b)$ for $p\in [1,\infty]$ endowed with its usual norm. We show that the objective function has not local extrema with the mentioned constraints for $p\in [1,4)$, and has (up to an additive constant) only a local maximizer for $p=\infty$, unlike the conclusion of the main result of the discussed paper where it is mentioned that there are (up to additive constants) two local minimizers and a local maximizer. We also show that the same conclusions are valid for the similar problem treated in the preprint by X. Lu and D.Y. Gao [On the extrema of a nonconvex functional with double-well potential in higher dimensions, arXiv:1607.03995]., 12 pages; in this version we added the forgotten condition $F(x) \ne 0$ for $x\in (a,b)$ on page 3

We show that a relative entropy condition recently shown by Leger and Vasseur to imply uniqueness and stable $L^2$ dependence on initial data of Lax 1- or $n$-shock solutions of an $n\times n$ system of hyperbolic conservation laws with convex entropy implies Lopatinski stability in the sense of Majda. This means in particular that Leger and Vasseur's relative entropy condition represents a considerable improvement over the standard entropy condition of decreasing shock strength and increasing entropy along forward Hugoniot curves, which, in a recent example exhibited by Barker, Freist\"uhler and Zumbrun, was shown to fail to imply Lopatinski stability, even for systems with convex entropy. This observation bears also on the parallel question of existence, at least for small $BV$ or $H^s$ perturbations, Comment: to appear in Proceedings of the AMS

In this paper, we discuss the validity of the result of Khuangsatung and Kangtunyakarn [Existence and convergence theorem for fixed point problem of various nonlinear mappings and variational inequality problems without some assumptions, Filomat 32(1) (2018) 305–309].

The abstract results and applications presented in “Controllability of fractional neutral stochastic functional differential systems, Z. Angew. Math. Phys. 65 (2014), no. 5, 941–959, are not correct. Moreover, the class of differential control problems studied in [1] is not H-controllable.

Congruences mod p for a prime p and partial products of the numbers 1,…, p − 1 are obtained. Bibliography: 2 titles.

The preeminent target of present study is to reveal the speed characteristic of ongoing outbreak COVID-19 due to novel coronavirus. On January 2020, the novel coronavirus infection (COVID-19) detected in India, and the total statistic of cases continuously increased to 7 128 268 cases including 109 285 deceases to October 2020, where 860 601 cases are active in India. In this study, we use the Hermite wavelets basis in order to solve the COVID-19 model with time- arbitrary Caputo derivative. The discussed framework is based upon Hermite wavelets. The operational matrix incorporated with the collocation scheme is used in order to transform arbitrary-order problem into algebraic equations. The corrector scheme is also used for solving the COVID-19 model for distinct value of arbitrary order. Also, authors have investigated the various behaviors of the arbitrary-order COVID-19 system and procured developments are matched with exiting developments by various techniques. The various illustrations of susceptible, exposed, infected, and recovered individuals are given for its behaviors at the various value of fractional order. In addition, the proposed model has been also supported by some numerical simulations and wavelet-based results.

Funder: John Simon Guggenheim Memorial Foundation; doi: https://doi.org/10.13039/100005851, Funder: Einstein Visiting Fellow Program, The Rayleigh–Bénard system with stress-free boundary conditions is shown to have a global attractor in each affine space where velocity has fixed spatial average. The physical problem is shown to be equivalent to one with periodic boundary conditions and certain symmetries. This enables a Gronwall estimate on enstrophy. That estimate is then used to bound the L 2 norm of the temperature gradient on the global attractor, which, in turn, is used to find a bounding region for the attractor in the enstrophy–palinstrophy plane. All final bounds are algebraic in the viscosity and thermal diffusivity, a significant improvement over previously established estimates. The sharpness of the bounds are tested with numerical simulations.

The first symptomatic infected individuals of coronavirus (Covid-19) was confirmed in December 2020 in the city of Wuhan, China. In India, the first reported case of Covid-19 was confirmed on 30 January 2020. Today, coronavirus has been spread out all over the world. In this manuscript, we studied the coronavirus epidemic model with a true data of India by using Predictor-Corrector scheme. For the proposed model of Covid-19, the numerical and graphical simulations are performed in a framework of the new generalised Caputo sense non-integer order derivative. We analysed the existence and uniqueness of solution of the given fractional model by the definition of Chebyshev norm, Banach space, Schauder's second fixed point theorem, Arzel's-Ascoli theorem, uniform boundedness, equicontinuity and Weissinger's fixed point theorem. A new analysis of the given model with the true data is given to analyse the dynamics of the model in fractional sense. Graphical simulations show the structure of the given classes of the non-linear model with respect to the time variable. We investigated that the mentioned method is copiously strong and smooth to implement on the systems of non-linear fractional differential equation systems. The stability results for the projected algorithm is also performed with the applications of some important lemmas. The present study gives the applicability of this new generalised version of Caputo type non-integer operator in mathematical epidemiology. We compared that the fractional order results are more credible to the integer order results.

We propose a complete-search algorithm for solving a class of non-convex, possibly infinite-dimensional, optimization problems to global optimality. We assume that the optimization variables are in a bounded subset of a Hilbert space, and we determine worst-case run-time bounds for the algorithm under certain regularity conditions of the cost functional and the constraint set. Because these run-time bounds are independent of the number of optimization variables and, in particular, are valid for optimization problems with infinitely many optimization variables, we prove that the algorithm converges to an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon $$\end{document}ε-suboptimal global solution within finite run-time for any given termination tolerance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon > 0$$\end{document}ε>0. Finally, we illustrate these results for a problem of calculus of variations.

Part of the work was done while the authors were visiting the South China University of Science and Technology (SUSTech), Shenzhen, in 2018, and we are grateful (in particular to Professor Cai Heng Li) for the hospitality that we received.The authors would like to thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and hospitality during the programme Groups, representations and applications: new perspectives (supported by EPSRC grant no.EP/R014604/1), where further work on this paper was undertaken. In particular we acknowledge a Simons Fellowship (Cameron) and a Kirk Distinguished Visiting Fellowship (Praeger) during this programme. Schneider thanks the Centre for the Mathematics of Symmetry and Computation of The University of Western Australia and Australian Research Council Discovery Grant DP160102323 for hosting his visit in 2017 and acknowledges the support of the CNPq projects Produtividade em Pesquisa (project no.: 308212/2019-3) and Universal (project no.:421624/2018-3). Diagonal groups are one of the classes of finite primitive permutation groups occurring in the conclusion of the O'Nan-Scott theorem. Several of the other classes have been described as the automorphism groups of geometric or combinatorial structures such as affine spaces or Cartesian decompositions, but such structures for diagonal groups have not been studied in general. The main purpose of this paper is to describe and characterise such structures, which we call diagonal semilattices. Unlike the diagonal groups in the O'Nan-Scott theorem, which are defined over finite characteristically simple groups, our construction works over arbitrary groups, finite or infinite. A diagonal semilattice depends on a dimension m and a group T. For m=2, it is a Latin square, the Cayley table of T, though in fact any Latin square satisfies our combinatorial axioms. However, for m≥3, the group T emerges naturally and uniquely from the axioms. (The situation somewhat resembles projective geometry, where projective planes exist in great profusion but higher-dimensional structures are coordinatised by an algebraic object, a division ring.) A diagonal semilattice is contained in the partition lattice on a set Ω, and we provide an introduction to the calculus of partitions. Many of the concepts and constructions come from experimental design in statistics. We also determine when a diagonal group can be primitive, or quasiprimitive (these conditions turn out to be equivalent for diagonal groups). Associated with the diagonal semilattice is a graph, the diagonal graph, which has the same automorphism group as the diagonal semilattice except in four small cases with m

Differential-algebraic equations (DAEs) are widely used for modeling of dynamical systems. In numerical analysis of DAEs, consistent initialization and index reduction are important preprocessing prior to numerical integration. Existing DAE solvers commonly adopt structural preprocessing methods based on combinatorial optimization. Unfortunately, the structural methods fail if the DAE has numerical or symbolic cancellations. For such DAEs, methods have been proposed to modify them to other DAEs to which the structural methods are applicable, based on the combinatorial relaxation technique. Existing modification methods, however, work only for a class of DAEs that are linear or close to linear. This paper presents two new modification methods for nonlinear DAEs: the substitution method and the augmentation method. Both methods are based on the combinatorial relaxation approach and are applicable to a large class of nonlinear DAEs. The substitution method symbolically solves equations for some derivatives based on the implicit function theorem and substitutes the solution back into the system. Instead of solving equations, the augmentation method modifies DAEs by appending new variables and equations. The augmentation method has advantages that the equation solving is not needed and the sparsity of DAEs is retained. It is shown in numerical experiments that both methods, especially the augmentation method, successfully modify high-index DAEs that the DAE solver in MATLAB cannot handle., Comment: A preliminary version of this paper is to appear in Proceedings of the 44th International Symposium on Symbolic and Algebraic Computation (ISSAC 2019), Beijing, China, July 2019

In a recent paper [JUKIĆ, D.: A necessary and sufficient criterion for the existence of the global minima of a continuous lower bounded function on a noncompact set, J. Comput. Appl. Math. 375 (2020)], a new existence level was introduced and then was used to obtain a necessary and sufficient criterion for the existence of the global minima of a continuous lower bounded function on a noncompact set. In this paper, we determined that existence level for the residual sum of squares of the power-law regression with an unknown location parameter, and so we obtained a necessary and sufficient condition which guarantee the existence of the least squares estimate.

Diffusion equations with Riemann–Liouville fractional derivatives are Volterra integro-partial differential equations with weakly singular kernels and present fundamental challenges for numerical computation. In this paper, we make a convergence analysis of the Schwarz waveform relaxation (SWR) algorithms with Robin transmission conditions (TCs) for these problems. We focus on deriving good choice of the parameter involved in the Robin TCs, at the continuous and fully discretized level. Particularly, at the space-time continuous level, we show that the derived Robin parameter is much better than the one predicted by the well-understood equioscillation principle. At the fully discretized level, the problem of determining a good Robin parameter is studied in the convolution quadrature framework, which permits us to precisely capture the effects of different temporal discretization methods on the convergence rate of the SWR algorithms. The results obtained in this paper will be preliminary preparations for our further study of the SWR algorithms for integro-partial differential equations.

In order to solve the complicated process and low efficiency and low accuracy of solving a class of matrix equations, this paper introduces the linear saturated system model neural network architecture to solve the bisymmetric solution of a class of matrix equations. Firstly, a class of matrix equations is constructed to determine the key problems of solving the equations. Secondly, the linear saturated system model neural network structure is constructed to determine the characteristic parameters in the process of bisymmetric solution. Then, the matrix equations is solved by using backpropagation neural network topology. Finally, the class normalization is realized by using the objective function of bisymmetric solution, and the bisymmetric solution of a class of matrix equations is realized. In order to verify the solving effect of the method in this paper, three indexes (accuracy, correction accuracy, and solving time) are designed in the experiment. The experimental results show that the proposed method can effectively reduce the solving time, can improve the accuracy and correction effect of the bisymmetric solution, and has high practicability.

PurposeIn this paper, the authors applied the empirical likelihood method, which was originally proposed by Owen, to the copula moment based estimation methods to take advantage of its properties, effectiveness, flexibility and reliability of the nonparametric methods, which have limiting chi-square distributions and may be used to obtain tests or confidence intervals. The authors derive an asymptotically normal estimator of the empirical likelihood based on copula moment estimation methods (ELCM). Finally numerical performance with a simulation experiment of ELCM estimator is studied and compared to the CM estimator, with a good result.Design/methodology/approachIn this paper we applied the empirical likelihood method which originally proposed by Owen, to the copula moment based estimation methods.FindingsWe derive an asymptotically normal estimator of the empirical likelihood based on copula moment estimation methods (ELCM). Finally numerical performance with a simulation experiment of ELCM estimator is studied and compared to the CM estimator, with a good result.Originality/valueIn this paper we applied the empirical likelihood method which originally proposed by Owen 1988, to the copula moment based estimation methods given by Brahimi and Necir 2012. We derive an new estimator of copula parameters and the asymptotic normality of the empirical likelihood based on copula moment estimation methods.

In recent years, the problem of estimating Hankel determinants has attracted the attention of many mathematicians. Their research have been focused mainly on deriving the bounds of H 2 , 2 {H_{2,2}} or H 3 , 1 {H_{3,1}} over different subclasses of 𝒮 {\mathcal{S}} . Only in a few papers third Hankel determinants for non-univalent functions were considered. In this paper, we consider two classes of analytic functions with real coefficients. The first one is the class 𝒯 {\mathcal{T}} of typically real functions. The second object of our interest is 𝒦 ℝ ⁢ ( i ) {\mathcal{K}_{\mathbb{R}}(i)} , the class of functions with real coefficients which are convex in the direction of the imaginary axis. In both classes, we find lower and upper bounds of the third Hankel determinant. The results are sharp.

In the present paper we describe the family of all closed convex plane curves of class $$C^1$$ C 1 which have circles as their isoptics. In the first part of the paper we give a certain characterization of all ellipses based on the notion of isoptic and we give a geometric characterization of curves whose orthoptics are circles. The second part of the paper contains considerations on curves which have circles as their isoptics and we show the form of support functions of all considered curves.

In this paper, we study a problem posed by Furstenberg on intersections between × 2 , × 3 \times 2, \times 3 invariant sets. We present here a direct geometrical counting argument to revisit a theorem of Wu and Shmerkin. This argument can be used to obtain further improvements. For example, we show that if A 2 , A 3 ⊂ [ 0 , 1 ] A_2,A_3\subset [0,1] are closed and × 2 , × 3 \times 2, \times 3 invariant respectively, assuming that dim ⁡ A 2 + dim ⁡ A 3 > 1 \dim A_2+\dim A_3>1 then A 2 ∩ ( u A 3 + v ) A_2\cap (uA_3+v) is sparse (defined in this paper) and has box dimension zero uniformly with respect to the real parameters u , v u,v such that u u and u − 1 u^{-1} are both bounded away from 0 0 .

In this paper, we develop the theory of a class of casual differential equations on a time scale. An existence theorem for casual differential equations on a time scale is given under mixed Lipschitz and compactness conditions by the fixed point theorem in Banach algebra due to Dhage. Some fundamental differential inequalities on a time scale are also presented which are utilized to investigate the existence of extremal solutions. The comparison principle on casual differential equations on a time scale is established. Our results in this paper extend and improve some well-known results.

This paper is a survey of results on countable families with natural degree spectra. These results were obtained by a modification of the methodology proposed by Wechner, who first found a family of sets with the spectrum consisting precisely of nonzero Turing degrees. Based on this method, many researchers obtained examples of families with other natural spectra. In addition, in this paper we extend these results and present new examples of natural spectra. In particular, we construct a family of finite sets with the spectrum consisting of exactly non-K-trivial degrees and also we find new sufficient conditions on $$ {\Delta}_2^0 $$ -degree a, which guarantees that the class {x : x ≰ a} is the degree spectrum of some family. Finally, we give a survey of our recent results on the degree spectra of α-families, where α is an arbitrary computable ordinal.

This paper presents some new formulas and properties of the generalized fractional derivative with non-singular kernel that covers various types of fractional derivatives such as the Caputo–Fabrizio fractional derivative, the Atangana–Baleanu fractional derivative, and the weighted Atangana–Baleanu fractional derivative. These new properties extend many recent results existing in the literature. Furthermore, the paper proposes some interesting inequalities that estimate the generalized fractional derivatives of some specific functions. These inequalities can be used to construct Lyapunov functions with the aim to study the global asymptotic stability of several fractional-order systems arising from diverse fields of science and engineering.

The paper focuses on the analysis of the Euler projection Galerkin finite element method (FEM) for the dynamics of magnetization in ferromagnetic materials, described by the Landau–Lifshitz equation with the point-wise constraint $|{\textbf{m}}|=1$. The method is based on a simple sphere projection that projects the numerical solution onto a unit sphere at each time step, and the method has been used in many areas in the past several decades. However, error analysis for the commonly used method has not been done since the classical energy approach cannot be applied directly. In this paper we present an optimal $\textbf{L}^2$ error analysis of the backward Euler sphere projection method by using quadratic or higher order finite elements under a time step condition $\tau =O(\epsilon _0 h)$ with some small $\epsilon _0>0$. The analysis is based on more precise estimates of the extra error caused by the sphere projection in both $\textbf{L}^2$ and $\textbf{H}^1$ norms, and the classical estimate of dual norm. Numerical experiment is provided to confirm our theoretical analysis.

The objective of this paper is two-fold. The first objective is to complete the former work of Bensoussan and Frehse [2] . One big limitation of this paper was the fact that they are systems of PDE. on a bounded domain. One can expect solutions to be bounded, since one looks for smooth solutions. This is a very important property for the development of the method. It is true also that solutions which exist in a bounded domain may fail to exist on R n , because of the lack of bounds. We give conditions so that the results of [2] can be extended to R n . The second objective is to consider the BSDE (Backward stochastic differential equations) version of the system of PDE. This is the objective of a more recent work of Xing and Zitkovic [8] . They consider systems of BSDE with quadratic growth, which is a well-known open problem in the BSDE literature. Since the BSDE are Markovian, the problem is equivalent to the analytic one. However, because of this motivation the analytic problem is in R n and not on a bounded domain. Xing and Zitkovic developed a probabilistic approach. The connection between the analytic problem and the BSDE is not apparent. Our objective is to show that the analytic approach can be completely translated into a probabilistic one. Nevertheless probabilistic concepts are also useful, after their conversion into the analytic framework. This is in particular true for the uniqueness result.

In this paper, we establish the global existence and uniqueness of the solution of the Cauchy problem of a one-dimensional compressible isentropic Euler system for a Chaplygin gas with large initial data in the space L loc 1 . The hypotheses on the initial data may be the least requirement to ensure the existence of a weak solution in the Lebesgue measurable sense. The novelty and also the essence of the difficulty of the problem lie in the fact that we have neither the requirement on the local boundedness of the density nor that which is bounded away from vacuum. We develop the previous results on this degenerate system. The method used is Lagrangian representation, the essence of which is characteristic analysis. The key point is to prove the existence of the Lagrangian representation and the absolute continuity of the potentials constructed with respect to the space and the time variables. We achieve this by finding a property of the fundamental theorem of calculus for Lebesgue integration, which is a sufficient and necessary condition for judging whether a monotone continuous function is absolutely continuous. The assumptions on the initial data in this paper are believed to also be necessary for ruling out the formation of Dirac singularity of density. The ideas and techniques developed here may be useful for other nonlinear problems involving similar difficulties.

The paper is devoted to discretization of integral norms of functions from a given finite dimensional subspace. Even though this problem is extremely important in applications, its systematic study has begun only recently. In this paper we obtain a conditional theorem for all integral norms $$L_q$$ , $$1\le q

Full-waveform inversion (FWI) is today a standard process for the inverse problem of seismic imaging. PDE-constrained optimization is used to determine unknown parameters in a wave equation that represent geophysical properties. The objective function measures the misfit between the observed data and the calculated synthetic data, and it has traditionally been the least-squares norm. In a sequence of papers, we introduced the Wasserstein metric from optimal transport as an alternative misfit function for mitigating the so-called cycle skipping, which is the trapping of the optimization process in local minima. In this paper, we first give a sharper theorem regarding the convexity of the Wasserstein metric as the objective function. We then focus on two new issues. One is the necessary normalization of turning seismic signals into probability measures such that the theory of optimal transport applies. The other, which is beyond cycle skipping, is the inversion for parameters below reflecting interfaces. For the first, we propose a class of normalizations and prove several favorable properties for this class. For the latter, we demonstrate that FWI using optimal transport can recover geophysical properties from domains where no seismic waves travel through. We finally illustrate these properties by the realistic application of imaging salt inclusions, which has been a significant challenge in exploration geophysics.

In this paper, the global exponential stability and periodicity are investigated for delayed neural network models with continuous coefficients and piecewise constant delay of generalized type. The sufficient condition for the existence and uniqueness of periodic solutions of the model is established by applying Banach’s fixed point theorem and the successive approximations method. By constructing suitable differential inequalities with generalized piecewise constant delay, some sufficient conditions for the global exponential stability of the model are obtained. Typical numerical examples with simulations are utilized to illustrate the validity and improvement in less conservatism of the theoretical results. This paper ends with a brief conclusion.

In this paper, we study flows of hypersurfaces in hyperbolic space, and apply them to prove geometric inequalities. In the first part of the paper, we consider volume preserving flows by a family of curvature functions including positive powers of $k$-th mean curvatures with $k=1,\cdots,n$, and positive powers of $p$-th power sums $S_p$ with $p>0$. We prove that if the initial hypersurface $M_0$ is smooth and closed and has positive sectional curvatures, then the solution $M_t$ of the flow has positive sectional curvature for any time $t>0$, exists for all time and converges to a geodesic sphere exponentially in the smooth topology. The convergence result can be used to show that certain Alexandrov-Fenchel quermassintegral inequalities, known previously for horospherically convex hypersurfaces, also hold under the weaker condition of positive sectional curvature. In the second part of this paper, we study curvature flows for strictly horospherically convex hypersurfaces in hyperbolic space with speed given by a smooth, symmetric, increasing and homogeneous degree one function $f$ of the shifted principal curvatures $\lambda_i=\kappa_i-1$, plus a global term chosen to impose a constraint on the quermassintegrals of the enclosed domain, where $f$ is assumed to satisfy a certain condition on the second derivatives. We prove that if the initial hypersurface is smooth, closed and strictly horospherically convex, then the solution of the flow exists for all time and converges to a geodesic sphere exponentially in the smooth topology. As applications of the convergence result, we prove a new rigidity theorem on smooth closed Weingarten hypersurfaces in hyperbolic space, and a new class of Alexandrov-Fenchel type inequalities for smooth horospherically convex hypersurfaces in hyperbolic space.

In this paper, we consider boundary-value problems for a new class of hybrid systems of differential equations whose coefficients contain the Dirac delta-function. Hybrid systems are systems that contain both ordinary and partial differential equations; such systems appear, for example, when equations of motion of mechanical systems of rigid bodies attached to a beam by elastic bonds are derived from the Hamilton–Ostrogradsky variational principle. We present examples that lead to such systems and introduce the notions of generalized solutions and eigenvalues of a boundary-value problem. We also compare results of numerical simulations based on methods proposed in this paper with results obtained by previously known methods and show that our approach is reliable and universal.

This is the first of a series of two papers dealing with local limit theorems in relatively hyperbolic groups. In this first paper, we prove rough estimates for the Green function. Along the way, we introduce the notion of relative automaticity which will be useful in both papers and we show that relatively hyperbolic groups are relatively automatic. We also define the notion of spectral positive recurrence for random walks on relatively hyperbolic groups. We then use our estimates for the Green function to prove that $p_n\asymp R^{-n}n^{-3/2}$ for spectrally positive-recurrent random walks, where $p_n$ is the probability of going back to the origin at time n and where R is the inverse of the spectral radius of the random walk.

Non-monotonicity height is an important index to describe the complexity of dynamics for piecewise monotone functions. Although it is used extensively in the theory of iterative roots, its calculation is still difficult especially in the infinite case. In this paper, by introducing the concept of spanning interval, we first present a sufficient condition for piecewise monotone functions to have height infinity and then an algorithm for finding the spanning intervals is given. We further investigate the density of all piecewise monotone functions with infinite and finite height, respectively, and the results indicate the instability of height. At the end of this paper, the variance of height under composition, especially for functions of height 1 and infinity, are also discussed.

PurposeIn this paper, Picard–S hybrid iterative process is defined, which is a hybrid of Picard and S-iterative process. This new iteration converges faster than all of Picard, Krasnoselskii, Mann, Ishikawa, S-iteration, Picard–Mann hybrid, Picard–Krasnoselskii hybrid and Picard–Ishikawa hybrid iterative processes for contraction mappings and to find the solution of delay differential equation, using this hybrid iteration also proved some results for Picard–S hybrid iterative process for nonexpansive mappings.Design/methodology/approachThis new iteration converges faster than all of Picard, Krasnoselskii, Mann, Ishikawa, S-iteration, Picard–Mann hybrid, Picard–Krasnoselskii hybrid, Picard–Ishikawa hybrid iterative processes for contraction mappings.FindingsShowed the fastest convergence of this new iteration and then other iteration defined in this paper. The author finds the solution of delay differential equation using this hybrid iteration. For new iteration, the author also proved a theorem for nonexpansive mapping.Originality/valueThis new iteration converges faster than all of Picard, Krasnoselskii, Mann, Ishikawa, S-iteration, Picard–Mann hybrid, Picard–Krasnoselskii hybrid, Picard–Ishikawa hybrid iterative processes for contraction mappings and to find the solution of delay differential equation, using this hybrid iteration also proved some results for Picard–S hybrid iterative process for nonexpansive mappings.

This is the first of two papers in which we investigate the properties of the displacement functions of automorphisms of free groups (more generally, free products) on Culler-Vogtmann Outer space and its simplicial bordification - the free splitting complex - with respect to the Lipschitz metric. The theory for irreducible automorphisms being well-developed, we concentrate on the reducible case. Since we deal with the bordification, we develop all the needed tools in the more general setting of deformation spaces, and their associated free splitting complexes. In the present paper we study the local properties of the displacement function. In particular, we study its convexity properties and the behaviour at bordification points, by geometrically characterising its continuity-points. We prove that the global-simplex-displacement spectrum of $Aut(F_n)$ is a well-ordered subset of $\mathbb R$, this being helpful for algorithmic purposes. We introduce a weaker notion of train tracks, which we call {\em partial train tracks} (which coincides with the usual one for irreducible automorphisms) and we prove that, for any automorphism, points of minimal displacement - minpoints - coincide with the marked metric graphs that support partial train tracks. We show that any automorphism, reducible or not, has a partial train track (hence a minpoint) either in the outer space or its bordification. We show that, given an automorphism, any of its invariant free factors is seen in a partial train track map. In a subsequent paper we will prove that level sets of the displacement functions are connected, and we will apply that result to solve certain decision problems., 50 pages. Originally part of arXiv:1703.09945 . We decided to split that paper following the recommendations of a referee. Updated subsequent to acceptance by Transactions of the American Mathematical Society

The classical result of L. Székelyhidi states that (under some assumptions) every solution of a general linear equation must be a polynomial function. It is known that Székelyhidi’s result may be generalized to equations where some occurrences of the unknown functions are multiplied by a linear combination of the variables. In this paper we study the equations where two such combinations appear. The simplest nontrivial example of such a case is given by the equation$$\begin{aligned} F(x + y) - F(x) - F(y) = yf(x) + xf(y) \end{aligned}$$F(x+y)-F(x)-F(y)=yf(x)+xf(y)considered by Fechner and Gselmann (Publ Math Debrecen 80(1–2):143–154, 2012). In the present paper we prove several results concerning the systematic approach to the generalizations of this equation.

In this paper we compute the rates of convergence of power series statistical convergence of sequences of positive linear operators. We also investigate some Korovkin type approximation properties of the $$q$$ -Meyer–Konig and Zeller operators and Durrmeyer variant of the $$q$$ -Meyer–Konig and Zeller operators via power series statistical convergence. We show that the approximation results obtained in this paper expand some previous approximation results of the corresponding operators.

The first aim of this paper is to investigate the growth of the entire function defined by the Laplace-Stieltjes transform converges on the whole complex plane. By introducing the concept of generalized order, we obtain two equivalence theorems of Laplace-Stieltjes transforms related to the generalized order, A * and λn. The second purpose of this paper is to study the problem on the approximation of this Laplace-Stieltjes transform. We also obtain some theorems about the generalized order, the error, and the coefficients of Laplace-Stieltjes transforms, which are generalization and improvement of the previous results.

In this paper we study multi-dimensional mean-field backward doubly stochastic differential equations (BDSDEs), that is, BDSDEs whose coefficients depend not only on the solution processes but also on their law. The first part of the paper is devoted to the comparison theorem for multi-dimensional mean-field BDSDEs with Lipschitz conditions. With the help of the comparison result for the Lipschitz case we prove the existence of a solution for multi-dimensional mean-field BDSDEs with an only continuous drift coefficient of linear growth, and we also extend the comparison theorem to such BDSDEs with a continuous coefficient.

A sequence $$(v_1,\ldots ,v_k)$$ of vertices in a graph G without isolated vertices is called a total dominating sequence if every vertex $$v_i$$ in the sequence totally dominates at least one vertex that was not totally dominated by $$\{v_1,\ldots , v_{i-1}\}$$ and $$\{v_1,\ldots ,v_k\}$$ is a total dominating set of G. The length of a shortest such sequence is the total domination number of G ( $$\gamma _{t}(G)$$ ), while the length of a longest such sequence is the Grundy total domination number of G ( $$\gamma _{gr}^t(G)$$ ). In this paper we study graphs with equal total and Grundy total domination numbers. We characterize bipartite graphs with both total and Grundy total dominations number equal to 4, and show that there is no connected chordal graph G with $$\gamma _{t}(G)=\gamma _{gr}^t(G)=4$$ . The main result of the paper is a characterization of bipartite graphs with $$\gamma _{t}(G)=\gamma _{gr}^t(G)=6$$ proved by establishing a surprising correspondence between the existence of such graphs and a classical but still open problem of the existence of certain finite projective planes.

This paper concerns on two types of integral boundary value problems of a nonlinear fractional differential system, i . e ., nonlocal strip integral boundary value problems and coupled integral boundary value problems. With the aid of the monotone iterative method combined with the upper and lower solutions, the existence of extremal system of solutions for the above two types of differential systems is investigated. In addition, a new comparison theorem for fractional differential system is also established, which is crucial for the proof of the main theorem of this paper. At the end, an example explaining how our studies can be used is also given.

A dynamical system $$ {\displaystyle \begin{array}{ll}{u}_{tt}-\Delta u-\nabla 1\mathrm{n}\;\rho \cdot \nabla u=0& in\kern0.6em {\mathrm{\mathbb{R}}}_{+}^3\times \left(0,T\right),\\ {}{\left.u\right|}_{t=0}={\left.{u}_t\right|}_{t=0}=0& in\kern0.6em \overline{{\mathrm{\mathbb{R}}}_{+}^3},\\ {}{\left.{u}_z\right|}_{z=0}=f& for\kern0.36em 0\le t\le T,\end{array}} $$ is under consideration, where ρ = ρ(x, y, z) is a smooth positive function; f = f(x, y, t) is a boundary control; u = uf (x, y, z, t) is a solution. With the system one associates a response operator R : f ↦ uf|z = 0. The inverse problem is to recover the function ρ via the response operator. A short representation of the local version of the BC-method, which recovers ρ via the data given on a part of the boundary, is provided. If ρ is constant, the forward problem is solved in explicit form. In the paper, the corresponding representations for the solutions and response operator are derived. A way to use them for testing the BC-algorithm, which solves the inverse problem, is outlined. The goal of the paper is to extend the circle of the BC-method users, who are interested in numerical realization of methods for solving inverse problems.

The study in this paper is made on the nonlinear fractional differential equation whose nonlinearity involves the explicit fractional order D0+β u(t). The corresponding Green's function is derived first, and then the completely continuous operator is proved. Besides, based on the Schauder's fixed point theorem and the Krasnosel'skii's fixed point theorem, the sufficient conditions for at least one or two existence of positive solutions are established. Furthermore, several other sufficient conditions for at least three, n or 2n-1 positive solutions are also obtained by applying the generalized AveryHenderson fixed point theorem and the Avery-Peterson fixed point theorem. Finally, several simulation examples are provided to illustrate the main results of the paper. In particularly, a novel efficient iterative method is employed for simulating the examples mentioned above, that is, the interesting point of this paper is that the approximation graphics for the solutions are given by using the iterative method.

In order to construct global solutions to two-dimensional (2D for short) Riemann problems for nonlinear hyperbolic systems of conservation laws, it is important to study various types of wave interactions. This paper deals with two types of wave interactions for a 2D nonlinear wave system with a nonconvex equation of state: Rarefaction wave interaction and shock-rarefaction composite wave interaction. In order to construct solutions to these wave interactions, the authors consider two types of Goursat problems, including standard Goursat problem and discontinuous Goursat problem, for a 2D self-similar nonlinear wave system. Global classical solutions to these Goursat problems are obtained by the method of characteristics. The solutions constructed in the paper may be used as building blocks of solutions of 2D Riemann problems.

The Monte Carlo method to solve the Cauchy problem for large systems of linear differential equations is proposed in this paper. Firstly, a quick overview of previously obtained results from applying the approach towards the Fredholm-type integral equations is made. In the main part of the paper, the method is applied towards a linear ODE system that is transformed into an equivalent system of the Volterra-type integral equations, which makes it possible to remove the limitations due to the conditions of convergence of the majorant series. The following key theorems are stated. Theorem 1 provides the necessary compliance conditions that should be imposed upon the transition propability and initial distribution densities that initiate the corresponding Markov chain, for which equality between the mathematical expectation of the estimate and the functional of interest would hold. Theorem 2 formulates the equation that governs the estimate’s variance. Theorem 3 states the Markov chain parameters that minimize the variance of the estimate of the functional. Proofs are given for all three theorems. In the practical part of this paper, the proposed method is used to solve a linear ODE system that describes a closed queueing system of ten conventional machines and seven conventional service persons. The solutions are obtained for systems with both constant and time-dependent matrices of coefficients, where the machine breakdown intensity is time dependent. In addition, the solutions obtained by the Monte Carlo and Runge–Kutta methods are compared. The results are presented in the corresponding tables.

This paper concerns the existence of transonic shock solutions to the 2-D steady compressible Euler system in an almost flat finite nozzle ( in the sense that it is a generic small perturbation of a flat one ), under physical boundary conditions proposed by Courant-Friedrichs in \cite{CourantFriedrichs1948}, in which the receiver pressure is prescribed at the exit of the nozzle. In the resulting free boundary problem, the location of the shock-front is one of the most desirable information one would like to determine. However, the location of the normal shock-front in a flat nozzle can be anywhere in the nozzle so that it provides little information on the possible location of the shock-front when the nozzle's boundary is perturbed. So one of the key difficulties in looking for transonic shock solutions is to determine the shock-front. To this end, a free boundary problem for the linearized Euler system will be proposed, whose solution will be taken as an initial approximation for the transonic shock solution. In this paper, a sufficient condition in terms of the geometry of the nozzle and the given exit pressure is derived which yields the existence of the solutions to the proposed free boundary problem. Once an initial approximation is obtained, a further nonlinear iteration could be constructed and proved to lead to a transonic shock solution., 53 pages

We establish a theory of $Q$-valued functions minimizing a suitable generalization of the Dirichlet integral. In a second paper the theory will be used to approximate efficiently area minimizing currents $\mathrm{mod}(p)$ when $p=2Q$, and to establish a first general partial regularity theorem for every $p$ in any dimension and codimension., 37 pages. First part of a two-papers work aimed at establishing a first general partial regularity theory for area minimizing currents modulo p, for any p and in any dimension and codimension. v3 is the final version, to appear on Comm. Pure Appl. Math

In this paper, we consider the operator sheaf $$ -\Delta +V+\varepsilon {\mathrm{\mathcal{L}}}_{\varepsilon}\left(\lambda \right)+{\lambda}^2 $$ in the space L2(ℝ2), where the real-valued potential V depends only on the first variable x1, e is a small positive parameter, λ is the spectral parameter, $$ {\mathrm{\mathcal{L}}}_{\varepsilon}\left(\lambda \right) $$ is a localized operator bounded with respect to the Laplacian −Δ, and the essential spectrum of this operator is independent of e and contains certain critical points defined as isolated eigenvalues of the operator $$ -\frac{d^2}{dx_1^2}+V\left({x}_1\right) $$ in L2(ℝ). The basic result obtained in this paper states that for small values of e, in neighborhoods of critical points mentioned, isolated eigenvalues of the sheaf considered arise. Sufficient conditions for the existence or absence of such eigenvalues are obtained. The number of arising eigenvalues is determined, and in the case where they exist, the first terms of their asymptotic expansions are found.

In this paper, we investigate a mixed fractional integral boundary value problem with p(t)-Laplacian operator. Firstly, we derive the Green function through the direct computation and obtain the properties of Green function. For $$p(t)\ne $$ constant, under the appropriate conditions of the nonlinear term, we establish the existence result of at least one positive solution of the above problem by means of the Leray–Schauder fixed point theorem. Meanwhile, we also obtain the positive extremal solutions and iterative schemes in view of applying a monotone iterative method. For $$p(t)=$$ constant, by using Guo–Krasnoselskii fixed point theorem, we study the existence of positive solutions of the above problem. These results enrich the ones in the existing literatures. Finally, some examples are included to demonstrate our main results in this paper and we give out an open problem.