function [d, T] = razdaljaK(K)
%[d, T] = razdaljaK(K) vrne razdaljo d in koordinate tock (v T) med 
%dvema kroznicama.

%pripravimo parametrizaciji obeh kroznic
p = @(t) [K(3,1)*cos(t) + K(1,1); K(3,1)*sin(t) + K(2,1)];
q = @(t) [K(3,2)*cos(t) + K(1,2); K(3,2)*sin(t) + K(2,2)];

%... in njuna odvoda
pdot = @(t) [-K(3,1)*sin(t); K(3,1)*cos(t)];
qdot = @(t) [-K(3,2)*sin(t); K(3,2)*cos(t)];

%pripravimo se gradient f
gradf = @(X) [2*pdot(X(1))'*(p(X(1)) - q(X(2))); -2*qdot(X(2))'*(p(X(1)) - q(X(2)))];

%pozenemo gradientno metodo
[X, k] = gradmet(gradf, 1e-3, [pi/2; pi], 1e-3, 10000);

%vstavimo parametre
T = [p(X(1)), q(X(2))];
d = norm(T(:, 1) - T(:, 2));