function y=arrhen(x);
% file name: arrhen.m
% the Arrhenius rate function
global gamma delta;
y=delta*exp(gamma-gamma./x);


