It seems to me that this formula is correct: L = (UrToaESUN*cos(theta))/PI; Watch the discussion here: