function [ MSD ] = MeanSquareDisplacement( data )
posX3=data(:,4);
posY3=data(:,5);
u=length(posY3);
uNo=data(end,2)+1;
uTime=floor(u/uNo);

startTime=1;
MSD=zeros(data(end,1),1);
counter=zeros(data(end,1),1);
for t1=startTime:uTime-1
    for t2=t1+1:uTime
        for pp=1:uNo
            time1=data(t1*uNo,1);
            time2=data(t2*uNo,1);
            if time2-time1>0
                MSD(time2-time1)=MSD(time2-time1)+((norm([posX3((t2-1)*uNo+pp);posY3((t2-1)*uNo+pp)]-[posX3((t1-1)*uNo+pp);posY3((t1-1)*uNo+pp)]))^2);
                counter(time2-time1)=counter(time2-time1)+1;
            end
        end
    end
end

MSD=MSD./counter;
loglog(MSD);
end

