#!/usr/bin/octave-cli
  
  filename = argv(){1};
  col = eval(argv(){2});
  mult = eval(argv(){3});
  
  data.freq = load(filename)(:,col).*mult;
  data.rate = 1;

  [ad, s, err, tau] = allan(data, 2.^[0:nextpow2(length(data.freq))-3]./data.rate, '', 0);

  tosave = vertcat(tau, ad, err)';
  save('-ascii', strcat(strsplit(filename, '.'){1}, '_ad.dat'), 'tosave');
  
  exit