Function py_pernarowski_simulate 

pub(crate) fn py_pernarowski_simulate<'py>(
    py: Python<'py>,
    v0: f64,
    w0: f64,
    z0: f64,
    alpha: f64,
    beta: f64,
    eps1: f64,
    eps2: f64,
    gamma: f64,
    dt: f64,
    v_threshold: f64,
    n_steps: usize,
    current: f64,
) -> (Bound<'py, PyArray1<f64>>, i64, f64, f64, f64)

N-step Pernarowski (1994) pancreatic beta-cell burster simulation.

Parity contract with sc_neurocore.neurons.models.pernarowski.PernarowskiNeuron.simulate: for the same parameters and constant input the returned v trace, upward-crossing spike count, and final (v, w, z) state are bit-identical to the Python RK4 reference (the cubic uses v.powi(3) = v*v*v, matching the Python v*v*v; no transcendental functions).