1 files changed, 56 insertions, 1 deletions

diff --git a/AufgabeFFP2.hs b/AufgabeFFP2.hs
index e66b7a8..7993c01 100644
--- a/AufgabeFFP2.hs
+++ b/AufgabeFFP2.hs
@@ -38,4 +38,59 @@ powFast :: Int -> Integer
 powFast 0 = 1
 powFast n = pof2s !! (n-1) + pof2s !! (n-1)
 
--------------------------------------------------------------------------------
\ No newline at end of file
+-------------------------------------------------------------------------------
+
+-- 3
+
+-- power series of N
+pofNs :: Integer -> [Integer]
+pofNs n = [1] ++ (map (n*) $ pofNs n)
+
+-- faculty function
+fac :: Integer -> Integer
+fac n = product [1..n]
+
+-- stream of faculties
+facGen :: [Integer]
+facGen = map (fac) [1..]
+
+-- function g with memoization (using hFast)
+fMT :: Int -> Int -> Float
+fMT z k = g z k hMT
+
+-- function g without memoization (uning hSlow)
+f :: Int -> Int -> Float
+f z k = g z k h
+
+-- actual function g (converts Int to Integer for more precision)
+g :: Int -> Int -> (Integer -> Integer -> Float) -> Float
+g z k h = sum $ map (h $ fromIntegral z) [1..(fromIntegral k)]
+        
+-- helper function h using mem-table for the power-series (z^i) and for faculty (i!)
+hMT :: Integer -> Integer -> Float
+hMT z i = (fromInteger $ pofNs z !! (fromInteger i)) / (fromInteger $ facGen !! ((fromIntegral i)-1))
+
+-- helper function h without memoization
+h :: Integer -> Integer -> Float
+h z i =  (fromInteger $ z^i) / (fromInteger $ fac i)
+
+-------------------------------------------------------------------------------
+
+-- 4
+
+-- gets the digits of an integer as a list
+digits :: Integer -> [Integer]
+digits x
+        | x<=0 = []
+        | otherwise = (digits $ x `div` 10)++[x `mod` 10]
+
+-- calculates the goedel-number for the given integer
+-- returns 0 for non-positive numbers
+gz :: Integer -> Integer
+gz n
+        | n<=0 = 0
+        | otherwise = product $ zipWith (^) primes (digits n)
+
+-- goedel-number generator
+gzs :: [Integer]
+gzs = map (gz) [1..]