I have solved a matrix ODE with state u
, and want to plot the solution sol
. But I want to scale the time axis from second to year. Can I do this for all variables, something like:
y2s = 365*24*60*60
plot(sol,vars=((t,u)->(t/y2s,u),0,:), xlim=(0,sol.t[end]/y2s))
or do I need to loop through all variables and do it individually as in :
plot()
for i in 1:length(u0)
plot!(sol,vars=((t,u)->(t/y2s,u),0,i), xlim=(0,sol.t[end]/y2s))
end
?