diff --git a/Maths/SimpsonRule.py b/Maths/SimpsonRule.py
index 51b5ed1e4..091c86c17 100644
--- a/Maths/SimpsonRule.py
+++ b/Maths/SimpsonRule.py
@@ -8,6 +8,8 @@ method 2:
 "Simpson Rule"
 
 '''
+from __future__ import print_function
+
 
 def method_2(boundary, steps):
 # "Simpson Rule"
@@ -41,7 +43,7 @@ def main():
 	steps = 10.0		#define number of steps or resolution
 	boundary = [a, b]	#define boundary of integration
 	y = method_2(boundary, steps)
-	print 'y = {0}'.format(y)
+	print('y = {0}'.format(y))
 
 if __name__ == '__main__':
         main()
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()