Hyers–Ulam stability for linear quaternion-valued differential equations with constant coefficient