По приведенному Вами тексту я ничего не могу понять. Но сведение двух обыкновенных дифференциальных уравнений второго порядка к системе четырех уравнений первого порядка - стандартный материал. (См. любой учебник, например, Эльсгольц "Обыкновенные дифференциальные уравнения и вариационное исчисление".) Если вопрос по дифференциальным уравнениям, то набирайте формулы в
![$\TeX$ $\TeX$](https://dxdy-02.korotkov.co.uk/f/1/f/0/1f0597fbd477973ca54be5f188c6438f82.png)
.
Если Вы спрашиваете о реализации в Matlab, то фрагмент слишком короткий, чтобы что-то содержательное говорить. (Но такое я ни разу не видел.) Приводите более объёмные исходные тексты. Посмотреть можно, например, в примере 14 (движение тела, брошенного с начальной скоростью
![$v_0$ $v_0$](https://dxdy-04.korotkov.co.uk/f/7/5/1/751613a1a4da78db7647a339cbf261c382.png)
под углом
![$\alpha$ $\alpha$](https://dxdy-01.korotkov.co.uk/f/c/7/4/c745b9b57c145ec5577b82542b2df54682.png)
к горизонту в предположении, что сопротивление воздуха пропорционально квадрату скорости) Доля П.Г. "Использование MATLAB. Решение дифференциальных уравнений." (В Сети был доступен pdf).