  function mu = fEnd(t)
% Let RS be the right tangent point on the Sun as
% seen from Mars at time t. 
% Let E and M be the location of Earth and Mars at time t.
% mu is the sine of angle RS-M-E.

  E = Earth(t);
  M = Mars(t);
  S = [0;0];
  rSun = .695500;
 
  [LS,RS] = TanPoints(M,S,rSun);
  
  u = RS-M; u = u/norm(u);      % Mars-to-RS vector
  v = E-M; v = v/norm(v);       % Mars-to-LE vector
  mu = u(1)*v(2) - u(2)*v(1);   % Cross-prod of u and v