Multivariate Mittag-Leffler Solution for a Forced Fractional-Order Harmonic Oscillator

The harmonic oscillator is a fundamental physical-mathematical system that allows for the description of a variety of models in many fields of physics. Utilizing fractional derivatives instead of traditional derivatives enables the modeling of a more diverse array of behaviors. Furthermore, if the effect of the fractional derivative is applied to each of the terms of the differential equation, this will involve greater complexity in the description of the analytical solutions of the fractional differential equation. In this work, by using the Laplace method, the solutions to the multiple-term forced fractional harmonic oscillator are presented, described through multivariate Mittag-Leffler functions. Additionally, the cases of damped and undamped free fractional harmonic oscillators are addressed. Finally, through simulations, the effect of the fractional non-integer derivative is demonstrated, and the consistency of the result is verified when recovering the integer case.