clear all; clc;
%t = linspace(1, 3, 2000);
a = 2;u = -1;  %Mudar na 
figure
hold;
xlabel('x_1');ylabel('x_2');
x2(1) = 0 ;x1(1) = 0;
c1 = x2(1) + u/a; c2 = x1(1) - c1/a;
k = 1;
for t = 1:0.001:3
%     x1 = @(t) (c1*exp(-a*t) + u*t)/a + c2;
%     x2 = @(t)  c1*exp(-a*t) + u/a;
    k = k+1;
    x1(k) = (c1*exp(-a*t) + u*t)/a + c2;
    x2(k) = c1*exp(-a*t) + u/a;
end

plot(x1(1:k), x2(1:k))
grid on
%axis equal