which is the difference formula for the third derivative in terms of the function at
the four points situated at unequal intervals. The corresponding truncation error
can be obtained by differentiating equation (19). Putting
1 2 3
h = h = h = h in above
expressions, we obtained the results for equal sub-intervals as obtained by Singh
and Thorpe [9].