function theta = phase(beta,xi)
% This function returns the pahse angle, theta, associated with the 
% parameters:
% beta - the frequency ratio
% xi - the damping ratio
% Written by Dr Colin Caprani - www.colincaprani.com

theta = atan2((2*xi.*beta),(1-beta.^2)); % refers to complex plane