Исследование двумерного распределения частиц при диффузии с поверхности

[Gloky]

Рассмотрим диффузию атомов примеси в кристалл по междоузельному механизму. Пусть междоузельные положения образуют квадратную решетку. Первоначально все они не заняты, а атомы находятся на левой границе.
В дальнейшем каждый атом может перепрыгнуть с вероятностью 1/4 в соседнее положение (если оно не занято). Граничные атомы могут прыгать только вправо, причем их место сразу после прыжка занимают новые примесные атомы (модель неограниченного источника примеси).
Смоделировать процесс диффузии атомов, показывая на экране их новые положения через 1 шаг (в данном случае это время, за которое будут перебраны все атомы, находящиеся в расчетной области). На верхней и нижней применять периодические граничные условия. По оси y взять не менее 100 узлов решетки. Число узлов по оси x взять таким, чтобы за заданное время диффузии атомы не достигали правой границы.
Вот такая задача. Помогите реализовать в Delphi. Хотя бы скажите с чего начать, очень прошу. Спасибо.

[evg_m]

По оси y взять не менее 100 узлов решетки. Число узлов по оси x взять таким, чтобы за заданное время диффузии атомы не достигали правой границы.

1. "нарисовать" пустую решетку и ось X ось Y.
2. нарисовать "атом" в координатах (X1, X2)
3. переместить "атом" на одну позицию влево/вправо/вверх/вниз.
4. после выполнения п.1-3 прочитать задание снова.

[Gloky]

Между 3 и 4 пунктами нужно добавить еще один атом на левую границу.