Factoring and decomposing a class of linear functional systems

Within a constructive homological algebra approach, we study the factorization and decomposition problems for a class of linear functional (determined, over-determined, under-determined) systems. Using the concept of Ore algebras of functional operators (e.g., ordinary/partial differential operators, shift operators, time-delay operators), we first concentrate on the computation of morphisms from a finitely presented left module M over an Ore algebra to another one M′, where M (resp., M′) is a module intrinsically associated with the linear functional system Ry = 0 (resp., R′z = 0). These morphisms define applications sending solutions of the system R′z = 0 to solutions of R y = 0. We explicitly characterize the kernel, image, cokernel and coimage of a general morphism. We then show that the existence of a non-injective endomorphism of the module M is equivalent to the existence of a non-trivial factorization R = R₂R₁ of the system matrix R. The corresponding system can then be integrated “in cascade”. Under certain conditions, we also show that the system Ry = 0 is equivalent to a system R′z = 0, where R′ is a block-triangular matrix of the same size as R. We show that the existence of idempotents of the endomorphism ring of the module M allows us to reduce the integration of the system Ry = 0 to the integration of two independent systems R₁y₁ = 0 and R₂y₂ = 0. Furthermore, we prove that, under certain conditions, idempotents provide decompositions of the system Ry = 0, i.e., they allow us to compute an equivalent system R′z = 0, where R′ is a block-diagonal matrix of the same size as R. Applications of these results in mathematical physics and control theory are given. Finally, the different algorithms of the paper are implemented in the Maple package Morphisms based on the library oremodules.     

Meromorphic Functions That Share Fixed-Points

In this paper, we prove the following result: Let f(z) and g(z) be two nonconstant meromorphic(entire) functions, n ≥ 11(n ≥ 6) a positive integer. If fⁿ(z)f′(z) and gⁿ(z)g′(z) have the same fixed-points, then either f(z) = c₁e^cz2, g(z) = c₂e^− cz2, where c₁, c₂, and c are three constants satisfying 4(c₁c₂)^n + 1c² = −1, or f(z) ≡ tg(z) for a constant t such that t^n + 1 = 1.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

New types of 3-D systems of quadratic differential equations with chaotic dynamics based on Ricker discrete population model

The wide class of 3-D autonomous systems of quadratic differential equations, in each of which either there is a couple of coexisting limit cycles or there is a couple of coexisting chaotic attractors, is found. In the second case the couple consists of either Lorentz-type attractor and another attractor of a new type or two Lorentz-type attractors. It is shown that the chaotic behavior of any system of the indicated class can be described by the Ricker discrete population model: z_i+1 = z_i exp(r − z_i), r > 0, z_i > 0, i = 0, 1, … . The values of parameters, at which in the 3-D system appears either the couple of limit cycles or the couple of chaotic attractors, or only one limit cycle, or only one sphere-shaped chaotic attractor, are indicated. Examples are given.     

Stable periodic traveling waves for a predator-prey model with non-constant death rate and delay

In this paper we will consider a predator-prey model with a non-constant death rate and distributed delay, described by a partial integro-differential system. The main goal of this work is to prove that the partial integro-differential system has periodic orbitally asymptotically stable solutions in the form of periodic traveling waves; i.e. N(x, t) = N(σt − μ · x), P(x, t) = P(σt − μ · x), where σ > 0 is the angular frequency and μ is the vector number of the plane wave, which propagates in the direction of the vector μ with speed c = σ/∥μ∥; and N(x, t) and P(x, t) are the spatial population densities of the prey and the predator species, respectively. In order to achieve our goal we will use singular perturbation’s techniques.     

Duality for a Class of Multiobjective Control Problems

In this paper, a class of multiobjective control problems is considered, where the objective and constraint functions involved are f(t, x(t), ?(t), y(t), z(t)) with x(t) ∈ Rⁿ, y(t) ∈ Rⁿ, and z(t) ∈ R^m, where x(t) and z(t) are the control variables and y(t) is the state variable. Under the assumption of invexity and its generalization, duality theorems are proved through a parametric approach to related properly efficient solutions of the primal and dual problems.     

Analysis of a semigroup approach in the inverse problem of identifying an unknown parameters

