ROMBERG ALGORITHM

代码说明：

*   *   To approximate I = integral ( ( f(x) dx ) ) from a to b:*   *   INPUT:   endpoints a, b; integer n.**   OUTPUT:  an array R. ( R(2,n) is the approximation to I. )**   R is computed by rows; only 2 rows saved in storage