def fcts(self): h = sympy.Matrix([r**2 * z]) # Create an isotropic sigma vector Sig = sympy.Matrix([200/(sympy.pi)*(r*z)**2]) return h, Sig