Function velocity 

pub fn velocity(
    amplitude: f64,
    angular_frequency: f64,
    time: f64,
    phase: f64,
) -> Option<f64>

Computes velocity for simple harmonic motion using v(t) = -Aω * sin(ωt + φ).