Константы используемые в программе :

aq

= 1

- правая граница области G

b

= 1

- левая граница области G

N

= 8

- колличество точек разбиения отрезка

[0,

a

]

M

= 8

- колличество точек разбиения отрезка [0,

b

]

h

1 =

aq

/

N

- шаг сетки по X

h

2 =

b

/

M

- шаг сетки по Y

Переменные :

u

0

- значения сеточной функции U

на k

-ом шаге

u

1

- значения сеточной функции U

на (

k

+1)

-ом шаге

a

- массив коэффициентов шаблона

Описание процедур :

procedure

Prt

(

u

:

masa

)

- печать результата

function

ff

(

x

1,

x

2:

real

):

real

- возвращает значение функции

f

в узле (

x

1,

x

2

)

procedure

Koef

- задаёт значения коэффициентов

Действие :

Берётся начальое приближение u

0

и с учётом краевых условий ведётся вычисление с i

=2 .

N

,

j

=2 .

M

. На каждом итерационном шаге получаем u

1

по u

0

. По достижении заданной точности eps

>0

вычисления прекращаются. И все элементы матрицы A, которые лежат ниже главной диагонали получают итерационный шаг (

k

+1

)

, а те элементы которые лежат выше главной диагонали (исключая главную диагональ) получают итерационный шаг k

.