ВУЗ:
Составители:
89
REM M - максимальное число иттераций
DIM dX(N) AS DOUBLE, B(N) AS DOUBLE
DIM A(N, N) AS DOUBLE
M = 500: S = 0
begin:
REM Формирование матрицы частных производных
CALL SystemUrawnenie(N, F(), X())
FOR i = 1 TO N
B(i) = -F(i)
NEXT i
FOR j = 1 TO N
X = X(j): h = E * ABS(X)
X(j) = X + h
CALL SystemUrawnenie(N, F(), X())
FOR i = 1 TO N
A(i, j) = (F(i) + B(i)) / h
NEXT i
X(j) = X
NEXT j
REM Подсчет количества итераций
S = S + 1
IF S = M + 1 THEN PRINT "Сделаномаксимальноечислоитераций"; S:
GOTO finis
REM Обращениематрицы
CALL Rmatr(N, A(), det)
IF det = 0 THEN GOTO finis
REM Определение приращений dX(i) и новых значений X(i)
FOR i = 1 TO N
dX(i) = 0
FOR j = 1 TO N
dX(i) = dX(i) + A(i, j) * B(j)
NEXT j
X(i) = X(i) + dX(i)
NEXT i
REM Проверка условий сходимости
R = 0
FOR i = 1 TO N
Nitro PDF Trial
www.nitropdf.com
Страницы
- « первая
- ‹ предыдущая
- …
- 87
- 88
- 89
- 90
- 91
- …
- следующая ›
- последняя »