fn main() {
use bijli::charge::{self, PointCharge};
use bijli::field;
use bijli::maxwell;
use bijli::wave;
let e = field::electric_field_point_charge(
1e-6, [0.0, 0.0, 0.0],
[1.0, 0.0, 0.0],
)
.unwrap();
println!("E field at 1m from 1μC: {:.2} V/m", e.magnitude());
let q1 = PointCharge::new(1e-6, 1.0, [0.0, 0.0, 0.0], [0.0; 3]).unwrap();
let q2 = PointCharge::new(-1e-6, 1.0, [0.5, 0.0, 0.0], [0.0; 3]).unwrap();
let f = charge::coulomb_force(&q1, &q2).unwrap();
println!("Coulomb force: {:.4} N", f.magnitude());
let c = maxwell::wave_speed(field::EPSILON_0, field::MU_0).unwrap();
println!("Speed of light: {:.0} m/s", c);
let intensity = wave::plane_wave_intensity(1000.0);
println!("Plane wave intensity (E₀=1kV/m): {:.2} W/m²", intensity);
}