Написать программу, которая рисует на экране кривую Гильберта (использовать рекурсию)

using System;
using System.Drawing;
using System.Windows.Forms;
 
namespace КриваяГильберта
{
    public partial class Form1 : Form
    {
        private Graphics ex;
        private int CurX;
        private int CurY;
        private int z = 10;
 
        public Form1()
        {
            InitializeComponent();
            this.ClientSize = new System.Drawing.Size(450, 450);
        }
 
        private void Form1_Paint(object sender, PaintEventArgs e)
        {
            ex = e.Graphics;
            int n = 2;
            int d;
            d = 8 * z * 2;
            f(d + 5 * z, d + 5 * z + 30, n);
        }
 
        private void SetPos(int x, int y)
        {
            CurX = x;
            CurY = y;
        }
 
        private void mylineTo(bool c, int x, int y)
        {
            if (CurX == 0 || CurY == 0)
                throw new Exception();
            ex.DrawLine(c ? Pens.Black : Pens.WhiteSmoke, CurX, CurY, x, y);
            CurX = x;
            CurY = y;
        }
 
        private void sq4(int x, int y)
        {
            int[] d = { -1, 0, 1, 0, -1 };
            SetPos(x - z, y);
            for (int i = 0; i < 4; ++i)
                mylineTo(true, x + z * d[i + 1], y + z * d[i]);
        }
 
        private void sqc(int x, int y)
        {
            int[] dx = { -1, 1, 1, -1 };
            int[] dy = { -1, -1, 1, 1 };
            int[,] k = new int[2, 3] { { 2 * z, z, 0 }, { z, 2 * z, 3 * z } };
            SetPos(x - 3 * z, y);
            for (int i = 0; i < 4; ++i)
                for (int j = 0; j < 3; ++j)
                    mylineTo(j == 1 ? false : true, x + dx[i] * k[i % 2, j], y + dy[i] * k[1 - i % 2, j]);
        }
 
        private void f(int x, int y, int n)
        {
            int d = 4 * z * (1 << (n - 1));
            int[] dx = { -1, -1, 1, 1 };
            int[] dy = { -1, 1, -1, 1 };
            for (int i = 0; i < 4; ++i)
                if (n > 0)
                    f(x + d * dx[i], y + d * dy[i], n - 1);
                else
                    sq4(x + 2 * z * dx[i], y + 2 * z * dy[i]);
            sqc(x, y);
        }
    }
}