find the center of a black hole and then divide by zero in rust

main.rs
const G: f64 = 6.67430e-11; // gravitational constant in m^3 kg^-1 s^-2
const M: f64 = 1.0e6 * 1.989e30; // mass of black hole in kg

// assume we have the position and velocity vectors of a nearby star
let r_star = Vector3::new(1.0, 0.0, 0.0); // in units of AU
let v_star = Vector3::new(0.0, 10.0, 0.0); // in units of AU/yr

// calculate the specific angular momentum of the star
let h_star = r_star.cross(&v_star); // in units of AU^2/yr
let h_scalar = h_star.magnitude(); // in units of AU^2/yr

// calculate the semi-major axis of the star's orbit
let a_star = 1.0 / (2.0 / r_star.magnitude() - v_star.dot(&v_star) / (G * M)); // in units of AU

// calculate the location of the black hole's center
let r_bh = (h_scalar.powi(2) / (G * M)).sqrt(); // in units of AU

// divide by zero approximation
let r_divzero = r_bh / 0.0f64; // returns f64::INFINITY

println!("The center of the black hole is located at ({}, {}, {})", r_bh.x, r_bh.y, r_bh.z);
954 chars
22 lines