Программа с заполнением пикселями по заданной форме (моделирование кристаллизации слитков)

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
 
namespace WindowsFormsApplication1
{
    public partial class Form1 : Form
    { 
        Graphics gr;
        Bitmap myBitmap;
        int N = 10; //начальное количество центров кристаллизации
        Random rnd;//генератор случайных чисел для выбора места и цвета зародышей
        Color c1;//Цвет зародыша, характеризующий его состояние
        int[,] Mat = { { 1, 1, 1 }, { 1, 1, 1 }, { 1, 1, 1 } };//Матрица проверки соседей
        
       //gr = this.panel1.CreateGraphics();
        //gr.Clear(Color.White);
        double R, kolvo;
        public Form1()
        {
            InitializeComponent();
        }
 
       
 
        private void button3_Click(object sender, EventArgs e)
        {
 
            gr = this.panel1.CreateGraphics();
        gr.Clear(Color.White);
        myBitmap = new Bitmap(this.panel1.Width, this.panel1.Height);
        for (int i = 0; i < myBitmap.Width; i++)
        {
        for (int k = 0; k < myBitmap.Height; k++)
        myBitmap.SetPixel(i, k, Color.White);
        }
        c1 = myBitmap.GetPixel(0, 0);
        }
        private Point p()
        {
        return new Point(rnd.Next( 0, this.panel1.Width-1), rnd.Next(0, this.panel1.Height-1));
        }
 
        private void button1_Click(object sender, EventArgs e)
        {
           
          R=(double)numericUpDown1.Value;
          kolvo=(double)numericUpDown2.Value;
         rnd = new Random();
         for (int i = 0; i < N; i++)
         {
            
            Point p1 = p();
             Color c = myBitmap.GetPixel(p1.X, p1.Y);
           
               //if ((c == c1) & (p1.X >= 0) & ((Math.Pow((p1.X - 100), 2) + Math.Pow((p1.Y - 90), 2))) == R * R)// Думаю здесь надо задать условия или ограничения
             if ((c == c1) & ((int)Math.Sqrt((Math.Pow(R, 2) - Math.Pow((p1.Y - 90), 2))) == (p1.X-100)))
                 //if ((c == c1) && ((p1.X * p1.X + p1.Y * p1.Y) == R * R))
                 {
                     myBitmap.SetPixel(p1.X, p1.Y, Color.FromArgb(rnd.Next(255), rnd.Next(255), rnd.Next(255)));
                 }
                 else
                 {
                     i--;
                 }
             
         
         }
 
 
        
        }
 
        private void button2_Click(object sender, EventArgs e)
        {
        int K = (this.panel1.Width-1) * (this.panel1.Height-1) - N;
        while (K > 0)
        {
        Point p1 = p();
        Color col;
        if (Rost(p1, out col))// Или здесь задать условия чтобы пиксели распределялись по правой части полуокружности
        {
        myBitmap.SetPixel(p1.X, p1.Y, col);
        K--;
        gr.DrawImage(myBitmap, 0, 0, myBitmap.Width, myBitmap.Height);
         }
         }
        }
 
        
        private bool Rost(Point p, out Color col)
{
int j0,i0,jk,ik;
bool result = false;
Color c = myBitmap.GetPixel(p.X, p.Y);
col = c;
if (c == c1)
{
if ((p.X == 0))
j0 = 0;
else
j0 = -1;
if ((p.Y == 0))
i0 = 0;
else
i0 = -1;
if ((p.X == myBitmap.Width))
jk = 0;
else
jk = 1;
if ((p.Y == myBitmap.Height))
ik = 0;
else
ik = 1;
for (int i = i0; i <= ik; i++)
{
for (int j = j0; j <= jk; j++)
{
Color c2 = myBitmap.GetPixel(p.X + j, p.Y + i);
if ((c2 != c1) & (Mat[i + 1, j + 1] == 1))
{
result = true;
col = c2;
break;
}
}
if (result) break;
}
}
return result;
}
 
        private void Form1_Load(object sender, EventArgs e)
        {
        
        }
 
 
 
    }
}