A numerical procedure of generating the matrices of any finite element by its digital nomenclature code dncm [1,2] is proposed. Digit d is the dimension, i.e. the number of arguments x,y, ... used in the interpolation polynomial Y(x, ...) of the element; n is the number of nodes of the element; c is the number of coordinates per node, or the number of derivatives of the field variable Y per node; m is the vectorization multiplier showing how many polynomials are interpolated simultaneously. The dncm code itself only specifies the interpolation polynomials of the finite element. The physical properties of the element should additionally be given by a variation functional which employs, e.g., strain energy. In most cases, the form of the functional can be also identified by dncm code.
