clear all

n = 10000

tic
a = rand(n, n);
b = rand(n, 1);
toc

tic
x = a \ b;
toc
norm(a * x - b)

tic
x = inv(a) * b;
toc
norm(a * x - b)
