diff --git a/hfractal.cabal b/hfractal.cabal
--- a/hfractal.cabal
+++ b/hfractal.cabal
@@ -1,5 +1,5 @@
 name:                hfractal
-version:             0.4.2.1
+version:             0.4.2.2
 cabal-version:       >= 1.2
 synopsis:            OpenGL fractal renderer
 description:         An OpenGL fractal browser with multicore support and the capability to output high quality png images.
@@ -16,6 +16,6 @@
 Executable hfractal
     main-is:            Hfractal.hs
     build-depends:      base >=3 && <5, array, gd < 3000.3.0 && >= 3000.2.0, data-accessor >=0.2, data-accessor-template >=0.2, OpenGL >= 2.3, OpenGLRaw < 1.1.0.0, GLUT >= 2.2.1.0, colour >=2.3.1, containers >= 0.2
-    other-Modules:      FracComp, FracColour, FracImg, FracState, Bindings
+    other-Modules:      FracComp, FracColour, FracImg, FracState, Bindings, PointComp
     ghc-options:        -O2 -threaded -fvia-C -optc-O3 -optc-ffast-math
     Hs-Source-Dirs:     src
diff --git a/src/PointComp.hs b/src/PointComp.hs
new file mode 100644
--- /dev/null
+++ b/src/PointComp.hs
@@ -0,0 +1,17 @@
+{-# LANGUAGE BangPatterns #-}
+module PointComp
+	where
+{- Module for computing fractal functions at an indivdual point -}
+
+-- Number of iterations to escape
+mandPoint :: Double -> Double -> Int -> Double
+mandPoint cx cy mi | inCardiod cx cy = 0.0
+				    | otherwise      = go 0 0.0 0.0 where
+	inCardiod cx cy = q * (q + (cx - 0.25)) < 0.25 * cy * cy where
+		q = (cx - 0.25) * (cx - 0.25) + cy*cy
+	go !n !x !y | x2 + y2 > 4.0    = 1.0 - logBase 2 (0.5 * logBase 2 (x2 + y2)) + fromIntegral n
+				| n > mi		   = 0.0
+				| otherwise		   = go (n+1) (x2 - y2 + cx) (2.0*x*y + cy) where
+		!x2 = x*x 
+		!y2 = y*y  --This CSE saves a few cycles
+