diff --git a/Maths/TrapezoidalRule.py b/Maths/TrapezoidalRule.py
index 2ad857390..52310c1ed 100644
--- a/Maths/TrapezoidalRule.py
+++ b/Maths/TrapezoidalRule.py
@@ -7,6 +7,7 @@ method 1:
 "extended trapezoidal rule"
 
 '''
+from __future__ import print_function
 
 def method_1(boundary, steps):
 # "extended trapezoidal rule"
@@ -39,7 +40,7 @@ def main():
 	steps = 10.0		#define number of steps or resolution	
 	boundary = [a, b]	#define boundary of integration
 	y = method_1(boundary, steps)
-	print 'y = {0}'.format(y)
+	print('y = {0}'.format(y))
 
 if __name__ == '__main__':
 	main()