t=-6:0.1:6;
y=exp(-t/3).*(-cos(3*t)+(5/9.)*sin(3*t));
w=sqrt(125/81.)*exp(-t/3);
plot(t,y,t,w,t,-w);
%set(axes,'FontUnits',16);
title('exp(-t/3)*(-cos(3*t)+(5/9)*sin(3*t))','FontSize',24);
xlabel('t','FontSize',24);
ylabel('y(t)','FontSize',24);
grid on;
t_size=size(t);
t_ne=t_size(2);
y_text_pos = round(t_ne*.17);
w_text_pos = round(t_ne*.17);
m_text_pos = round(t_ne*.05);
text(t(y_text_pos),y(y_text_pos),'\leftarrow y(t) (solution)',...
     'HorizontalAlignment','left','FontSize',24);
text(t(w_text_pos),w(w_text_pos),'\leftarrow 1.24226*exp(-t/3)',...
     'HorizontalAlignment','left','FontSize',24);
text(t(m_text_pos),-w(m_text_pos),'\leftarrow -1.24226*exp(-t/3)',...
     'HorizontalAlignment','left','FontSize',24);
