How can I pull stress from a symModel in XY fashion?

What I mean is I am running the following example:

openExample('pde/StressConcentrationInPlateWithCircularHoleExample')

I'll keep it simple: can someone explain how I extract the XY Data?

I am trying to write a function that will receive an XY position to test, and return the relevant values at that position in the solution plot, such as for Von Mises stress at that XY position. However the model doesn't seem to have a readily available vector or matrix of corresponding XY values, and 'XYdata' looks like a special call to model, so not necessarily stored in the format useful to what I am doing.

"At [4.6783, 10.4534], what is the stress in the X-direction?" etc.