
        format none;

        function xdot = f (x, t)
            # Zmena hodnoty x za jednotku casu
            xdot = zeros (8,1);
            xdot(2) = 0; # fixed, Medium 
xdot(1) = 0; # not entered !! ; # reaction, [CO2] 
xdot(3) = 0; # not entered !! ; # reaction, [HCO3.-] 
xdot(4) = 0; # not entered !! ; # reaction, [CO3.2-] 
xdot(5) = 0; # not entered !! ; # reaction, [H.+] 
xdot(6) = 0; # not entered !! ; # reaction, [OH.-] 
xdot(7) = 0; # not entered !! ; # reaction, CP 
xdot(8) = 0; # not entered !! ; # reaction, TA 

        endfunction

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

        # default values
        x0 = zeros (8,1);
        x0(1) = 0; # fixed specie, [CO2]
x0(2) = 1; # fixed compartment, Medium
x0(3) = 0; # fixed specie, [HCO3.-]
x0(4) = 0; # fixed specie, [CO3.2-]
x0(6) = 0; # fixed specie, [OH.-]
x0(7) = 0; # fixed specie, CP
x0(8) = 0; # fixed specie, TA
x0(5) = KW/[H+]; # assignment specie, [H.+]


        #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]
        