Correct option is (b) N1=c(1-ξ)(1-η)(1+ξ+η)
The explanation: The shape function is the function which interpolates the solution between the discrete values obtained at the mesh nodes. Therefore, appropriate functions have to be used and, as already mentioned, low order polynomials are typically chosen as shape functions. In this work linear shape functions are used.