WIP julia

Julia.fs

@@ -6,12 +6,15 @@ open System.Drawing
 
 let generatePixel x y height width zoom rMax maxIter = 
     let c = MandelBrot.getComplexValueForPixel (float x) (float y) (float height) (float width) zoom
+    let z = Complex(-0.8, 0.156)
 
     let rec findI i z =
         if i >= maxIter then i
         else if Complex.Abs(z) >= (float rMax) then i
-        //else findI (i+1) (z*z+c)
-        else findI (i+1) (Complex.Add(c, (Complex.Multiply(z, z))))
+        else
+            let newRe = c.Real*c.Real - c.Imaginary*c.Imaginary
+            let newIm = 2.0 * c.Real * c.Imaginary
+            findI (i+1) (Complex(newRe, newIm)+z)
 
     MandelBrot.colorMap (findI 0 (Complex(0.0, 0.0))) maxIter
 

MandelBrot.fs

@@ -66,8 +66,7 @@ let generatePixel x y height width zoom rMax maxIter =
     let rec findI i z =
         if i >= maxIter then i
         else if Complex.Abs(z) >= (float rMax) then i
-        //else findI (i+1) (z*z+c)
-        else findI (i+1) (Complex.Add(c, (Complex.Multiply(z, z))))
+        else findI (i+1) (c+z*z)
 
     colorMap (findI 0 (Complex(0.0, 0.0))) maxIter
 

Program.fs

@@ -6,8 +6,10 @@ open System
 let main argv =
     printfn "Mandelbrot!"
 
-    let result = MandelBrot.generate 800 600 150.0 32 570
+    let resultMandelbrot = MandelBrot.generate 800 600 150.0 32 570
+    resultMandelbrot.Save("mandelbrot.bmp")
 
-    result.Save("mandelbrot.bmp")
+    let resultJulia = Julia.generate 800 600 150.0 32 570
+    resultJulia.Save("julia.bmp")
 
     0 // return an integer exit code