hitr.m

function [vx,vy,v]=hitrost(x,y,tdif)
        xdif=diff(x);
        ydif=diff(y); 
        vx=xdif/tdif;
        vy=ydif/tdif;    

        vx=[0;vx];
        vy=[0;vy];
        v=sqrt(vx.*vx+vy.*vy);
end