pub fn strinit(q: f64) -> Result<(Vec3, Vec3), RocheError>Expand description
strinit sets a particle just inside the L1 point with the correct velocity as given in Lubow and Shu.
Arguments:
q: mass ratio = M2/M1
Returns:
- start position
- start velocity