function [x1,x2] = kv_rovnice(a,b,c)
% kv_rovnice - reseni kvadraticke rovnice ax^2+bx+c=0
% [x1,x2] = kv_rovnice(a,b,c)
% a,b,c ... koeficienty kvadratikce rovnice (a<>0)
% vystupy:
% x1,x2 ... koreny kvadraticke rovnice
% priklady volani: [x,y]=kv_rovnice(1,-2,1)
%                  [x,y]=kv_rovnice(1,0,1)
%                  [x,y]=kv_rovnice(2,-5,-3)

if nargin<3 % kontrola poctu vstupu - musi byt 3
    error('Zadejte TRI vstupy!')
end
if abs(a) < 2*eps % a==0? Kontrolujeme pomoci konstanty 'eps'
    error('Zadejte NEnulovy koeficient u kvadratickeho clenu (1. vstup)!')
end
D = b^2-4*a*c; % diskriminant
x1 = (-b+sqrt(D))/(2*a); % nebo (-b+sqrt(D))/2/a
x2 = (-b-sqrt(D))/(2*a);