Заполнение окружности рандомными пикселями

double[] cwidth = new double[(int)E];
            double[] cheight = new double[(int)E];
            var rnd = new Random();
            for (int z = 0; z < E - 1; z++)
            {
                cwidth[z] = rnd.Next((int)(W - RK), (int)(W + RK));
                cheight[z] = rnd.Next((int)(H - RK), (int)(H + RK));
                this.g.SetPixel((int)cwidth[z], (int)cheight[z], button8.BackColor);
                pictureBox1.Image = this.g;
            }

double r = RK * Math.Sqrt(rnd.NextDouble());
double th = rnd.NextDouble() * 2 * Math.PI;
 
int x = (int) (W + r * Math.Cos(th));
int y = (int) (H + r * Math.Sin(th));

this.g.SetPixel(x, y, button8.BackColor);

pictureBox1.Image = this.g;

for (int i = 0; i < 100; i++)
{
    double r = RK * Math.Sqrt(rnd.NextDouble());
    double th = rnd.NextDouble() * 2 * Math.PI;
 
    int x = (int) (W + r * Math.Cos(th));
    int y = (int) (H + r * Math.Sin(th));
 
    this.g.SetPixel(x, y, button8.BackColor);
}
pictureBox1.Image = this.g;

    public static partial class StaticMethods
    {
        private static readonly Random random = new Random();
        public static (int x, int y)[] GetRandomPoints(int x0, int y0, int x1, int y1, int count, Func<int, int, bool> func)
        {
            (int x, int y)[] points = new (int x, int y)[count];
            if (x0 > x1)
                (x0, x1) = (x1, x0);
            if (y0 > y1)
                (y0, y1) = (y1, y0);
            x1++;
            y1++;
 
            for (int i = 0; i < points.Length; i++)
            {
                int x = random.Next(x0, x1);
                int y = random.Next(y0, y1);
                if (func(x, y))
                {
                    points[i] = (x, y);
                }
                else
                {
                    i--;
                }
            }
 
            return points;
        }
 
        public static (int x, int y)[] GetRandomPointsOfRectabgle(int x0, int y0, int x1, int y1, int count)
            => GetRandomPoints(x0, y0, x1, y1, count, (x, y) => true);
 
        public static (int x, int y)[] GetRandomPointsOfCircle(int xCenter, int yCenter, int radius, int count)
            => GetRandomPoints(xCenter - radius, yCenter - radius, xCenter + radius, yCenter + radius, count,
                (x, y) =>
                {
                    int dx = xCenter - x;
                    int dy = yCenter - y;
                    return (dx * dx + dy * dy) <= radius * radius;
                });
 
    }