Using a quadratic segment fixes the discontinuity. The resulting friction curve seems to be somewhere between constant and parabolic pressure distributions.
It would be interesting to parametrize the piecewise distribution, be able to adjust the length of each segment, to see how it affects the friction curve shape. But that is quite a bit of work. So I'll stick with the parabolic distribution for now, will be implementing fy (lateral force) next.