New numerical method for ordinary differential equations: Newton polynomial

Abstract

An error estimate of optimal order is established for the correspondingnumerical solutions in a scaled residual norm. In addition, a mathematical convergenceis established in a weak L2topology for the new numerical method. Numerical resultsare reported to demonstrate the efficiency of the primal–dual weak Galerkin method aswell as the accuracy of the numerical approximations