Pow() with x -> 0.0 smoothing
function softpow input Real x; input Real e; input Real r; output Real y; end softpow;