A = construct(3,4, @(i,j) i+j);

expected = [ 2, 3, 4, 5;
	     3, 4, 5, 6;
	     4, 5, 6, 7];

unit_test_equals("construct a_ij == i+j", ...
		 expected, ...
		 A);
