function [A] = zero_matrix(N);
  A = zeros(N,N);
endfunction;