22.07.2014 Views

Lecture 4 - Computer Science Department - University of Kentucky

Lecture 4 - Computer Science Department - University of Kentucky

Lecture 4 - Computer Science Department - University of Kentucky

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

On the other hand, define<br />

We expand F(a+2h) as<br />

Basic Simpson’s Rule<br />

F (<br />

x<br />

)<br />

= x f<br />

(<br />

t<br />

)<br />

dt<br />

a<br />

∫<br />

5<br />

2 4<br />

3 2<br />

4<br />

(4) 2<br />

5<br />

(5)<br />

F ( a<br />

+ 2 h<br />

) =<br />

F ( a<br />

)<br />

+<br />

2<br />

hF '( a<br />

)<br />

+<br />

2<br />

h F<br />

"( a<br />

)<br />

+<br />

h F<br />

"'( a<br />

)<br />

+<br />

h<br />

f<br />

( a<br />

)<br />

+<br />

h F<br />

( a<br />

)<br />

+L<br />

3 3 5!<br />

Note that F’=f, F(a)=0, F”=f’, F”’=f”, we have<br />

∫ + 5<br />

a 2h<br />

2 4 3 2 4 2 5 (4)<br />

f(<br />

x)<br />

dx=<br />

2hf(<br />

a)<br />

+ 2h<br />

f'(<br />

a)<br />

+ h f"(<br />

a)<br />

+ h f"'(<br />

a)<br />

+ h f ( a)<br />

+L<br />

a<br />

3 3 5⋅4!<br />

From the previous page, we have<br />

h<br />

[ f ( a)<br />

+ 4 f ( a + h)<br />

+ f ( a + 2h)]<br />

=<br />

3<br />

2 4 3<br />

2hf<br />

( a)<br />

+ 2h<br />

f '(<br />

a)<br />

+ h f "( a)<br />

3<br />

+ 2 4 20 5 (4<br />

h f "'(<br />

a)<br />

+ h f<br />

) ( a)<br />

+L<br />

3 3⋅<br />

4!<br />

36

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!