Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Efficient Numerical Solution of the Modified Mild-Slope Equation

An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.

Numerical Solution of Constrained Mechanical System Motions Equations and Inverse Problems of Dynamics

In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.

ON THE NUMERICAL SOLUTION OF QUASILINEAR WAVE EQUATION WITH STRONG DISSIPATIVE TERM

The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.

NUMERICAL SOLUTIONS OF DISCONTINUOUS BOUNDARY VALUE PROBLEMS FOR GENERAL ELLIPTIC COMPLEX EQUATIONS OF FIRST ORDER

In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.

On numerical stationary distribution of overdamped Langevin equation in harmonic system

Efficient numerical algorithm for stochastic differential equation has been an important object in the research of statistical physics and mathematics for a long time.In this work we study the highly accurate numerical algorithm for the overdamped Langevin equation.In particular,our interest is in the behaviour of the numerical schemes for solving the overdamped Langevin equation in the harmonic system.Based on the large friction limit of the underdamped Langevin dynamic scheme,three algorithms for overdamped Langevin equation are obtained.We derive the explicit expression of the stationary distribution of each algorithm by analysing the discrete time trajectory for both one-dimensional case and multi-dimensional case.The accuracy of the stationary distribution of each algorithm is illustrated by comparing with the exact Boltzmann distribution.Our results demonstrate that the“BAOA-limit”algorithm generates an accurate distribution of the harmonic system in a canonical ensemble,within a stable range of time interval.The other algorithms do not produce the exact distribution of the harmonic system.

The θ－Methods in Numerical Solution of Systems of Differential Equations with Two Delay Terms

This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.

New Soliton Wave Solutions to a Nonlinear Equation Arising in Plasma Physics

The extraction of traveling wave solutions for nonlinear evolution equations is a challenge in various mathematics,physics,and engineering disciplines.This article intends to analyze several traveling wave solutions for themodified regularized long-wave(MRLW)equation using several approaches,namely,the generalized algebraic method,the Jacobian elliptic functions technique,and the improved Q-expansion strategy.We successfully obtain analytical solutions consisting of rational,trigonometric,and hyperbolic structures.The adaptive moving mesh technique is applied to approximate the numerical solution of the proposed equation.The adaptive moving mesh method evenly distributes the points on the high error areas.This method perfectly and strongly reduces the error.We compare the constructed exact and numerical results to ensure the reliability and validity of the methods used.To better understand the considered equation’s physical meaning,we present some 2D and 3D figures.The exact and numerical approaches are efficient,powerful,and versatile for establishing novel bright,dark,bell-kink-type,and periodic traveling wave solutions for nonlinear PDEs.

Existence and Numerical Solution for a Coupled System of Multi-term Fractional Differential Equations

An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.

The Numerical Solutions of Systems of Nonlinear Integral Equations with the Spline Functions

The main goal of this work is to develop an effective technique for solving nonlinear systems of Volterra integral equations. The main tools are the cardinal spline functions on small compact supports. We solve a system of algebra equations to approximate the solution of the system of integral equations. Since the matrix for the algebraic system is nearly triangular, It is relatively painless to solve for the unknowns and an approximation of the original solution with high precision is accomplished. In order to enhance the accuracy, several cardinal splines are employed in the paper. Our schemes were compared with other techniques proposed in recent papers and the advantage of our method was exhibited with several numerical examples.

Existence and Uniqueness of 2π-Periodic Solution About Duffing Equation and Its Numerical Method

A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large convergence domain is given. Finally, a numerical example is given.

A ONE-STEP EXPLICIT FORMULA FOR THE NUMERICAL SOLUTION OF STIFF ORDINARY DIFFERENTIAL EQUATION

In this paper, a new one-step explicit method of fourth order is derived. The new method is proved to be A-stable and L-stable, and it gives exact results when applied to the test equation y’=λy with Re(λ)【0, Also several numerical examples are included.

ON FICTITIOUS DOMAIN METHOD FOR THE NUMERICAL SOLUTION TO HEAT CONDUCTION EQUATION WITH DERIVATIVE BOUNDARY CONDITIONS

This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.

On L ∞ Stability and Convergence of Fictitious Domain Method for the Numerical Solution to Parabolic Differential Equation with Derivative Boundary Conditions

This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.

Approximate Formulation and Numerical Solution for Hypersingular Boundary Integral Equations in Plane Elasticity

Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.

Numerical Solution for Fractional Partial Differential Equation with Bernstein Polynomials

A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational matrix of Bernstein polynomials is derived. With the operational matrix, the equation is transformed into the products of several dependent matrixes which can also be regarded as the system of linear equations after dispersing the variable. By solving the linear equations, the numerical solutions are acquired. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples are provided to show that the method is computationally efficient.

Lie group analysis, numerical and non-traveling wave solutions for the (2+1)-dimensional diffusion–advection equation with variable coefficients

In this paper, the variable-coefficient diffusion-advection （DA） equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended （Gl/G）-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.

Numerical Complexiton Solutions of Complex KdV Equation

In this paper, we directly extend the applications of the Adomian decomposition method to investigate the complex KdV equation. By choosing different forms of wave functions as the initial values, three new types of realistic numerical solutions： numerical positon, negaton solution, and particularly the numerical analytical complexiton solution are obtained, which can rapidly converge to the exact ones obtained by Lou et al. Numerical simulation figures are used to illustrate the efficiency and accuracy of the proposed method.

Numerical Solutions of a New Type of Fractional Coupled Nonlinear Equations

In this paper, we investigate a new type of fractional coupled nonlinear equations. By introducing the fractional derivative that satisfies the Caputo＇s definition, we directly extend the applications of the Adomian decomposition method to the new system. As a result, with the aid of Maple, the realistic and convergent rapidly series solutions are obtained with easily computable components. Two famous fractional coupled examples： KdV and mKdV equations, are used to illustrate the efficiency and accuracy of the proposed method.

Numerical Solutions of a Class of Nonlinear Evolution Equations with Nonlinear Term of Any Order

In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt＋auxx ＋ bu ＋ cu^p＋ du^2p-1=0, which contains some important famous equations. When setting the initial conditions in different forms, some new generalized numerical solutions： numerical hyperbolic solutions, numerical doubly periodic solutions are obtained. The numerical solutions are compared with exact solutions. The scheme is tested by choosing different values of p, positive and negative, integer and fraction, to illustrate the efficiency of the ADM method and the generalization of the solutions.