class ElementTriArgyris(ElementH2): nodal_dofs = 6 facet_dofs = 1 dim = 2 maxdeg = 5 def gdof(self, u, du, ddu, v, e, n, i): if i < 18: j = i % 6 k = int(i/6) if j == 0: return u(*v[k]) elif j == 1: return du[0](*v[k]) elif j == 2: return du[1](*v[k]…

""" Author: kinnala Visualise the Argyris basis functions. """ from skfem import * import matplotlib.pyplot as plt m = MeshTri.init_sqsymmetric() map = MappingAffine(m) e = ElementTriArgyris() ib = InteriorBasis(m, e, map, 5)…