[t,x]=ode45(@reg_perm_sin,[0 30],[5 1 8]);

z=x'; y=z(1,:)+z(2,:); 
plot(t,y);
hold on;
u=sin(t); 
plot(t,u,'g');

u=sqrt(2)*sin(t-pi/4); 
plot(t,u,'r');

[t,x]=ode45(@reg_perm_sin,[0 30],[0 0 0]);

z=x'; y=z(1,:)+z(2,:); 
plot(t,y,'magenta');

legend('y=x_1+x_2,x(0)=(5,1,8)','sin(t)(ingresso)','sqrt(2)sin(t-pi/4)(regime permanente)','y=x_1+x_2,x(0)=(0,0,0)');