function [x,y] = hyperbola(f1,f2,tdoa)
% f1 = souradnice jedne BTS (x,y) kartezske
% f2 = souradnice druhe stanice BTS (x,y) kartezske
% tdoa = rozdíl doby pøíchodu signálu vzdálenost v m

    % Vektor mezi ohnisky
    delta = f1-f2;
    
    % pøevod do polárních souøadnic
    [phi,r] = cart2pol(delta(1),delta(2));    

    % vzdálenost mezi ohniskem osou hyperboly     
    rd = (r+tdoa)/2;        
    
    % parametry hyperboly
    a = (r/2)-rd;
    e = r/(2*a);
    c = a*e;
    b = sqrt(c^2-a^2);
    
    % Výpoèet jednotlivých bodù hyperboly z parametrické rovnice.
    x = zeros(1,499);
    y = zeros(1,499);
    for k = 1:499                
        mu = (k-249)/50;      
        x(k) = a*cosh(mu);
        y(k) = b*sinh(mu);      
        x(k) = x(k)-c;
        y(k) = y(k);                
        [phi2,r2] = cart2pol(real(x(k)), real(y(k)));
        [x(k), y(k)] = pol2cart(phi2+phi, r2);                            
    end
    x = x+f1(1);
    y = y+f1(2);    
end

