It is a minor variation of the example given by Scott for some data set y, cosy, y2, sqrty, siny, logy (doesn't matter the values, could be random)
plot_1 = plot([y cosy], title = "Data y", xlims = (0,10), ylims = (-0.1,1.1), grid = true, xlabel = "Iteration", ylabel = "y & cos(y)" ); plot_2 = plot(y2, title = "Data y Squared", xlims = (0,10), ylims = (-0.1,0.5), grid = false, xlabel = "Iteration", ylabel = "y^2", legend = false, ); plot_3 = plot( [sqrty siny], title = "Square Root of y", xlims = (0,10), ylims = (-0.1,0.5), grid = false, xlabel = "Iteration", ylabel = "y^2" ); plot_4 = plot(logy, title = "Log of y", xlims = (0,10), ylims = (-10,0), grid = false, xlabel = "Iteration", ylabel = "y^2" ); plot_5 = plot(expy, title = "Exp of y", xlims = (0,10), ylims = (1,2), grid = false, xlabel = "Iteration", ylabel = "y^2" ); lay = @layout [a grid(2,2); b] plot( plot_1, plot_2, plot_3, plot_4, plot_5, plot_5, layout = lay )