% Testprogramm zu SubMat

n = 5

for j=0:n-1
  a(j+1) = (j==0)/4 + (3+2*cos(2*pi*j/n))/4/n;
end

A = SubMat(a);

load Aref

if norm(A-Aref)<.01
  disp('Ihr Programm arbeitet korrekt')
else
  disp('Ihr Programm arbeitet nicht korrekt')
end