1function [n, dxi, deta] = plate_fem_shape(xi, eta)
2% PLATE_FEM_SHAPE Bilinear shape functions of the four-node element.
3%
4% [n, dxi, deta] = plate_fem_shape(xi, eta) at a point of the reference
5% square [-1, 1]^2: the four functions, ordered anticlockwise from the
6% corner (-1, -1), and their derivatives with respect to xi and eta.
8n = 0.25 * [(1 - xi) * (1 - eta), (1 + xi) * (1 - eta), (1 + xi) * (1 + eta), (1 - xi) * (1 + eta)];
9dxi = 0.25 * [-(1 - eta), 1 - eta, 1 + eta, -(1 + eta)];
10deta = 0.25 * [-(1 - xi), -(1 + xi), 1 + xi, 1 - xi];
11end