A bit late to the party... but there is no need for external packages. approx actually does achieve what you want. What you want is to interpolate only between the min and max values, so simply do:

approx(range(s), c(0,1), s)$y

or equivalently,

approxfun(range(s), c(0,1))(s)