d/dα ln[(gamma*10^12 e^(α (3.64568×10^6)^β (-2^(-4.7619×10^-6 b β)) s^(-β)))/s]