%ind_sha(4)%

shape function for boundary points