There is an excellent number theoretical BASIC called UBASIC. It can handle arbitrary precision numbers and so it is ideally suited for approximations and number theoretic manipulations like factoring large numbers. If you cannot find it on the net, you can download it here.
I have ported a MINIMAX VIA ERROR FLUFFING ALGORITHM by Steven A. Ruzinsky to UBASIC. This program demonstrated the determination of polynomial coefficients to a minimum maximum absolute error criterion. The UBASIC version, by virtue of the larger precision, is more practical and useful.
Using UBASIC I have made
approximate functions that are known analytically,
by calculating the true Chebychev coefficients directly from its
definition. As everybody knows, a Chebychev approximation is a
polynomial with the property that positive and negative deviations
alternate and all the deviations are equal. This program makes those
deviations equal to a certain precision specified by the user. 20 bits
is easily obtainable if the order of the approximation is not too
large. If you can't run it you may want to look at the
NOTE: Chebychev approximations are no longer the best choice for very high precision approximation. Use rational approximations, i.e. quotients of polynomials, instead.
For the testing of function approximations I have made a toolset that generates tables of functions with a given interval and a given precision. Again this requires UBASIC. The forth program used for testing the floating point of tForth is supplied as an example. The tables with extension .ftb are interpreted as programs in Forth. These tables are examples. In a real test tables are much denser and occupy many megabytes. You have to create them yourself using UBASIC.
Go to the home page of Albert van der Horst