From feb43795b2be471912d34c6f40570a43f9b91b17 Mon Sep 17 00:00:00 2001
From: rlar <rlar>
Date: Tue, 8 Oct 2013 18:38:43 +0200
Subject: [PATCH] numparm, fix pow() implementation

---
 src/frontend/numparam/xpressn.c | 7 +------
 1 file changed, 1 insertion(+), 6 deletions(-)

diff --git a/src/frontend/numparam/xpressn.c b/src/frontend/numparam/xpressn.c
index b62965f68..350da5b49 100644
--- a/src/frontend/numparam/xpressn.c
+++ b/src/frontend/numparam/xpressn.c
@@ -1038,7 +1038,6 @@ operate(char op, double x, double y)
     /* bug:   x:=x op y or simply x:=y for empty op?  No error signalling! */
     double u = 1.0;
     double z = 0.0;
-    double epsi = 1e-30;
     double t;
 
     switch (op)
@@ -1060,11 +1059,7 @@ operate(char op, double x, double y)
         x = x / y;
         break;
     case '^':                   /* power */
-        t = absf(x);
-        if (t < epsi)
-            x = z;
-        else
-            x = exp(y * ln(t));
+        x = pow(fabs(x), y);
         break;
     case '&':                   /* && */
         if (y < x)