Lecture 4 - Computer Science Department - University of Kentucky
Lecture 4 - Computer Science Department - University of Kentucky
Lecture 4 - Computer Science Department - University of Kentucky
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