clc
echo on
[t,y]=ode45('hmotadif', [0:0.1:50], [5;0]);
plot (t, y);
grid;
title ('prùbìh výchylky a rychlosti homoty');
xlabel ('doba simulace');
ylabel ('výchylka, rychlost');