function y = U(THETA, PHI)
%  **Example**
%  y = sin(THETA).*(sin(PHI)).^2;
y = besselj(1,pi/10*sin(THETA)).*besselj(1,pi/10*sin(THETA));