Lab fenòmens: afegir informe i ultims canvis al codi
Change-Id: I7f01665996295625fc390f9382dd0ba1c9882bbc
diff --git a/quad10/fenomens/lab/p4/cv_via_derivative.cc b/quad10/fenomens/lab/p4/cv_via_derivative.cc
new file mode 100644
index 0000000..c5308a2
--- /dev/null
+++ b/quad10/fenomens/lab/p4/cv_via_derivative.cc
@@ -0,0 +1,29 @@
+#include <iostream>
+#include <vector>
+using namespace std;
+
+// Fitxer que calcula a partir d'un fitxer de resultats del programa mc-2 els
+// valors de d<e>/dT amb derivació numèrica.
+
+struct Row {
+ long double T;
+ long double E;
+};
+
+int main() {
+ vector<Row> data;
+ long double t; // Trash
+ long double T;
+ long double E;
+ while (cin >> t >> T >> E >> t >> t >> t >> t >> t >> t >> t >> t) {
+ Row r;
+ r.T = T;
+ r.E = E;
+ data.push_back(r);
+ }
+
+ for (int i = 0; i < (int)data.size() - 1; ++i) {
+ long double derivada = (data[i + 1].E - data[i].E)/(data[i + 1].T - data[i].T);
+ cout << data[i].T << " " << derivada << endl;
+ }
+}