function [cas,cas_premica]=risi_brah(T1,T2,st_tock,brisi);
%RISI_BRAH narise brahistohrono
%RISI_BRAH(T1,T2,st_tock)
%T1=[x_1;y1]; T_2=[x_2;y_2]
%st_tock je stevilo razdelitev intervala [0,theta*]
%brisi je 1, ce zelimo pobrisati sliko
%vrnemo tudi cas potovanja po brahistohroni ter cas potovanja po premici
%med danima tockama

if brisi==1
    clf;
end

%naredimo premik tock:
b=T2(1)-T1(1);
B=T2(2)-T1(2);
[k,theta0]=isci_theta_k(b,B);
theta=linspace(0,theta0,st_tock);
x=1/2*k^2*(theta-sin(theta))+T1(1);
y=-1/2*k^2*(1-cos(theta))+T1(2);
hold on;
plot(x,y,'b');
plot([T1(1),T2(1)],[T1(2),T2(2)],'ro');
grid;
axis equal;
hold off;
g = 9.81;
cas = sqrt(1/(2*g))*theta0*k; % čas po cikloidi

d = sqrt(b^2+B^2);
a = g/sqrt(2);
cas_premica = sqrt(2*d/a); % čas po premici