многократное повторение реализаций функций импульсного отклика

Есть часть кода по вычислению функций импульсного отклика в нелинейной модели(файл прилагается). В этот код входят случайные шоки. В результате выполнения этого кода имеем функции отклика в виде векторов (ТI,1). ТI=200 у нас. Необходимо реализовать 500 повторений этого кода при разных случайных шоках (они будут разными автоматически), чтобы получилась матрица функций отклика (200, 500). Помогите, пожалуйста.