🔧 Adjusts constants

This commit is contained in:
Daniel Svitan
2025-04-18 22:23:09 +02:00
parent 620bd8be51
commit 0e43e6433c
2 changed files with 71 additions and 4 deletions

View File

@@ -4,7 +4,7 @@ use physical_constants;
use std::{thread, time};
fn main() {
let f = 100_f64; // Hz
let f = 20_f64; // Hz
let time_delta = 1_f64 / f; // s
let speed_dial = 1_f64;
let pause = time::Duration::from_millis((time_delta * 1000_f64 / speed_dial) as u64);
@@ -18,21 +18,27 @@ fn main() {
let Fg = mass * physical_constants::STANDARD_ACCELERATION_OF_GRAVITY; // N
let mut Fd = 0_f64;
let mut a; // m/s^2
let mut a = 0_f64; // m/s^2
let mut v = 0_f64; // m/s
let mut d = 0_f64; // m
let mut elapsed = 0_f64;
let roof = 3_f64;
loop {
print!("t = {:.3} s Fg = {:.3} N Fd = {:.3} N a = {:.3} m/s^2 v = {:.3} m/s d = {:.3} m\n", elapsed, Fg, Fd, a, v, d);
a = (Fg - Fd) / mass;
v += a * time_delta;
d += v * time_delta;
print!("t = {:.3} s Fg = {:.3} N Fd = {:.3} N a = {:.3} m/s^2 v = {:.3} m/s d = {:.3} m\r", elapsed, Fg, Fd, a, v, d);
Fd = p * Cd * A * v.powf(2_f64) / 2_f64;
elapsed += time_delta;
if elapsed >= roof {
break;
}
thread::sleep(pause);
}
}