Thank your for your response. But the key issue is there are no plots with labels in the given axis:
julia> axislegend(ax)
ERROR: There are no plots with labels in the given axis that can be put in the legend. Supply labels to plotting functions like `plot(args...; label = "My label")`
I wonder if I can add labels to each axis through ax?