Here is the output for Julia 1.5.2.
Now I have a new question: how do I adjust the position of the axes labels so that they do not overlap with the tick labels?
Also, How do I put z axis label to the left of the plot? And how do I rotate it by 90 degrees?
(I am curious if there are the corresponding attributes in Plots.jl. If there are not, I guess I should better stick with PyPlot.jl)