feat: add awful natuurkunde code

Author: mmvanheusden

Cargo.toml

@@ -30,6 +30,7 @@ edition = "2021"
 
 [dependencies]
 colored = "2.1.0"
+plotters = "0.3.5"
 rand = "0.8.5"
 
 
@@ -88,3 +89,7 @@ path = "Experiments/Minion/Minion.rs"
 [[bin]]
 name = "ErrorHandling"
 path = "Experiments/ErrorHandling/ErrorHandling.rs"
+
+[[bin]]
+name = "NatuurkundeModellen"
+path = "Experiments/NatuurkundeModellen/NatuurkundeModellen.rs"

Experiments/NatuurkundeModellen/NatuurkundeModellen.rs

@@ -0,0 +1,24 @@
+#![allow(non_snake_case)]
+mod Util;
+
+/// This is my awful code for school. Please don't look at it
+fn main() {
+    let mut punten: Vec<(f64, f64)> = vec![]; // We store the x,y points here
+    const HOOGTE: f32 = 1.0e3; // The height the kogel is dropped at. in metres
+    let mut DELTA_HOOGTE: f32 = 1.0; // the difference in height
+    const GRAVITATIECONSTANTE_NL: f32 = 9.81; // Gravitatieconstante
+    const MASSA: f32 = 2.0; // mass of kogel in kg
+
+    while DELTA_HOOGTE >= 0.0 {
+        println!(
+            "x: {DELTA_HOOGTE}, y: {}",
+            Util::kogelSnelheidBerekenen(HOOGTE, DELTA_HOOGTE, MASSA, GRAVITATIECONSTANTE_NL)
+        );
+        punten.push((
+            DELTA_HOOGTE as f64,
+            Util::kogelSnelheidBerekenen(HOOGTE, DELTA_HOOGTE, MASSA, GRAVITATIECONSTANTE_NL),
+        ));
+        DELTA_HOOGTE -= 0.001; // each step is 0.001.
+    }
+    Util::plot(punten).expect("ERROR!");
+}

Experiments/NatuurkundeModellen/Util.rs

@@ -0,0 +1,39 @@
+use plotters::prelude::*;
+
+/// Creates a plot.
+pub fn plot(lijnen: Vec<(f64, f64)>) -> Result<(), Box<dyn std::error::Error>> {
+    let root = BitMapBackend::new("plot.png", (640, 480)).into_drawing_area();
+    root.fill(&WHITE)?;
+    let mut chart = ChartBuilder::on(&root)
+        .caption("Kogelsnelheid", ("sans-serif", 30).into_font())
+        .x_label_area_size(40)
+        .y_label_area_size(40)
+        .build_cartesian_2d(0.0..1.0, 0.0..4.5)?;
+
+    chart
+        .configure_mesh()
+        .x_desc("Hoogteverschil in m")
+        .y_desc("Snelheid in m/s^2")
+        .draw()?;
+
+    chart
+        .draw_series(LineSeries::new(lijnen, &RED))?
+        .legend(|(x, y)| PathElement::new(vec![(x, y), (x + 20, y)], &RED));
+
+    root.present()?;
+
+    Ok(())
+}
+
+/// Calculates kogel snelheid
+pub fn kogelSnelheidBerekenen(
+    mut hoogte: f32,
+    deltaH: f32,
+    massa: f32,
+    valversnelling: f32,
+) -> f64 {
+    hoogte -= deltaH;
+    let zwaarte_energie = massa * valversnelling * hoogte;
+    let kinetische_energie = 19620.0 - zwaarte_energie;
+    (2f32 * kinetische_energie / massa).sqrt() as f64
+}