May I ask if Mathematica provides a trapezoidal method instruction for solving differential equations? I checked the reference materials and for NDSolve, it seems that only explicit Euler, implicit Euler, midpoint method, etc. are provided, but it seems that the trapezoidal method is not provided...