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

// *** Funkce definujici system - kruhovy test
function xdot = ctest(t,x)
  xdot = [ x(2); -x(1) ];

// take lze pouzit ekvivalentni maticovy zapis:
//  xdot = [  0 1  ;
//           -1 0  ]  *  x;

endfunction

t = 0:0.1:20; // kdy chceme vysledne 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,:));