diff --git a/p5/numerical_integration/main.m b/p5/numerical_integration/main.m
new file mode 100644
index 0000000..e8fee82
--- /dev/null
+++ b/p5/numerical_integration/main.m
@@ -0,0 +1,19 @@
+T0 = 77.36
+T1 = 293.15
+
+CuA = 17.72891;
+CuB = 28.09870;
+CuC = -31.25289;
+CuD = 13.97243;
+CuE = 0.068611;
+f = @(t) CuA + CuB*t + CuC*t^2 + CuD*t^3 + CuE/(t^2);
+
+AlA = 28.08920;
+AlB = -5.414849;
+AlC = 8.560423;
+AlD = 3.427370;
+AlE = -0.277375;
+g = @(t) AlA + AlB*t + AlC*t^2 + AlD*t^3 + AlE/(t^2);
+
+quad(f, T0, T1)
+quad(g, T0, T1)
diff --git a/p5/numerical_integration/main2.m b/p5/numerical_integration/main2.m
new file mode 100644
index 0000000..cc89db4
--- /dev/null
+++ b/p5/numerical_integration/main2.m
@@ -0,0 +1,12 @@
+# 0 -> 6.25, 50 -> 16.1
+# x -> 6.25 + ((16.1 - 6.25)/50)*x
+# 27.36 -> 11.63992
+
+# 0 -> 24.0, 50 -> 24.5
+# x -> 24 + ((24.5 - 24)/50)*x
+# 43.15 -> 24.4315
+
+x = [77.36, 100, 150, 200, 250, 293.15];
+y = [11.63992, 16.1, 20.5, 22.8, 24.0, 24.4315];
+
+trapz(x, y)/63.5463