function [D] = actualdat(A,f,sigma)
save('sigma.txt','sigma','-ascii')
M = modtimedat(A,f)
randn('state',0)
D = M + sigma*randn(100,1)
D;
save('Ddata.txt', 'D', '-ascii')