diff --git a/stubs/PIL/ImageMath.pyi b/stubs/PIL/ImageMath.pyi
new file mode 100644
--- /dev/null
+++ b/stubs/PIL/ImageMath.pyi
@@ -0,0 +1,60 @@
+# Stubs for PIL.ImageMath (Python 3)
+#
+# NOTE: This dynamically typed stub was automatically generated by stubgen.
+
+import __builtin__
+from typing import Any, Optional
+
+builtins = __builtin__
+VERBOSE: int
+
+class _Operand:
+    im: Any = ...
+    def __init__(self, im: Any) -> None: ...
+    def apply(self, op: Any, im1: Any, im2: Optional[Any] = ..., mode: Optional[Any] = ...): ...
+    def __bool__(self): ...
+    __nonzero__: Any = ...
+    def __abs__(self): ...
+    def __pos__(self): ...
+    def __neg__(self): ...
+    def __add__(self, other: Any): ...
+    def __radd__(self, other: Any): ...
+    def __sub__(self, other: Any): ...
+    def __rsub__(self, other: Any): ...
+    def __mul__(self, other: Any): ...
+    def __rmul__(self, other: Any): ...
+    def __truediv__(self, other: Any): ...
+    def __rtruediv__(self, other: Any): ...
+    def __mod__(self, other: Any): ...
+    def __rmod__(self, other: Any): ...
+    def __pow__(self, other: Any): ...
+    def __rpow__(self, other: Any): ...
+    __div__: Any = ...
+    __rdiv__: Any = ...
+    def __invert__(self): ...
+    def __and__(self, other: Any): ...
+    def __rand__(self, other: Any): ...
+    def __or__(self, other: Any): ...
+    def __ror__(self, other: Any): ...
+    def __xor__(self, other: Any): ...
+    def __rxor__(self, other: Any): ...
+    def __lshift__(self, other: Any): ...
+    def __rshift__(self, other: Any): ...
+    def __eq__(self, other: Any): ...
+    def __ne__(self, other: Any): ...
+    def __lt__(self, other: Any): ...
+    def __le__(self, other: Any): ...
+    def __gt__(self, other: Any): ...
+    def __ge__(self, other: Any): ...
+
+def imagemath_int(self): ...
+def imagemath_float(self): ...
+def imagemath_equal(self, other: Any): ...
+def imagemath_notequal(self, other: Any): ...
+def imagemath_min(self, other: Any): ...
+def imagemath_max(self, other: Any): ...
+def imagemath_convert(self, mode: Any): ...
+
+ops: Any
+
+def eval(expression: Any, _dict: Any = ..., **kw: Any): ...