I figure out a simple solution:

theta_east = full( theta(east) );
disp(sprintf('East Center: %05d (%f, %f) pi',east,theta_east/pi,phi(east)/pi));

It works on Matlab 2013a.