function p = presek_premic(T1,T2,s1,s2)
% PRESEK_PREMIC vrne koordinati presecisca premic p =[x; y].
%
% Vhodni podatki:
% T1 in s1 sta tocka in smerni vektor, ki podajata prvo premico.
% T2 in s2 sta tocka in smerni vektor, ki podajata drugo premico.
%
% Vektorji in koordinate točk so podani kot stolpci.
p = T1+linsolve([s1(1) , -s2(1) ; s1(2) , -s2(2)], T2-T1)(1).*s1
end
