Numerical Solutions to Partial Differential Equations
|
|
- Osborn McDowell
- 6 years ago
- Views:
Transcription
1 Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University
2 The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme The Crank-Nicolson scheme [ U m+1 U m = 1 U m +1 2U m + U 1 m τ 2 h 2 + Um+1 +1 ] 2Um+1 + U m+1 1 h 2 ; (1 + µ)u m+1 = (1 µ)u m + µ ( ) U 1 m + U+1 m + U m Um+1 +1 [ 1 T m+ 1 2 = 1 12 u ttt (x, η)τ 2 + u xxxx (ξ, t m+ 1 )h ]; 2 2 [ ] [ ] 2 λ k = 1 2 µ sin 2 kπ x 2 / µ sin 2 kπ x 2, thus unconditionally L 2 stable. 3 The computational cost is about twice that of the explicit scheme if solved by the Thompson method. 4 The maximum principle holds for µ 1. 5 Convergence rate is O(τ 2 + h 2 ), efficient if τ = O(h). 2 / 37
3 The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme The θ-scheme (0 < θ < 1, θ 1/2) U m+1 τ U m = (1 θ) Um +1 2Um + U m 1 h 2 (1+2µθ)U m+1 = (1 2µ(1 θ))u m +µ(1 θ) ( U m 1 + U m +1 + θ Um Um+1 + U m+1 1 h 2 ; ) +µθ (U m Um+1 +1 ). 1 T m+ 1 2 = O(τ 2 + h 4 ), if θ = [ [ 12µ ; = O(τ + h2 ), otherwise. ], thus ] / 2 λ k = 1 4(1 θ) µ sin 2 kπ x θ µ sin 2 kπ x 2 L 2 stable for 2µ (1 2θ) 1 (unconditional for θ 1/2). 3 The maximum principle holds for 2µ(1 θ) 1. 3 / 37
4 The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme The maximum principle and L stability and convergence Remark 1: For a finite difference scheme, L 2 stability conditions are generally weaker than L stability conditions. Remark 2: The maximum principle is only a sufficient condition for L stability. A finite difference scheme can be L stable, but does not satisfy the maximum principle. Remark 3: Under certain circumstances, for parabolic problems, L 2 stability and convergence can lead to L stability and convergence (see page 45-46, conclusion 2.6 for an example). 4 / 37
5 The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme The maximum principle and L stability and convergence Remark 4: In applications, it is safe to use L 2 stability if the data involved is smooth with high frequencies modes decay fast, otherwise, use the maximum principle conditions until it is such. Remark 5: In general applications, it should be safe to use the Crank-Nicolson scheme with µ = 1/2 initially, since then, for the high frequency modes, λ k π2 4 ( N k N )2, and the high frequency modes will decay fast in limited steps, and it will be safe to use arbitrary µ afterwards. Notice also that, λ N2 = 1 3 and 0 λ k < λ N2 if N 2 < k N. 5 / 37
6 Direct Difference Discretization The Variable-coefficient Linear Heat Equation Consider the variable-coefficient linear heat equation: u t = a(x, t)u xx + f (x, t), with the homogeneous Dirichlet boundary condition, where a(x, t) a 0 > 0 and f (x, t) are given real functions. The 1st order forward explicit scheme: U m+1 U m = a m U+1 m 2Um + U 1 m τ h 2, where a m = a(x, t m ), or equivalently U m+1 = (1 2µ m )U m + µ m ( U m +1 + U 1 m ), where µ m = a m τ/h 2 is the grid ratio. 6 / 37
7 Explicit Scheme for Variable-coefficient Heat Equation 1 The truncation operator of the explicit scheme is T (x, t) = [ +t δ 2 x ] t a(x, t) [ ( x) 2 t a(x, t) x 2 ]. 2 The local truncation error is Tu(x, t) = 1 2 u tt(x, η)τ a(x, t) 1 12 u xxxx(ξ, t)h 2 = O(τ + h 2 ). 3 Denote µ m = a m τ/h 2, the error equation is ( ) e m+1 = (1 2µ m )e m + µ m e+1 m + em 1 τt m. 4 The maximum principle holds on Ω tmax = [0, 1] [0, t max ], if µ(x, t) a(x, t) τ h 2 1 2, (x, t) Ω t max. 5 Variable-coefficient schemes do not have Fourier modes solutions, so, the Fourier analysis method does not apply. However, necessary conditions for L 2 stability can be derived by local Fourier analysis.
8 Direct Difference Discretization θ-scheme for Variable-coefficient Heat Equation [ 1 U m+1 = U m + µ m θδxu 2 m+1 + (1 θ)δxu 2 m µ m ] + τf m, where = µ m+θ, if θ = 0, 1; and µ m = µ m+ 1 2, if 0 < θ < 1. 2 The local truncation error is T m+ = { O(τ 2 + h 2 ), if θ = 1 2, O(τ + h 2 ), if θ The maximum principle holds on Ω tmax = [0, 1] [0, t max ], if µ(x, t) (1 θ) a(x, t) τ h 2 (1 θ) 1 2, (x, t) Ω t max. 8 / 37
9 θ-scheme for u t = (a(x, t)u x ) x in Conservation Form xr 1 The integral form of the equation x l [u(x, t a ) u(x, t b )] dx = ta t b [a(x r, t)u x (x r, t) a(x l, t)u x (x l, t)] dt, 0 x l < x r, t a > t b 0, 2 Take the control volume x l = x 1, x r = x 2 + 1, t b = t m, 2 t a = t m+1, and apply middle point quadrature rules; 3 The corresponding θ-scheme [ ( U m+1 = U m +θ µ m+ U m+1 +1 Um [ + (1 θ) µ m ( U m +1 U m ) ( ) ] µ m+ U m+1 1 U m ) ( µ m+ U m 1 U m ) ] 1. 2
10 Direct Difference Discretization Explicit Scheme for u t = a(u)u xx (with a( ) a 0 > 0) 1 The 1st order forward explicit difference scheme is given as U m+1 = U m + µ a(u m ) ( U+1 m 2U m + U 1) m, where µ = τ/h 2 is the coefficient independent grid ratio. 2 The maximum principle holds, if [ ] µ max a(u m ) 1 (x,t m) Ω tmax 2. 3 Substituting the real solution u into the scheme, we have u m+1 = u m + µ a(u m ) ( u+1 m 2u m + u 1 m ) + τ T m. where the local truncation error is again T m = O(τ + h 2 ). 10 / 37
11 Error Equation of Explicit Scheme in Nonlinear Case 4 Taylor expanding a(u m ) at U m leads to a(u m ) = a(u m ) + a (η m )(u m U m ) = a(u m ) a (η m ) e m. 5 The error equation e m+1 = e m + µ a(u m ) ( e+1 m 2e m + e 1 m ) τ T m + µ a (η m ) e m ( u m +1 2u m + u 1 m ). Nonlinear θ-scheme and its consistency and stability can also be established in a similar way. Note, because of the extra term in the error equation, the stability and consistency are not sufficient to guarantee convergence. Under certain additional conditions, the difference schemes for the nonlinear equations do converge, but the errors can often be significantly larger than in the linear case (see Exercise 2.16).
12 Difference Schemes Based on Semi-discretization Semi-discrete Methods of Parabolic Equations The idea of semi-discrete methods (or the method of lines) is to discretize the equation L(u) = f u t as if it is an elliptic equation, i.e. for a fixed t, we 1 introduce a spatial grid J Ω on Ω; 2 introduce spatial grid functions U h (t) = {U (t)} JΩ and f h = {f (t)} JΩ ; 3 approximate the differential operator L by a corresponding difference operator L h. 12 / 37
13 Difference Schemes Based on Semi-discretization Semi-discrete Methods of Parabolic Equations The procedure transforms an initial-boundary value problem of the parabolic partial differential equation into an initial problem of a system of ordinary differential equations on the grid nodes J Ω : U (t) = L h (U h (t)) + f (t), J Ω, U (0) = u(x, 0), J Ω. 13 / 37
14 Difference Schemes Based on Semi-discretization Convergence Analysis of Semi-discrete Methods For linear problems, the convergence analysis of semi-discrete methods is similar as that of the finite difference methods, i.e. 1 analyze the consistency of the approximation by calculating the local truncation error via Taylor series expansions; 2 analyze the stability of the semi-discrete system of the ordinary differential equations (the uniformly well-posedness of the semi-discrete problem). 14 / 37
15 Difference Schemes Based on Semi-discretization An Example of the Semi-discrete Method Consider the 1D heat equation on (0, 1) with homogeneous Dirichlet boundary condition. Let h = 1/N, on the spatial grid x = /N, = 0, 1,, N; substitute 2 / x 2 by δ 2 x/h 2. The corresponding initial problem of a system of ordinary differential equations: U (t) = N 2 [U +1 (t) 2 U (t) + U 1 (t)], t > 0, 1 N 1, U 0 (t) = U N (t) = 0, t > 0, U (0) = u(x, 0), 1 N / 37
16 Difference Schemes Based on Semi-discretization An Example of the Semi-discrete Method The initial problem of a system of ordinary differential equations: U (t) = N 2 [U +1 (t) 2 U (t) + U 1 (t)], t > 0, 1 N 1, The matrix form of the equation is U 0 (t) = U N (t) = 0, t > 0, U (0) = u(x, 0), 1 N 1. U h (t) = N2 A h U h (t), where A h is a tridiagonal, positive definite and symmetric matrix with diagonal elements 2, subdiagonal and superdiagonal nonzero elements 1. The eigenvalues and corresponding eigenvectors of A h are λ k = 4 sin 2 kπ 2N, V k = (V k ) = (sin kπ N ), 1 k N / 37
17 Difference Schemes Based on Semi-discretization Semi-discrete Method + ODE Solver Fully Discrete Scheme Discretize the ODE system obtained by the semi-discretization by 1 forward Euler method leads to the forward explicit scheme; 2 backward Euler method leads to the backward implicit scheme; 3 a linear combination of the two yields the θ-scheme; Note that the ODE system obtained is generally a stiff system, special care should be taken in choosing ODE solvers. Recommended solvers including Gear and Runge-Kutta methods. 17 / 37
18 General Boundary Conditions for 1D Problems The 3rd Type Boundary Condition for 1D Problems We consider at x = 0 the 3rd type boundary condition: u x (0, t) = α(t)u(0, t) + g 0 (t), α(t) 0, t > 0. The simplest finite difference approximation of the boundary condition: U1 m Um 0 = α m U0 m + g0 m, m 1, x or equivalently U m 0 = β m U m 1 β m g m 0 x, m 1, where β m = α m x. 18 / 37
19 General Boundary Conditions for 1D Problems The Effective Difference Scheme on the Node (x 1, t m ) By the Taylor series expansion, the local truncation error of the numerical boundary condition is ] u0 m = T m b = [ +x x x [ ] 1 m 2 xu xx + = O( x). 0 Since U0 m will be used to calculate Um 1, the truncation error at the boundary will spread into the interior nodes. 19 / 37
20 General Boundary Conditions for 1D Problems The Effective Difference Scheme on the Node (x 1, t m ) To see how this will affect the local truncation error there, we substitute the boundary condition U0 m = βm U1 m βm g0 m x into the 1st order forward explicit scheme to obtain the effective difference scheme on the node (x 1, t m ) U m+1 1 U m 1 t = Um 2 2Um 1 + (βm U m 1 βm g m 0 x) ( x) / 37
21 General Boundary Conditions for 1D Problems Truncation Error of the Effective Difference Scheme Remember g0 m = u x(0, t m ) α m u0 m and βm = (1 + α m x) 1, we have u2 m 2um 1 + (βm u1 m βm g0 m x) ( x) 2 u xx (x 1, t m ) [ u m = 2 2u1 m + ] [ um 0 ( x) 2 u xx (x 1, t m ) + βm u m 1 u m ] 0 u x (x 0, t m ) x x [ ] [ ] δ 2 = x u m1 + βm +x u0 m = ( x) 2 2 x [ ( x)2 u xxxx + ] m x 1 + βm x x x = 1 2 βm u xx (x 0, t m ) + O( x) = O(1). [ 1 2 xu xx + ] m 0 21 / 37
22 General Boundary Conditions for 1D Problems Truncation Error of the Effective Difference Scheme Thus the local truncation error T [ 1 m u1 m+1 u1 m τ = [ u m+1 1 u m 1 τ of the effective scheme is ] um 2 2um 1 + (βm u m 1 βm g m 0 h) h 2 um 2 2um 1 + um 0 h 2 (u t u xx )(x 1, t m ) [u t u xx ](x 1, t m ) ] = (βm u1 m βm g0 mh) um 0 [ h τu tt 1 ] m [ ] 12 h2 u xxxx + βm 1 m 1 h 2 hu xx + = O(1). 0 In other words, we have T m 1 = ˆT m 1 βm h T m b. 22 / 37
23 General Boundary Conditions for 1D Problems L Stability of the Effective Difference Scheme Compare the standard 1st order forward explicit scheme on (, m) U m+1 = (1 2µ)U m + µ(u m +1 + U m 1), > 1, and the effective scheme on (1, m) U m+1 1 = (1 µ(2 β m ))U m 1 + µu m 2 µβ m g m 0 h, we see that the conditions for the maximum principle to hold are 2µ 1 and (2 β m )µ 1. Since 0 < β m < 1, the condition for the maximum principle to hold is still µ 1/2. 23 / 37
24 Convergence Rate of the Effective Difference Scheme By taking a properly defined comparison function and applying the maximum principle, we can prove that the error of the overall scheme obtained by combining the 1st order forward explicit scheme with the 1st order forward approximation boundary condition is O(τ + h) (see Exercise 2.18). The overall convergence rate is eopardized by the numerical boundary condition, which has a local truncation error O(h). To decrease the local truncation error, especially to cope with the finite volume approximations, which create conservative schemes for conservative equations, we must consider other kinds of numerical boundary conditions.
25 General Boundary Conditions for 1D Problems A Numerical Boundary Condition with a Ghost Node 1 set h = x = ( N 1 2) 1, and let x = ( 1 2 )h, = 0, 1,, N, where x 0 = h/2 is a ghost node; 2 boundary condition: U m 1 Um 0 x = 1 2 αm (U m 1 + U m 0 ) + g m 0, or equivalently U m 0 = ξ m U m 1 η m g m 0 x, where ξ m = αm x αm x and ηm = αm x ; 3 local truncation error of the numerical boundary condition T m b = O(h2 ); 25 / 37
26 General Boundary Conditions for 1D Problems A Numerical Boundary Condition with a Ghost Node 4 the effective scheme: U m+1 1 U m 1 τ = Um 2 (2 ξm )U m 1 ηm g m 0 h h 2 ; 5 local truncation error of the effective scheme: T m 1 = ˆT m 1 ηm h T m b = O(τ + h); 6 the condition for the maximum principle: µ 1/2; 7 overall convergence rate can be shown to be O(τ + h 2 ). 26 / 37
27 General Boundary Conditions for 1D Problems Another Numerical Boundary Condition with a Ghost Node 1 set h = x = N 1, and let x = h, = 1, 0,, N, where x 1 = h is a ghost node, [x 1, x 0 ] is a ghost cell; 2 boundary condition: U1 m Um 1 2 x = α m U m 0 + g m 0, or equivalently U m 1 = U m 1 2 α m U m 0 x 2 g m 0 x; 3 local truncation error of the numerical boundary condition T m b = O(h2 ); 27 / 37
28 General Boundary Conditions for 1D Problems Another Numerical Boundary Condition with a Ghost Node 4 the effective scheme: U0 m+1 U0 m = 2 Um 1 2 (1 + αm h) U0 m 2 g 0 mh τ h 2 ; 5 local truncation error of the effective scheme: T m 0 = ˆT m 0 2 h T m b = O(τ + h); 6 the condition for the maximum principle: (1 + α m h)µ 1/2. 7 overall convergence rate can be shown to be O(τ + h 2 ). 28 / 37
29 Dissipation and Conservation Dissipation and Stability of Difference Schemes For the heat equation u t = u xx defined on(0, 1), 1 analytical Fourier mode solutions e k2 π 2t e ikπx decay fast; 2 discrete Fourier modes λ m k eikπh ; 3 for 1st order forward explicit scheme, λ k = 1 4 µ sin 2 kπh 2, h = 1 N, k = 1, 2,, N, all Fourier modes decay if µ < 1/2; 4 for any N 1, λ N = 1 4 µ = 1, if µ = 1/2, the corresponding Fourier mode does not decay (nor grow). 29 / 37
30 Dissipation and Conservation Dissipation and Stability of Difference Schemes 5 lim h 0 λ m k ek2 π 2 mτ = lim h 0 (1 k 2 π 2 τ) m e k2 π 2 mτ = 1, 0 < mτ t max and any fixed k, convergence of mode k; 6 for any N 1, λ N = 1 4 µ < 1, if µ > 1/2, unstable; 7 a necessary and sufficient condition for L 2 stability: µ 1/2. 30 / 37
31 Dissipation and Conservation Local L 2 Stability of Variable-coefficient Difference Schemes For variable-coefficient heat equation u t = a(x)u xx defined on (0, 1), we can apply the Fourier method to study the necessary stability conditions of a scheme by analyzing its local dissipation property. 1 discrete Fourier modes λ m k eikπh ; 2 for the forward explicit scheme U m+1 = U m + µ a(ξ) ( U m +1 2U m + U m 1), = 1, 2,, N. λ k = 1 4 µ a(ξ) sin 2 kπh 2, h = 1 N, k = 1, 2,, N, all Fourier modes decay, if µ max ξ [0,1] a(ξ) < 1/2; 3 for any N 1, if µ max ξ [0,1] a(ξ) > 1/2, then λ N < 1, unstable; 31 / 37
32 Parabolic Equation in Conservation Form Parabolic equations with the differential operator L being of the divergence form can be written into a conservative integral form: xr x l u(x, t a ) dx = xr x l In a heat flow problem, u(x, t b ) dx+ ta t b [a u x (x r, t) a u x (x l, t)] dt, 0 x l < x r 1, t a > t b 0. 1 h([x l, x r ]; t) x r x l u(x, t) dx is the total heat of the system on the interval (x l, x r ) at time t; 2 f (u) a u x (x, t) is the heat flux on x at time t; 3 u x (0, t) = g 0 (t), u x (1, t) = g 1 (t) (the Neumann boundary condition) given inflow heat rate at the boundary; Remark: The equation describes the local heat conservation (or balance) property of a source free heat flow problem.
33 Dissipation and Conservation Finite Volume Method and Conservative Schemes We would like to establish a scheme which preserves the property. 1 U m : the heat density of the discrete system on the grid cell (x 1, x ) at time t m ; 2 2 H([x l, x r ]; t m+k ) r h: the total heat of the discrete system on the interval (x l, x r ) at time t m+k ; = l U m+k 3 F (U m+1, U m, U m+1 1, Um 1 ): the numerical heat flux on x 1 2 the time period (t m, t m+1 ); in 33 / 37
34 Dissipation and Conservation Finite Volume Method and Conservative Schemes 4 a general conservative finite volume scheme U m+1 = U m τ [ ] F (U m+1 +1 h, Um +1, U m+1, U m ) F (U m+1, U m, U m+1 1, Um 1), 5 local heat conservation of the discrete system r U m+k h = = l +F (U m+l l +1, Um+l 1 l +1 r k [ U m h+ F (U m+l r +1, Um+l 1 r +1, U m+l r, U m+l 1 r ) = l l=1 ], U m+l, U m+l 1 l ) τ, 0 l < r < N, m 0, k > 0. l Remark: To preserve the global conservation of the system, the initial and boundary data must also be properly handled. 34 / 37
35 Dissipation and Conservation Numerical Initial and Boundary Data for Conservative Schemes 1 a grid with ghost nodes x = ( 1 2 )h, 0 N, h = 1/(N 1), t m = mτ, m 0, τ = µ(h 2 ); 2 F (U m+1, U m, U m+1 1, Um 1 ) = a Um U m 1 h 3 numerical Neumann boundary conditions: U1 m Um 0 = ḡ m UN m 0, Um N 1 = ḡ1 m, h h where ḡ m 0 = 1 τ tm+1 t m g 0 (t) dt, ḡ m 1 = 1 τ explicit scheme; tm+1 t m g 1 (t) dt, 4 global numerical boundary flux equals global boundary flux [( m τ a Ul N ) ( Ul N 1 a Ul 1 ) ] Ul tm+1 0 = a(g 1 (t) g 0 (t))dt; h h l= / 37
36 Dissipation and Conservation Numerical Initial and Boundary Data for Conservative Schemes 5 numerical initial data: U 0 = 1 h x the initial heat of the system: H([0, 1]; 0) = N 1 =1 x 1 2 u(x, 0) dx, = 1,, N 1; U 0 h = N 1 =1 7 the heat of the system at t m+1 : H([0, 1]; t m+1 ) = h([0, 1]; 0) + x+ 1 2 x 1 2 u(x, 0) dx = h([0, 1]; 0); tm+1 0 a (g 1 (t) g 0 (t)) dt = h([0, 1]; t m+1 ), m / 37
37 SK 2µ18, 19; þåš 2 Thank You!
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
More informationFinite difference method for heat equation
Finite difference method for heat equation 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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationImplicit Scheme for the Heat Equation
Implicit Scheme for the Heat Equation Implicit scheme for the one-dimensional heat equation Once again we consider the one-dimensional heat equation where we seek a u(x, t) satisfying u t = νu xx + f(x,
More informationIntroduction to standard and non-standard Numerical Methods
Introduction to standard and non-standard Numerical Methods Dr. Mountaga LAM AMS : African Mathematic School 2018 May 23, 2018 One-step methods Runge-Kutta Methods Nonstandard Finite Difference Scheme
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
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 informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
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 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 informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
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 informationNumerical Integration of Linear and Nonlinear Wave Equations
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations, Theses, and Student Research Papers in Mathematics Mathematics, Department of 12-2004 Numerical Integration
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationTMA4212 Numerical solution of partial differential equations with finite difference methods
TMA4 Numerical solution of partial differential equations with finite difference methods Brynjulf Owren January 3, 07 Translated and amended by E. Celledoni Preface The preparation of this note started
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationChapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0
Chapter 4. MAC Scheme 4.1. MAC Scheme and the Staggered Grid. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes equation in the velocity-pressure formulation: (4.1.1)
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationWorksheet 8 Sample Solutions
Technische Universität München WS 2016/17 Lehrstuhl für Informatik V Scientific Computing Univ.-Prof. Dr. M. Bader 19.12.2016/21.12.2016 M.Sc. S. Seckler, M.Sc. D. Jarema Worksheet 8 Sample Solutions Ordinary
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 informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few
More informationNumerical PDE. DYKwak. September 6, 2007
Numerical PDE DYKwak September 6, 2007 ii Contents 1 Finite Difference Method 7 1.1 2nd order linear p.d.e. in two variables................ 7 1.2 Finite difference operator....................... 8 1.2.1
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationSelected HW Solutions
Selected HW Solutions HW1 1 & See web page notes Derivative Approximations. For example: df f i+1 f i 1 = dx h i 1 f i + hf i + h h f i + h3 6 f i + f i + h 6 f i + 3 a realmax 17 1.7014 10 38 b realmin
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
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 informationNotes on Numerical Analysis
Notes on Numerical Analysis Alejandro Cantarero This set of notes covers topics that most commonly show up on the Numerical Analysis qualifying exam in the Mathematics department at UCLA. Each section
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
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 informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More information