Скорость объекта который падает свободно вследствие земного притяжения

Перепишем уравнение в следующем виде:
dvdt=-g+kv2m
Используемые команды:
@ - создание анонимной функции
Ode45 - Решение нежестких дифференциальных уравнений
Plot – создание графика
Код матлаб:
clear;clc %ochistka sredi
m=5; g=9.81; k=0.05; %zadanie constant
dvdt=@(t,v)-g+k*v^2/m; %zapicb uravneniya s prisvaivaniem anonimnoi fucnkcii
[t, v]=ode45(dvdt,[0:0.1:15],0); %integrirovanie diff uravneniya v diapazone t s nachalnim usloviem v0
plot(t,v) % sozdanie grafika
xlabel('Time (s)') %podpisb x
ylabel('Velocity (m/s)') %podpisb y
Результат выполнения кода:
Рис