function beta = calc_beta_simple(params,V,z)
% params: [kf, zf, kr, zr];
R = 8.314;
T = 310;
F = 96485;

beta = params(2)*exp(z*F*V/R/T);

end