We consider the fast solution of large, piecewise smooth minimization problems as typically arising from the nite element discretization of porous media ow. For lack of smoothness, usual Newton multigrid methods cannot be applied. We propose a new approach based on a combination of convex minization with constrained Newton linearization. No regularization is involved. We show global convergence...

The collocation method based on cubic B-spline, is developed to approximate the solution of second kind nonlinear Fredholm integral equations. First of all, we collocate the solution by B-spline collocation method then the Newton-Cotes formula use to approximate the integrand. Convergence analysis has been investigated and proved that the quadrature rule is third order convergent. The presented...

Restarting GMRES, a linear solver frequently used in numerical schemes, is known to suffer from stagnation. In this paper, a simple strategy is proposed to detect and avoid stagnation, without modifying the standard GMRES code. Numerical tests with the modified GMRES(m), GMRESH(m) procedure, alone and as part of an inexact Newton procedure with several choices for the forcing term, demonstrate ...

This paper investigates the global chaos synchronization of identical hyperchaotic Newton-Leipnik systems (Ghosh and Bhattacharya, 2010) by sliding mode control. The stability results derived in this paper for the complete synchronization of identical hyperchaotic Newton-Leipnik chaotic systems are established using Lyapunov stability theory. Since the Lyapunov exponents are not required for th...

In the Laurent expansion ζ(s, a) = 1 s− 1 + ∞ ∑ k=0 (−1)γk(a) k! (s− 1) , 0 < a ≤ 1, of the Riemann-Hurwitz zeta function, the coefficients γk(a) are known as Stieltjes, or generalized Euler, constants. [When a = 1, ζ(s, 1) = ζ(s) (the Riemann zeta function), and γk(1) = γk.] We present a new approach to high-precision approximation of γk(a). Plots of our results reveal much structure in the gr...

A function and its first two derivatives are estimated by convolutions with well-chosen non-differentiable kernels. The convolutions are in turn approximated by Newton–Cotes integration techniques with the aid of a polynomial interpolation based on an arbitrary finite set of points. Precise numerical results are obtained with far fewer points than that in classic SPH, and error bounds are deriv...

We study the kernels of the remainder term Rn,s(f) of GaussTurán quadrature formulas ∫ 1 −1 f(t)w(t) dt = n ∑

We study the dynamical behavior of the discontinuous Galerkin finite element method for initial value problems in ordinary differential equations. We make two different assumptions which guarantee that the continuous problem defines a dissipative dynamical system. We show that, under certain conditions, the discontinuous Galerkin approximation also defines a dissipative dynamical system and we ...

in this paper, a novel hybrid method based on learning algorithmof fuzzy neural network and newton-cotesmethods with positive coefficient for the solution of linear fredholm integro-differential equation of the second kindwith fuzzy initial value is presented. here neural network isconsidered as a part of large field called neural computing orsoft computing. we propose alearning algorithm from ...