function sparsetest(nx)

% kode som sammenligner tidsbruken ved full lagring og lagring i
% sparse-formatet i Matlab ved løsning av Ax=b
  
n = nx^2;
A = zeros(n,n);
b = ones(n,1);

% Genererer den glisne A-matrisa
for j=1:n
  A(j,j) = 4;
end
for j=1:n-1
  A(j+1,j) = -1;
  A(j,j+1) = -1;
end
for j=nx+1:n
  A(j,j-nx) = -1;
  A(j-nx,j) = -1;
end
%%%

figure(1)
spy(A)
Ainv = inv(A);
figure(2)
spy(Ainv)

disp('x=inv(A)*b')
tic
x = inv(A)*b;
toc

disp('x=A\b - A full')
tic
x = A\b;
toc

As = sparse(A);
disp('x=A\b - A sparse')
tic
x = As\b;
toc
