SIMULINK!!!

function[sys, x0] = sistem(t,x,u,flags) % funkcija -
    if(flag == 0)  % inicijalizacija
        sys = [3;0;1;1;0;0] %brPromStanja, nista, brIzlaza, brUlaza, 0, veza ul i iz
        x0 = [0;1;-1] % vektor pocetnih uslova
    else if(flag == 1)  % racunanje izvoda
            sys(1) = x(2) %%%%
            sys(2) = x(3)%%%%%   izvodi
            sys(3) = ......%%%
        else if(flag == 3)  % racunanje izlaza
                sys = [x(1)] % gde je izlaz
            else
                sys = []
            end