I've modeled a state space system and now want to plot how the state variables change in time. This is how I've set it up: hx = phi_t * B hy = C * phi_t * B hx1 = []; hx2 = []; hx3 = []; hx4 = []; hx5 = []; for i = 0:0.05:10 hx1 = [hx1 subs(hx(1), t, i)]; hx2 = [hx2 subs(hx(2), t, i)]; hx3 = [hx3 subs(hx(3), t, i)]; hx4 = [hx4 subs(hx(4), t, i)]; hx5 = [hx5 subs(hx(5), t, i)]; end hold on title('Impulse Response: State Vectors') plot(hx1) plot(hx2) plot(hx3) plot(hx4) plot(hx5) hold off I'd like to know how I can show the names of the variables in the figure this code gives out, so I can know which is which (x1,x2,x3,x4 and x5). How can I do that?