function R = rot_y(alpha)
	c = cos(alpha);
	s = sin(alpha);
%	R = [ c 0 -s; 0 1 0; s 0 c];
	R = [ c 0 s; 0 1 0; -s 0 c];