Matlab: решение связанных ОДУ с изменяющимися во времени параметрами

Я пытаюсь решить набор ODE в Matlab (я думаю, ODE45 - правильный инструмент), которые имеют изменяющиеся во времени параметры. Вот иллюстрация:

x' = x/y+z*a(t)
y' = y*x+a(t)/y+b(t)
z' = x*y + c(t)/z + d

У меня есть данные для a(t), b(t) и c(t), указанные в таблице поиска, подготовленной заранее. d — это то, что нужно вычислять из внешних данных на каждой итерации.

Я занимаюсь решением уравнений движения системы. Как я могу это сделать? Спасибо!