Constant GAMMA 

const GAMMA: f64 = 2.67522128e8;

Gyromagnetic ratio of the proton, in rad/s/T: the reference implementation’s value (getconstants.m), not crate::constants. Keep it.