This article presents a semigroup approach to the mathematical analysis of the inverse parameter problems of identifying the unknown parameters p(t) and q in the linear parabolic equation u_t(x, t)  = u_xx + qu_x(x, t) + p(t)u(x, t), with Dirichlet boundary conditions u(0, t) = ψ₀, u(1, t) = ψ₁. The main purpose of this paper is to investigate the distinguishability of the input-output mapping Φ[·]:P→H^1,2[0,T], via semigroup theory. In this paper, it is shown that if the nullspace of the semigroup T(t) consists of only zero function, then the input-output mapping Φ[·] has the distinguishability property. It is also shown that the types of the boundary conditions and the region on which the problem is defined play an important role in the distinguishability property of the mapping. Moreover, under the light of the measured output data u_x(0, t) = f(t) the unknown parameter p(t) at (x, t) = (0, 0) and the unknown coefficient q are determined via the input data. Furthermore, it is shown that measured output data f(t) can be determined analytically by an integral representation. Hence the input-output mapping Φ[·]:P→H^1,2[0,T] is given explicitly interms of the semigroup.     

An adaptive Huber method for non-linear systems of weakly singular second kind Volterra integral equations

Numerical methods for systems of weakly singular Volterra integral equations are rarely considered in the literature, especially if the equations involve non-linear dependencies between unknowns and their integrals. In the present work an adaptive Huber method for such systems is proposed, by extending the method previously formulated for single weakly singular second kind Volterra equations. The method is tested on example systems of integral equations involving integrals with kernels K(t, τ) = (t − τ)^−1/2, K(t, τ) = exp[−λ(t − τ)](t − τ)^−1/2 (where λ > 0), and K(t, τ) = 1. The magnitude of the errors, and practical accuracy orders, observed for IE systems, are comparable to those for single IEs. In cases when the solution vector is not differentiable at t = 0, the estimation of errors at t = 0 is found somewhat less reliable for IE systems, than it was for single IEs. The stability of the IE systems solved appears to be sufficient, in practice, for the numerical stability of the method.     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

The improved linear multistep methods for differential equations with piecewise continuous arguments

This paper deals with the convergence of the linear multistep methods for the equation x′(t) = ax(t) + a₀x([t]). Numerical experiments demonstrate that the 2-step Adams-Bashforth method is only of order p = 0 when applied to the given equation. An improved linear multistep methods is constructed. It is proved that these methods preserve their original convergence order for ordinary differential equations (ODEs) and some numerical experiments are given.     

Distributive equations of implications based on nilpotent triangular norms

In this paper, we explore the distributive equations of implications, both independently and along with other equations. In detail, we consider three classes of equations. (1) By means of the section of I, we give out the sufficient and necessary conditions of solutions for the distributive equation of implication I(x, T(y, z)) = T(I(x, y), (x, z)) based on a nilpotent triangular norm T and an unknown function I, which indicates that there are no continuous solutions satisfying the boundary conditions of implications. Under the assumptions that I is continuous except the vertical section I(0, y), y ∈ [0, 1), we get its complete characterizations. (2) We prove that there are no solutions for the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, I(y, z)) = I(T(x, y), z). (3) We obtain the sufficient and necessary conditions on T and I to be solutions of the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, y) = I(N(y), N(x)).     

Properties of non-simultaneous blow-up in heat equations coupled via different localized sources

This paper deals with u_t = Δu + u^m(x, t)e^pv(0,t), v_t = Δv + u^q(0, t)e^nv(x,t), subject to homogeneous Dirichlet boundary conditions. The complete classification on non-simultaneous and simultaneous blow-up is obtained by four sufficient and necessary conditions. It is interesting that, in some exponent region, large initial data u₀(v₀) leads to the blow-up of u(v), and in some betweenness, simultaneous blow-up occurs. For all of the nonnegative exponents, we find that u(v) blows up only at a single point if m > 1(n > 0), while u(v) blows up everywhere for 0 ? m ? 1 (n = 0). Moreover, blow-up rates are considered for both non-simultaneous and simultaneous blow-up solutions.     

Elliptic complex numbers with dual multiplication

