Solution of Fourth Order Obstacle Problems Using Quintic B-Splines
|
|
- Claribel Atkins
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol. 6, 202, no. 94, Solution of Fourth Order Obstacle Problems Using Quintic B-Splines Shahid S. Siddiqi, Ghazala Akram and Kalsoom Arshad Department of Mathematics University of the Punjab, Lahore, Pakistan Abstract In this paper, a numerical method is developed for solving a system of fourth order boundary value problem associated with obstacle, contact and unilateral problems. It is shown that the Quintic B-spline method is very effective tool and yields better results. Numerical examples are presented to illustrate the applicability of the new method. Keywords: Quintic B-Spline; Boundary Value Problems; Variational Inequalities; Obstacle Problems; Unilateral Problems Introduction For studying the contact, unilateral, obstacle and equilibrium problems arising in different branches of pure and applied sciences, variational inequality theory has become an effective and powerful tool. Variational inequality theory has proved to be immensely useful in the study of many branches of mathematical and engineering sciences. Penalty methods and projection methods have been developed for the solution of general variational inequalities(see[-4, 6, 8-0, 2, 4, 9]). The projection methods are not supposed to be suitable as in projection methods, the projection is needed, which is difficult to be obtained. The penalty methods are inefficient as in these methods, instability is created. However, the general variational inequalities can be characterized by a system of differential equations using the penalty function technique, if the obstacle function is known. This technique was used by Lewy and Stampacchia [8] to study the regularity of the solution of variational inequalities. The main advantage of this technique is its simple applicability in solving obstacle and
2 4652 S. S. Siddiqi, Gh. Akram and K. Arshad unilateral problems. This technique has been used for solving fourth order system of differential equations associated with obstacle and unilateral problems by finite difference and quintic spline methods [, 2, 9]. Noor and Tirmizi [0] solved the system of second order boundary value problems using pade approximation. Al-Said [4] developed the method for the solution of system of second order boundary value problems using cubic spline. It is claimed that the method can be considered as an improvement for the cubic spline method developed in [3]. Gao and Chi [6] solved a system of third-order boundary value problems associated with third-order obstacle problems using the quartic B-splines and the method is claimed to be of second order. Usmani [5], developed the method of the solution of fourth order boundary value problem, considering it to be the problem of bending a rectangular clamped beam of length l resting on an elastic foundation. The vertical deflection w of the beam satisfies the system [ ( )] L + K D w = D } q(x), (. ) w(0) = w(l) =w (0) = w (l) =0, where L d4 dx 4, D is the flexural rigidity of the beam, K is the spring constant of the elastic foundation and the load q(x) acts vertically downwards per unit length of the beam. Usmani [6], developed the discrete methods for the solution of fourth order linear special case boundary value problem, similar to the problem (.) with the change in boundary conditions in terms of second order derivatives instead of first order derivatives. Twizell and Tirmizi [5] developed and analysed a sixth-order method for the numerical solution of the linear fourth-order boundary value problem for which the boundary conditions are given in terms of functional values and (i) first-order derivatives (the clamped-clamped beam problem) or (ii) second-order derivatives (the simple-simple beam problem). Papamichael and Worsey [3] developed the cubic spline method for the solution of problem (.). Daele et al. [] developed second order method for the solution of fourth-order boundary value problem using non polynomial spline which is mixture of mixed spline function consisting of cubic polynomial function and a trigonometric function. Siddiqi and Ghazala [7, 8] developed the non polynomial spline technique for the solution of linear special case fifth and eighth order boundary value problems, respectively. Siddiqi et al. [20] developed the method for the solution of linear special case fifth order boundary value problems using the non polynomial sextic spline. The method is proved to be fifth order convergent. Siddiqi et al. [2] presented the method for the solution of sixth order boundary value problems using quintic spline. Ghazala and Siddiqi [7] solved the sixth order linear special case boundary value prob-
3 Fourth order obstacle problems 4653 lems using non polynomial spline. The following system of fourth order boundary value problems, is considered f(x), a x c, y (4) (x) = f(x)+y(x)g(x)+r, c x d, f(x), d x b, (.2 ) along with the boundary conditions y(a) =y(b) =α 0, y (a) =y (b) = α, y(c) =y(d) =α 2, y (c) =y (d) = α 3, } (.3 ) where r and α i, i =0,, 2, 3 are finite real constants and the functions f(x) and g(x) are continuous on [a, b] and [c, d], respectively. Such type of systems arise in connection with contact, obstacle and unilateral problems. In this paper, quintic B-spline function is used to develop a technique for the solution of system (.2). In Section 2, a class of contact problems in elasticity is considered and is formulated in terms of variational inequalities. The quintic B-spline method for the solution of system (.2) is derived in Section 3. In Section 4, two examples are considered for the usefulness of the method developed. 2 Formulation Khan et al. [] considered the linear fourth order boundary value problem describing the equilibrium configuration of an elastic beam, pulled at the ends and lying over an elastic obstacle of constant height /4 and unit rigidity of the type y (4) f(x), y ψ(x), (y (4) f(x)) (y(x) ψ(x)) = 0, y(0) = y() = y (0) = y () = 0, on Ω = [0, ], (2. ) where f is a given force acting on the beam string and ψ(x) is the elastic obstacle. To study the problem (2.) via the variational inequality formulation, the set K can be defined as K = {v : v H 2 0(Ω) : v ψ,on Ω}, (2.2 ) which is a closed convex set in H 2 0 (Ω), where H2 0 (Ω) is a Sobolev space [9, 2, 4], which is basically a Hilbert space.
4 4654 S. S. Siddiqi, Gh. Akram and K. Arshad The Kikuchi and Oden technique [2], shows that the energy functional associated with the obstacle problem (2.) is where I[v] = = 0 0 {d 4 v/dx 4 2f(x)} v(x)dx, v H 2 0(Ω), (d 2 v/dx 2 ) 2 dx 2 0 f(x)v(x)dx, = a(v, v) 2 <f,v>, (2.3 ) a(u, v) = and <f,v> = 0 0 ( d 2 u/dx 2)( d 2 v/dx 2) dx (2.4 ) f(x)v(x)dx. (2.5 ) It can be proved that the form a(u, v) defined by (2.4) is bilinear, positive and symmetric. Moreover, the functional f defined by (2.5) is linear and continuous. The minimum y of the functional I[v] defined by (2.3) on the closed convex set K in H0 2 (Ω) can be characterized by the following variational inequality [9, 2, 4] a(y, v y) <f,v y>, v K. (2.6 ) The obstacle problem (2.) is, thus, equivalent to solving the variational inequality problem (2.6). The equivalence has been used to study the existence of a unique solution of (2.), see for example [9, 4]. Following Lewy and Stampacchia [8], the problem (2.6) can be written as y (4) + μ(y ψ)(y ψ) =f, 0 <x<, (2.7 ) y(0) = y() = 0, y (0) = y () = ɛ, (2.8 ) where ɛ is a small positive constant, ψ is the obstacle function and μ(t) is the penalty function defined by μ(t) = { 4, t 0, 0, t < 0. (2.9 ) Equation (2.6) describes the equilibrium configuration of an elastic beam, pulled at the ends and lying over an elastic obstacle of constant height /4. Since the obstacle function ψ is known, it is possible to find the exact solution of the problem in the interval /4 x 3/4. Assuming that the obstacle function ψ is defined by { /4, 0 x /4, 3/4 x, ψ(x) = /4, /4 x 3/4. (2.0 )
5 Fourth order obstacle problems 4655 From equation (2.6) (2.9), the following system of equations can be obtained as { f, 0 x /4, y (4) = 4y + f, /4 x 3/4, 3/4 x, (2. ) with the following boundary conditions y(0) = y(/4) = y(3/4) = y() = 0, y (0) = y (/4) = y (3/4) = y () = ɛ, (2.2 ) such that y and y are continuous at x =/4 and 3/4, (see []). 3 Quintic B-Spline Method Extending the technique developed by Gao and Chi [6], the boundary value problem (2.) is transformed into the following initial value problems { y (4) 0, 0 x /4, 3/4 x, = 4y, /4 x 3/4, (3. ) with initial values y (0) = 0, y () (0) = 0, y (2) (0) = 0, y (3) (0) = and { y (4) f, 0 x /4, 3/4 x, 2 = 4y 2 + f, /4 x 3/4, with initial values (3.2 ) y 2 (0) = 0, y () 2 (0) = 0, y (2) 2 (0) = 0, y (3) 2 (0) = 0. It may be noted that the exact solution of the problem (2.) is y(x) =y 2 (x) (y 2 ()/y ()) y (x). (3.3 ) To determine the quintic B-spline solution, the interval [0, ] is divided into n equal subintervals with knots x i = ih, i =0,, 2,..., n, h =/n, where n is a multiple of 4. Let Ω n = {x 5,x 4,x 3,x 2,x,x 0,x..., x n } (3.4 )
6 4656 S. S. Siddiqi, Gh. Akram and K. Arshad be the set of extended equally spaced knots with x i = ih. It may be mentioned that the quintic B-spline s(x) = n i= 5 γ i B i,5 (x) (3.5 ) is a polynomial of degree five in each subinterval [x i, x i+ ]. It may also be noted that the quintic B-spline s, by definition, C 4 [0, ]. Let s (x) = n i= 5 α ib i,5 (x) and s 2 (x) = n i= 5 β ib i,5 (x) be the approximate solution to y (x) and y 2 (x), respectively. Moreover, s (x) and s 2 (x) are determined such that they must satisfy the initial problems (3.) and (3.2) respectively. i.e. s (4) (x i )= 0, 0 x i /4, 4s (x i ), /4 x i 3/4, 0, 3/4 x i, alongwith the conditions s (0)=0, s () (0)=0, s (2) (0)=0, s (3) (0) =, and s (4) 2 (x i )= f(x i ), 0 x i /4, f(x i ) 4s 2 (x i )+, /4 x i 3/4, f(x i ), 3/4 x i, (3.6 ) (3.7 ) along with the conditions s 2 (0)=0, s () 2 (0)=0, s (2) 2 (0)=0, s (3) 2 (0) = 0. The first five coefficient α 5, α 4, α 3,α 2, α and β 5, β 4, β 3, β 2, β of s (x) and s 2 (x), x [x 0, xn ] can be determined by applying the initial 4 conditions which leads to the following s (4) (x i )=/h 4 [α i 5 3α i 4 +6α i 3 3α i 2 + α i ], (3.8 ) and s (x i )=/20[α i 5 +26α i 4 +66α i 3 +26α i 2 + α i ], (3.9 ) s (4) 2 (x i )=/h 4 [β i 5 3β i 4 +6β i 3 3β i 2 + β i ] (3.0 ) s 2 (x i )=/20[β i 5 +26β i 4 +66β i 3 +26β i 2 + β i ]. (3. ) From Eqns. (3.5) (3.0), it can be written as α i = α i 5 +4α i 4 6α i 3 +4α i 2, α i 5 ( 4/h 4 +26/30)/(/h 4 +/30)α i 4 (6/h 4 +66/30)/(/h 4 +/30)α i 3 ( 4/h 4 +26/30)/(/h 4 +/30)α i 2, i n/4, (3.2 ) n/4 i 3n/4, α i 5 +4α i 4 6α i 3 +4α i 2, 3n/4 x i n,
7 Fourth order obstacle problems 4657 β i = β i 5 +4β i 4 6β i 3 +4β i 2 + h 4 f(x i ), x i n/4, β i 5 ( 4/h 4 +26/30)/(/h 4 +/30)β i 4 (6/h 4 +66/30)/(/h 4 +/30)β i 3 ( 4/h 4 +26/30)/(/h 4 (3.3 ) +/30)β i 2 +(f(x i )+)/(/h 4 +/30), n/4 x i 3n/4, β i 5 +4β i 4 6β i 3 +4β i 2 + h 4 f(x i ), 3n/4 x i n, 4 Error Bound Considering the quintic splines s (x) and s 2 (x) to be the approximate solutions to y (x) and y 2 (x) respectively. Since s (x) and s 2 (x) are defined piecewise, on [x 0, x n ], therefore error bound of s (x) over any of the three subintervals [x 0,x n/4 ], [x n/4,x 3n/4 ] and [x 3n/4,x n ] will help in calculating the error bound over [x 0, x n ]. The error bound of s (x) over [x 0, x n/4 ] can be calculated as under: Consider the following initial value problem y (x) =f, x [x 0,x n/4 ] (4. ) along with the following conditions s (0) = y (0), s () (0) = y () (0), s (2) (0) = y (2) (0), s (3) (0) = y (3) (0), s (4) (0) = y (4) (0), (4.2 ) let the error term corresponding to each knot can be denoted by then e (x i )=s (x i ) y (x i ), i =0,,..., n 4, (4.3 ) e (x 0 )=s (x 0 ) y (x 0 )=0, e () (x 0 )=s () (x 0 ) y () e (2) (x 0 )=s (2) (x 0 ) y (2) (x 0 )=0, e (3) (x 0 )=s (3) (x 0 ) y (3) e (4) (x 0 )=s (4) (x 0 ) y (4) Expanding s (x) and y (x) as (x 0 )=0. (x 0 )=0, (x 0 )=0, (4.4 ) s (x ) = s(x 0 )+hs () (x 0 )+s (2) (x 0 )h 2 /2+s (3) (x 0 )h 3 /6+s (4) (x 0 )h 4 /4! + O(h 5 ) y (x ) = y(x 0 )+hy () (x 0 )+y (2) (x 0 )h 2 /2+y (3) (x 0 )h 3 /6+s (4) (x 0 )h 4 /4! + O(h 5 )
8 4658 S. S. Siddiqi, Gh. Akram and K. Arshad The order of error terms can be expressed as e (x ) = s(x ) y(x )=O(h 5 ) e () (x ) = s () (x ) y () (x )=O(h 4 ) e (2) (x ) = s (2) (x ) y (2) (x )=O(h 3 ) e (3) (x ) = s (3) (x ) y (3) (x )=O(h 2 ) and e (4) (x ) = s (4) (x ) y (4) (x )=O(h). Similarly, expanding s (x) and y (x) using the Taylor s series, give s (x 2 ) = s(x )+hs () (x )+s (2) (x )h 2 /2+s (3) (x )h 3 /6+s (4) (x )h 4 /4! + O(h 5 ) y (x 2 ) = y(x )+hy () (x )+y (2) (x )h 2 /2+y (3) (x )h 3 /6+s (4) (x )h 4 /4! + O(h 5 ) The order of error terms can be expressed as In general, it can be written as e (x 2 ) = e(x )+O(h 5 ), e () (x 2 ) = e () (x )+O(h 4 ), e (2) (x 2 ) = e (2) (x )+O(h 3 ), e (3) (x 2 ) = e (3) (x )+O(h 2 ), and e (4) (x 2 ) = e (4) (x )+O(h). which shows that e (x i+ ) = e (x i )+O(h 5 ), (4.5 ) and e () (x i+ ) = e () (x i )+O(h 4 ). (4.6 ) e (x n ) = n 4 4 O(h5 )=O(h 4 ) and e () (x n ) = n 4 4 O(h4 )=O(h 3 ). (4.7 ) Similarly, the error term of s 2 can be determined to be of O(h 4 ), which shows that the error bounds of s and s 2 over the remaining two subintervals are also of O(h 4 ). Thus from Eq. (3.3) it can be concluded that quintic B-spline solution of BVP (2.) is second order convergent method.
9 Fourth order obstacle problems Numerical Results In this section, two examples are considered to illustrate the quintic B-spline method. Example When f = 0, the following problem can be considered as 0, 0 x /4, y (4) (x) = 4y(x), /4 x 3/4, (5. ) 0, 3/4 x, with the conditions y(0) = 0, y () (0) = 0, y() = 0, y () () = 0. The analytic solution of the problem (5.) is y(x) = 0, 0 x /4, exp(x)( cos(x) sin(x) + exp( x)( cos(x) sin(x), /4 x 3/4, 0, 3/4 x. (5.2 ) Table : Maximum absolute errors for the problem (5.) in y i. h Maximum absolute errors Example 2 For f =, the following BV P is considered as, 0 x /4, y (4) (x) = 4y(x)+2, /4 x 3/4,, 3/4 x, (5.3 ) with the conditions y(0) = 0, y () (0) = 0, y() = 0, y () () = 0.
10 4660 S. S. Siddiqi, Gh. Akram and K. Arshad The analytic solution of the problem (5.2) is x 4 /24 x 3 /48 + x 2 /384, 0 x /4, exp(x)( cos(x) sin(x) y(x) = (5.4 ) + exp( x)( cos(x) sin(x), /4 x 3/4, x 4 /24 7x 3 / x 2 /384 7/64x +3/28, 3/4 x. Conclusion. A new B-spline method for solving a system of fourth order Table 2: Maximum absolute errors for the problem (5.2) in y i. h Maximum absolute errors problems is developed. The present method enables us to approximate the solution as well as its first, second, third, derivative at every point of range of integration. It has been observed that the results obtained from the method developed in the paper is acceptable. References [] A. Khan, M. A. Noor and T. Aziz, Parametric Quintic-Spline Approach to the Solution of a System of Fourth-Order Boundary-Value Problems, Journal of Optimization Theory and Applications, Vol 22, No. 2, pp , [2] E.A. Al-Said and M. A. Noor, Computational Methods for Fourth-Order Obstacle Boundary-Value Problems, Communications in Applied Nonlinear Analysis, Vol 2, pp , 995. [3] E.A. Al-Said, Spline Methods for Solving System of Second-Order Boundary-Value Problems, Int. J. Comput. Math. Vol 70, pp , 999. [4] E.A. Al-Said, The Use of Cubic Splines in the Numerical Solution of a System of Second-Order Boundary Value Problems, Computers and Mathematics with Applications, Vol 42, pp , 200.
11 Fourth order obstacle problems 466 [5] E.H. Twizell and S.I.A. Tirmizi, A Sixth-Order Multiderivative Method For Two Beam Problems, International Journal for Numerical Methods in Engineering, Vol 23, pp , 986. [6] Feng Gao and Chun-Mei Chi, Solving Third-Order Obstacle Problems with Quartic B-Splines, Applied Mathematics and Computation, Vol 80, pp , [7] Ghazala Akram and Shahid S. Siddiqi, Solution of Sixth Order Boundary Value Problems using Non-Polynomial Spline Technique, Applied Mathematics and Computation 8 (2006) [8] H. Lewy and G. Stampacchia, On the Regularity of the Solution of the Variational Inequalities, Communications in Pure and Applied Mathematics, Vol 22, pp , 969. [9] J. Crank, Free and Moving Boundary-Value Problems, Clarendon Press, Oxford, UK, 984. [0] M. A. Noor and S. I. A. Tirmizi, Highly Accurate Methods for Solving Unilateral Problems, Punjab Univ. Journal of Mathematics, Vol. 22, pp. 9-30, [] M. Van Daele, G. Vanden Berghe and H. De Meyer, A Smooth Approximation for the Solution of a Fourth-Order Boundary Value Problem based on Nonpolynomial Splines, Journal of Computational and Applied Mathematics, Vol 5, No. 3, pp , 994. [2] N. Kikuchi and T. J. Oden, Contact Problems in Elasticity, SIAM, Philadelphia, Pennsylvania 988. [3] N. Papamichael and A. J. Worsey, A Cubic Spline Method For The Solution of a Linear Fourth-Order Two-Point Boundary Value Problem, J. Comptut. Appl. Math. Vol. 7, pp , 98. [4] R. Glowinski, J. L. Lions and R. Tremolieres, Numerical Analysis of Variational Inequalities, North-Holland, Amsterdam, Holland, 98. [5] Riaz A. Usmani, Discrete Variable Methods for a Boundary Value Problem with Engineering Applications, Mathematics of Computation, Vol. 32 No. 44, pp , 978. [6] Riaz A. Usmani, Discrete Methods for Boundary Value Problems with Applications in Plate Deflection Theory, ZAMP, Vol. 30, pp , 979.
12 4662 S. S. Siddiqi, Gh. Akram and K. Arshad [7] Shahid S. Siddiqi and Ghazala Akram, Solution of Fifth-Order Boundary- Value Problems using Non Polynomial Spline Technique, Applied Mathematics and Computation, Vol. 75, No. 2, pp , [8] Shahid S. Siddiqi and Ghazala Akram, Solution of Eighth-Order Boundary-Value Problems using Non Polynomial Spline Technique, International Journal of Computer Mathematics, Vol. 84,No. 3, pp , [9] Shahid S. Siddiqi and Ghazala Akram, Solution of tenth order boundary value problems using non polynomial spline technique, Applied Mathematics and Computation, Vol. 90, pp , [20] Shahid S. Siddiqi, Ghazala Akram and Salman Amin Malik, Non Polynomial Sextic Spline method for the Solution along with the convergence of linear Special Case Fifth Order Two-point boundary value problems, Applied Mathematics and Computation, Vol. 90, pp , [2] Shahid S. Siddiqi, Ghazala Akram and Saima Nazeer, Quintic Spline Solutions of Linear Sixth Order Boundary Value Problems, Applied Mathematics and Computation, Vol. 89, pp , Received: April, 202
QUINTIC SPLINE SOLUTIONS OF FOURTH ORDER BOUNDARY-VALUE PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Number, Pages 0 c 2008 Institute for Scientific Computing Information QUINTIC SPLINE SOLUTIONS OF FOURTH ORDER BOUNDARY-VALUE PROBLEMS
More informationSolution of a Fourth Order Singularly Perturbed Boundary Value Problem Using Quintic Spline
International Mathematical Forum, Vol. 7, 202, no. 44, 279-290 Solution of a Fourth Order Singularly Perturbed Boundary Value Problem Using Quintic Spline Ghazala Akram and Nadia Amin Department of Mathematics
More informationSolution of Fourth Order Boundary Value Problems by Numerical Algorithms Based on Nonpolynomial Quintic Splines
Journal of Numerical Mathematics and Stochastics, 4(1) : 13-25, 2012 http://www.jnmas.org/jnmas4-2.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 Solution of Fourth Order Boundary Value Problems
More informationNon-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems
Non-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems Jalil Rashidinia Reza Jalilian Abstract In this paper we use quintic non-polynomial spline functions to develop numerical
More informationList of Publications of Dr. Ghazala Akram
List of Publications List of Publications of Dr. Ghazala Akram Published 1). End Conditions for Interpolatory Septic Ghazala Akram and Shahid S. Siddiqi, International Journal of Computer Mathematics,
More information1. Introduction We consider the following self-adjoint singularly perturbed boundary value problem: ɛu + p(x)u = q(x), p(x) > 0,
TWMS J. Pure Appl. Math. V., N., 00, pp. 6-5 NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR THE SOLUTION OF SINGULARLY-PERTURBED BOUNDARY VALUE PROBLEMS J. RASHIDINIA, R. MOHAMMADI Abstract. We consider the
More informationSolution of Seventh Order Boundary Value Problem by Differential Transformation Method
World Applied Sciences Journal 16 (11): 1521-1526, 212 ISSN 1818-4952 IDOSI Publications, 212 Solution of Seventh Order Boundary Value Problem by Differential Transformation Method Shahid S. Siddiqi, Ghazala
More informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationVariation of Parameters Method for Solving Fifth-Order. Boundary Value Problems
Applied Mathematics & Information Sciences 2(2) (28), 135 141 An International Journal c 28 Dixie W Publishing Corporation, U. S. A. Variation of Parameters Method for Solving Fifth-Order Boundary Value
More informationSolving Linear Sixth-Order Boundary Value Problems by Using Hyperbolic Uniform Spline Method
International Journal of Mathematical Modelling & Computations Vol. 03, No. 03, 013, 169-180 Solving Linear Sixth-Order Boundary Value Problems by Using Hyperbolic Uniform Spline Method J. Dabounou a,
More informationTransactions on Modelling and Simulation vol 6, 1993 WIT Press, ISSN X
Auxiliary principle technique for a class of nonlinear variational inequalities M.A. Noor, E.A. Al-Said Mathematics Department, College of Science, King Saud University, PO Box 2455, Riyadh 11451, Saudia
More informationNumerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets
Applied Mathematical Sciences, Vol. 5, 20, no. 63, 33-346 Numerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets Fazal-i-Haq Department of Maths/Stats/CS Khyber Pukhtoon Khwa Agricultral
More informationVARIATIONAL ITERATION HOMOTOPY PERTURBATION METHOD FOR THE SOLUTION OF SEVENTH ORDER BOUNDARY VALUE PROBLEMS
VARIATIONAL ITERATION HOMOTOPY PERTURBATION METHOD FOR THE SOLUTION OF SEVENTH ORDER BOUNDARY VALUE PROBLEMS SHAHID S. SIDDIQI 1, MUZAMMAL IFTIKHAR 2 arxiv:131.2915v1 [math.na] 1 Oct 213 Abstract. The
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationChapter 1 Numerical approximation of data : interpolation, least squares method
Chapter 1 Numerical approximation of data : interpolation, least squares method I. Motivation 1 Approximation of functions Evaluation of a function Which functions (f : R R) can be effectively evaluated
More informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More informationSolving Analytically Singular Sixth-Order Boundary Value Problems
Int. Journal of Math. Analysis, Vol. 3, 009, no. 39, 1945-195 Solving Analytically Singular Sixth-Order Boundary Value Problems Fazhan Geng 1 Department of Mathematics, Changshu Institute of Technology
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationEXP-FUNCTION METHOD FOR SOLVING HIGHER-ORDER BOUNDARY VALUE PROBLEMS
Bulletin of the Institute of Mathematics Academia Sinica (New Series) Vol. 4 (2009), No. 2, pp. 219-234 EXP-FUNCTION METHOD FOR SOLVING HIGHER-ORDER BOUNDARY VALUE PROBLEMS BY SYED TAUSEEF MOHYUD-DIN,
More informationGauss-Seidel Type Algorithms for a Class of Variational Inequalities
Filomat 32:2 2018, 395 407 https://doi.org/10.2298/fil1802395n Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Gauss-Seidel Type
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationarxiv: v1 [math.na] 20 May 2013
A numerical method based on the reproducing kernel Hilbert space method for the solution of fifth-order boundary-value problems arxiv:1305.4445v1 [math.na] 0 May 013 Mustafa Inc, Ali Akgül and Mehdi Dehghan
More informationTaylor Series. Math114. March 1, Department of Mathematics, University of Kentucky. Math114 Lecture 18 1/ 13
Taylor Series Math114 Department of Mathematics, University of Kentucky March 1, 2017 Math114 Lecture 18 1/ 13 Given a function, can we find a power series representation? Math114 Lecture 18 2/ 13 Given
More informationMulti-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media
Multi-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media Mohamed A. Hajji Abstract In this paper we propose a numerical scheme based on finite differences for the
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationComputational Non-Polynomial Spline Function for Solving Fractional Bagely-Torvik Equation
Math. Sci. Lett. 6, No. 1, 83-87 (2017) 83 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.18576/msl/060113 Computational Non-Polynomial Spline Function for Solving Fractional
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationSecond Order ODEs. CSCC51H- Numerical Approx, Int and ODEs p.130/177
Second Order ODEs Often physical or biological systems are best described by second or higher-order ODEs. That is, second or higher order derivatives appear in the mathematical model of the system. For
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More informationG-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS
Journal of Al-Nahrain University Vol.0(), December, 00, pp.8- Science G-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS Osama H.
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationFusion Higher -Order Parallel Splitting Methods for. Parabolic Partial Differential Equations
International Mathematical Forum, Vol. 7, 0, no. 3, 567 580 Fusion Higher -Order Parallel Splitting Methods for Parabolic Partial Differential Equations M. A. Rehman Department of Mathematics, University
More informationA parallel algorithm for the heat equation with derivative boundary conditions
International Mathematical Forum, 2, 2007, no. 12, 565-574 A parallel algorithm for the heat equation with derivative boundary conditions M. Akram University College of Information Technology University
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 informationContinuum mechanism: Plates
Observations of plate tectonics imply that the thin near-surface rocks, that constitute the lithosphere, are rigid, and therefore behave elastically on geological time scales. From the observed bending,
More informationRECURSIVE DIFFERENTIATION METHOD FOR BOUNDARY VALUE PROBLEMS: APPLICATION TO ANALYSIS OF A BEAM-COLUMN ON AN ELASTIC FOUNDATION
Journal of Theoretical and Applied Mechanics, Sofia, 2014, vol. 44, No. 2, pp. 57 70 RECURSIVE DIFFERENTIATION METHOD FOR BOUNDARY VALUE PROBLEMS: APPLICATION TO ANALYSIS OF A BEAM-COLUMN ON AN ELASTIC
More informationDifferential Equations Class Notes
Differential Equations Class Notes Dan Wysocki Spring 213 Contents 1 Introduction 2 2 Classification of Differential Equations 6 2.1 Linear vs. Non-Linear.................................. 7 2.2 Seperable
More informationNumerical solution of fourth order parabolic partial dierential equation using parametric septic splines
Hacettepe Journal of Mathematics and Statistics Volume 5 20, 07 082 Numerical solution of fourth order parabolic partial dierential equation using parametric septic splines Arshad Khan and Talat Sultana
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
More informationInterpolation Theory
Numerical Analysis Massoud Malek Interpolation Theory The concept of interpolation is to select a function P (x) from a given class of functions in such a way that the graph of y P (x) passes through the
More informationYou can learn more about the services offered by the teaching center by visiting
MAC 232 Exam 3 Review Spring 209 This review, produced by the Broward Teaching Center, contains a collection of questions which are representative of the type you may encounter on the exam. Other resources
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationAdvanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x
. Define f n, g n : [, ] R by f n (x) = Advanced Calculus Math 27B, Winter 25 Solutions: Final nx2 + n 2 x, g n(x) = n2 x 2 + n 2 x. 2 Show that the sequences (f n ), (g n ) converge pointwise on [, ],
More informationData Analysis-I. Interpolation. Soon-Hyung Yook. December 4, Soon-Hyung Yook Data Analysis-I December 4, / 1
Data Analysis-I Interpolation Soon-Hyung Yook December 4, 2015 Soon-Hyung Yook Data Analysis-I December 4, 2015 1 / 1 Table of Contents Soon-Hyung Yook Data Analysis-I December 4, 2015 2 / 1 Introduction
More informationAP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:
WORKSHEET: Series, Taylor Series AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period: 1 Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The
More informationLecture 34: Recall Defn: The n-th Taylor polynomial for a function f at a is: n f j (a) j! + f n (a)
Lecture 34: Recall Defn: The n-th Taylor polynomial for a function f at a is: n f j (a) P n (x) = (x a) j. j! j=0 = f(a)+(f (a))(x a)+(1/2)(f (a))(x a) 2 +(1/3!)(f (a))(x a) 3 +... + f n (a) (x a) n n!
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 informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationNumerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity
ANZIAM J. 46 (E) ppc46 C438, 005 C46 Numerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity Aliki D. Muradova (Received 9 November 004, revised April 005) Abstract
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationVARIATION OF PARAMETERS METHOD FOR SOLVING SIXTH-ORDER BOUNDARY VALUE PROBLEMS
Commun. Korean Math. Soc. 24 (29), No. 4, pp. 65 615 DOI 1.4134/CKMS.29.24.4.65 VARIATION OF PARAMETERS METHOD FOR SOLVING SIXTH-ORDER BOUNDARY VALUE PROBLEMS Syed Tauseef Mohyud-Din, Muhammad Aslam Noor,
More informationSection 0.2 & 0.3 Worksheet. Types of Functions
MATH 1142 NAME Section 0.2 & 0.3 Worksheet Types of Functions Now that we have discussed what functions are and some of their characteristics, we will explore different types of functions. Section 0.2
More informationA Short Essay on Variational Calculus
A Short Essay on Variational Calculus Keonwook Kang, Chris Weinberger and Wei Cai Department of Mechanical Engineering, Stanford University Stanford, CA 94305-4040 May 3, 2006 Contents 1 Definition of
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationComparison of Optimal Homotopy Asymptotic Method with Homotopy Perturbation Method of Twelfth Order Boundary Value Problems
Abstract Comparison of Optimal Homotopy Asymptotic Method with Homotopy Perturbation Method of Twelfth Order Boundary Value Problems MukeshGrover grover.mukesh@yahoo.com Department of Mathematics G.Z.S.C.E.T
More informationP-stable Exponentially fitted Obrechkoff methods
P-stable Exponentially fitted Obrechkoff methods Marnix Van Daele, G. Vanden Berghe Marnix.VanDaele@UGent.be Vakgroep Toegepaste Wiskunde en Informatica Universiteit Gent SciCADE 7 p. /4 Outline Introduction
More informationMathematical Economics: Lecture 2
Mathematical Economics: Lecture 2 Yu Ren WISE, Xiamen University September 25, 2012 Outline 1 Number Line The number line, origin (Figure 2.1 Page 11) Number Line Interval (a, b) = {x R 1 : a < x < b}
More informationand verify that it satisfies the differential equation:
MOTIVATION: Chapter One: Basic and Review Why study differential equations? Suppose we know how a certain quantity changes with time (for example, the temperature of coffee in a cup, the number of people
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More informationPolynomial Functions
Polynomial Functions Equations and Graphs Characteristics The Factor Theorem The Remainder Theorem http://www.purplemath.com/modules/polyends5.htm 1 A cross-section of a honeycomb has a pattern with one
More informationOn Newton-type methods with cubic convergence
Journal of Computational and Applied Mathematics 176 (2005) 425 432 www.elsevier.com/locate/cam On Newton-type methods with cubic convergence H.H.H. Homeier a,b, a Science + Computing Ag, IT Services Muenchen,
More informationNumerical Solution of Non-Linear Biharmonic. Equation for Elasto-Plastic Bending Plate
Applied Mathematical Sciences, Vol. 9, 5, no. 6, 769-78 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ams.5.575 Numerical Solution of Non-Linear Biharmonic Equation for Elasto-Plastic Bending Plate
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
More informationQuintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems
Stud. Univ. Babeş-Bolyai Math. 63(2018, No. 1, 141 151 DOI: 10.24193/subbmath.2018.1.09 Quintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems Ram
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
More informationCalculus I (Math 241) (In Progress)
Calculus I (Math 241) (In Progress) The following is a collection of Calculus I (Math 241) problems. Students may expect that their final exam is comprised, more or less, of one problem from each section,
More informationOn an iterative algorithm for variational inequalities in. Banach space
MATHEMATICAL COMMUNICATIONS 95 Math. Commun. 16(2011), 95 104. On an iterative algorithm for variational inequalities in Banach spaces Yonghong Yao 1, Muhammad Aslam Noor 2,, Khalida Inayat Noor 3 and
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More informationThe accurate numerical solution of the Schrödinger equation with an explicitly time-dependent Hamiltonian
The accurate numerical solution of the Schrödinger equation with an explicitly time-dependent Hamiltonian Marnix Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationAn elegant operational matrix based on harmonic numbers: Effective solutions for linear and nonlinear fourth-order two point boundary value problems
ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 2016, Vol. 21, No. 4, 448 464 http://dx.doi.org/10.15388/na.2016.4.2 An elegant operational matrix based on harmonic numbers: Effective solutions
More informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationNonlinear Thermo- Mechanics of Plates and Shallow Shells
Nonlinear Thermo- Mechanics of Plates and Shallow Shells Payam Khazaeinejad 1, Asif S. Usmani 1, Omar Laghrouche 1 IIE, School of Engineering, The University of Edinburgh IIE, School of the Built Environment,
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationMATH ASSIGNMENT 07 SOLUTIONS. 8.1 Following is census data showing the population of the US between 1900 and 2000:
MATH4414.01 ASSIGNMENT 07 SOLUTIONS 8.1 Following is census data showing the population of the US between 1900 and 2000: Years after 1900 Population in millions 0 76.0 20 105.7 40 131.7 60 179.3 80 226.5
More informationMath 581 Problem Set 5 Solutions
Math 581 Problem Set 5 Solutions 1. Show that the set { 2, 2 + i, 3 i} is linearly independent over Q. Proof: Suppose there exists a 0, a 1, and a 2 in Q so that a 0 2 + a1 ( 2 + i) + a 2 ( 3 i) = 0. Then
More informationNumerical solution of General Rosenau-RLW Equation using Quintic B-splines Collocation Method
Available online at www.ispacs.com/cna Volume 2012, Year 2012 Article ID cna-00129, 16 pages doi:10.5899/2012/cna-00129 Research Article Numerical solution of General Rosenau-RLW Equation using Quintic
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationBoundary-Value Problems for Ordinary Differential Equations
Chapter 11 Boundary-Value Problems for Ordinary Differential Equations 11.1 Introduction The differential equations in Chapter 5 are of first order and have one initial condition to satisfy. Later in the
More informationAhmed Abbas Mizeal and Mudhir A. Abdul Hussain
ARCHIVUM MATHEMATICUM (BRNO) Tomus 8 (01) 7 37 TWO-MODE BIFURCATION IN SOLUTION OF A PERTURBED NONLINEAR FOURTH ORDER DIFFERENTIAL EQUATION Ahmed Abbas Mizeal and Mudhir A. Abdul Hussain Abstract. In this
More informationAlgebraic Rational Cubic Spline with Constrained Control
Global Journal of Mathematical Sciences: Theory and Practical ISSN 0974-3200 Volume 6, Number 2 (2014), pp. 101-109 International Research Publication House http://www.irphouse.com Algebraic Rational Cubic
More information