Is there any way to establish roughly the global discretization error for the ode15s solver?

Hello, I am working with a large set of linear DAEs in the form M(t)*y'(t) = A*y + B, where M is a singular mass matrix, A is a regular matrix (the DAEs describe a large electrical circuit) and B is a non-zero vector. I use the ode15s solver to achieve a numerical solution and I want to compare the results for the different values of the entries - coefficients in the matrixes M, A and B. The question is whether I can neglect the discretization errors of the method itself with respect to the changes of the entries. Is there any way to establish, however roughly, the global discretization error of the ode15s for this problem?