posp1.m

function [ax,ay,a]=posp1(x,y,tdif)
        [vx,vy,v]=hitr(x,y,tdif)
        xdif=diff(vx);
        ydif=diff(vy);
        ax=xdif/tdif;
        ay=ydif/tdif;    

        ax=[0;ax];
        ay=[0;ay];
        a=sqrt(ax.*ax+ay.*ay);
end