Use Richardson Extrapolation to compute the function

320 Views Asked by At

A regular n sided polygon inscribed in a circle of diameter one has the circumference as

cn = 2 n Sin( Π / n )

put h=1/n

How can i show whether c(h) = c1/h satisfies the assumption for repeated Richardson extrapolation

Also i need to calculate c2 , c3 , c6 (which i can do ) and perform repeated Richardson Extrapolation

My knowledge in Richardson Extrapolation is very bad and i have got no idea how to do it

So please help

1

There are 1 best solutions below

4
On

Just use Taylor's formula: $$ C(h)=\frac 2h \sin (\pi h) = \frac 2h \left((\pi h)- \frac{(\pi h)^3}{3!} + O(h^5) \right) = 2 \pi + \frac{2 \pi^3}{3!} h^2 + O(h^4) $$

In general, when you are approximating a quantity $C^*$ by a formula like $$ C^* = C(h) + a_0 h^{k_0} + O(h^{k_1}), $$

the extrapolation based on $C(h)$ and $C(h/t)$ is given by $$ R(h,t)=\dfrac{t^{k_0} C(h/t)-C(h)}{t^{k_0}-1} $$

and in fact $$ C^* = R(h,t)+O(h^{k_1}) $$

In your case, the extrapolation based on $c_3, c_6$ is given by $$ C^* \approx \dfrac{2^2 c_6-c_3}{2^2-1}=\frac 43 c_6-\frac 13 c_3 $$

You can see the advantage of this procedure just by computing the errors: $$ |2 \pi -C(1/3)|\approx 1.09, \quad |2 \pi -C(1/6)|\approx 0.28, \quad |2\pi - R(1/3,2)| \approx 0.015 $$