function printtime(t1)
t1=diff(t1);
n=length(t1);
disp(['cpu time (seconds)  ' num2str(t1(n))])