function y = U(THETA, PHI)
%  **Example**
%  y = sin(THETA).*(sin(PHI)).^2;
y = sin(THETA).*(sin(PHI)).^2;