c=[0:3 250:255]*257;
cc=(c*255*2+1)/(2^16-1)/2;
ccc=bitshift(floor((c*255*2+1)/(2^16-1)),-1);
[c'/256 cc' ccc' ccc'-floor(c'/256)]
