Public Function Gamma_process(ByVal a, ByVal B, ByVal N, ByVal start, ByVal fin, Optional Gamma0 As Double)
'Simulate a Gamma process
    Dim vec() As Double
    Dim i As Long
    Dim dt As Double
    
    ReDim vec(0 To N)
    vec(0) = Gamma0
    dt = (fin - start) / N
    
    For i = 1 To N
        vec(i) = vec(i - 1) + Gamma(a * dt, B)
    Next i
    Gamma_process = vec
End Function