Numerical solution of the well-known Bagley-Torvik equation is considered. The fractional-order derivative in the equation is converted, approximately, to ordinary-order derivatives up to second order. Approximated Bagley-Torvik equation is obtained using finite number of terms from the infinite series of integer-order derivatives expansion for the Riemann–Liouville fractional derivative. The Bagley-Torvik equation is a second-order differential equation with constant coefficients. The derived equation, by considering only the first three terms from the infinite series to become a second-order ordinary differential equation with variable coefficients, is numerically solved after it is transformed into a system of first-order ordinary differential equations. The approximation of fractional-order derivative and the order of the truncated error are illustrated through some examples. Comparison between our result and exact analytical solution are made by considering an example with known analytical solution to show the preciseness of our proposed approach.
Numerical Solution of the Bagley-Torvik Equation Using the Integer-Order Derivatives Expansion
