Я сразу не заметил, а потом не было времени. В уравнениях
post1182021.html#p1182021 ошибка (проектируя на оси координат силу пропорциональную квадрату скорости квадраты компонент скоростей не получить). Если ТС открыл Доля, то уже на это обратил внимание. И я нагуглил (это было не сложно)
страницу П.Г. Доля. Там есть ещё много материала по matlab.