I see. Yes the layout option produces one “plot” with several subplots. The other suggestion would then be to try to specify vars. Looking at the source code of the package you are using the plot recipes are in
Each recipe has a vars option. If it is not specified then it produces one series for each function, which without layout makes one plot, and with layout gets split out in to different subplots. However, with vars the behaviour is different. So I would try
plot(sol,vars=(0,3))
to plot the third function against time (the code seems single 0 out for the time).