function [cps, de]=readcirc (filename);
% Syntax: [cps, de]=readcirc (filename);
%
% Simple function to read an CIRC file containing 12 coeffs and energy.
% skips 12 byte header. Divides the data into 12 ceps. coeffs and delta
% energy.


ff=fopen (filename,'r');
dummy=fread (ff, 12, 'char');
dummy
aux=fread (ff, [13 inf], 'float');
cps=aux (2:13, :);
de=aux (1, :);
fclose (ff);