Compute values at certain points

[s.h.tong]

I want to compute the value of u at several certain points. Where can I add the positions? Thanks

[julienG (Julien)]

cout << u(x1,y1) <<endl;