Function quad4_dshape 

pub fn quad4_dshape(xi: f64, eta: f64) -> [[f64; 2]; 4]

[dN_i/dxi, dN_i/deta] for the quadrilateral.