function makeimage(C, colors)

% C: (un-normalized) modulus matrix.
% colors: number of colors.

if (nargin == 1)
	colors = 256;
end

a = min(min(C));
C = (C-a) * ((colors-1)/(max(max(C)) - a));
C = C + 1;
C = round(C);
colormap(gray(colors));
image(C);
%Title = sprintf('s: %f, c: %f, step: %f', s, c, steps(1));
%title(Title);
