An Extended Numerical Method by Stancu Polynomials for Solution of Integro-Differential Equations Arising in Oscillating Magnetic Fields

In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.

Linear Functional Equations and Twisted Polynomials

A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view a number of algorithms.

Nonlinear Algebraic Equations Solved by an Optimal Splitting-Linearizing Iterative Method

How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).

A Collocation Technique via Pell-Lucas Polynomials to Solve Fractional Differential EquationModel for HIV/AIDS with Treatment Compartment

In this study,a numerical method based on the Pell-Lucas polynomials(PLPs)is developed to solve the fractional order HIV/AIDS epidemic model with a treatment compartment.The HIV/AIDS mathematical model with a treatment compartment is divided into five classes,namely,susceptible patients(S),HIV-positive individuals(I),individuals with full-blown AIDS but not receiving ARV treatment(A),individuals being treated(T),and individuals who have changed their sexual habits sufficiently(R).According to the method,by utilizing the PLPs and the collocation points,we convert the fractional order HIV/AIDS epidemic model with a treatment compartment into a nonlinear system of the algebraic equations.Also,the error analysis is presented for the Pell-Lucas approximation method.The aim of this study is to observe the behavior of five populations after 200 days when drug treatment is applied to HIV-infectious and full-blown AIDS people.To demonstrate the usefulness of this method,the applications are made on the numerical example with the help of MATLAB.In addition,four cases of the fractional order derivative(p=1,p=0.95,p=0.9,p=0.85)are examined in the range[0,200].Owing to applications,we figured out that the outcomes have quite decent errors.Also,we understand that the errors decrease when the value of N increases.The figures in this study are created in MATLAB.The outcomes indicate that the presented method is reasonably sufficient and correct.

THE GROWTH OF SOLUTIONS TO HIGHER ORDER DIFFERENTIAL EQUATIONS WITH EXPONENTIAL POLYNOMIALS AS ITS COEFFICIENTS

By looking at the situation when the coefficients Pj(z)(j=1,2,…,n-1)(or most of them) are exponential polynomials,we investigate the fact that all nontrivial solutions to higher order differential equations f((n))+Pn-1(z)f((n-1))+…+P0(z)f=0 are of infinite order.An exponential polynomial coefficient plays a key role in these results.

Operator Methods and SU(1,1) Symmetry in the Theory of Jacobi and of Ultraspherical Polynomials

Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting properties. It leads to an alternative definition of the Ultraspherical polynomials by a fixed integral operator in application to powers of the variable u in an analogous way as it is possible for Hermite polynomials. From this follows a generating function which is apparently known only for the Legendre and Chebyshev polynomials as their special case. Furthermore, we show that the Ultraspherical polynomials form a realization of the SU(1,1) Lie algebra with lowering and raising operators which we explicitly determine. By reordering of multiplication and differentiation operators we derive new operator identities for the whole set of Jacobi polynomials which may be applied to arbitrary functions and provide then function identities. In this way we derive a new “convolution identity” for Jacobi polynomials and compare it with a known convolution identity of different structure for Gegenbauer polynomials. In short form we establish the connection of Jacobi polynomials and their related orthonormalized functions to the eigensolution of the Schrödinger equation to Pöschl-Teller potentials.

HERMITE MATRIX POLYNOMIALS AND SECOND ORDER MATRIX DIFFERENTIAL EQUATIONS

In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.

New Implementation of Legendre Polynomials for Solving Partial Differential Equations

In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.

Chebyshev Polynomials for Solving a Class of Singular Integral Equations

This paper is devoted to studying the approximate solution of singular integral equations by means of Chebyshev polynomials. Some examples are presented to illustrate the method.

Numerical Solutions of Volterra Equations Using Galerkin Method with Certain Orthogonal Polynomials

This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations of the first and second kind respectively using orthogonal polynomials as trial functions which are constructed in the interval [-1,1] with respect to the weight function w(x)=1+x2. The efficiency of the proposed method is tested on several numerical examples and compared with the analytic solutions available in the literature.

The Degenerate Form of the Adomian Polynomials in the Power Series Method for Nonlinear Ordinary Differential Equations

