x=input('data? ');
y=importdata(x);
t=y.data(:,1);
mx=y.data(:,2);
my=y.data(:,3);
mz=y.data(:,4);
plot(t,mx, 'p')

Fs = length(t)/t(end);
analyze1(t,Fs)

p=find(t>(24.485*.9));
q=t(t>(24.485*.9));
