Function py_wilson_hr_simulate 

pub(crate) fn py_wilson_hr_simulate<'py>(
    py: Python<'py>,
    v0: f64,
    r0: f64,
    tau_r: f64,
    v_peak: f64,
    dt: f64,
    n_steps: usize,
    current: f64,
) -> (Bound<'py, PyArray1<f64>>, i64, f64, f64)

N-step Wilson (1999) polynomial cortical-neuron simulation.

Parity contract with sc_neurocore.neurons.models.wilson_hr.WilsonHRNeuron.simulate: for the same parameters and constant input the returned v trace (already hard-reset to -0.7 on spiking steps), spike count, and final (v, r) state are bit-identical to the Python RK4 reference (the right-hand side is exact polynomial arithmetic — no transcendental functions).