From 05861223c8e1d6aadec146a47a59140695380680 Mon Sep 17 00:00:00 2001
From: Ramakrishnan Muthukrishnan <vu3rdd@gmail.com>
Date: Tue, 4 May 2010 23:38:04 +0530
Subject: [PATCH] solution for exercise 1.31

---
 src/sicp/ex1_31.clj | 43 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)
 create mode 100644 src/sicp/ex1_31.clj

diff --git a/src/sicp/ex1_31.clj b/src/sicp/ex1_31.clj
new file mode 100644
index 0000000..a1c8dd7
--- /dev/null
+++ b/src/sicp/ex1_31.clj
@@ -0,0 +1,43 @@
+(ns sicp.ex1_31
+  (:use [clojure.contrib test-is]
+	[sicp utils]))
+
+(defn product [a b]
+  (if (> a b)
+    1
+    (* a (product (+ a 1) b))))
+
+(defn prod [term a next b]
+  (if (> a b)
+    1
+    (* (term a)
+       (prod term (next a) next b))))
+
+;; factorial
+;; fact (n) = n * factorial (n - 1)
+
+(defn factorial [n]
+  (prod identity 1 inc n))
+
+(deftest test-factorial-of-6
+  (is (factorial 6) 720))
+
+(deftest test-factorial-of-10
+  (is (factorial 10) 3628800))
+
+(defn prod-pi [a b]
+  (prod #(/ (* % (+ % 2)) (square (+ % 1.0))) a #(+ % 2) b))
+
+(deftest test-pi-2-digits-accuracy
+  (is (> (Math/abs (- 3.14159 (* 4 (prod-pi 2 1000)))) 0.001) true))
+
+(defn iter-prod [term a next b result]
+  (if (> a b)
+    result
+    (iter-prod term (next a) next b (* a result))))
+
+(defn iprod [term a next b]
+  (iter-prod term a next b 1))
+
+(deftest test-iprod-for-prod-of-1-to-10
+  (is (iprod identity 1 inc 10) (prod identity 1 inc 10)))
-- 
2.45.2