Is there a way to know all dependent variables in across the whole domain at a particular step when using bvp4c?

The derivative function for bvp4c, in my understanding, acts like an ODE solver (eg. ode45). It will pass in the current position in the domain and the current best guess of what the dependent variables are at that position. My system of PDEs contains integrals across the domain (requiring the dependent variables at every place in the domain) and was wondering if bvp4c has the capability to handle this? If not are there suggestions out there of what could potentially handle this better? I was thinking potentially a non-linear solver like fsolve with finite difference matrices to handle spatial derivatives.