format none;

        function xdot = f (x, t)
            # Zmena hodnoty x za jednotku casu
            xdot = zeros (0,1);
            
        endfunction

        # time scale
        t = linspace(1, 100, 200);

        # default values
        x0 = zeros (0,1);
        

        #lsode_options ( 'absolute tolerance', max (50*eps, 0.5e-28) );
        lsode_options ( 'relative tolerance', max (50*eps, 0.5e-28) ); 

        y = lsode ("f", x0, t);
        
        for i = 1:length(t)
            
        endfor
        
        [t',y]