--- theano/scalar/basic_scipy.py.orig	2019-01-15 14:13:57.000000000 -0700
+++ theano/scalar/basic_scipy.py	2019-08-05 13:12:02.587477090 -0600
@@ -596,7 +596,7 @@ class GammaIncC(BinaryScalarOp):
         if node.inputs[0].type in float_types:
             dtype = 'npy_' + node.outputs[0].dtype
             return """%(z)s =
-                (%(dtype)s) gammaQ(%(k)s, %(x)s);""" % locals()
+                (%(dtype)s) GammaQ(%(k)s, %(x)s);""" % locals()
         raise NotImplementedError('only floatingpoint is implemented')
 
     def __eq__(self, other):