function y = g_eta(phi,theta,W,B)

y = [0 0 W-B -4820*cos(theta)*sin(phi) -4820*sin(theta)*cos(theta)*cos(phi) 0]';