function y=nsd(x, y)
% nsd nejvetsi spolecny delitel cisel x, y
if (x<=0 | y<= 0),
	error('argumenty musi byt kladne');
end
X=x;
Y=y;
while X ~= Y
	if X > Y
		X=X-Y;
	else
		Y=Y-X;
	end
end
y=X;