%usage C = fracdim(membed, logimax, mindelay, timeseries);
%
% membed=3 	The The maximum embedding dimension.
% logimax=7	The number of boxes that the spaces is divided into.
%		7 seems to work pretty well.
% mindelay=1	The minimum spacing between points used to compute C(r).
% timeseries  	The time series for fracdim to process

function fracdim(membed, logimax, mindelay, timeseries)

localpath = pwd;
'fracdim reading data...'

of = fopen('tseries','w');
fprintf(of, '%.6f\n', timeseries);
fclose(of);

'fracdim computing C(r)...'
command = sprintf('fracdim membed=%d logimax=%d mindelay=%d file=%s', membed, logimax, mindelay, 'tseries');
unix(command);

'To load data into matlab, please use loaddat.m'