bode phase incorrect for larger delay

[albheim]

The phase plot does some funky business here, shows up around a delay of 0.8 and grows with more delay. Numerical problem, or just hard to simulate system?

using Plots, ControlSystems
sys = delay(2) / (tf("s") + 1)
bodeplot(sys)

[baggepinnen]

I've noticed this too, when the frequency grid is not dense enough, the jump in phase between two frequency points gets over 360 deg and then unwrap does not work properly