function[output] = koreny(a, b, c)
    D = b^2 - 4*a*c;
    if D < 0  
        disp("dva komplexne sdruzene koreny")
        output = [(- b - sqrt(D)) / 2*a, (- b + sqrt(D)) / 2*a];
    elseif D == 0
        disp("jeden dvounasobny koren")
        output = - b / 2*a;
    else
        disp("dva realne koreny")
        output = [(- b - sqrt(D)) / 2*a, (- b + sqrt(D)) / 2*a];
    end
end