solution of equations by divided difference interpolaton
Divided diference interpolation can be used to obtain approximate solutions to equations and to invert functions numerically. The idea is that, given an equation which we want to solve for , we first take several numbers and compute as . Then we compute the divided differences of the ’s regarded as functions of the ’s and form the divided difference series. Substituting in this series provides an approximation to .
To illustrate how this works, we will examine the transcendental equation . We note that and , so there will be a solution between and , likely closer to than . Therefore, as our values of the ’s, we shall take , , , , , , . We now tabulate for those values:
Next, we form a divided difference table of the ’s as a function of the ’s:
From this table, we form the series
Substituting for , we obtain . Given that
this answer is correct to all decimal places.
In the presentation above, we tacitly assumed that there was a solution to our equation and focussed our attention on finding that answer numerically. To complete the treatment we will now show that there indeed exists a unique solution to the equation in the interval .
Since depends continuously on , it follows that there exists such that .
As for uniqueness, note that the derivative of is . When , we have , or . Hence, is a strictly increaing function of , so there can be at most one such that .
|Title||solution of equations by divided difference interpolaton|
|Date of creation||2013-03-22 16:49:22|
|Last modified on||2013-03-22 16:49:22|
|Last modified by||rspuzio (6075)|