// both iterate2d and iterate3d in one file for i = 1:10, X = A*X; N = norm(X); X = X/N; if n == 2, printf("N = %f, X = [%f; %f], i = %d\n", N, X(1,1), X(2,1), i); else printf("N = %f, X = [%f; %f; %f], i = %d\n", N, X(1,1), X(2,1), X(3,1), i); end; end;