function y = smoothstep(slope, x)
% function that maps [0 1] to [0 1] with parametrized slope 0.5
% used for non uniform quantization

z = exp(4*slope*(x-.5));
y = z ./ (1+z);

