Rotating boson stars using finite differences and global Newton methods