Based on this:
http://www.math.unipd.it/~alvise/CS_2008/APPROSSIMAZIONE_2009/MFILES/Remez.pdf
Please find attached remez.mcd and remeztest (simple test harness).
Not exhaustively tested but seems to work for me although I do get some errors with high order polynomial approximation.
Parameters:
F: function to be approximated
a: low limit of x range over which F is to be approximated
b: high limit ditto
n: order of polynomial
convergence limit: e.g. if 0.01 will stop iterating when error in approximation is within 1% of error in previous iteration
returns approximating polynomial coefficients c0, c1, ..., cn, maximum error in approximation
Any comments, improvements welcome.
Dave W.
