Added example exponential code.

1 files changed, 23 insertions, 0 deletions

diff --git a/ccc2012/exp.py b/ccc2012/exp.py
new file mode 100644
index 0000000..39b5035
--- /dev/null
+++ b/ccc2012/exp.py
@@ -0,0 +1,23 @@
+from math import factorial
+
+def exp_right(x):
+    ret = 1
+    negative = (x < 0)
+
+    xa = abs(x)
+
+    for i in range(1, 128):
+        ret = ret + xa**float(i) / float(factorial(i))
+
+    if negative:
+        ret = 1.0/ret
+
+    return ret
+
+def exp_wrong(x):
+    ret = 1
+
+    for i in range(1, 128):
+        ret = ret + x**float(i) / float(factorial(i))
+
+    return ret