The paper describes different approaches to generalize the trapezoidal method to fractional differential equations. We analyze the main theoretical properties and we discuss computational aspects to implement efficient algorithms. Numerical experiments are provided to illustrate potential and limitations of the different methods under investigation. [ABSTRACT FROM AUTHOR]