function ret=anglebetween( a,b,c )
  c1=roots([a b c])(1);
  c2=roots([a b c])(2);
  theta=2*atan( imag(c1)/real(c1)/pi );
  
  ret=theta;
endfunction
