I have some input parameters:
alpha[r]:= (650250000 r^3 (0.9989 - 0.449505 (1 - Tanh[(-1450 + r)/1160])))/22468879468420441 beta[r]:= (650250000 r^2 (1/1000 + 0.