Is there a way to evaluate a solution of a PDE at one of the domain boundaries? E.g. if u is the solution that I’ve just obtained, is there a simple way to then obtain u(gamma) where gamma is a boundary of the domain on which u is defined? So then u(gamma) should only contain the values of u on gamma.

Hope this is clear, please say if it isn’t and thanks!

Hi Frederic, thanks for your answer! This seems close, although I need the function evaluated at a specific border, rather than the entire boundary. Apologies, I used the word ‘boundary’ in my question when I actually meant to say ‘border’. I believe that Chris’ answer above should work in my code though!