Investigated is a number system in which the square of a basis number: (w)², and the square of its additive inverse: (−w)², are not equal. Termed W space, a vector space over the reals, this number system will be introduced by restating defining relations for complex space C, then changing a defining conjugacy relation from conj(z) + z = 0 in the complexes to conj(z) + z = 1 for W space. This change produces a dual-represented vector space consisting of two dual, isomorphic fields, which are unified under one “context-sensitive” multiplication. Fundamental algebraic and geometric properties will be investigated. W space can be interpreted as a generalization of the complexes but is characterized by an interacting duality which seems to produce two of everything: two representations, two multiplications, two norm values, and two solutions to a linear equation. W space will be compared to a previous suggestion of a similar algebra, and then possible applications will be offered, including a W space fractal.     

Analysis of nonlinear integral equations with Erdélyi-Kober fractional operator

This paper initiates the investigation of nonlinear integral equations with Erdélyi-Kober fractional operator. Existence and uniqueness results of solutions in a closed ball are obtained by using a directly computational method and Schauder fixed point theorem via a weakly singular integral inequality due to Ma and Pec?ari? [20]. Meanwhile, three certain solutions sets Y_K,σ, Y_1,λ and Y_1,1, which tending to zero at an appropriate rate t^−ν, 0 < ν = σ (or λ or 1) as t → +∞, are constructed and local stability results of solutions are obtained based on these sets respectively under some suitable conditions. Two examples are given to illustrate the results.     

Eigenvalue intervals for higher order problems with singular nonlinearities

We study eigenvalue intervals for higher order problems with multipoint boundary conditions. We consider the case when the nonlinearity may be singular at t = 0 or t = 1. Our approach is based on topological methods and cover both sublinear and superlinear cases. We also study the continuous dependence of solutions on functional parameters.     

Generalized Jordan algebras

We study commutative algebras which are generalizations of Jordan algebras. The associator is defined as usual by (x, y, z) = (x y)z − x(y z). The Jordan identity is (x², y, x) = 0. In the three generalizations given below, t, β, and γare scalars. ((x x)y)x + t((x x)x)y = 0, ((x x)x)(y x) − (((x x)x)y)x = 0, β((x x)y)x + γ((x x)x)y − (β + γ)((y x)x)x = 0. We show that with the exception of a few values of the parameters, the first implies both the second and the third. The first is equivalent to the combination of ((x x)x)x = 0 and the third. We give examples to show that our results are in some reasonable sense, the best possible.     

Well-posedness of a class of perturbed optimization problems in Banach spaces

Let X be a Banach space and Z a nonempty subset of X. Let J:Z→R be a lower semicontinuous function bounded from below and p?1. This paper is concerned with the perturbed optimization problem of finding z₀∈Z such that ‖x−z₀p^‖+J(z₀)=infz∈Z{‖x−zp^‖+J(z)}, which is denoted by minJ(x,Z). The notions of the J-strictly convex with respect to Z and of the Kadec with respect to Z are introduced and used in the present paper. It is proved that if X is a Kadec Banach space with respect to Z and Z is a closed relatively boundedly weakly compact subset, then the set of all x∈X for which every minimizing sequence of the problem minJ(x,Z) has a converging subsequence is a dense Gδ-subset of X?Z₀, where Z₀ is the set of all points z∈Z such that z is a solution of the problem minJ(z,Z). If additionally p>1 and X is J-strictly convex with respect to Z, then the set of all x∈X for which the problem minJ(x,Z) is well-posed is a dense Gδ-subset of X?Z₀.     

Eigenvalues and the One-Dimensional p-Laplacian

We consider the boundary value problem (?_p(u′))′ + λF(t, u) = 0, with p > 1, t ∈ (0, 1), u(0) = u(1) = 0, and with λ > 0. The value of λ is chosen so that the boundary value problem has a positive solution. In addition, we derive an explicit interval for λ such that, for any λ in this interval, the existence of a positive solution to the boundary value problem is guaranteed. In addition, the existence of two positive solutions for λ in an appropriate interval is also discussed.     

Landau’s theorem for functions with logharmonic Laplacian

In this paper, we show the existence of Landau constant for functions with logharmonic Laplacian of the form F(z) = ∣z∣²L(z) + K(z), ∣z∣ < 1, where L is logharmonic and K is harmonic. Moreover, the problem of minimizing the area is solved