function [n time gflops xc] = LinpackMC;

% Derek O'Connor.   Jan 29 2008.  www.derekroconnor.net

disp('Starting Linpack Benchmark');
pause;
maxNumCompThreads(8);
n = [20 15 10 5 4 3 2 1]*10^3;
nvals = length(n);
time = zeros(nvals,1);
gflops = zeros(nvals,1);
for k = 1:nvals
    n(k)
    A = randn(n(k),n(k));
    x = ones(n(k),1);
    b = A*x;
    tic;
    A\b;
    time(k) = toc;
    gflops(k) = 2*n(k)^3/3 + 2*n(k)^2;
    clear A x b
end;
n = n';