rpim shape function

I am solving a simple solid mechanics problem by meshless method.  use RPIM to  calculate shape function.

when i calculate the derivatives of shape function, i found the derivative is not close to zero at the compute point, when the point is on the boundry of problem domain. 

 

My question is, should the derivative be zero if the node is on the boundry?