Метод конечных разностей, опрос по задаче МКР

[nihi]

Люди, если кто программировал МКР, помогите! ЕСть тепловая задча только вместо T - температуры у меня C- концентрация. Пытаюсь реализовать на делфи явную разностную схему по http://ikt.muctr.ru/html2/7/lek7_4_3.html Хэлпаните плс!
ВОт кусок про-мы:

delta_fi:=2*pi/20;
  C0:=C0+C0*sin(delta_fi);//концентрация в соответствии с нач. условиями
  tau_na_dva:=sqr(delta_fi)/2*D;//определяем расчетный шаг по времени, руководствуясь условием устойчивости
  Rc:=D/sqr(r1+h_u);


  //{определяем расчетный шаг сетки по пространственной координате}
  hx:=delta_fi;


  //определяем концентрацию в начальный момент времени
   for i:=1 to Nx do
      for j:=1 to Ny do

      C[i,j]:=C0;

      time:=0;
      while time<t_end do

         begin
           time:=time+tau_na_dva;

         for i:=2 to Nx-1 do
          for j:=2 to Ny-1 do

         end;

          CC[i,j]:=C[i,j]+Rc*tau_na_dva/sqr(hx)*(C[i,j-1]-2*C[i,j]+C[i,j+1]);

           begin

          for i:=1 to Nx do
           for j:=1 to Ny do
           CC[1,j]:=0;CC[Nx,j]:=Cz;CC[i,1]:=0;CC[i,Ny]:=Cz;

           time:=time+tau_na_dva;
       end;


      begin

        for i:=1 to Nx do
          for j:=1 to Ny do

       writeln(CC[i,j]:8:5);
      end;
   readln;readln;






end.

В делфи не силен, прошу строго не судить

[s-andriano]

Попытайтесь:
1. Правильно "лесенкой" расставить операторы.
2. Написать комментарии - зачем нужны те или иные строки программы.
Вероятнее всего, это поможет Вам самостоятельно найти ошибки.

Вообще же ошибки в программе бывают двух видов:
1. Синтаксические - эти ошибки Вам покажет компилятор.
2. Логические - когда формально правильная программа делает не то, что от нее требуется. Чтобы найти такие ошибки, нужно знать, что должна делать программа. Без комментариев это установить практически невозможно, можно только догадываться. Но мы ведь здесь собрались не для того, чтобы упражняться в телепатии?