1 files changed, 41 insertions, 0 deletions
diff --git a/mdk-stage1/dietlibc/libm/poly.c b/mdk-stage1/dietlibc/libm/poly.c
new file mode 100644
index 000000000..cdcfb8c5c
--- /dev/null
+++ b/mdk-stage1/dietlibc/libm/poly.c
@@ -0,0 +1,41 @@
+/*--------------------------------------------------------------------------*
+
+Name		__poly - generates a polynomial from arguments
+
+Usage		double  __poly ( double x, int n, const double* c );
+
+Prototype in	math.h
+
+Description	__poly generates a polynomial in x, of degree n, with
+		coefficients c[0], c[1], ..., c[n]. For example, if n=4,
+		the generated polynomial is
+
+			c[4]*x^4 + c[3]*x^3 + c[2]*x^2 + c[1]*x + c[0]
+
+		The polynomial is calculated using Horner's method:
+
+			polynom = (..((x*c[n] + c[n-1])*x + c[n-2])..)*x + c[0]
+
+Return value	__poly returns the value of the polynomial as evaluated for
+		the given x.
+		A range error occurs if the result exceeds double range.
+
+*---------------------------------------------------------------------------*/
+
+#include <stdio.h>
+#include "dietlibm.h"
+
+double  __poly ( double x, size_t n, const double* c) 
+{
+    long double ret;
+    size_t      i;
+    
+    i   = n;
+    c  += n;
+    ret = 0;
+    do
+        ret = ret * x + *c--;
+    while ( i-- );
+    
+    return ret;
+}