Function acceleration 

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

Computes acceleration for simple harmonic motion using a(t) = -Aω² * cos(ωt + φ).