Расчитать момент инерции системы - VB

Формулировка задачи:

Дана система точечных N масс c координатами . Центр масс системы определяется по формуле: Момента инерции системы относительно оси Z, проходящей через центр масс системы: Составить программу расчета момента инерции системы относительно оcи Z, проходящей через центр масс системы. Исходные параметры N5

Решение задачи: «Расчитать момент инерции системы»

textual
                            Листинг программы 
                        

                            Option Base 1
Private Sub Command1_Click()
Dim N As Double, m() As Double, x() As Double, y() As Double, z() As Double
Dim i As Integer
Dim X1 As Double, Y1 As Double, Iz As Double, a As Double, b As Double, c As Double
N = Val(InputBox("Введите N (от 1 до 5)"))
If N < 1 Or N > 5 Then
    MsgBox "Неверное значение N", 0, ""
    Exit Sub
End If
ReDim m(N): ReDim x(N): ReDim y(N): ReDim z(N)
Randomize
a = 0: b = 0
For i = 1 To N
    m(i) = Rnd: x(i) = Rnd: y(i) = Rnd
    a = a + m(i) * x(i)
    c = c + m(i) * y(i)
    b = b + m(i)
    Debug.Print i & "." & " m = " & m(i) & " x = " & x(i) & " y = " & y(i)
Next i
X1 = a / b
Y1 = c / b
Debug.Print " X = " & X1 & " Y = " & Y1
Iz = 0
For i = 1 To N
    Iz = Iz + m(i) * Abs((X1 - x(i)) ^ 2 + (Y1 - y(i)) ^ 2)
Next i
Debug.Print "  Iz = "; Iz
End Sub
                        