Join := function(b)

	local nrows, ncolumns, F, v;
	
	nrows := NumberOfRows(b);
	ncolumns := NumberOfColumns(b);
	F := Parent(b[1][1]);
	v := Matrix(F, 1, nrows * ncolumns, []);
	V := VectorSpace(F, nrows * ncolumns);

	for i in [1..nrows] do	
		for j in [1..ncolumns] do
			v[1][(i-1)*ncolumns + j] := b[i][j];
		end for;
	end for;
	
	return V!v;
end function;