diff --git a/nifty5/linearization.py b/nifty5/linearization.py
index 17675aade38d5507fa958823658de680d574b805..566c6801be7d8de9a3962a5f80804157552b2cf9 100644
--- a/nifty5/linearization.py
+++ b/nifty5/linearization.py
@@ -1,3 +1,4 @@
+
 from __future__ import absolute_import, division, print_function
 
 import numpy as np
@@ -104,7 +105,7 @@ class Linearization(object):
     def __pow__(self, power):
         if not np.isscalar(power):
             return NotImplemented
-        return self.new(self._val**power, (self._jac**(power-1)).scale(power))
+        return self.new(self._val**power, makeOp(self._val**(power-1)).scale(power)(self._jac))
 
     def inverse(self):
         return self.new(1./self._val, makeOp(-1./(self._val**2))(self._jac))