25c with l 3 and 327 with k 3 ie b1 c2 b2 c2

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 1) 3 3 Y2 = yn 1 + h f (tn 1 + h Y1) + f (tn Y2)] 2 3 yn = yn 1 + h 3f (tn 1 + h Y1) + f (tn Y2)]: 4 3 We can check by constructing a Taylor's series that this method is indeed third order. Hairer et al. 16], Section II.2, additionally show that our necessary conditions for thirdorder accuracy are also su cient in this case. The computation of Y1 can be recognized as the backward Euler method for one-third of the time step h. The computation of Y2 and yn are not recognizable in terms of simple quadrature rules. Since the method is third-order, its local error is O(h4). We can also construct an SDIRK method by insisting that a11 = a22 . Enforcing this condition and using the previous relations gives two methods having the tableau ; ; ; ; ; ; 1; 1;2 1=2 where 0 1=2 1 = 1 (1 p ): 2 3 p The method with = (1 + 1= 3)=2 is A-stable while the other method has a bounded stability region. Thus, this would be the method of choice. 28 Let us conclude this Section by noting a relationship between implicit R...
View Full Document

This document was uploaded on 03/16/2014 for the course CSCI 6820 at Rensselaer Polytechnic Institute.

Ask a homework question - tutors are online