Mine

function R = RT_2D(B, s, e)
	// 2D --> 1D
	// *** Funktioniert noch nicht richtig ***
	m = (s + e - 1) / 2
	sizor = size(B,1)

	if e - s == 0 then
		R = RT(B, 1, size(B,1))
	else
		Y = [];
		Y(1:sizor, s:m) = B(1:sizor, s:m) + Y(1:sizor, m+1:e);
		Y(1:sizor, m+1:e) = abs(B(1:sizor, s:m) - B(1:sizor, m+1:e))

		RT_2D(B, s, m);
		RT_2D(B, m+1, e);

		R = Y;
	end;

endfunction