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)