// Priklad pro SciLab
// Kodovani: ASCII
// Dalsi informace viz https://help.scilab.org/
// Pouziti: scilab -f "scilab-test.txt"

// Funkce definujici system - kruhovy test
function xdot = ctest(t,x)
  xdot = [ x(2); -x(1) ];
//
// nebo muzete pouzit ekvivalentni maticovy zapis:
//  xdot = [  0 1  ;
//           -1 0  ]  *  x;
//
endfunction

t = 0:0.1:20; // kdy chceme vysledné hodnoty

// Reseni obyc. dif. rovnic
x = ode([1;0], 0.0, t, ctest);
//      ^^^^^ pocatecni  podminky

// Vizualizace vysledku
//xbasc();
plot2d(t,x(1,:));
plot2d(t,x(2,:));

//xinit(); // nove okno
//xbasc();
//fchamp(ctest, 0, -1.2:0.1:1.2, -1.2:0.1:1.2);
//plot2d(x(1,:),x(2,:));