pub fn exponential_solution<F: IntegrateFloat>( amplitude: F, decay_rate: F, ) -> ExponentialSolution<F>