Smooth saturation function
function smoothSat input Real x; input Real xmin "Lower bound of range where y = x"; input Real xmax "Upper bound of range where y = x"; input Real dxmin "Width of lower smoothing range"; input Real dxmax = dxmin "Width of upper smoothing range"; output Real y; end smoothSat;