Small bug in calculating taper corrected

For r < r0:
exp((r-r0)**p) < 1 if p is odd,
               > 1 if p is even.
The second case is not intended (should be <1) and is corrected.