function [x,y,dx,dy,time,M,centroids]=pyanalyze(mainout,N)
n=0;
centroids=[];
while n<N
    display(['Image Number ' num2str(n)])
    if n==0
        I=pyimportsingle(mainout,n);
        cent=centroid_image(I);
    else
        I=pyimportsingle(mainout,n);
        cent=centroid_image(I,cent);
    end
    centroidn=cent.current_centroids;
    centroids=cat(3,centroids,centroidn);
        
n=n+1;
end
M=pyimportM(mainout,N);
time=0.5*(M(:,2)-M(:,1))+M(:,1)-M(1,1);

x=centroids(:,1,:);
y=centroids(:,2,:);
i=1;
dx=[];
dy=[];
while i<(size(centroids,3)+1)
dxi=centroids(:,1,i)-centroids(:,1,1);
dyi=centroids(:,2,i)-centroids(:,2,1);
dx=cat(3,dx,dxi);
dy=cat(3,dy,dyi);
i=i+1;
end