Как я могу использовать ODE45 для решения связанных дифференциальных уравнений, зависящих от частоты?

Question

Уравнение 1

Уравнение 2

FREQ = 1:1000; P = 1i * (2 * pi * FREQ);

Постоянные значения и уравнения переменных состояния

Для этой конкретной задачи функция форсирования P зависит от частоты.Это будет двигать массы.Как я могу реализовать ODE45 для решения уравнений с частотой, являющейся изменяющимися условиями?(Не обращайте внимания на значение P на изображении 3.)