samer@4: function x=clamp(a,b,x)
samer@4: % clamp - Clamp values into given (closed) interval
samer@4: %
samer@4: % clamp :: real, real, [Size] -> [Size].
samer@4: 
samer@4: x(x<a)=a;
samer@4: x(x>b)=b;
samer@4: %y=min(b,max(a,x));