In this paper, we propose a new variation of the Adomian polynomials, which we call the degenerate Adomian polynomials, for the power series solutions of nonlinear ordinary differential equations with nonseparable nonlinearities. We establish efficient algorithms for the degenerate Adomian polynomials. Next we compare the results by the Adomian decomposition method using the classic Adomian polynomials with the results by the Rach-Adomian-Meyers modified decomposition method incorporating the degenerate Adomian polynomials, which itself has been shown to be a confluence of the Adomian decomposition method and the power series method. Convergence acceleration techniques including the diagonal Pade approximants are considered, and new numeric algorithms for the multistage decomposition are deduced using the degenerate Adomian polynomials. Our new technique provides a significant advantage for automated calculations when computing the power series form of the solution for nonlinear ordinary differential equations. Several expository examples are investigated to demonstrate its reliability and efficiency.

SOlvaBILITY OF HIGHER INDEX TIME-VARYING LINEAR DIFFERENTIAL-ALGEBRAIC EQUATIONS

Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.

Reconstructed Elzaki Transform Method for Delay Differential Equations with Mamadu-Njoseh Polynomials

One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transform method (ETM)) for computing delay differential equations (DDEs). Here, a reconstructed Elzaki transform method (RETM) is proposed for the solution of DDEs where Mamadu-Njoseh polynomials are applied as basis functions in the approximation of the analytic solution. Using this strategy, a numerical illustration as in Ref.[1] is provided to the RETM as a basis for comparison to guarantee accuracy and consistency of the method. All numerical computations were performed with MAPLE 18 software.

A NEW ALGORITHM FOR SOLVING DIFFERENTIAL/ALGEBRAIC EQUATIONS OF MULTIBODY SYSTEM DYNAMICS

The second order Euler-Lagrange equations are transformed to a set of first order differential/algebraic equations, which are then transformed to state equations by using local parameterization. The corresponding discretization method is presented, and the results can be used to implementation of various numerical integration methods. A numerical example is presented finally.

Efficient Numerical Methods for Solving Differential Algebraic Equations

This research aims to solve Differential Algebraic Equation (DAE) problems in their original form, wherein both the differential and algebraic equations remain. The Newton or Newton-Broyden technique along with some integrators such as the Runge-Kutta method is coupled together to solve the problems. Experiments show that the method developed in this paper is efficient, as it demonstrates that implementation of the method is not difficult, and such method is able to provide approximate solutions with ease within some desired accuracy standards.

Application of Bernstein polynomials for solving Fredholm integro-differential-difference equations

In this paper,the Bernstein polynomials method is proposed for the numerical solution of Fredholm integro-differential-difference equation with variable coefficients and mixed conditions.This method is using a simple computational manner to obtain a quite acceptable approximate solution.The main characteristic behind this method lies in the fact that,on the one hand,the problem will be reduced to a system of algebraic equations.On the other hand,the efficiency and accuracy of the Bernstein polynomials method for solving these equations are high.The existence and uniqueness of the solution have been proved.Moreover,an estimation of the error bound for this method will be shown by preparing some theorems.Finally,some numerical experiments are presented to show the excellent behavior and high accuracy of this algorithm in comparison with some other well-known methods.

ON A REGULARIZATION OF INDEX 2 DIFFERENTIAL-ALGEBRAIC EQUATIONS WITH PROPERLY STATED LEADING TERM

In this article, linear regular index 2 DAEs A（t）[D（t）x（t）]＇ ＋ B（t）x（t） = q（t） are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.

NONNEGATIVITY OF SOLUTIONS OF NONLINEAR FRACTIONAL DIFFERENTIAL-ALGEBRAIC EQUATIONS

Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.

Finite Element Orthogonal Collocation Approach for Time Fractional Telegraph Equation with Mamadu-Njoseh Polynomials

Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.

On the Numerical Solution of Singular Integral Equation with Degenerate Kernel Using Laguerre Polynomials

In this paper, we derive a simple and efficient matrix formulation using Laguerre polynomials to solve the singular integral equation with degenerate kernel. This method is based on replacement of the unknown function by truncated series of well known Laguerre expansion of functions. This leads to a system of algebraic equations with Laguerre coefficients. Thus, by solving the matrix equation, the coefficients are obtained. Some numerical examples are included to demonstrate the validity and applicability of the proposed method.