pub fn tri3_shape(xi: f64, eta: f64) -> [f64; 3]
N_i(xi, eta) for the triangle: 1 - xi - eta, xi, eta.
N_i(xi, eta)
1 - xi - eta
xi
eta