Скачиваний:
23
Добавлен:
07.01.2014
Размер:
1.38 Кб
Скачать
Option Base 1
Sub main()
Dim n!(2, 2), n_obr!(2, 2), x(2), b(2), f(2)
e = 0.01
For i = 1 To 2
x(i) = Cells(i, 4)
Next i

While p < 1
Call om(n, x)
Call obr(n, n_obr)
Call schet(f, x)
Call umnoj(f, b, n_obr)
Call proverka(x, b, e, p)
g = g + 1
For i = 1 To 2
Cells(5 + i, g) = x(i)
Next i
Wend

End Sub
Sub om(n, x)
n(1, 1) = -x(1) / (0.16 - x(1) ^ 2) ^ 0.5 +1
n(1, 2) = 1
n(2, 1) = 1.2 * 2.72 ^ x(1) - 1
n(2, 2) = -1
End Sub
Sub proverka(x, b, e, p)
For i = 1 To 2
s = s + b(i) ^ 2
x(i) = x(i) - b(i)
Next i
If Sqr(s) < e Then p = 1
End Sub
Sub schet(f, x)
f(1) = (0.16 - x(1) ^ 2) ^ 0.5 + x(2)
f(2) = 1.2 * 2.72 ^ x(1) - x(2) - 1.5
End Sub
Sub umnoj(f, b, n_obr)
s = 0
For i = 1 To 2
For j = 1 To 2
s = s + n_obr(i, j) * f(j)
Next j
b(i) = s
Next i
End Sub
Sub obr(n, n_obr)
Dim a!(2, 4)

For i = 1 To 2
For j = 1 To 4
If j < 3 Then a(i, j) = n(i, j) Else: If i = j - 2 Then a(i, j) = 1 Else a(i, j) = 0


Next j
Next i

For i = 1 To 2
s = a(i, i)
For j = 1 To 4
a(i, j) = a(i, j) / s
Next j
Next i

s = a(2, 1)
For j = 1 To 4
a(2, j) = a(2, j) - s * a(1, j)
Next j

For i = 1 To 2
s = a(i, i)
For j = 1 To 4
a(i, j) = a(i, j) / s
Next j
Next i

s = a(1, 2)
For j = 1 To 4
a(1, j) = a(1, j) - s * a(2, j)
Next j

For i = 1 To 2
For j = 1 To 2
n_obr(i, j) = a(i, j + 2)
Next j
Next i

End Sub

Соседние файлы в папке Решение СНУ - Парилов - 2008
  • #
    07.01.20141.38 Кб23222.txt
  • #
    07.01.201415.87 Кб22Книга1.xls
  • #
    07.01.201430.21 Кб25Ньютон-Рафсон.xls