%ISCI_THETA doloci theta^*, resitev g(theta)=0
%k=ISCI_THETA(b,B);
function k=isci_theta_k(b,B);
  f=inline('1-cos(theta)+B/b*(theta-sin(theta))','theta','b','B');
  opt=optimset('fzero');
  k=fzero(f,2*pi,opt,b,B);