solution of higher order differential equations matrix method