u1 = [1;2];  u2 = [-1 ;3] ; z = zeros(2,1);




X = [z , u1, u1+u2, u2, z]
X =
     0     1     0    -1     0
     0     2     5     3     0



plot(X(1,:), X(2,:),'linewidth',2)
hold on
plot([-1 1], [ 0 0])
plot([-1 0], [ 1 0])
close

plot(X(1,:), X(2,:),'linewidth',2)
A = [2 -3; 2 3]
A =
     2    -3
     2     3
Y = A*X
Y =
     0    -4   -15   -11     0
     0     8    15     7     0
X
X =
     0     1     0    -1     0
     0     2     5     3     0
plot(Y(1,:), Y(2,:),'r--')
hold on
plot(Y(1,:), Y(2,:),'r--')
plot(X(1,:), X(2,:),'linewidth',2)



A = [2 2 -4; -1 0 2; 1 3 -1]
A =
     2     2    -4
    -1     0     2
     1     3    -1


rrechl0([A, eye(3)])
{Undefined function 'rrechl0' for input arguments of type 'double'.} 
addpath ../ECHLN/
rrechl0([A, eye(3)])
A =
    1.0000    1.0000   -2.0000    0.5000         0         0
         0    1.0000         0    0.5000    1.0000         0
         0    2.0000    1.0000   -0.5000         0    1.0000
A =
    1.0000         0   -2.5000    0.7500         0   -0.5000
         0    1.0000    0.5000   -0.2500         0    0.5000
         0         0   -0.5000    0.7500    1.0000   -0.5000
A =
    1.0000         0         0   -3.0000   -5.0000    2.0000
         0    1.0000         0    0.5000    1.0000         0
         0         0    1.0000   -1.5000   -2.0000    1.0000
ans =
    1.0000         0         0   -3.0000   -5.0000    2.0000
         0    1.0000         0    0.5000    1.0000         0
         0         0    1.0000   -1.5000   -2.0000    1.0000

diary off