The fractional oscillation equation with two fractional derivative terms in the sense of Caputo, where the orders $$\alpha$$
α
and $$\beta$$
β
satisfy $$1<\alpha \le 2$$
1
<
α
≤
2
and $$0<\beta \le 1$$
0
<
β
≤
1
, is investigated, and the unit step response and two initial value responses are obtained in two forms by using different methods of inverse Laplace transform. The first method yields series solutions with the nonnegative powers of t, which converge fast for small t. The second method is our emphasis, where the complex path integral formula of the inverse Laplace transform is used. In order to determine singularities of integrand we first seek for the roots of the characteristic equation, which is a transcendental equation with four parameters, two coefficients and two noninteger power exponents. The existence conditions and properties of the roots on the principal Riemann surface are given. Based on the results on the characteristic equation, we derive these responses as a sum of a classical exponentially damped oscillation, which vanishes in an indicated case, and an infinite integral of the Laplace type, which converges fast for large t, with a steady component in the unit step response. Asymptotic behaviors of solutions for large t are derived as algebraic decays in negative power laws characterized by the orders $$\alpha$$
α
and $$\beta$$
β
. The fractional system exhibits a magical transition from oscillation to monotonic decay in negative power law.