Order of bvp5c - Documentation or Implementation error?

I'm writing an overview about the use of MATLAB's bvp5c boundary value problem solver.

In the literature (for example [1], pg 36), it is stated that the Lobatto IIIa methods are order 2s-2 for an s stage implementation.

At the bottom of the pages of MATLAB's documentation for bvp5c [2] and bvp4c, it states that they implement the Lobatto IIIa 4 and 3 stage methods respectively.

I would expect, therefore, that bvp5c is 6th order and bvp4c is 4th order. However, bvp5c is stated in the documentation as 5th order.

Does anyone have any idea why or if the documentation just has an error?

Thanks

[1] Geometric Numerical Integration: Hairer et al. [2] bvp5c Documentation

[Torsten]

At least for IVPs, order of global error = order of local error - 1. Maybe that's similar for BVPs, but I'm not sure.

Best wishes

Torsten.