Code matches with the math spec. Can you verify that truth data is correct?
Definition: The Beta angle, β, is defined as the angle between the orbit normal vector, and the vector from the celestial body to the sun.
ˆh = r⊕ × v⊕ / |r⊕ × v⊕k|
ˆrs⊕ = rs⊕ / |s⊕|
β = sin−1(ˆh · ˆrs⊕) (3.102)
• r⊕: Position vector of spacecraft with respect to celestial body, in the EarthMJ2000Eq system.
• v⊕: Velocity vector of spacecraft with respect to celestial body, in the EarthMJ2000Eq system.
• rs⊕: Position vector from celestial body, to the sun.