function [x,ifail] = getcorrelation

% Read data from file
fid = fopen('rawdata.dat');
at = zeros(10,91);
[at,count] = fscanf(fid, '%f', [10 91]);
a = at';
clear at

% Using the data in A, form a correlation matrix
[xbar, std, v, x, ifail] = g02bx(a);
clear xbar std v;
disp(['Ifail ' int2str(ifail)])

end