Phase angle

`theta = angle(z)`

`angle`

takes a complex number *z* = *x* +
*i**y* and uses the `atan2`

function to compute the angle between
the positive *x*-axis and a ray from the origin to the point
(*x*,*y*) in the *xy*-plane.