An eighth order frozen Jacobian iterative method for solving nonlinear IVPs and BVPs
|
|
- Janel Brooks
- 5 years ago
- Views:
Transcription
1 Available online at J. Math. Computer Sci., 7 7, Research Article Journal Homepage: An eighth order frozen Jacobian iterative method for solving nonlinear IVPs and BVPs Dina Abdullah Alrehaili a,, Dalal Adnan Al-Maturi a, Eman Salem Al-Aidarous a, Fayyaz Ahmad b,c,d a Department of Mathematics, Faculty of Science, King Abdulaziz University, P. O. Box 83, Jeddah 589, Saudi Arabia. b Dipartimento di Scienza e Alta Tecnologia, Universita dell Insubria, Via Valleggio, Como, Italy. c Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Eduard Maristany, Barcelona 89, Spain. d UCERD Islamabad, Pakistan. Abstract A frozen Jacobian iterative method is proposed for solving systems of nonlinear equations. In particular, we are interested in solving the systems of nonlinear equations associated with initial value problems IVPs and boundary value problems BVPs. In a single instance of the proposed iterative method DEDF, we evaluate two Jacobians, one inversion of the Jacobian and four function evaluations. The direct inversion of the Jacobian is computationally expensive, so, for a moderate size, LU factorization is a good direct method to solve the linear system. We employed the LU factorization of the Jacobian to avoid the direct inversion. The convergence order of the proposed iterative method is at least eight, and it is nine for some particular classes of problems. The discretization of IVPs and BVPs is employed by using Jacobi-Gauss-Lobatto collocation J-GL-C method. A comparison of J-GL-C methods is presented in order to choose best collocation method. The validity, accuracy and the efficiency of our DEDF are shown by solving eleven IVPs and BVPs problems. c 7 All rights reserved. Keywords: Frozen Jacobian iterative methods, systems of nonlinear equations, nonlinear initial-boundary value problems, Jacobi-Gauss-Lobatto quadrature, collocation method. MSC: 65H, 65N.. Introduction The closed form solutions of nonlinear problems are not always available. It means we need an iterative method to solve them. We are interested in iterative numerical methods for solving systems of nonlinear equations associated with nonlinear IVPs and BVPs. At the first place, we find some efficient and accurate discretization method to approximate the nonlinear differential equations. Once, we translated the continuous problem into a discrete nonlinear problem. We proceed by finding some efficient iterative solvers for the associated systems of nonlinear equations. The pseudospectral collocation methods offer excellent accuracy. Doha et al. [7] employed J-GL-C method for the discretization of nonlinear + Schrödinger for spatial dimensions and got a system of Corresponding author addresses: dalrehaili@kau.edu.sa Dina Abdullah Alrehaili, dmaturi@kau.edu.sa Dalal Adnan Al-Maturi, ealaidarous@kau.edu.sa Eman Salem Al-Aidarous, fahmad@uninsubria.it,fayyaz.ahmad@upc.edu Fayyaz Ahmad doi:.436/jmcs Received 7-4-3
2 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, ordinary differential equations. They used implicit Runge-Kutta method to solve the system of ordinary differential equations and obtained highly accurate numerical solutions, for four different kinds of nonlinear + Schrödinger equations. In another article Bhrawy et al. [5] solved the nonlinear reaction-diffusion equations by using J-GL-C method. In the solution of the complex generalized Zakharov system of equation [4], the J-GL-C method is also the method of discretization. The further application of pseudospectral collocation techniques for solving nonlinear IVPs and BVPs can be found in [6, 8, 9, 7] and references therein. The J-GL-C method is a parametric pseudospectral collocation method. By selecting different values of the parameters, we can get different pseudospectral collocation methods. The Legendre, Chebyshev, and Gegenbauer pseudospectral collocation methods are special cases of J-GL-C method [6]. The Jacobi polynomials are the eigenfunctions of a singular Strum-Liouville problem [6] x σ x + θ φ + θ + φ + xσ x + nn + θ + φ + σx =. The following recurrence relation produces the Jacobi polynomials J θ,φ k+ x = a θ,φ k b θ,φ k J θ,φ k x c θ,φ k J θ,φ k J θ,φ x =, J θ,φ x = θ + φ + x + φ θ, x, k, where a θ,φ k + θ + φ + k + θ + φ + k =, k + k + θ + φ + b θ,φ k = c θ,φ k = θ φ k + θ + φ + k + k + θ + φ k + θ + φ, k + θk + φk + θ + φ + k + k + θ + φ k + θ + φ. The p-th derivative of k-degree Jacobi polynomial J θ,φ k can be computed as D p J θ,φ Γj + θ + φ + p + k x = p Γj + θ + φ + Jθ+p,φ+p k p x. The Jacobi polynomials are orthogonal over the domain [, ] with respect to weight function + x θ x φ. The J-GL-C method is attractive because the numerical differentiation matrix is easy to construct for approximating differential operators of different orders. The construction of such differentiation matrices can be found in [4]. Assume Q is the Jacobi-Gauss-Lobatto differentiation matrix of the first order derivative operator over the domain [, ].Then a derivative of order p can be approximated over the interval [a, b] as follows d p p dx p b a Q. The discretization of IVPs and BVPs gives a system of nonlinear algebraic equations Fy =, where y = [y, y, y 3,, y n ] T. The classical iterative method for solving system of nonlinear equations is the Newton method [, 8], which can be written as NR = { y = initial guess, y n+ = y n F y n Fy n, where det F y n. Many researchers [ 3,,, 3, 5, 9, ] have proposed the frozen Jacobian multi-step iterative method for solving the system of nonlinear equations. The frozen Jacobian iterative
3 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, methods are computationally efficient, because inversion of the Jacobian is too expensive. For moderate size system of nonlinear equations, it is good idea to use the LU factorization for solving the system of linear equations. The benefit of using the LU factorization of frozen Jacobian is apparent because in DEDF method we use eight systems of linear equations with fixed augmented matrix, which is a frozen Jacobian. In fact, we solve eight lower, and eight upper triangular systems and the solution of the triangular system is computationally economical. In the DEDF method, we get an increment in the convergence order of one per solving a system of linear equations. It is a good idea to design higher order frozen Jacobian iterative methods, which provide us rapid convergence in the vicinity of the root.. Frozen Jacobian iterative method A new iterative method DEDF can be described as y = initial guess, F y φ = Fy, y = y φ, Convergence order 8 F y φ = Fy, Function evaluations = 4 y = y φ, Jacobian evaluations = F y φ 3 = Fy, LU-factorization = Matrix-vector y 3 = y α φ 3, DEDF = y 3 = y α φ 3, multiplications = 4 F y φ 4 = Fy 3, Vector-vector F y φ 5 = F y 3 φ 4, multiplications = 8 F y φ 6 = F y 3 φ 5, Number of lower and F y φ 7 = F y 3 φ 6, upper triangular systems = 8 F y φ 8 = F y 3 φ 7, y 4 = y β φ 3 β φ 4 β 3 φ 5 β 4 φ 6 β 5 φ 7 β 6 φ 8, y = y 4, where α = , α = 4α 3, 6 α 5 β = 4 α 3 α 3, β = /3 96 α 3 56 α + 6 α 659 α 3 7 α + 8 α 3 4 α 3, β 3 = /8 6 α 35 α + 46 α 3 7 α + 8 α 3, β 4 = 3/6 α 6 α + 7 α 3 7 α + 8 α 3, β 5 = /8 96 α 79 α + 84 α 3 7 α + 8 α 3, β 6 = /3 8 α 48 α + 69 α 3 7 α + 8 α 3.
4 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, The DEDF is an efficient iterative method, since it requires only one inversion of the Jacobian regarding LU factorization. In summary, the DEDF employs four functions, two Jacobian evaluations, and eight matrix-vector multiplications. 3. Convergence analysis We used the symbolic algebra of Maple software to deal with symbolic computations. The Fréchet differentiability condition on the system of nonlinear equations is the essential because it ensure the linearization of the nonlinear function. The existence of the following limits Fy + h Fy Ah lim =, h h ensures the Fréchet differentiability. The linear operator A is called the first order Fréchet derivative and we denote it by F y. The higher order Fréchet derivatives can be computed recursively where u is a vector independent from y. F y = Jacobian Fy, F j yu j = Jacobian F j yv j, j, Theorem 3.. Let F : Γ R n R n be a sufficiently Fréchet differentiable function on an open convex neighborhood Γ of y R n with Fy = and detf y, where F y denotes the Fréchet derivative of Fy. Let A = F y and A j = j! F y F j y, for j, where F j y denotes j-order Fréchet derivative of Fy. Then, with an initial guess in the neighborhood of y, the sequence {y m } generated by DEDF converges to y with local order of convergence at least eight and error e 4 = Le 8 + O e 9, p times {}}{ where e = y y, e p = e, e,..., e, and L = 6A 3 A 3A α4 + 48A A 3 A 4 α4 + 6A A 3A 3 α4 48A 3A 5 α4 86A7 α3 + 56A3 A 3A α3 6α 3 A A 3 A 4 56α3 A A 3A 3 + 6α3 A 3A 5 + 8α A7 7α A3 A 3A + 36α A A 3 A 4 + 7α A A 3A 3 36α A 3A 5 44α A 7 + 4α A 3 A 3A 64α A A 3 A 4 4α A A 3A α A 3 A A7 8A3 A 3A + 7A A 3 A 4 + 8A A 3A 3 7A 3A 5 / α 3α α +, 8 times {}}{ is a 8-linear function, i.e., L L R n, R n, R n,, R n with Le 8 R n. Proof. We define the error at the n-th step e n = y n y. To complete the convergence proof, we performed the detailed computations by using Maple and details are provided below in sequence. F y = I A e + 3A 3 + 4A e + 4A 4 + 6A 3 A + 6A A 3 8A 3 e 3 + 5A 5 A A 3 A 3 A A A 3 A + 8A 4 A + 9A 3 + 8A A 4 + 6A 4 e 4 + 6A 6 6A A 4 3A 5 6A 4 A 8A 3 A 6A A 4 A 8A 3 A A 3 8A A 3 + 4A 3A 3 + 4A A 3 A + 4A A 3A + A 5 A + A 4 A 3 + A 3 A 4 + A A 5 + 4A 3 A 3 e O e 9 A, Fy =A e + A e + A 3e 3 + A 4e 4 + A 5e 5 + A 6e 6 + A 7e 7 + A 8e 8 e + O 9,
5 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, φ =e A e + A 3 + A e 3 + 3A 4 + 4A A 3 + 3A 3 A 4A 3 e 4 + 4A 5 + 6A A 4 + 4A 4 A 6A A 3 A + 6A 3 + 8A4 6A 3A 8A A 3 e 5 + 5A 6 A 3 A A 3 8A A 4 A + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A + A A 3A + A A 3 A + A 3A 3 9A 3 A A A 3 6A5 A A 4 8A 4 A + 6A3 A 3 e O e 9, e =A e A + 3 A e 3 + 3A 4 4A A 3 3A 3 A + 4A 3 e 4 + 4A 5 6A A 4 4A 4 A + 6A A 3 A 6A 3 8A4 + 6A 3A + 8A A 3 e 5 + 5A 6 + A 3 A A 3 + 8A A 4 A 8A A 5 9A 3 A 4 8A 4 A 3 5A 5 A A A 3A A A 3 A A 3A 3 + 9A 3 A + A A 3 + 6A5 + A A 4 + 8A 4 A 6A 3 A 3 e O e 9, Fy =A A e A + 3 A e 3 + 3A 4 4A A 3 3A 3 A + 5A 3 e 4 + 4A 5 6A A 4 6A 3 4A 4A A 4 + 8A A 3 A + 6A 3 A + A A 3 e 5 + 5A 6 + 8A 4 A + A 3A A 3 4A 3 A 3 9A A 3A 9A A 3 A A 3A 3 + 9A 3 A + A A 4 A + 6A A 3 + 8A5 8A A 5 9A 3 A 4 8A 4 A 3 5A 5 A + 5A A 4 e O e 9, φ =A e + 4A + A 3 e 3 3A A 3 8A A 3 6A 3 A e 4 + 4A 5 A A 4 A 3 8A 4A 38A 4 + A A 3 A + 8A 3 A + 6A A 3 e 5 + 5A 6 + 4A 4 A + 36A 3A A 3 76A 3 A 3 59A A 3A 55A A 3 A 5A 3A 3 + 7A 3 A + 7A A 4 A + 4A A 3 + 4A5 6A A 5 8A 3 A 4 6A 4 A 3 A 5 A + 39A A 4 e O e 9, e =e 3 A + 9A 3 + 4A A 3 + 3A 3 A e 4 3A A A 4 + 6A 3 + 4A 4A 4A A 3 A A 3 A 8A A 3 e A A 3A + 43A A 3 A + 38A 3A 3 6A 4A 8A 3 A 9A A 4 A 4A 3 A A 3 8A A 3 + 6A3 A 3 88A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 7A A 4 e O e 9, Fy =A A e3 + 9A 3 + 4A A 3 + 3A 3 A e 4 3A A A 4 + 6A 3 + 4A 4A 4A A 3 A A 3 A 8A A 3 e A A 3A + 43A A 3 A + 38A 3A 3 6A 4A 8A 3 A 9A A 4 A 4A 3 A A 3 8A A 3 + 6A3 A 3 84A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 7A A 4 e O e 9, φ 3 =A e3 + 3A 3 + 4A A 3 + 3A 3 A e 4 6A + A 4 + 4A 4 A A A 3 A + 56A 4 + 6A 3 8A 3A 6A A 3 e 5 + 7A A 4 A 36A 3 A A 3 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A + 87A A 3A + 79A A 3 A + 77A 3A 3 4A 4A 7A 3 A 4A A 3 + A3 A 3 96A 5 39A A 4 e O e 9, e 3 =4A 3 e4 + 6A 4 + 6A A 3 A + 6A 3 A + 8A A 3 e 5 + 4A A 3A 36A A 3 A 39A 3A 3 + 8A 4 A + 9A 3 A + 8A A 4 A + A 3 A A 3 + A A 3 5A3 A 3 + 8A 5 + A A 4 e 6 + 8A A 3 7A A 3 A A 3 78A 3 A A 3 + 6A 4 A A 3 + 8A A A 4 A 3 + 8A A 3 A 4 + 8A 3 A A 4 + 6A 4 A 3 78A 3 A 4 356A 6 + A 5A + A 4A 3 A + A 3 A 4 A + A A 5 A + 68A 3 A 3A + 5A A 3A + 48A A 3 A A 3A 4 + 6A A 5 5A 4 A 3 54A 3 A 48A A 4 A 54A A 3 A 6A 3 A A 3 A 54A A 4A e 7 e + + O 9,
6 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Fy 3 =A 4A 3 e4 + 6A 4 + 6A A 3 A + 6A 3 A + 8A A 3 e 5 + 4A A 3A 36A A 3 A 39A 3A 3 + 8A 4 A + 9A 3 A + 8A A 4 A + A 3 A A 3 + A A 3 5A3 A 3 + 8A 5 + A A 4 e 6 + 8A A 3 7A A 3 A A 3 78A 3 A A 3 + 6A 4 A A 3 + 8A A A 4 A 3 + 8A A 3 A 4 + 8A 3 A A 4 + 6A 4 A 3 78A 3 A 4 356A 6 + A 5A + A 4A 3 A + A 3 A 4 A + A A 5 A + 68A 3 A 3A + 5A A 3A + 48A A 3 A A 3A 4 + 6A A 5 5A 4 A 3 54A 3 A 48A A 4 A 54A A 3 A 6A 3 A A 3 A 54A A 4A e 7 e + + O 9, e 3 = A α A e 3 3α + A 3 4α A A 3 3α A 3 A 9A 3 + 4A A 3 + 3A 3 A e 4 + α A A 3 A + 8α A 3 A 6α A A 4 6α A 3 4α A 4 A 56α A 4 + 6α A A 3 4A A 3 A A 3 A + 6A A 4 + 6A 3 + 4A 4A + 3A 4 8A A 3 e A A 3A + 43A A 3 A + 38A 3A 3 6A 4A 8A 3 A 9A A 4 A 8A A 3 4A 3A A 3 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 39α A A 4 + 6A 3 A 3 88A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 87α A A 3A 79α A A 3 A 77α A 3 A 3 + 4α A 4 A + 7α A 3 A + 7α A A 4 A + 4α A A α A 3 A A 3 α A 3 A α A 5 8α A A 5 7A A 4 e O e 9, Fy 3 =A A α A e 3 3α + A 3 4α A A 3 3α A 3 A 9A 3 + 4A A 3 + 3A 3 A e 4 + α A A 3 A + 8α A 3 A 6α A A 4 6α A 3 4α A 4 A 56α A 4 + 6α A A 3 4A A 3 A A 3 A + 6A A 4 + 6A 3 + 4A 4A + 3A 4 8A A 3 e A A 3A + 43A A 3 A + 38A 3A 3 6A 4A 8A 3 A 9A A 4 A 8A A 3 4A 3A A 3 87α A A 3A 79α A A 3 A 77α A 3 A 3 + 4α A 4 A + 7α A 3 A + 7α A A 4 A + 4α A A α A 3 A A 3 α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 39α A A 4 + 6A 3 A 3 84A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 7A A 4 + 4α A5 e O e 9, F y 3 =A I + 4α A 3 + 4A3 e 3 + 8α A A 3 + 6α A 4 6α A A 3 A + 8A A 3 8A 4 + 6A A 3 A e 4 + α A A 4 + 5α A 3 A 3 α A 5 α A A α A A 3 A + 4α A A 3A 8α A A 4 A + A A 4 36A 3 A 3 + 6A 5 + A A 3 4A A 3 A 8A A 3A + 8A A 4 A e A A 3 48A A 3 A A 3 + 6A A 4 A 3 + 8A A 3 A 4 + α A 3A 4 36A A 3 A 38A A 4A 3A A 4 A + A4 A A 3 A 3A + 86A A 3A + 76A A 3 A 3 + A 3 A 4 76A6 + 6A A 5 + A A 5 A 54A 3 A 4 + 8α A A 3 + 7α A A 3 A A 3 6α A A 4 A 3 8α A A 3 A α A A 3 A + 54α A A 4A + 48α A A 4 A 4α A 4 A 3 74α A 3 A 3A 58α A A 3A 54α A A 3 A 3 4α A 3 A α A 6 6α A A 5 α A A 5 A + 78α A 3 A 4 e O e 9, φ 4 = A α A e 3 + 4α A A 3 3α A 3 A + 7α A 3 + 4A A 3 + 3A 3 A 3A 3 e 4 + 6α A A 3 A 6α A A 4 4α A 4 A + 4α A 3 A 6α A 3 9α A α A A 3 A A 3 A + 6A A 4 + 4A 4 A 8A 3 A + 6A A4 6A A 3 e A A 3A + 79A A 3 A + 77A 3A 3 4A 4A 7A 3 A 7A A 4 A 4A A 3 36A 3A A 3 39α A A 3A 7α A A 3 A 8α A 3 A 3 + 3α A 4 A + 36α A 3 A + 35α A A 4 A + 5α A A α A 3 A A 3 8α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3
7 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, α A 5 A + 5α A A 4 + A 3 A 3 96A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 39A A 4 + 4α A5 e O e 9, φ 5 = A α A e 3 α + A 3 4α A A 3 3α A 3 A 7A 3 + 4A A 3 + 3A 3 A e 4 + 3α A A 3 A + 3α A 3 A 6α A A 4 6α A 3 4α A 4 A 3α A 4 + 4α A A 3 6A A 3 A 4A 3 A + 6A A 4 + 6A 3 + 4A 4A + 9A 4 34A A 3 e A A 3A + 7A A 3 A + 8A 3A 3 + 8α α A 5 3A 4A 36A 3 A 35A A 4 A 5A A 3 48A 3A A 3 3α A A 3A 87α A A 3 A 9α A 3 A 3 + 4α A 4 A + 45α A 3 A + 43α A A 4 A + 64α A A 3 + 6α A 3 A A 3 64α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 63α A A 4 + 8A 3 A 3 368A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 8α A 5 5A A 4 + 4α A5 e O e 9, φ 6 = A α A e 3 5α + A 3 4α A A 3 3α A 3 A A 3 + 4A A 3 + 3A 3 A e α A A 3 A + 36α A 3 A 6α A A 4 6α A 3 4α A 4 A 8α A 4 + 5α A A 3 3A A 3 A 3A 3 A + 6A A 4 + 6A 3 + 4A 4A + 3A 4 4A A 3 e 5 + 3A A 3A + 87A A 3 A + 9A 3A 3 + 6α α A 5 4A 4A 45A 3 A 43A A 4 A 64A A 3 6A 3A A 3 79α A A 3A 59α A A 3 A 66α A 3 A α A 4 A + 54α A 3 A + 5α A A 4 A + 76α A A 3 + 7α A 3 A A 3 364α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 75α A A A 3 A 3 64A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 6α A 5 63A A 4 + 4α A5 e O e 9, φ 7 = A α A e 3 9α + A 3 4α A A 3 3α A 3 A 5A 3 + 4A A 3 + 3A 3 A e α A A 3 A + 4α A 3 A 6α A A 4 6α A 3 4α A 4 A 4α A α A A 3 38A A 3 A 36A 3 A + 6A A 4 + 6A 3 + 4A 4A + 8A 4 5A A 3 e A A 3A + 59A A 3 A + 66A 3A 3 + 4α α A 5 48A 4A 54A 3 A 5A A 4 A 76A A 3 7A 3A A 3 367α A A 3A 343α A A 3 A 353α A 3 A α A 4 A + 63α A 3 A + 59α A A 4 A + 88α A A α A 3 A A 3 48α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 87α A A A 3 A 3 98A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 4α A 5 75A A 4 + 4α A5 e O e 9, φ 8 = A α A e 3 33α + A 3 4α A A 3 3α A 3 A 9A 3 + 4A A 3 + 3A 3 A e 4 + 5α A A 3 A + 48α A 3 A 6α A A 4 6α A 3 4α A 4 A 36α A α A A 3 44A A 3 A 4A 3 A + 6A A 4 + 6A 3 + 4A 4A + 4A 4 58A A 3 e A A 3A + 343A A 3 A + 353A 3A 3 + 3α α A 5 56A 4A 63A 3 A 59A A 4 A 88A A 3 84A 3A A 3 467α A A 3A 439α A A 3 A 45α A 3 A α A 4 A + 7α A 3 A + 67α A A 4 A + α A A α A 3 A A 3 6α A 3 A α A 5 8α A A 5 9α A 3 A 4 8α A 4 A 3 5α A 5 A + 99α A A A 3 A 3 45A 5 + 8A A 5 + 9A 3 A 4 + 8A 4 A 3 + 5A 5 A 3α A 5 87A A 4 + 4α A5 e O e 9, e 4 = 6A 3 A 3A α4 + 48A A 3 A 4 α4 + 6A A 3A 3 α4 48A 3A 5 α4 86A7 α3 + 56A3 A 3A α3 6α 3 A A 3 A 4 56α3 A A 3A 3 + 6α3 A 3A 5 + 8α A7 7α A3 A 3A + 36α A A 3 A 4 + 7α A A 3A 3 36α A 3A 5 44α A 7 + 4α A 3 A 3A 64α A A 3 A 4 4α A A 3A 3
8 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, α A 3 A A7 8A3 A 3A + 7A A 3 A 4 + 8A A 3A 3 7A 3A 5 / α 3 α α + e 8 + O e Numerical simulations Numerical simulations are the direct way for checking the validity, accuracy and efficiency of an iterative method. In this section, we solved eleven initial and boundary value problems. The discretizations are performed by using four pseudospectral collocation methods. The details of absolute errors are tabulated for different parameters on different pseudospectral collocation methods. The verification of theoretical convergence order is crucial, and we adopt the following definition of computational convergence order COC COC = log Fx k+ / Fx k log Fx k / Fx k. 4.. Verification of computational convergence order First we consider the following system of nonlinear equations Fx = [F x, F x, F 3 x, F 4 x] T =, F x = x x 3 + x 4 x + x 3 =, F x = x x 3 + x 4 x + x 3 =, F 3 x = x x + x 4 x + x =, F 4 x = x x + x 3 x + x =. 4. Table : MFA : verification of convergence order for the problem 4.. Iters Fx k COC.5e e e e Table shows that COC is nine which is higher than our claimed convergence order. The verification of theoretical convergence order is hard to verify by solving the system of nonlinear equations associated with IVPs and BVPs. 4.. Solving the initial and boundary value problems The temporal and spatial discretizations of all IVPs and BVPs are made by employing different J-GL-C methods. The discrete problems are the systems of nonlinear algebraic equations and we used our proposed iterative method DEDF to solve them. In the majority of the solved problems, the vector is the initial guess. The initial and boundary conditions are embedded in the zero vector, and this inclusion makes the initial guess non-smooth. In most of the numerical simulations, the non-smooth initial guess works well. But in some numerical simulations of schrödinger equations, we get divergence and then we make the initial guess smooth to get the convergence.
9 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Lane-Emden equation x t + t x t + xt p =, x =, x =. 4. The discrete Lane-Emden equation 4. is obtained by using J-GL-C method Fx = S xx x + diag S x x + x p =, t F x = S xx + diag S x + p diag x p, t where diag means a diagonal matrix, S x = /b Q x, S xx = /b Q x, x = [x, x,, x n ] T and [, b] is the domain of the problem. Figure depicts the numerical solutions of different Lane-Emden equations and initial guess. Table shows infinity norm of residue function for different values of index p of Lane-Emden equation. We performed three iterations Iters to show the convergence of our proposed iterative method DEDF. Table : Norm of residue of system of nonlinear equation associated with the problem 4. over the domain [,3], number of grid points = 5, θ, φ = /, /. Iters \ p Fx k 6.58e-3 3.4e-.39e- 3.5e e-35.43e e-7 4.7e e e e-6 3.7e-7.8 initial guess index= index=3 index=4 index= Figure : Solution of the Lane-Emden equation 4. with different indices, number of gird points= Bratu-problem x t + α e xt =, x =, x =, 4.3 where α is a parameter. The associated system of nonlinear equations in 4.3, by using J-GL-C method,
10 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, can be written as Fx = S xx x + α e x =, F x = S xx + α diage x, where S xx = Q x, x = [x, x,, x n ] T. The Bratu problem has a parameter α and we solved it against the several values of parameter α. Figure shows the solution of the Bratu problem for different values of α and norm of residues of associated system of nonlinear equations are listed in Table 3. We can see that the DEDF iterative method leads to a rapid convergence. Table 3: Norm of residue of the Bratu problem 4.3 over the domain [,], number of grid points 5, θ, φ = /, /. Iters \ α 3 Fx k 6.e e-7 6.5e-7 -.7e e-59.5e e-54 5.e-45.e initial guess α= α= α= Figure : Plot of the solution of Bratu problem 4.3 with different parameter values, number of gird points Frank-Kamenetzkii problem Consider the problem where α is a parameter. Hence, we have x t + t x t + α e xt =, x =, x =, 4.4 Fx = S xx x + diag S x x + α e x =, t F x = S xx + diag S x + α diag e x, t where x T = [x, x,, x n ] T. The Frank-Kamenetzkii problem also has a parameter α. The simulated results are depicted in Table 4 and Figure 3.
11 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Table 4: Norm of residue of the problem 4.4 over the domain [,], number of grid points 5, θ, φ = /, /. Iters \ α...3 Fx k.96e-7 5.8e-7 8.3e-7.9e-6-3.e-6.85e-6.3e e e-448.8e-435.e-46.4e initial guess α= α=. α=. α= Figure 3: Plot of the solution of Frank-Kamenetzkii problem 4.4 with different parameter values, number of gird points nonlinear Schrödinger equations Consider the problem where i t ψx, t + xx ψx, t + γ ψx, t ψx, t δ Rx, t ψx, t =, x, t D x D t, with the initial-boundary conditions D x = [a x, b x ], D t = [, t f ], ψa x, t = η t, ψb x, t = η t, ψx, = ζx. The complex function ψx, t can be written as ψx, t = ψ x, t + i ψ x, t. The Schrödinger equation in the real function ψ x, t and ψ x, t is t xx ψ δrx, t ψ + γ ψ + ψ ψ =, 4.5 xx t ψ with initial-boundary conditions ψ a x, t + i ψ a x, t = η t + i η, ψ x, + i ψ x, = ζ x + i ζ x. ψ ψ ψ b x, t + i ψ b x, t = η t + i η t, The discretization of 4.5 by using J-GL-C method leads to Fψ = S t I x I t S xx ψ δ O diagr ψ I t S xx S t I x ψ diagr O ψ O diag ψ + γ + ψ diag ψ + ψ =, 4.6 ψ O ψ
12 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, F ψ = S t I x I t S xx I t S xx δ O diagr + γ diagψ ψ 3 diagψ, S t I x diagr O 3 diagψ diagψ ψ where ψ j = [ψ j,, ψ j,,, ψ j,n t,, ψ jn x,, ψ jn x,,, ψ jn x,n t ]T for j =,, ψ = [ψ T ψt ] T, R = [R,, R,,, R,nt,, R nx,, R nx,,, R nx,n t ] T, is the Kronecker product, is the element-wise multiplication between two vectors, O is the zeros matrix, I is the identity matrix, S t = Q t, S xx = Q x and Q is J-GL-C operational matrix. Further the system of nonlinear t f b x a x equations 4.6 can be written as follows Fψ = Aψ + gψ p, where A = S t I x I t S xx I t S xx δ O diagr, S t I x diagr O gψ = γ diagψ ψ 3 diagψ, 3 diagψ diagψ ψ and p is the vector to incorporate the initial-boundary conditions. A list of four nonlinear Schrödinger equations is given in Table 5, with their corresponding analytical solutions. The numerically computed solutions of four nonlinear Schrödinger equations are tabulated in Tables 6, 7, 8 and 9. In most of the nonlinear Schrödinger equations, the Chebyshev collocation method of first kind shows better accuracy. The numerical solution and absolute errors are visualized in Figures 4, 5, 6 and 7. When the initial guess is non-smooth, due to the inclusion of initial and boundary conditions, the DEDF diverges. Consequently we make the initial guess smooth by applying the following iteration ψ = A gψ p, single time for the problem and two times for the problems and 4. Table 5: Nonlinear Schrödinger equations. Problem Schrödinger equation Analytical solution i t ψ + xx ψ + ψ ψ = ψx, t = e i x+t i t ψ + xx ψ ψ ψ = ψx, t = e i x 3t 3 i t ψ + xxψ ψ ψ cosx ψ = ψx, t = e 3 i t sinx 4 i t ψ + xx ψ + ψ ψ δψ = ψx, t = e i x+ δt x 3.9 ψ ψ.5.5 Absolute error ψ ψ temporal axis spatial axis index Figure 4: Plot of the solution of nonlinear Schrödinger equation problem, θ, φ = /, /, n x = 5, n t = 3.
13 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, x ψ ψ ψ Absolute error Ψ..4.6 temporal axis.8.5 spatial axis index Figure 5: Plot of the solution of nonlinear Schrödinger equation problem, θ, φ = /, /, n x = 5, n t = x 4 ψ ψ 3 ψ ψ Absolute error.5.5 temporal axis.5 spatial axis index Figure 6: Plot of the solution of nonlinear Schrödinger equation problem 3, θ, φ = /, /, n x = 5, n t = 3. 9 x 4 ψ 8 ψ ψ Absolute error temporal axis..5 ψ spatial axis index Figure 7: Plot of the solution of nonlinear Schrödinger equation problem 4, θ, φ = /, /, n x = 5, n t = 3.
14 ux,t D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, x temporal axis. 5 spatial axis Figure 8: Plot of the solution of nonlinear Klein Gorden equation, θ, φ = /, /, n x =, n t = 3. 6 x 5 4 Absolute error index Figure 9: Plot of absolute error for the nonlinear -D wave equation, θ, φ = /, /, n x =, n y =, n t =. 3.5 x Absolute error index Figure : Plot of absolute error for the nonlinear 3-D wave equation, θ, φ = /, /, n x = 7, n y = 7, n z = 7, n t =.
15 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, x Absolute error index Figure : Plot of absolute error for the nonlinear 3-D Poisson equation, θ, φ = /, /, n x =, n y =, n z =. Table 6: Problem in Table 5 over the domain x, t [, ] [, ], number of iterations=, initial guess is ψ = ψ =. Iters \ Grid size ψ ψ analytical 4.33e- 7.e- 7.33e- Legendre polynomials.85e e e-6 3 θ, φ =,.85e-4.6e-7.86e e-4 3.7e-9.57e e e- 5.9e-3 ψ ψ analytical 7.e- divergent 7.7e- θ, φ = /,.59e e e-3.75e e-3 9.7e e-3.45e-3 ψ ψ analytical 4.83e- 7.6e- 7.84e- Chebyshev polynomials of second kind.4e e e-6 3 θ, φ = /, /.4e e-9 5.5e e-3.9e-9 3.e e-3.9e-9 5.9e-3 ψ ψ analytical 3.6e- 6.88e- 7.5e- Chebyshev polynomials of first kind.5e-3 7.e e-6 3 θ, φ = /, /.5e-3.e-9 8.5e e-3 4.9e- 3.54e e-3 4.9e- 8.68e-4
16 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Table 7: Problem in Table 5 over the domain x, t [, ] [, ], number of iterations=, initial guess is ψ = ψ =. Iters \ Grid size ψ ψ analytical 4.33e- 7.e- 7.33e- Legendre polynomials.85e e e-6 3 θ, φ =,.85e-4.6e-7.86e e-4 3.7e-9.57e e e- 5.9e-3 ψ ψ analytical 7.e- divergent 7.7e- θ, φ = /,.59e-3 divergent 7.39e e-3 divergent.75e e-3 divergent 9.7e e-3 divergent.45e-3 ψ ψ analytical 4.83e- 7.6e- 7.84e- Chebyshev polynomials of second kind.4e e e-6 3 θ, φ = /, /.4e e-9 5.5e e-3.9e-9 3.e e-3.9e-9 5.9e-3 ψ ψ analytical 3.6e- 6.88e- 7.5e- Chebyshev polynomials of first kind.5e-3 7.e e-6 3 θ, φ = /, /.5e-3.e-9 8.5e e-3 4.9e- 3.54e e-3 4.9e- 8.68e Klein Gorden equation We consider the problem u tt x, t c u xx x, t + k u γ ux, t 3 =, x, t D x D t, 4.7 where with the initial-boundary conditions D x = [a x, b x ], D t = [, t f ], ua x, t = ρ t, ub x, t = ρ t, ux, = ξ x, u t x, = ξ x. Consequently we find, Fu = S tt u c S xx u + k u γ u 3 =, F u = S tt c S xx + k I 3 γ diag u,
17 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, where u = [u,, u,,, u,nt,, u nx,, u nx,,, u nx,n t ] T, S tt = Q t, S xx = Q x. t f b x a x k The analytical solution of 4.7 is ux, t = δ sech κx v t, where κ = c v and δ = k γ. The parameters with their numerical values are c =, γ =, v =.5 and k =.5. Table shows that the Chebyshev collocation method of first kind exhibits has the best accuracy in the numerical solution for finer grid. The numerical solution with achieved best absolute error is plotted in Figure 8. Table 8: Problem 3 in Table 5 over the domain x, t [, ] [, ], number of iterations=, initial guess is ψ = ψ =. Iters \ Grid size ψ ψ analytical 4.33e- 7.e- 7.33e- Legendre polynomials.85e e e-6 3 θ, φ =,.85e-4.6e-7.86e e-4 3.7e-9.57e e e- 5.9e-3 ψ ψ analytical 7.e- divergent 7.7e- θ, φ = /,.59e e e-3.75e e-3 9.7e e-3.45e-3 ψ ψ analytical 4.83e- 7.6e- 7.84e- Chebyshev polynomials of second kind.4e e e-6 3 θ, φ = /, /.4e e-9 5.5e e-3.9e-9 3.e e-3.9e-9 5.9e-3 ψ ψ analytical 3.6e- 6.88e- 7.5e- Chebyshev polynomials of first kind.5e-3 7.e e-6 3 θ, φ = /, /.5e-3.e-9 8.5e e-3 4.9e- 3.54e e-3 4.9e- 8.68e-4
18 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Table 9: Problem 4 in Table 6 over the domain x, t [, ] [, ], number of iterations=, initial guess is ψ = ψ =. Iters \ Grid size ψ ψ analytical 4.33e- 7.e- 7.33e- Legendre polynomials.85e e e-6 3 θ, φ =,.85e-4.6e-7.86e e-4 3.7e-9.57e e e- 5.9e-3 ψ ψ analytical 7.e- divergent 7.7e- θ, φ = /,.59e e e-3.75e e-3 9.7e e-3.45e-3 ψ ψ analytical 4.83e- 7.6e- 7.84e- Chebyshev polynomials of second kind.4e e e-6 3 θ, φ = /, /.4e e-9 5.5e e-3.9e-9 3.e e-3.9e-9 5.9e-3 ψ ψ analytical 3.6e- 6.88e- 7.5e- Chebyshev polynomials of first kind.5e-3 7.e e-6 3 θ, φ = /, /.5e-3.e-9 8.5e e-3 4.9e- 3.54e e-3 4.9e- 8.68e D nonlinear wave equation We consider the problem u tt x, y, t c u xx x, y, t + u yy x, y, t + fu = pt, x, y, x, y, t D x D y D t, 4.8 where with the initial-boundary conditions D x = [a x, b x ], D y = [a y, b y ], D t = [, t f ], ua x, y, t = ρ y, t, ux, a y, t = ρ x, t, ux, y, = ζx, y. ub x, y, t = ρ y, t, ux, b y, t = ρ x, t, As a consequence we have Fu = S tt I x I y u c I t S xx I y + I t I x S yy u + fu p =, F u = S tt I x I y c I t S xx I y + I t I x S yy + diagf u, where u = [u,,, u,,,, u,,nt,, u nx,n y,, u nx,n y,, u nx,n y,n z ] T, S tt = Q t, t f
19 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, S xx = Q x and S yy = Q y. We assume e t sinx + y as the analytical solution of b x a x b y a y 4.8 with c = and fu = u 3. The Chebyshev collocation method of first kind over the grid achieve best accuracy. Table and Figure 9 represent the results of numerical simulations for the -D nonlinear wave equation. Table : Klein Gorden equation 4.7 over the domain x, t [, ] [, ], number of iterations=, initial guess is u =. Iters \ Grid size u u analytical 8.54e-3.43e-5 3.9e-6 Legendre polynomials 8.53e-3.4e e-9 3 θ, φ =, 8.53e-3.4e e-9 u u analytical 7.67e e-5 3.4e-6 θ, φ = /, 7.67e e-5.45e e e-5.45e-9 u u analytical 7.85e-3.9e-5 3.3e-6 Chebyshev polynomials of second kind 7.85e-3.9e-5.5e-9 3 θ, φ = /, / 7.85e-3.9e-5.5e-9 u u analytical 9.8e-3.58e-5 3.3e-6 Chebyshev polynomials of first kind 9.8e-3.58e e- 3 θ, φ = /, / 9.8e-3.58e-5 3.6e- Table : -D nonlinear wave equation 4.8 over the domain x, y, t [, ] [, ] [, ], number of iterations=, initial guess is u =. Iters \ Grid size 5 5 u u analytical.e-3 divergent.73e- Legendre polynomials.e e- 3 θ, φ =,.e-3.48e- u u analytical.3e-3.45e-7.54e-3 θ, φ = /,.3e-3.95e-7.69e e-3.95e-7.69e- u u analytical.7e-3.e-7 4.4e- Chebyshev polynomials of second kind.7e-3.74e e-6 3 θ, φ = /, /.7e-3.74e-7.57e- u u analytical 7.58e-4 8.3e e-5 Chebyshev polynomials of first kind 7.58e-4.e e- 3 θ, φ = /, / 7.58e-4.e e-
20 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, D nonlinear wave equation We consider the problem u tt x, y, z, t c u xx x, y, z, t + u yy x, y, z, t + u zz x, y, z, t + fu = px, y, z, where with the initial-boundary conditions Therefor we find x, y, z, t D x D y D z D t, D x = [a x, b x ], D y = [a y, b y ], D z = [a z, b z ], D t = [, t f ], ua x, y, z, t = ρ y, z, t, ux, a y, z, t = ρ x, z, t, ux, y, a z, t = ρ 3 x, y, t, ux, y, z, = ζx, y, z. ub x, y, z, t = ρ y, z, t, ux, b y, z, t = ρ x, z, t, ux, y, b z, t = ρ 3 x, y, t, Fu = S tt I x I y I z u c I t S xx I y I z + I t I x S yy I z + I t I x I y S zz u + fu p =, F u = S tt I x I y I z c I t S xx I y I z + I t I x S yy I z + I t I x I y S zz + diagf u, where u = [u,,,, u,,,,, u,,,nt,, u nx,n y,n z,, u nx,n y,n z,,, u nx,n y,n z,n t ] T, S tt = Q t, S xx = Q x, S yy = Q y and S zz = Q z. We assume t f b x a x b y a y b z a z e t sinx + y + z as the solution of 4.9 with c =.8 and fu = u. The Legendre collocation method produces the best numerical results. The absolute error is displayed with respect to different girds in Table. The numerical solution and absolute error can be visualized in Figure D nonlinear Poisson equation We consider the problem where u xx x, y, z + u yy x, y, z + u zz x, y, z + fu = px, y, z, x, y, z D x D y D z, 4. with the initial-boundary conditions Consequently we have D x = [a x, b x ], D y = [a y, b y ], D z = [a z, b z ], ua x, y, z = ρ y, z, ux, a y, z = ρ x, z, ux, y, a z = ρ 3 x, y, ub x, y, z = ρ y, z, ux, b y, z = ρ x, z, ux, y, b z = ρ 3 x, y. Fu = S xx I y I z + I x S yy I z + I x I y S zz u + fu p =, F u = S xx I y I z + I x S yy I z + I x I y S zz + diagf 4. u, where u = [u,,, u,,,, u,,nz,, u nx,n y,, u nx,n y,, u nx,n y,n z ] T, S xx = Q x, b x a x S yy = Q y and S zz = Q z. We assume sinx + y + z as the solution of 4. with b y a y b z a z fu = u 4. In the case of 3-D nonlinear Poisson equation, the Legendre collocation method offers best accuracy over different sizes grids. Figure depicts the absolute error in the numerically computed solution. 4.9
21 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, Table : 3-D nonlinear wave equation 4.9 over the domain x, y, z, t [, ] [, ] [, ] [, ], number of iterations=, initial guess is u =. Iters \ Grid size u u analytical 5.54e-6 7.5e e-9 Legendre polynomials 5.54e e-8 3.9e-9 3 θ, φ =, 5.54e e-8 3.9e-9 u u analytical.4e-5.e-6.45e- θ, φ = /,.4e-5.e-6.6e e-5.e-6.6e-7 u u analytical.8e-5.76e-6 9.5e-7 Chebyshev polynomials.8e-5.76e e-8 3 of second kind.8e-5.76e e-8 u u analytical.9e-5.5e- 3.4e-8 Chebyshev polynomials of first kind.9e-5.8e-6 3.7e-8 3 θ, φ = /, /.9e-5.8e-6 3.7e-8 Table 3: 3-D nonlinear Poisson equation 4. over the domain x, y, z [, ] [, ] [, ], number of iterations=, initial guess is u =. Iters \ Grid size u u analytical 6.83e-8 7.5e-8 7.9e-8 Legendre polynomials.95e-.55e e-5 3 θ, φ =,.95e-.54e-4 7.e-5 u u analytical 7.49e-8 8.7e-8 8.3e-8 θ, φ = /,.33e e e e e e-5 u u analytical 6.93e-8 7.6e e-8 Chebyshev polynomials 9.59e- 6.78e-3.75e-4 3 of second kind 9.59e- 6.78e-3.8e-4 u u analytical 6.65e-8 7.4e e-8 Chebyshev polynomials of first kind 5.63e- 3.6e e-5 3 θ, φ = /, / 5.63e- 3.6e e-5 5. Conclusions We consider higher order frozen Jacobian iterative methods which are computationally efficient. We compute once the frozen Jacobian LU factors. We use them to solve eight lower and upper triangular
22 D. A. Alrehaili, et al., J. Math. Computer Sci., 7 7, systems per iteration. The numerical solution of lower and the upper triangular system is computational economic that makes the entire method computationally efficient. We can see that in most of the simulations, we achieved good accuracy in almost three iterations. Different discretizations collocation methods are used that are derived from J-GL-C method, by changing the parameters. We conclude that in the majority of the cases, the Chebyshev collocation method of first kind gives best accuracy for different sizes grids. Eleven IVPs and BVPs are solved to show the validity, accuracy and efficiency of our higher order iterative method DEDF. Acknowledgment The third author was funded by King Abdulaziz University during scientific communication year 7-8. References [] F. Ahmad, E. Tohidi, J. A. Carrasco, A parameterized multi-step Newton method for solving systems of nonlinear equations, Numer. Algorithms, 7 6, [] F. Ahmad, E. Tohidi, M. Z. Ullah, J. A. Carrasco, Higher order multi-step Jarratt-like method for solving systems of nonlinear equations: application to PDEs and ODEs, Comput. Math. Appl., 7 5, [3] E. S. Alaidarous, M. Z. Ullah, F. Ahmad, A. S. Al-Fhaid, An efficient higher-order quasilinearization method for solving nonlinear BVPs, J. Appl. Math., 3 3, pages. [4] A. H. Bhrawy, An efficient Jacobi pseudospectral approximation for nonlinear complex generalized Zakharov system, Appl. Math. Comput., 47 4, [5] A. H. Bhrawy, E. H. Doha, M. A. Abdelkawy, R. A. Van Gorder, Jacobi-Gauss-Lobatto collocation method for solving nonlinear reaction-diffusion equations subject to Dirichlet boundary conditions, Appl. Math. Model., 4 6, [6] M. Dehghan, F. Fakhar-Izadi, The spectral collocation method with three different bases for solving a nonlinear partial differential equation arising in modeling of nonlinear waves, Math. Comput. Modelling, 53, [7] E. H. Doha, A. H. Bhrawy, M. A. Abdelkawy, R. A. Van Gorder, Jacobi-Gauss-Lobatto collocation method for the numerical solution of + nonlinear Schrdinger equations, J. Comput. Phys., 6 4, [8] E. H. Doha, A. H. Bhrawy, S. S. Ezz-Eldien, Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations, Appl. Math. Model., 35, [9] E. H. Doha, A. H. Bhrawy, R. M. Hafez, On shifted Jacobi spectral method for high-order multi-point boundary value problems, Commun. Nonlinear Sci. Numer. Simul., 7, [] H. Montazeri, F. Soleymani, S. Shateyi, S. S. Motsa, On a new method for computing the numerical solution of systems of nonlinear equations, J. Appl. Math.,, 5 pages. [] J. M. Ortega, W. C. Rheinboldt, Iterative solution of nonlinear equations in several variables, Academic Press, New York-London, 97. [] U. Qasim, Z. Ali, F. Ahmad, F. Ahmad, S. Serra-Capizzano, M. Z. Ullah, M. Asma, Constructing frozen Jacobian iterative methods for solving systems of nonlinear equations, associated with ODEs and PDEs using the homotopy method, Algorithms Basel, 9 6, 7 pages. [3] S. Qasim, Z. Ali, F. Ahmad, S. Serra-Capizzano, M. Z. Ullah, A. Mahmood, Solving systems of nonlinear equations when the nonlinearity is expensive, Comput. Math. Appl., 7 6, [4] J. Shen, T. Tang, L.-L. Wang, Spectral methods, Algorithms, analysis and applications, Springer Series in Computational Mathematics, Springer, Heidelberg,. [5] F. Soleymani, T. Lotfi, P. Bakhtiari, A multi-step class of iterative methods for nonlinear systems, Optim. Lett., 8 4, 5. [6] G. Szegö, Orthogonal Polynomials, American Mathematical Society Colloquium Publications, American Mathematical Society, New York, 939. [7] E. Tohidi, S. Lotfi Noghabi, An efficient Legendre pseudospectral method for solving nonlinear quasi bang-bang optimal control problems, J. Appl. Math. Stat. Inform., 8, [8] J. F. Traub, Iterative methods for the solution of equations, Prentice-Hall Series in Automatic Computation Prentice-Hall, Inc., Englewood Cliffs, N.J., 964. [9] M. Z. Ullah, S. Serra-Capizzano, F. Ahmad, An efficient multi-step iterative method for computing the numerical solution of systems of nonlinear equations associated with ODEs, Appl. Math. Comput., 5 5, [] M. Z. Ullah, F. Soleymani, A. S. Al-Fhaid, Numerical solution of nonlinear systems by a general class of iterative methods with application to nonlinear PDEs, Numer. Algorithms, 67 4, 3 4.
Multi-step frozen Jacobian iterative scheme for solving IVPs and BVPs based on higher order Fréchet derivatives
Punjab University Journal of Mathematics ISSN 116-2526 Vol. 491217 pp. 125 138 Multi-step frozen Jacobian iterative scheme for solving IVPs and BVPs based on higher order Fréchet derivatives Iqra Ilyas
More informationA higher order frozen Jacobian iterative method for solving Hamilton-Jacobi equations
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (26), 62 6227 Research Article A higher order frozen Jacobian iterative method for solving Hamilton-Jacobi equations Ebraheem O. Alzahrani a,
More informationA Preconditioned Iterative Method for Solving Systems of Nonlinear Equations Having Unknown Multiplicity
Article A Preconditioned Iterative Method for Solving Systems of Nonlinear Equations Having Unknown Multiplicity Fayyaz Ahmad 1,2,3, *, Toseef Akhter Bhutta 4, Umar Shoaib 4, Malik Zaka Ullah 1,5, Ali
More informationA Family of Iterative Methods for Solving Systems of Nonlinear Equations Having Unknown Multiplicity
Article A Family of Iterative Methods for Solving Systems of Nonlinear Equations Having Unknown Multiplicity Fayyaz Ahmad 1,2, *, S Serra-Capizzano 1,3, Malik Zaka Ullah 1,4 and A S Al-Fhaid 4 Received:
More informationNumerical Iterative Methods For Nonlinear Problems
UNIVERSITÀ DELL INSUBRIA DIPARTIMENTO DI SCIENZA E ALTA TECNOLOGIA DOCTORAL THESIS Numerical Iterative Methods For Nonlinear Problems Author: Malik Zaka Ullah Supervisor: Prof. Stefano Serra-Capizzano
More informationMulti-step derivative-free preconditioned Newton method for solving systems of nonlinear equations
Multi-step derivative-free preconditioned Newton method for solving systems of nonlinear equations Fayyaz Ahmad Abstract Preconditioning of systems of nonlinear equations modifies the associated Jacobian
More informationFour-Point Optimal Sixteenth-Order Iterative Method for Solving Nonlinear Equations
Four-Point Optimal Sixteenth-Order Iterative Method for Solving Nonlinear Equations Malik Zaka Ullah a, A.S. Al-Fhaid a, Fayyaz Ahmad b, a Department of Mathematics, King Abdulaziz University, Jeddah 21589,
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationTwo improved classes of Broyden s methods for solving nonlinear systems of equations
Available online at www.isr-publications.com/jmcs J. Math. Computer Sci., 17 (2017), 22 31 Research Article Journal Homepage: www.tjmcs.com - www.isr-publications.com/jmcs Two improved classes of Broyden
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
More informationRELIABLE TREATMENT FOR SOLVING BOUNDARY VALUE PROBLEMS OF PANTOGRAPH DELAY DIFFERENTIAL EQUATION
(c) 216 217 Rom. Rep. Phys. (for accepted papers only) RELIABLE TREATMENT FOR SOLVING BOUNDARY VALUE PROBLEMS OF PANTOGRAPH DELAY DIFFERENTIAL EQUATION ABDUL-MAJID WAZWAZ 1,a, MUHAMMAD ASIF ZAHOOR RAJA
More informationNEW NUMERICAL APPROXIMATIONS FOR SPACE-TIME FRACTIONAL BURGERS EQUATIONS VIA A LEGENDRE SPECTRAL-COLLOCATION METHOD
Romanian Reports in Physics, Vol. 67, No. 2, P. 340 349, 2015 NEW NUMERICAL APPROXIMATIONS FOR SPACE-TIME FRACTIONAL BURGERS EQUATIONS VIA A LEGENDRE SPECTRAL-COLLOCATION METHOD A.H. BHRAWY 1,2, M.A. ZAKY
More informationSixth Order Newton-Type Method For Solving System Of Nonlinear Equations And Its Applications
Applied Mathematics E-Notes, 17(017), 1-30 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Sixth Order Newton-Type Method For Solving System Of Nonlinear Equations
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationOptimal derivative-free root finding methods based on the Hermite interpolation
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (2016), 4427 4435 Research Article Optimal derivative-free root finding methods based on the Hermite interpolation Nusrat Yasmin, Fiza Zafar,
More informationNUMERICAL SIMULATION OF TIME VARIABLE FRACTIONAL ORDER MOBILE-IMMOBILE ADVECTION-DISPERSION MODEL
Romanian Reports in Physics, Vol. 67, No. 3, P. 773 791, 2015 NUMERICAL SIMULATION OF TIME VARIABLE FRACTIONAL ORDER MOBILE-IMMOBILE ADVECTION-DISPERSION MODEL M.A. ABDELKAWY 1,a, M.A. ZAKY 2,b, A.H. BHRAWY
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationNUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 37 (2017) 161-174 NUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD Mahua Jahan Rupa 1 and Md. Shafiqul Islam 2
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationChebyshev finite difference method for solving a mathematical model arising in wastewater treatment plants
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 448-455 Chebyshev finite difference method for solving a mathematical model arising in wastewater treatment
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationSolving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted Legendre Functions
International Journal of Mathematical Modelling & Computations Vol. 5, No. 3, Summer 215, 219-23 Solving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted
More informationULTRASPHERICAL WAVELETS METHOD FOR SOLVING LANE-EMDEN TYPE EQUATIONS
ULTRASPHERICAL WAVELETS METHOD FOR SOLVING LANE-EMDEN TYPE EQUATIONS Y. H. YOUSSRI, W. M. ABD-ELHAMEED,, E. H. DOHA Department of Mathematics, Faculty of Science, Cairo University, Giza 63, Egypt E-mail:
More informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationLecture 24: Starting to put it all together #3... More 2-Point Boundary value problems
Lecture 24: Starting to put it all together #3... More 2-Point Boundary value problems Outline 1) Our basic example again: -u'' + u = f(x); u(0)=α, u(l)=β 2) Solution of 2-point Boundary value problems
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction 1 Introduction In this module, we develop solution techniques for numerically solving ordinary
More informationPower Series and Analytic Function
Dr Mansoor Alshehri King Saud University MATH204-Differential Equations Center of Excellence in Learning and Teaching 1 / 21 Some Reviews of Power Series Differentiation and Integration of a Power Series
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationRATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages 72 83 c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR
More informationA new family of four-step fifteenth-order root-finding methods with high efficiency index
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 3, No. 1, 2015, pp. 51-58 A new family of four-step fifteenth-order root-finding methods with high efficiency index Tahereh
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationNUMERICAL SOLUTION FOR CLASS OF ONE DIMENSIONAL PARABOLIC PARTIAL INTEGRO DIFFERENTIAL EQUATIONS VIA LEGENDRE SPECTRAL COLLOCATION METHOD
Journal of Fractional Calculus and Applications, Vol. 5(3S) No., pp. 1-11. (6th. Symp. Frac. Calc. Appl. August, 01). ISSN: 090-5858. http://fcag-egypt.com/journals/jfca/ NUMERICAL SOLUTION FOR CLASS OF
More informationExact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol 6, 212, no 122, 697-618 Exact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method Abdelhalim Ebaid 1 and Mona D Aljoufi
More informationSolving high-order partial differential equations in unbounded domains by means of double exponential second kind Chebyshev approximation
Computational Methods for Differential quations http://cmde.tabrizu.ac.ir Vol. 3, No. 3, 015, pp. 147-16 Solving high-order partial differential equations in unbounded domains by means of double exponential
More informationBoundary Layer Problems and Applications of Spectral Methods
Boundary Layer Problems and Applications of Spectral Methods Sumi Oldman University of Massachusetts Dartmouth CSUMS Summer 2011 August 4, 2011 Abstract Boundary layer problems arise when thin boundary
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationA Chebyshev-Gauss-Radau Scheme For Nonlinear Hyperbolic System Of First Order
Appl. Math. Inf. Sci. 8, o., 535-544 (014) 535 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.1785/amis/08011 A Chebyshev-Gauss-Radau Scheme For onlinear Hyperbolic
More informationAn Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley Equation
Applied Mathematical Sciences, Vol. 11, 2017, no. 30, 1467-1479 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.7141 An Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley
More informationExact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Published by Canadian Center of Science and Education Exact Solutions for the Nonlinear +-Dimensional Davey-Stewartson Equation
More informationA robust uniform B-spline collocation method for solving the generalized PHI-four equation
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 11, Issue 1 (June 2016), pp. 364-376 Applications and Applied Mathematics: An International Journal (AAM) A robust uniform B-spline
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More informationResearch Article Hermite Wavelet Method for Fractional Delay Differential Equations
Difference Equations, Article ID 359093, 8 pages http://dx.doi.org/0.55/04/359093 Research Article Hermite Wavelet Method for Fractional Delay Differential Equations Umer Saeed and Mujeeb ur Rehman School
More informationResearch Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden Type Equations
Applied Mathematics Volume 2012, Article ID 103205, 13 pages doi:10.1155/2012/103205 Research Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationGeneralized Lagrange Jacobi Gauss-Lobatto (GLJGL) Collocation Method for Solving Linear and Nonlinear Fokker-Planck Equations
Commun. Theor. Phys. 69 (2018 519 531 Vol. 69, No. 5, May 1, 2018 Generalized Lagrange Jacobi Gauss-Lobatto (GLJGL Collocation Method for Solving Linear and Nonlinear Fokker-Planck Equations K. Parand,
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationAn Efficient Numerical Solution of Nonlinear Hunter Saxton Equation
Commun. Theor. Phys. 67 (2017) 483 492 Vol. 67, No. 5, May 1, 2017 An Efficient Numerical Solution of Nonlinear Hunter Saxton Equation Kourosh Parand 1,2, and Mehdi Delkhosh 1 1 Department of Computer
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationChebyshev Collocation Spectral Method for Solving the RLW Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.7(2009) No.2,pp.131-142 Chebyshev Collocation Spectral Method for Solving the RLW Equation A. H. A. Ali Mathematics
More informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationMath Assignment 11
Math 2280 - Assignment 11 Dylan Zwick Fall 2013 Section 8.1-2, 8, 13, 21, 25 Section 8.2-1, 7, 14, 17, 32 Section 8.3-1, 8, 15, 18, 24 1 Section 8.1 - Introduction and Review of Power Series 8.1.2 - Find
More informationThe automatic solution of PDEs using a global spectral method
The automatic solution of PDEs using a global spectral method 2014 CBMS-NSF Conference, 29th June 2014 Alex Townsend PhD student University of Oxford (with Sheehan Olver) Supervised by Nick Trefethen,
More informationA Posteriori Existence in Numerical Computations
Oxford University Computing Laboratory OXMOS: New Frontiers in the Mathematics of Solids OXPDE: Oxford Centre for Nonlinear PDE January, 2007 Introduction Non-linear problems are handled via the implicit
More informationLeast-squares Solutions of Linear Differential Equations
1 Least-squares Solutions of Linear Differential Equations Daniele Mortari dedicated to John Lee Junkins arxiv:1700837v1 [mathca] 5 Feb 017 Abstract This stu shows how to obtain least-squares solutions
More informationPartitioned Methods for Multifield Problems
Partitioned Methods for Multifield Problems Joachim Rang, 10.5.2017 10.5.2017 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 Contents Blockform of linear iteration schemes Examples 10.5.2017
More informationHigh-order Newton-type iterative methods with memory for solving nonlinear equations
MATHEMATICAL COMMUNICATIONS 9 Math. Commun. 9(4), 9 9 High-order Newton-type iterative methods with memory for solving nonlinear equations Xiaofeng Wang, and Tie Zhang School of Mathematics and Physics,
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationOrdinary Differential Equations
Ordinary Differential Equations Part 4 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Ordinary Differential Equations p. 1/23 Two-point Boundary Value Problems NRiC 17. BCs specified at two
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationResearch Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal Polynomial
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 4673, pages doi:.54//4673 Research Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationChebyshev-Halley s Method without Second Derivative of Eight-Order Convergence
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 4 2016, pp. 2987 2997 Research India Publications http://www.ripublication.com/gjpam.htm Chebyshev-Halley s Method without
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationAn Efficient Numerical Scheme for Solving Fractional Optimal Control Problems. 1 Introduction
ISS 1749-3889 (print), 1749-3897 (online) International Journal of onlinear Science Vol.14(1) o.3,pp.87-96 An Efficient umerical Scheme for Solving Fractional Optimal Control Problems M. M. Khader, A.
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationA new kind of double Chebyshev polynomial approximation on unbounded domains
Koç and Kurnaz Boundary Value Problems 2013, 2013:10 R E S E A R C H Open Access A new kind of double Chebyshev polynomial approximation on unbounded domains Ayşe Betül Koç * and Aydın Kurnaz * Correspondence:
More informationA Note on Integral Transforms and Partial Differential Equations
Applied Mathematical Sciences, Vol 4, 200, no 3, 09-8 A Note on Integral Transforms and Partial Differential Equations Adem Kılıçman and Hassan Eltayeb Department of Mathematics and Institute for Mathematical
More informationResearch Article New Algorithm for the Numerical Solutions of Nonlinear Third-Order Differential Equations Using Jacobi-Gauss Collocation Method
Mathematical Problems in Engineering Volume 20, Article ID 83728, 4 pages doi:0.55/20/83728 Research Article New Algorithm for the Numerical Solutions of Nonlinear hird-order Differential Equations Using
More informationS.F. Xu (Department of Mathematics, Peking University, Beijing)
Journal of Computational Mathematics, Vol.14, No.1, 1996, 23 31. A SMALLEST SINGULAR VALUE METHOD FOR SOLVING INVERSE EIGENVALUE PROBLEMS 1) S.F. Xu (Department of Mathematics, Peking University, Beijing)
More informationA CRANK NICHOLSON METHOD WITH FIVE POINTS TO SOLVE FISHER S EQUATION
International J. of Math. Sci. & Engg. Appls. (IJMSEA ISSN 097-9424, Vol. 10 No. I (April, 2016, pp. 177-189 A CRANK NICHOLSON METHOD WITH FIVE POINTS TO SOLVE FISHER S EQUATION MUSA ADAM AIGO Umm Al-Qura
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationA PSEUDOSPECTRAL METHOD FOR SOLVING THE TIME-FRACTIONAL GENERALIZED HIROTA SATSUMA COUPLED KORTEWEG DE VRIES SYSTEM
A PSEUDOSPECTRAL METHOD FOR SOLVING THE TIME-FRACTIONAL GENERALIZED HIROTA SATSUMA COUPLED KORTEWEG DE VRIES SYSTEM M. A. SAKER,2, S. S. EZZ-ELDIEN 3, A. H. BHRAWY 4,5, Department of Basic Science, Modern
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationThe Derivation of Interpolants for Nonlinear Two-Point Boundary Value Problems
European Society of Computational Methods in Sciences and Engineering ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics JNAIAM) vol. 1, no. 1, 2006, pp. 49-58 ISSN 1790 8140 The
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More information(Received 05 August 2013, accepted 15 July 2014)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.18(2014) No.1,pp.71-77 Spectral Collocation Method for the Numerical Solution of the Gardner and Huxley Equations
More information