Это очень легко сделать в Maxima , вы можете вставить туда свою функцию:
y(t) := (m*v0*sin(theta)/gamma + g*m^2/gamma^2) * (1 - exp(-gamma/m*(t-t0))) - g*m/gamma * (t - t0) + y0;
v(t) := diff(y(t),t);
display2d:false;
print(v(t));
Попробуйте онлайн!
Вывод
(((m*sin(theta)*v0)/gamma+(g*m^2)/gamma^2)*gamma*%e^-(((t-t0)*gamma)/m))/m -(g*m)/gamma
можно почти напрямую вставить в MATLAB, вам просто нужно заменить %e^
на exp()
:
v(t) = @(t)(((m*sin(theta)*v0)/gamma+(g*m^2)/gamma^2)*gamma*exp(-(((t-t0)*gamma)/m)))/m -(g*m)/gamma