Создание фрактала Ньютона

import javafx.fxml.FXML;
import javafx.fxml.Initializable;
import javafx.scene.canvas.Canvas;
import javafx.scene.canvas.GraphicsContext;
import javafx.scene.layout.StackPane;
import javafx.scene.paint.Color;
 
import java.net.URL;
import java.util.ResourceBundle;
 
import static java.lang.Math.*;
 
public class MainController implements Initializable {
    @FXML private StackPane stackPane;
    @FXML private Canvas canvas;
    private GraphicsContext graphicsContext;
    private double x0 = -1.0;
    private double xn = 1.0;
    private double y0 = -1.0;
    private double yn = 1.0;
 
    public void initialize(URL location, ResourceBundle resources) {
        graphicsContext = canvas.getGraphicsContext2D();
        stackPane.setStyle("-fx-background-color: aqua");
        try {
            draw();
        } catch (Exception e) {
            e.printStackTrace();
        }
    }
 
    private class Complex2 {
        public double x;
        public double y;
    }
 
    private void draw() throws CloneNotSupportedException {
        int iter = 50;
        final double max = 1e+6;
        final double min = 1e-6;
        final double tol = 0.005;
        int n;
        double p = 0;
        // начало экранных координат, смещено влево и вверх
        long xc = round(canvas.getWidth() / 2);
        long yc = round(canvas.getHeight() / 2);
        // масштаб для перевода координат в экранные
//        double mx = xc;
//        double my = yc / -y0;
        double mx = canvas.getWidth() / 2;
        double my = canvas.getHeight() / 2;
 
        Complex2 z = new Complex2();
        Complex2 z1;
        Complex2 delta;
        int col1;
        for(double y = y0; y < yn; y = y + tol) {
            for(double x = x0; x < xn; x = x + tol) {
 
                n = 0;
                z.x = x;
                z.y = y;
                z1 = z;
                delta = z;
                while(n <= iter && abs(z.x + z.y) < max && abs(delta.x + delta.y) > min) {
                    z.x += (3 * pow(z.x, 4) + 1) / (4 * pow(z.x, 3));
                    z.y += (3 * pow(z.y, 4) + 1) / (4 * pow(z.y, 3));
                    delta.x += abs(z1.x - z.x);
                    delta.y += abs(z1.y - z.y);
                    z1 = z;
                    n++;
                }
                col1 = (n * 256) / iter;
                if(col1 <= 255) {
                    graphicsContext.setFill(Color.grayRgb(col1));
                }
                graphicsContext.fillOval((x + n * y) + mx, (y + n * x) + my, 1, 1);
            }
        }
    }
}