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