pub type FnTime = fn(t: f64) -> f64;
Defines a function of time that returns f64 (e.g., to calculate boundary condition values)