Метод последовательной линейной верхней релаксации решения сеточных уравнений

Лабораторная работа № 4
Метод последовательной
линейной верхней релаксации
решения сеточных уравнений

Цель работы: ознакомиться с итерационным методом решения сеточных уравнений на компьютере.
Приборы и принадлежности: компьютер.

1)

procedure Button1OnClick;

  const n = 10;

        h = 1/n;

        pi =3.141592654;

        epsilon = 1e-3;

        //gamma = 2/(1+sqrt(sin(pi/2/n*(2-sin(pi/2/n)))));

        gamma=1.9;

   var T,Tx: array [0..n] of real;

       aa,bb,cc,ff,delta : real;

       i,iter : integer;

begin

    for i:=0 to n-1 do T[i]:= 110;   {ввод исходных данных}

        T[n]:=strtofloat(Edit1.Text);

        Tx[0]:=strtofloat(Edit2.Text);

        Tx[0]:=T[n];

        Tx[n]:=T[0];

        aa := 1;     {рабочий блок}

        bb := -2;

        cc := 1;

        ff := 0;

        iter:=1;              {релаксация}

   repeat

        delta:=0;

      for i:=1 to n-1 do

         begin

          Tx[i]:= gamma/bb*(ff-aa*Tx[i-1]-

          cc*T[i+1])+(1-gamma)*T[i];

            if abs(1-T[i]/Tx[i]) > delta

            then delta:=abs(1-T[i]/Tx[i]);

          end;

           T := Tx;

             iter := iter + 1

   until (delta < epsilon) {or (iter > 100)};

   Edit3.Text:= inttostr(iter);

        for i:=0 to n do               {вывод результата}

        begin

        TextBox1.Add(inttostr(i)+' .  '+floattostr(T[i]));

        end;

end;

4)