function res = Lab1()
  U = [1 0];
  V = [0 1];
  T = 2 * pi / 3;
  t = linspace(0, T, 1000);
  x_u = lsode("F", U, t);
  x_v = lsode("F", V, t);
  C = [ transpose(x_u(1000, :)) transpose(x_v(1000, :))];
  m = 1000;
  C
  for i = 1:m
    C = C * C;
  endfor
  C

  function res = F(x,t)
    res(1) = 0.5 * x(1) + 0.1 * sin(2 * t / 3) * x(2);
    res(2) = e ^(sin(0.4 * t)) * x(1) + x(2) * 2;
  endfunction
endfunction