function y = pseudocolor(x) global L; L = 256; x = double(x); r = graytran(x); g = graytran(x - (L-1)/12); b = graytran(x - (L-1)/6); y(:,:,1) = r; y(:,:,2) = g; y(:,:,3) = b; y = uint8(y); function s = graytran(r) global L; s = round((L-1)*abs(sin(5*pi/(L-1)/2*r)));