Как расчитать функцию Лапласса для параметра Z

Private Sub Form_Load()
 
Dim WorkNorm(8, 8) As Double
Dim WorkOpt(8, 8) As Double
Dim WorkPes(8, 8) As Double
Dim OptOcenka As Double
Dim PesOcenka As Double
Dim Otkloneniye(8, 8) As Double
Dim SredKvadOtkl As Double
Dim Dispersia As Double
Dim NV(8, 8) As Double
Dim MatWait As Double
Dim Z As Double
Dim Ver As Double
Dim TDirect As Double
 
WorkNorm(1, 2) = 4: WorkNorm(1, 3) = 3: WorkNorm(1, 4) = 5: WorkNorm(2, 5) = 7: WorkNorm(3, 6) = 8: WorkNorm(4, 7) = 9
WorkNorm(2, 6) = 10: WorkNorm(4, 6) = 12: WorkNorm(6, 8) = 10: WorkNorm(5, 8) = 8: WorkNorm(1, 2) = 11
 
OptOcenka = 0.2
PesOcenka = 0.25
TDirect = 30
 
For i = 1 To 8
  For j = 1 To 8
  WorkOpt(i, j) = WorkNorm(i, j) - WorkNorm(i, j) * OptOcenka
  Next j
Next i
 
 
For i = 1 To 8
  For j = 1 To 8
  WorkPes(i, j) = WorkNorm(i, j) + WorkNorm(i, j) * PesOcenka
  Next j
Next i
 
For i = 1 To 8
  For j = 1 To 8
  Otkloneniye(i, j) = 0.028 * ((WorkPes(i, j) - WorkOpt(i, j)) ^ 2)
  Next j
Next i
 
SredKvadOtkl = 0
For i = 1 To 8
  For j = 1 To 8
  SredKvadOtkl = SredKvadOtkl + Otkloneniye(i, j)
  Next j
Next i
 
Dispersia = Sqr(SredKvadOtkl)
 
 
For i = 1 To 8
  For j = 1 To 8
  NV(i, j) = (WorkOpt(i, j) + 4 * WorkNorm(i, j) + WorkPes(i, j)) / 6
  Next j
Next i
 
MatWait = NV(1, 4) + NV(4, 6) + NV(6, 8)
 
Z = Round(((TDirect - MatWait) / Dispersia), 8)
 
Ver = Ôóíêöèÿ_Ëàïëàññà(Z)
 
MsgBox Z 'Ver
 
End Sub
 
 
Public Function Ôóíêöèÿ_Ëàïëàññà(a As Double) As Double
Dim h As Double
Dim n As Double
Dim Èíòåãðàëë As Double
Dim e As Double
Dim sum As Double, F0 As Double, Fn As Double
Dim pi As Double
Dim Êîýôôèöèåíò As Double
 
pi = 3.1415926535
e = 2.71828182845905
n = 10
h = Round(((a - 0) / n), 8)
 
For i = 1 To n - 1
sum = sum + e ^ (-0.5 * (i) ^ 2)
Next i
 
F0 = e ^ (-0.5 * 0 ^ 2)
Fn = e ^ (-0.5 * 1000 ^ 2)
 
Èíòåãðàëë = Round(((h * ((F0 + Fn) / 2 + sum) + (((a - 0) * h ^ 2) / (-12)))), 8)
 
Êîýôôèöèåíò = Round(((1 / Sqr(2 * pi))), 8)
Ôóíêöèÿ_Ëàïëàññà = Êîýôôèöèåíò * Èíòåãðàëë
 
End Function

Public Function Ôóíêöèÿ_Ëàïëàññà(a As Double) As Double
Dim h As Double
Dim n As Double
Dim Èíòåãðàëë As Double
Dim e As Double
Dim sum As Double, F0 As Double, Fn As Double
Dim pi As Double
Dim Êîýôôèöèåíò As Double
 
pi = 3.1415926535
e = 2.71828182845905
n = Round(a * 10, 5)
h = Round(((30 - a) / n), 8)
 
