function [r] = urci_rad(a,b)
% function [r] = urci_rad(a,b)
% Urceni radu vicekrokove metody dane koeficienty a,b
% Pokud delka vektoru b je o 1 vetsi nez delka a, jedna se o implicitni
% metodu


la=length(a);
lb=length(b);
if la==lb   % explicitni metoda
    b=[0 b];
    lb=lb+1;
end
if lb-la~=1
    error('Nespravne delky vstupnich hodnot!')
end
r=-1;
v=sum(a);   % vyraz pro testovani radu
if v==1
    r=0;
else
    return;
end
while abs(v-1)<1e-6
   r=r+1;
   j=-(0:la-1);
   v1=sum(j.^r.*a);
   j=-(-1:lb-2);
   v2=r*sum(j.^(r-1).*b);
   v=v1+v2;
end
r=r-1;
end

