def perform(self, node, inputs, output_storage): x = inputs[0] z = output_storage[0] z[0] = sp.polygamma(self.n, x)