For i = 1 To n - 1
sum = sum + e ^ (-0.5 * (i) ^ 2)
Next i
 
F0 = e ^ (-0.5 * 0 ^ 2)
Fn = e ^ (-0.5 * n ^ 2)
 
Èíòåãðàëë = Round(((h * ((F0 + Fn) / 2 + sum)) + (((a - 0) * h ^ 2) / (-12))), 8)
 
Êîýôôèöèåíò = Round(((1 / Sqr(2 * pi))), 8)
Ôóíêöèÿ_Ëàïëàññà = Round(Êîýôôèöèåíò * Èíòåãðàëë, 8)
 
End Function

Option Explicit
Const pi = 3.141592653589
Private Sub Form_Load()
    Dim x As Double, ox As Double, ov As Double, v As Double
    Me.Scale (-4, 1)-(4, -1): Me.AutoRedraw = True: Me.DrawMode = vbXorPen
    Me.Line (-4, 0)-Step(8, 0), vbYellow: Me.Line (0, -1)-Step(0, 2), vbYellow
    For x = -4 To 4 Step 0.5
        If x = -4 Then
            ox = x: ov = Laplas(x)
        Else
            v = Laplas(x)
            Me.Line (x, -1)-Step(0, 2), &H101010
            Me.Line (-4, x / 4)-Step(8, 0), &H101010
            Me.Line (x, 0.01)-Step(0, -0.02), vbWhite
            Me.Line (0.03, x / 4)-Step(-0.06, 0), vbWhite
            Me.Line (ox, ov)-(x, v), vbCyan
            ox = x: ov = v
        End If
    Next
End Sub
Private Function Simpson(a As Double, b As Double, e As Double) As Double
    Dim i As Long, n As Long, h As Double, ov As Double, v As Double, q As Double
    n = 1
    Do
        ov = v: n = n * 2: h = (b - a) / n
        v = Exp(-a * a / 2) + Exp(-b * b / 2)
        For i = 1 To n - 1
           q = a + i * h
           v = v + 2 * (1 + i Mod 2) * Exp(-q * q / 2)
           DoEvents
        Next
        v = v * (b - a) / (3 * n)
        DoEvents
    Loop Until Abs(v - ov) <= e
    Simpson = v
End Function
Private Function Laplas(x As Double, Optional e As Double = 0.00000001) As Double
    Laplas = Simpson(0, x, e) / Sqr(2 * pi)
End Function

Private Function Simpson(a As Double, b As Double, e As Double) As Double
    Dim i As Long, n As Long, h As Double, ov As Double, v As Double, q As Double
    n = 1
    Do
        ov = v: n = n * 2: h = (b - a) / n
        v = Exp(-a * a / 2) + Exp(-b * b / 2)
        For i = 1 To n - 1
           q = a + i * h
           v = v + 2 * (1 + i Mod 2) * Exp(-q * q / 2)
           DoEvents
        Next
        v = v * (b - a) / (3 * n)
        DoEvents
    Loop Until Abs(v - ov) <= e
    Simpson = v
End Function

Public Function Метод_Трапеций(a As Double, b As Double) As Double
Dim h As Double
Dim n As Double
Dim e As Double
Dim sum As Double, F0 As Double, Fn As Double, Fx As Double
 
 
e = 2.71828182845905
n = 1000
h = Round(((a - b) / n), 8)
 
sum = 0
For i = 1 To n - 1
  Fx = b + i * h
  sum = sum + e ^ (-Fx * Fx / 2)
Next i
 
F0 = e ^ (-0.5 * 0 ^ 2)
Fn = e ^ (-0.5 * n ^ 2)
 
Метод_Трапеций = Round(((h * ((F0 + Fn) / 2 + sum))), 8)
 
End Function
 
Public Function Лапласс(a As Double) As Double
Dim pi As Double
pi = 3.1415926535
 
  Лапласс = Round((Метод_Трапеций(a, 0) / Sqr(2 * pi)), 8)
 
End Function