Chapter 9 Collocation Methods

110 downloads 330 Views 572KB Size Report
Chapter 9. Collocation Methods. 9.1 Introduction. Let's continue the discussion of collocation methods with the second-order linear BVP. Ly = y. 00. + p(x)y. 0.
Chapter 9 Collocation Methods 9.1 Introduction Let's continue the discussion of collocation methods with the second-order linear BVP

Ly = y + p(x)y + q(x)y = r(x) 00

a J convergence at the nodes and collocation points is at a higher rate than implied by Theorem 9.3.1. This is the phenomenon of superconvergence that we illustrated in Section 9.1. Proof. As in Section 9.1, introduce the Green's function

e( ) =

Z

b

a

XZ

G( x)Le(x)dx =

J

j

=1

xj

xj

G( x)Le(x)dx:

;1

(9.3.14)

Let us assume that  2= (x 1 x ) and, following the logic introduced in Section 9.1, write j;

j

Y J

G( x)Le(x) = w(x) (x ;  ) j

ij

=1

The function w(x) involves the J th derivative of Le. Using this and (9.3.7b), we may expect w(x) to have p ; J bounded derivatives. Thus, expand w(x) in a Taylor's series of the form w(x) = P 1(x) + O(h ) p;J

p;J ;

i

where P 1(x) is a polynomial of degree p ; J ; 1. If the one-step method is accurate to order p then p;J ;

Z

xj

xj

Since

;1

P

p; J ;

J

j

=1

(x ;  )dx = 0: ij

Y J

j

we have

1(x)

Y

Z

xj

xj

;1

=1

(x ;  ) = O(h ) J

ij

G( x)Le(x)dx = O(h

P ;J i

i

)O(h )h J i

i

 2= (x

1

j;

x ): j

(9.3.15a)

The result (9.3.13a) is obtained by summing the above relation over the subintervals. When  2 (x 1 x ) then we are only able to show that

Z

j;

j

xj

xj

;1

G( x)Le(x)dx = O(h +1) J i

Summing (9.3.15a,b) yields (9.3.13b). Superconvergence occurs whenever p > J + 1. 24

 2 (x

1

j;

x ): j

(9.3.15b)

Bibliography

1] U.M. Ascher, R. Mattheij, and R. Russell. Numerical Solution of Boundary Value Problems for Ordinary Di erential Equations. SIAM, Philadelphia, second edition, 1995.

2] C. de Boor. A Practical Guide to Splines. Springer-Verlag, New York, 1978.

3] C. de Boor and B. Swartz. Collocation at gaussian points. SIAM J. Numer. Anal., 10:582{687, 1973.

4] J.E. Flaherty and W. Mathon. Collocation with polynomial and tension splines for singularly perturbed boundary value problems. SIAM J. Sci. Stat. Comput, 1:260{ 289, 1980.

5] E. Isaacson and H.B. Keller. Analysis of Numerical Methods. John Wiley and Sons, New York, 1966.

25