# here "x" is the distance from the center of the atom in Angstrom
# B is the Debye-Waller-Ott "temperature factor" (same units as a Protein Data Bank file)
# retun values of the Ee_ff(x,B) are electron density (in electrons per cubic Angstrom)

# keep exp() from blowing up at large distances
safexp(x) = (x>700 ? 1e300 : (x<-700 ? 0 : exp(x)))

# single electron (approaches a delta-function at B=0)
Ano_ff(x,B) = \
  +1.0*(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
H_c = 0.003038;
H_a1 = 0.493002;H_a2 = 0.322912;H_a3 = 0.140191;H_a4 = 0.040810;
H_b1 = 10.510900;H_b2 = 26.125700;H_b3 = 3.142360;H_b4 = 57.799698;
H_ff(x,B) = \
  +H_a1*(4*pi/(H_b1+B))**1.5*safexp(-4*pi**2/(H_b1+B)*x*x) \
  +H_a2*(4*pi/(H_b2+B))**1.5*safexp(-4*pi**2/(H_b2+B)*x*x) \
  +H_a3*(4*pi/(H_b3+B))**1.5*safexp(-4*pi**2/(H_b3+B)*x*x) \
  +H_a4*(4*pi/(H_b4+B))**1.5*safexp(-4*pi**2/(H_b4+B)*x*x) \
  +H_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=2
H_minus_1_c = 0.002389;
H_minus_1_a1 = 0.897661;H_minus_1_a2 = 0.565616;H_minus_1_a3 = 0.415815;H_minus_1_a4 = 0.116973;
H_minus_1_b1 = 53.136799;H_minus_1_b2 = 15.187000;H_minus_1_b3 = 186.575989;H_minus_1_b4 = 3.567090;
H_minus_1_ff(x,B) = \
  +H_minus_1_a1*(4*pi/(H_minus_1_b1+B))**1.5*safexp(-4*pi**2/(H_minus_1_b1+B)*x*x) \
  +H_minus_1_a2*(4*pi/(H_minus_1_b2+B))**1.5*safexp(-4*pi**2/(H_minus_1_b2+B)*x*x) \
  +H_minus_1_a3*(4*pi/(H_minus_1_b3+B))**1.5*safexp(-4*pi**2/(H_minus_1_b3+B)*x*x) \
  +H_minus_1_a4*(4*pi/(H_minus_1_b4+B))**1.5*safexp(-4*pi**2/(H_minus_1_b4+B)*x*x) \
  +H_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=2
He_c = 0.006400;
He_a1 = 0.873400;He_a2 = 0.630900;He_a3 = 0.311200;He_a4 = 0.178000;
He_b1 = 9.103700;He_b2 = 3.356800;He_b3 = 22.927601;He_b4 = 0.982100;
He_ff(x,B) = \
  +He_a1*(4*pi/(He_b1+B))**1.5*safexp(-4*pi**2/(He_b1+B)*x*x) \
  +He_a2*(4*pi/(He_b2+B))**1.5*safexp(-4*pi**2/(He_b2+B)*x*x) \
  +He_a3*(4*pi/(He_b3+B))**1.5*safexp(-4*pi**2/(He_b3+B)*x*x) \
  +He_a4*(4*pi/(He_b4+B))**1.5*safexp(-4*pi**2/(He_b4+B)*x*x) \
  +He_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=3
Li_c = 0.037700;
Li_a1 = 1.128200;Li_a2 = 0.750800;Li_a3 = 0.617500;Li_a4 = 0.465300;
Li_b1 = 3.954600;Li_b2 = 1.052400;Li_b3 = 85.390503;Li_b4 = 168.261002;
Li_ff(x,B) = \
  +Li_a1*(4*pi/(Li_b1+B))**1.5*safexp(-4*pi**2/(Li_b1+B)*x*x) \
  +Li_a2*(4*pi/(Li_b2+B))**1.5*safexp(-4*pi**2/(Li_b2+B)*x*x) \
  +Li_a3*(4*pi/(Li_b3+B))**1.5*safexp(-4*pi**2/(Li_b3+B)*x*x) \
  +Li_a4*(4*pi/(Li_b4+B))**1.5*safexp(-4*pi**2/(Li_b4+B)*x*x) \
  +Li_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=2
Li_plus_1_c = 0.016700;
Li_plus_1_a1 = 0.696800;Li_plus_1_a2 = 0.788800;Li_plus_1_a3 = 0.341400;Li_plus_1_a4 = 0.156300;
Li_plus_1_b1 = 4.623700;Li_plus_1_b2 = 1.955700;Li_plus_1_b3 = 0.631600;Li_plus_1_b4 = 10.095300;
Li_plus_1_ff(x,B) = \
  +Li_plus_1_a1*(4*pi/(Li_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Li_plus_1_b1+B)*x*x) \
  +Li_plus_1_a2*(4*pi/(Li_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Li_plus_1_b2+B)*x*x) \
  +Li_plus_1_a3*(4*pi/(Li_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Li_plus_1_b3+B)*x*x) \
  +Li_plus_1_a4*(4*pi/(Li_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Li_plus_1_b4+B)*x*x) \
  +Li_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=4
Be_c = 0.038500;
Be_a1 = 1.591900;Be_a2 = 1.127800;Be_a3 = 0.539100;Be_a4 = 0.702900;
Be_b1 = 43.642700;Be_b2 = 1.862300;Be_b3 = 103.483002;Be_b4 = 0.542000;
Be_ff(x,B) = \
  +Be_a1*(4*pi/(Be_b1+B))**1.5*safexp(-4*pi**2/(Be_b1+B)*x*x) \
  +Be_a2*(4*pi/(Be_b2+B))**1.5*safexp(-4*pi**2/(Be_b2+B)*x*x) \
  +Be_a3*(4*pi/(Be_b3+B))**1.5*safexp(-4*pi**2/(Be_b3+B)*x*x) \
  +Be_a4*(4*pi/(Be_b4+B))**1.5*safexp(-4*pi**2/(Be_b4+B)*x*x) \
  +Be_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=2
Be_plus_2_c = -6.109200;
Be_plus_2_a1 = 6.260300;Be_plus_2_a2 = 0.884900;Be_plus_2_a3 = 0.799300;Be_plus_2_a4 = 0.164700;
Be_plus_2_b1 = 0.002700;Be_plus_2_b2 = 0.831300;Be_plus_2_b3 = 2.275800;Be_plus_2_b4 = 5.114600;
Be_plus_2_ff(x,B) = \
  +Be_plus_2_a1*(4*pi/(Be_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Be_plus_2_b1+B)*x*x) \
  +Be_plus_2_a2*(4*pi/(Be_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Be_plus_2_b2+B)*x*x) \
  +Be_plus_2_a3*(4*pi/(Be_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Be_plus_2_b3+B)*x*x) \
  +Be_plus_2_a4*(4*pi/(Be_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Be_plus_2_b4+B)*x*x) \
  +Be_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=5
B_c = -0.193200;
B_a1 = 2.054500;B_a2 = 1.332600;B_a3 = 1.097900;B_a4 = 0.706800;
B_b1 = 23.218500;B_b2 = 1.021000;B_b3 = 60.349800;B_b4 = 0.140300;
B_ff(x,B) = \
  +B_a1*(4*pi/(B_b1+B))**1.5*safexp(-4*pi**2/(B_b1+B)*x*x) \
  +B_a2*(4*pi/(B_b2+B))**1.5*safexp(-4*pi**2/(B_b2+B)*x*x) \
  +B_a3*(4*pi/(B_b3+B))**1.5*safexp(-4*pi**2/(B_b3+B)*x*x) \
  +B_a4*(4*pi/(B_b4+B))**1.5*safexp(-4*pi**2/(B_b4+B)*x*x) \
  +B_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=6
C_c = 0.215600;
C_a1 = 2.310000;C_a2 = 1.020000;C_a3 = 1.588600;C_a4 = 0.865000;
C_b1 = 20.843899;C_b2 = 10.207500;C_b3 = 0.568700;C_b4 = 51.651199;
C_ff(x,B) = \
  +C_a1*(4*pi/(C_b1+B))**1.5*safexp(-4*pi**2/(C_b1+B)*x*x) \
  +C_a2*(4*pi/(C_b2+B))**1.5*safexp(-4*pi**2/(C_b2+B)*x*x) \
  +C_a3*(4*pi/(C_b3+B))**1.5*safexp(-4*pi**2/(C_b3+B)*x*x) \
  +C_a4*(4*pi/(C_b4+B))**1.5*safexp(-4*pi**2/(C_b4+B)*x*x) \
  +C_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=6
Cv_c = 0.286977;
Cv_a1 = 2.260690;Cv_a2 = 1.561650;Cv_a3 = 1.050750;Cv_a4 = 0.839259;
Cv_b1 = 22.690701;Cv_b2 = 0.656665;Cv_b3 = 9.756180;Cv_b4 = 55.594898;
Cv_ff(x,B) = \
  +Cv_a1*(4*pi/(Cv_b1+B))**1.5*safexp(-4*pi**2/(Cv_b1+B)*x*x) \
  +Cv_a2*(4*pi/(Cv_b2+B))**1.5*safexp(-4*pi**2/(Cv_b2+B)*x*x) \
  +Cv_a3*(4*pi/(Cv_b3+B))**1.5*safexp(-4*pi**2/(Cv_b3+B)*x*x) \
  +Cv_a4*(4*pi/(Cv_b4+B))**1.5*safexp(-4*pi**2/(Cv_b4+B)*x*x) \
  +Cv_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=7
N_c = -11.528999;
N_a1 = 12.212600;N_a2 = 3.132200;N_a3 = 2.012500;N_a4 = 1.166300;
N_b1 = 0.005700;N_b2 = 9.893300;N_b3 = 28.997499;N_b4 = 0.582600;
N_ff(x,B) = \
  +N_a1*(4*pi/(N_b1+B))**1.5*safexp(-4*pi**2/(N_b1+B)*x*x) \
  +N_a2*(4*pi/(N_b2+B))**1.5*safexp(-4*pi**2/(N_b2+B)*x*x) \
  +N_a3*(4*pi/(N_b3+B))**1.5*safexp(-4*pi**2/(N_b3+B)*x*x) \
  +N_a4*(4*pi/(N_b4+B))**1.5*safexp(-4*pi**2/(N_b4+B)*x*x) \
  +N_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=8
O_c = 0.250800;
O_a1 = 3.048500;O_a2 = 2.286800;O_a3 = 1.546300;O_a4 = 0.867000;
O_b1 = 13.277100;O_b2 = 5.701100;O_b3 = 0.323900;O_b4 = 32.908897;
O_ff(x,B) = \
  +O_a1*(4*pi/(O_b1+B))**1.5*safexp(-4*pi**2/(O_b1+B)*x*x) \
  +O_a2*(4*pi/(O_b2+B))**1.5*safexp(-4*pi**2/(O_b2+B)*x*x) \
  +O_a3*(4*pi/(O_b3+B))**1.5*safexp(-4*pi**2/(O_b3+B)*x*x) \
  +O_a4*(4*pi/(O_b4+B))**1.5*safexp(-4*pi**2/(O_b4+B)*x*x) \
  +O_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=9
O_minus_1_c = 21.941200;
O_minus_1_a1 = 4.191600;O_minus_1_a2 = 1.639690;O_minus_1_a3 = 1.526730;O_minus_1_a4 = -20.306999;
O_minus_1_b1 = 12.857300;O_minus_1_b2 = 4.172360;O_minus_1_b3 = 47.017899;O_minus_1_b4 = -0.014040;
O_minus_1_ff(x,B) = \
  +O_minus_1_a1*(4*pi/(O_minus_1_b1+B))**1.5*safexp(-4*pi**2/(O_minus_1_b1+B)*x*x) \
  +O_minus_1_a2*(4*pi/(O_minus_1_b2+B))**1.5*safexp(-4*pi**2/(O_minus_1_b2+B)*x*x) \
  +O_minus_1_a3*(4*pi/(O_minus_1_b3+B))**1.5*safexp(-4*pi**2/(O_minus_1_b3+B)*x*x) \
  +O_minus_1_a4*(4*pi/(O_minus_1_b4+B))**1.5*safexp(-4*pi**2/(O_minus_1_b4+B)*x*x) \
  +O_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=9
F_c = 0.277600;
F_a1 = 3.539200;F_a2 = 2.641200;F_a3 = 1.517000;F_a4 = 1.024300;
F_b1 = 10.282499;F_b2 = 4.294400;F_b3 = 0.261500;F_b4 = 26.147600;
F_ff(x,B) = \
  +F_a1*(4*pi/(F_b1+B))**1.5*safexp(-4*pi**2/(F_b1+B)*x*x) \
  +F_a2*(4*pi/(F_b2+B))**1.5*safexp(-4*pi**2/(F_b2+B)*x*x) \
  +F_a3*(4*pi/(F_b3+B))**1.5*safexp(-4*pi**2/(F_b3+B)*x*x) \
  +F_a4*(4*pi/(F_b4+B))**1.5*safexp(-4*pi**2/(F_b4+B)*x*x) \
  +F_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
F_minus_1_c = 0.653396;
F_minus_1_a1 = 3.632200;F_minus_1_a2 = 3.510570;F_minus_1_a3 = 1.260640;F_minus_1_a4 = 0.940706;
F_minus_1_b1 = 5.277560;F_minus_1_b2 = 14.735300;F_minus_1_b3 = 0.442258;F_minus_1_b4 = 47.343700;
F_minus_1_ff(x,B) = \
  +F_minus_1_a1*(4*pi/(F_minus_1_b1+B))**1.5*safexp(-4*pi**2/(F_minus_1_b1+B)*x*x) \
  +F_minus_1_a2*(4*pi/(F_minus_1_b2+B))**1.5*safexp(-4*pi**2/(F_minus_1_b2+B)*x*x) \
  +F_minus_1_a3*(4*pi/(F_minus_1_b3+B))**1.5*safexp(-4*pi**2/(F_minus_1_b3+B)*x*x) \
  +F_minus_1_a4*(4*pi/(F_minus_1_b4+B))**1.5*safexp(-4*pi**2/(F_minus_1_b4+B)*x*x) \
  +F_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
Ne_c = 0.351500;
Ne_a1 = 3.955300;Ne_a2 = 3.112500;Ne_a3 = 1.454600;Ne_a4 = 1.125100;
Ne_b1 = 8.404200;Ne_b2 = 3.426200;Ne_b3 = 0.230600;Ne_b4 = 21.718399;
Ne_ff(x,B) = \
  +Ne_a1*(4*pi/(Ne_b1+B))**1.5*safexp(-4*pi**2/(Ne_b1+B)*x*x) \
  +Ne_a2*(4*pi/(Ne_b2+B))**1.5*safexp(-4*pi**2/(Ne_b2+B)*x*x) \
  +Ne_a3*(4*pi/(Ne_b3+B))**1.5*safexp(-4*pi**2/(Ne_b3+B)*x*x) \
  +Ne_a4*(4*pi/(Ne_b4+B))**1.5*safexp(-4*pi**2/(Ne_b4+B)*x*x) \
  +Ne_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=11
Na_c = 0.676000;
Na_a1 = 4.762600;Na_a2 = 3.173600;Na_a3 = 1.267400;Na_a4 = 1.112800;
Na_b1 = 3.285000;Na_b2 = 8.842199;Na_b3 = 0.313600;Na_b4 = 129.423996;
Na_ff(x,B) = \
  +Na_a1*(4*pi/(Na_b1+B))**1.5*safexp(-4*pi**2/(Na_b1+B)*x*x) \
  +Na_a2*(4*pi/(Na_b2+B))**1.5*safexp(-4*pi**2/(Na_b2+B)*x*x) \
  +Na_a3*(4*pi/(Na_b3+B))**1.5*safexp(-4*pi**2/(Na_b3+B)*x*x) \
  +Na_a4*(4*pi/(Na_b4+B))**1.5*safexp(-4*pi**2/(Na_b4+B)*x*x) \
  +Na_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
Na_plus_1_c = 0.404000;
Na_plus_1_a1 = 3.256500;Na_plus_1_a2 = 3.936200;Na_plus_1_a3 = 1.399800;Na_plus_1_a4 = 1.003200;
Na_plus_1_b1 = 2.667100;Na_plus_1_b2 = 6.115300;Na_plus_1_b3 = 0.200100;Na_plus_1_b4 = 14.039000;
Na_plus_1_ff(x,B) = \
  +Na_plus_1_a1*(4*pi/(Na_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Na_plus_1_b1+B)*x*x) \
  +Na_plus_1_a2*(4*pi/(Na_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Na_plus_1_b2+B)*x*x) \
  +Na_plus_1_a3*(4*pi/(Na_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Na_plus_1_b3+B)*x*x) \
  +Na_plus_1_a4*(4*pi/(Na_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Na_plus_1_b4+B)*x*x) \
  +Na_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=12
Mg_c = 0.858400;
Mg_a1 = 5.420400;Mg_a2 = 2.173500;Mg_a3 = 1.226900;Mg_a4 = 2.307300;
Mg_b1 = 2.827500;Mg_b2 = 79.261101;Mg_b3 = 0.380800;Mg_b4 = 7.193700;
Mg_ff(x,B) = \
  +Mg_a1*(4*pi/(Mg_b1+B))**1.5*safexp(-4*pi**2/(Mg_b1+B)*x*x) \
  +Mg_a2*(4*pi/(Mg_b2+B))**1.5*safexp(-4*pi**2/(Mg_b2+B)*x*x) \
  +Mg_a3*(4*pi/(Mg_b3+B))**1.5*safexp(-4*pi**2/(Mg_b3+B)*x*x) \
  +Mg_a4*(4*pi/(Mg_b4+B))**1.5*safexp(-4*pi**2/(Mg_b4+B)*x*x) \
  +Mg_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
Mg_plus_2_c = 0.485300;
Mg_plus_2_a1 = 3.498800;Mg_plus_2_a2 = 3.837800;Mg_plus_2_a3 = 1.328400;Mg_plus_2_a4 = 0.849700;
Mg_plus_2_b1 = 2.167600;Mg_plus_2_b2 = 4.754200;Mg_plus_2_b3 = 0.185000;Mg_plus_2_b4 = 10.141100;
Mg_plus_2_ff(x,B) = \
  +Mg_plus_2_a1*(4*pi/(Mg_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Mg_plus_2_b1+B)*x*x) \
  +Mg_plus_2_a2*(4*pi/(Mg_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Mg_plus_2_b2+B)*x*x) \
  +Mg_plus_2_a3*(4*pi/(Mg_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Mg_plus_2_b3+B)*x*x) \
  +Mg_plus_2_a4*(4*pi/(Mg_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Mg_plus_2_b4+B)*x*x) \
  +Mg_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=13
Al_c = 1.115100;
Al_a1 = 6.420200;Al_a2 = 1.900200;Al_a3 = 1.593600;Al_a4 = 1.964600;
Al_b1 = 3.038700;Al_b2 = 0.742600;Al_b3 = 31.547199;Al_b4 = 85.088600;
Al_ff(x,B) = \
  +Al_a1*(4*pi/(Al_b1+B))**1.5*safexp(-4*pi**2/(Al_b1+B)*x*x) \
  +Al_a2*(4*pi/(Al_b2+B))**1.5*safexp(-4*pi**2/(Al_b2+B)*x*x) \
  +Al_a3*(4*pi/(Al_b3+B))**1.5*safexp(-4*pi**2/(Al_b3+B)*x*x) \
  +Al_a4*(4*pi/(Al_b4+B))**1.5*safexp(-4*pi**2/(Al_b4+B)*x*x) \
  +Al_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
Al_plus_3_c = 0.706786;
Al_plus_3_a1 = 4.174480;Al_plus_3_a2 = 3.387600;Al_plus_3_a3 = 1.202960;Al_plus_3_a4 = 0.528137;
Al_plus_3_b1 = 1.938160;Al_plus_3_b2 = 4.145530;Al_plus_3_b3 = 0.228753;Al_plus_3_b4 = 8.285240;
Al_plus_3_ff(x,B) = \
  +Al_plus_3_a1*(4*pi/(Al_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Al_plus_3_b1+B)*x*x) \
  +Al_plus_3_a2*(4*pi/(Al_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Al_plus_3_b2+B)*x*x) \
  +Al_plus_3_a3*(4*pi/(Al_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Al_plus_3_b3+B)*x*x) \
  +Al_plus_3_a4*(4*pi/(Al_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Al_plus_3_b4+B)*x*x) \
  +Al_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=14
Si_c = 1.140700;
Si_a1 = 6.291500;Si_a2 = 3.035300;Si_a3 = 1.989100;Si_a4 = 1.541000;
Si_b1 = 2.438600;Si_b2 = 32.333698;Si_b3 = 0.678500;Si_b4 = 81.693695;
Si_ff(x,B) = \
  +Si_a1*(4*pi/(Si_b1+B))**1.5*safexp(-4*pi**2/(Si_b1+B)*x*x) \
  +Si_a2*(4*pi/(Si_b2+B))**1.5*safexp(-4*pi**2/(Si_b2+B)*x*x) \
  +Si_a3*(4*pi/(Si_b3+B))**1.5*safexp(-4*pi**2/(Si_b3+B)*x*x) \
  +Si_a4*(4*pi/(Si_b4+B))**1.5*safexp(-4*pi**2/(Si_b4+B)*x*x) \
  +Si_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=14
Siv_c = 1.247070;
Siv_a1 = 5.662690;Siv_a2 = 3.071640;Siv_a3 = 2.624460;Siv_a4 = 1.393200;
Siv_b1 = 2.665200;Siv_b2 = 38.663399;Siv_b3 = 0.916946;Siv_b4 = 93.545799;
Siv_ff(x,B) = \
  +Siv_a1*(4*pi/(Siv_b1+B))**1.5*safexp(-4*pi**2/(Siv_b1+B)*x*x) \
  +Siv_a2*(4*pi/(Siv_b2+B))**1.5*safexp(-4*pi**2/(Siv_b2+B)*x*x) \
  +Siv_a3*(4*pi/(Siv_b3+B))**1.5*safexp(-4*pi**2/(Siv_b3+B)*x*x) \
  +Siv_a4*(4*pi/(Siv_b4+B))**1.5*safexp(-4*pi**2/(Siv_b4+B)*x*x) \
  +Siv_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=10
Si_plus_4_c = 0.746297;
Si_plus_4_a1 = 4.439180;Si_plus_4_a2 = 3.203450;Si_plus_4_a3 = 1.194530;Si_plus_4_a4 = 0.416530;
Si_plus_4_b1 = 1.641670;Si_plus_4_b2 = 3.437570;Si_plus_4_b3 = 0.214900;Si_plus_4_b4 = 6.653650;
Si_plus_4_ff(x,B) = \
  +Si_plus_4_a1*(4*pi/(Si_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Si_plus_4_b1+B)*x*x) \
  +Si_plus_4_a2*(4*pi/(Si_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Si_plus_4_b2+B)*x*x) \
  +Si_plus_4_a3*(4*pi/(Si_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Si_plus_4_b3+B)*x*x) \
  +Si_plus_4_a4*(4*pi/(Si_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Si_plus_4_b4+B)*x*x) \
  +Si_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=15
P_c = 1.114900;
P_a1 = 6.434500;P_a2 = 4.179100;P_a3 = 1.780000;P_a4 = 1.490800;
P_b1 = 1.906700;P_b2 = 27.157000;P_b3 = 0.526000;P_b4 = 68.164497;
P_ff(x,B) = \
  +P_a1*(4*pi/(P_b1+B))**1.5*safexp(-4*pi**2/(P_b1+B)*x*x) \
  +P_a2*(4*pi/(P_b2+B))**1.5*safexp(-4*pi**2/(P_b2+B)*x*x) \
  +P_a3*(4*pi/(P_b3+B))**1.5*safexp(-4*pi**2/(P_b3+B)*x*x) \
  +P_a4*(4*pi/(P_b4+B))**1.5*safexp(-4*pi**2/(P_b4+B)*x*x) \
  +P_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=16
S_c = 0.866900;
S_a1 = 6.905300;S_a2 = 5.203400;S_a3 = 1.437900;S_a4 = 1.586300;
S_b1 = 1.467900;S_b2 = 22.215099;S_b3 = 0.253600;S_b4 = 56.172001;
S_ff(x,B) = \
  +S_a1*(4*pi/(S_b1+B))**1.5*safexp(-4*pi**2/(S_b1+B)*x*x) \
  +S_a2*(4*pi/(S_b2+B))**1.5*safexp(-4*pi**2/(S_b2+B)*x*x) \
  +S_a3*(4*pi/(S_b3+B))**1.5*safexp(-4*pi**2/(S_b3+B)*x*x) \
  +S_a4*(4*pi/(S_b4+B))**1.5*safexp(-4*pi**2/(S_b4+B)*x*x) \
  +S_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=17
Cl_c = -9.557400;
Cl_a1 = 11.460400;Cl_a2 = 7.196400;Cl_a3 = 6.255600;Cl_a4 = 1.645500;
Cl_b1 = 0.010400;Cl_b2 = 1.166200;Cl_b3 = 18.519400;Cl_b4 = 47.778400;
Cl_ff(x,B) = \
  +Cl_a1*(4*pi/(Cl_b1+B))**1.5*safexp(-4*pi**2/(Cl_b1+B)*x*x) \
  +Cl_a2*(4*pi/(Cl_b2+B))**1.5*safexp(-4*pi**2/(Cl_b2+B)*x*x) \
  +Cl_a3*(4*pi/(Cl_b3+B))**1.5*safexp(-4*pi**2/(Cl_b3+B)*x*x) \
  +Cl_a4*(4*pi/(Cl_b4+B))**1.5*safexp(-4*pi**2/(Cl_b4+B)*x*x) \
  +Cl_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
Cl_minus_1_c = -16.378000;
Cl_minus_1_a1 = 18.291500;Cl_minus_1_a2 = 7.208400;Cl_minus_1_a3 = 6.533700;Cl_minus_1_a4 = 2.338600;
Cl_minus_1_b1 = 0.006600;Cl_minus_1_b2 = 1.171700;Cl_minus_1_b3 = 19.542400;Cl_minus_1_b4 = 60.448601;
Cl_minus_1_ff(x,B) = \
  +Cl_minus_1_a1*(4*pi/(Cl_minus_1_b1+B))**1.5*safexp(-4*pi**2/(Cl_minus_1_b1+B)*x*x) \
  +Cl_minus_1_a2*(4*pi/(Cl_minus_1_b2+B))**1.5*safexp(-4*pi**2/(Cl_minus_1_b2+B)*x*x) \
  +Cl_minus_1_a3*(4*pi/(Cl_minus_1_b3+B))**1.5*safexp(-4*pi**2/(Cl_minus_1_b3+B)*x*x) \
  +Cl_minus_1_a4*(4*pi/(Cl_minus_1_b4+B))**1.5*safexp(-4*pi**2/(Cl_minus_1_b4+B)*x*x) \
  +Cl_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
Ar_c = 1.444500;
Ar_a1 = 7.484500;Ar_a2 = 6.772300;Ar_a3 = 0.653900;Ar_a4 = 1.644200;
Ar_b1 = 0.907200;Ar_b2 = 14.840700;Ar_b3 = 43.898300;Ar_b4 = 33.392899;
Ar_ff(x,B) = \
  +Ar_a1*(4*pi/(Ar_b1+B))**1.5*safexp(-4*pi**2/(Ar_b1+B)*x*x) \
  +Ar_a2*(4*pi/(Ar_b2+B))**1.5*safexp(-4*pi**2/(Ar_b2+B)*x*x) \
  +Ar_a3*(4*pi/(Ar_b3+B))**1.5*safexp(-4*pi**2/(Ar_b3+B)*x*x) \
  +Ar_a4*(4*pi/(Ar_b4+B))**1.5*safexp(-4*pi**2/(Ar_b4+B)*x*x) \
  +Ar_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=19
K_c = 1.422800;
K_a1 = 8.218599;K_a2 = 7.439800;K_a3 = 1.051900;K_a4 = 0.865900;
K_b1 = 12.794900;K_b2 = 0.774800;K_b3 = 213.186996;K_b4 = 41.684097;
K_ff(x,B) = \
  +K_a1*(4*pi/(K_b1+B))**1.5*safexp(-4*pi**2/(K_b1+B)*x*x) \
  +K_a2*(4*pi/(K_b2+B))**1.5*safexp(-4*pi**2/(K_b2+B)*x*x) \
  +K_a3*(4*pi/(K_b3+B))**1.5*safexp(-4*pi**2/(K_b3+B)*x*x) \
  +K_a4*(4*pi/(K_b4+B))**1.5*safexp(-4*pi**2/(K_b4+B)*x*x) \
  +K_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
K_plus_1_c = -4.997800;
K_plus_1_a1 = 7.957800;K_plus_1_a2 = 7.491700;K_plus_1_a3 = 6.359000;K_plus_1_a4 = 1.191500;
K_plus_1_b1 = 12.633100;K_plus_1_b2 = 0.767400;K_plus_1_b3 = -0.002000;K_plus_1_b4 = 31.912800;
K_plus_1_ff(x,B) = \
  +K_plus_1_a1*(4*pi/(K_plus_1_b1+B))**1.5*safexp(-4*pi**2/(K_plus_1_b1+B)*x*x) \
  +K_plus_1_a2*(4*pi/(K_plus_1_b2+B))**1.5*safexp(-4*pi**2/(K_plus_1_b2+B)*x*x) \
  +K_plus_1_a3*(4*pi/(K_plus_1_b3+B))**1.5*safexp(-4*pi**2/(K_plus_1_b3+B)*x*x) \
  +K_plus_1_a4*(4*pi/(K_plus_1_b4+B))**1.5*safexp(-4*pi**2/(K_plus_1_b4+B)*x*x) \
  +K_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=20
Ca_c = 1.375100;
Ca_a1 = 8.626600;Ca_a2 = 7.387300;Ca_a3 = 1.589900;Ca_a4 = 1.021100;
Ca_b1 = 10.442100;Ca_b2 = 0.659900;Ca_b3 = 85.748398;Ca_b4 = 178.436996;
Ca_ff(x,B) = \
  +Ca_a1*(4*pi/(Ca_b1+B))**1.5*safexp(-4*pi**2/(Ca_b1+B)*x*x) \
  +Ca_a2*(4*pi/(Ca_b2+B))**1.5*safexp(-4*pi**2/(Ca_b2+B)*x*x) \
  +Ca_a3*(4*pi/(Ca_b3+B))**1.5*safexp(-4*pi**2/(Ca_b3+B)*x*x) \
  +Ca_a4*(4*pi/(Ca_b4+B))**1.5*safexp(-4*pi**2/(Ca_b4+B)*x*x) \
  +Ca_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
Ca_plus_2_c = -14.875000;
Ca_plus_2_a1 = 15.634800;Ca_plus_2_a2 = 7.951800;Ca_plus_2_a3 = 8.437200;Ca_plus_2_a4 = 0.853700;
Ca_plus_2_b1 = -0.007400;Ca_plus_2_b2 = 0.608900;Ca_plus_2_b3 = 10.311600;Ca_plus_2_b4 = 25.990499;
Ca_plus_2_ff(x,B) = \
  +Ca_plus_2_a1*(4*pi/(Ca_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ca_plus_2_b1+B)*x*x) \
  +Ca_plus_2_a2*(4*pi/(Ca_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ca_plus_2_b2+B)*x*x) \
  +Ca_plus_2_a3*(4*pi/(Ca_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ca_plus_2_b3+B)*x*x) \
  +Ca_plus_2_a4*(4*pi/(Ca_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ca_plus_2_b4+B)*x*x) \
  +Ca_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=21
Sc_c = 1.332900;
Sc_a1 = 9.189000;Sc_a2 = 7.367900;Sc_a3 = 1.640900;Sc_a4 = 1.468000;
Sc_b1 = 9.021299;Sc_b2 = 0.572900;Sc_b3 = 136.108002;Sc_b4 = 51.353100;
Sc_ff(x,B) = \
  +Sc_a1*(4*pi/(Sc_b1+B))**1.5*safexp(-4*pi**2/(Sc_b1+B)*x*x) \
  +Sc_a2*(4*pi/(Sc_b2+B))**1.5*safexp(-4*pi**2/(Sc_b2+B)*x*x) \
  +Sc_a3*(4*pi/(Sc_b3+B))**1.5*safexp(-4*pi**2/(Sc_b3+B)*x*x) \
  +Sc_a4*(4*pi/(Sc_b4+B))**1.5*safexp(-4*pi**2/(Sc_b4+B)*x*x) \
  +Sc_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
Sc_plus_3_c = -6.666700;
Sc_plus_3_a1 = 13.400800;Sc_plus_3_a2 = 8.027300;Sc_plus_3_a3 = 1.659430;Sc_plus_3_a4 = 1.579360;
Sc_plus_3_b1 = 0.298540;Sc_plus_3_b2 = 7.962900;Sc_plus_3_b3 = -0.286040;Sc_plus_3_b4 = 16.066200;
Sc_plus_3_ff(x,B) = \
  +Sc_plus_3_a1*(4*pi/(Sc_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Sc_plus_3_b1+B)*x*x) \
  +Sc_plus_3_a2*(4*pi/(Sc_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Sc_plus_3_b2+B)*x*x) \
  +Sc_plus_3_a3*(4*pi/(Sc_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Sc_plus_3_b3+B)*x*x) \
  +Sc_plus_3_a4*(4*pi/(Sc_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Sc_plus_3_b4+B)*x*x) \
  +Sc_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=22
Ti_c = 1.280700;
Ti_a1 = 9.759500;Ti_a2 = 7.355800;Ti_a3 = 1.699100;Ti_a4 = 1.902100;
Ti_b1 = 7.850800;Ti_b2 = 0.500000;Ti_b3 = 35.633801;Ti_b4 = 116.104996;
Ti_ff(x,B) = \
  +Ti_a1*(4*pi/(Ti_b1+B))**1.5*safexp(-4*pi**2/(Ti_b1+B)*x*x) \
  +Ti_a2*(4*pi/(Ti_b2+B))**1.5*safexp(-4*pi**2/(Ti_b2+B)*x*x) \
  +Ti_a3*(4*pi/(Ti_b3+B))**1.5*safexp(-4*pi**2/(Ti_b3+B)*x*x) \
  +Ti_a4*(4*pi/(Ti_b4+B))**1.5*safexp(-4*pi**2/(Ti_b4+B)*x*x) \
  +Ti_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=20
Ti_plus_2_c = 0.897155;
Ti_plus_2_a1 = 9.114230;Ti_plus_2_a2 = 7.621740;Ti_plus_2_a3 = 2.279300;Ti_plus_2_a4 = 0.087899;
Ti_plus_2_b1 = 7.524300;Ti_plus_2_b2 = 0.457585;Ti_plus_2_b3 = 19.536100;Ti_plus_2_b4 = 61.655800;
Ti_plus_2_ff(x,B) = \
  +Ti_plus_2_a1*(4*pi/(Ti_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ti_plus_2_b1+B)*x*x) \
  +Ti_plus_2_a2*(4*pi/(Ti_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ti_plus_2_b2+B)*x*x) \
  +Ti_plus_2_a3*(4*pi/(Ti_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ti_plus_2_b3+B)*x*x) \
  +Ti_plus_2_a4*(4*pi/(Ti_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ti_plus_2_b4+B)*x*x) \
  +Ti_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=19
Ti_plus_3_c = -14.652000;
Ti_plus_3_a1 = 17.734400;Ti_plus_3_a2 = 8.738160;Ti_plus_3_a3 = 5.256910;Ti_plus_3_a4 = 1.921340;
Ti_plus_3_b1 = 0.220610;Ti_plus_3_b2 = 7.047160;Ti_plus_3_b3 = -0.157620;Ti_plus_3_b4 = 15.976800;
Ti_plus_3_ff(x,B) = \
  +Ti_plus_3_a1*(4*pi/(Ti_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ti_plus_3_b1+B)*x*x) \
  +Ti_plus_3_a2*(4*pi/(Ti_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ti_plus_3_b2+B)*x*x) \
  +Ti_plus_3_a3*(4*pi/(Ti_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ti_plus_3_b3+B)*x*x) \
  +Ti_plus_3_a4*(4*pi/(Ti_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ti_plus_3_b4+B)*x*x) \
  +Ti_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
Ti_plus_4_c = -13.280000;
Ti_plus_4_a1 = 19.511400;Ti_plus_4_a2 = 8.234730;Ti_plus_4_a3 = 2.013410;Ti_plus_4_a4 = 1.520800;
Ti_plus_4_b1 = 0.178847;Ti_plus_4_b2 = 6.670180;Ti_plus_4_b3 = -0.292630;Ti_plus_4_b4 = 12.946400;
Ti_plus_4_ff(x,B) = \
  +Ti_plus_4_a1*(4*pi/(Ti_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Ti_plus_4_b1+B)*x*x) \
  +Ti_plus_4_a2*(4*pi/(Ti_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Ti_plus_4_b2+B)*x*x) \
  +Ti_plus_4_a3*(4*pi/(Ti_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Ti_plus_4_b3+B)*x*x) \
  +Ti_plus_4_a4*(4*pi/(Ti_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Ti_plus_4_b4+B)*x*x) \
  +Ti_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=23
V_c = 1.219900;
V_a1 = 10.297100;V_a2 = 7.351100;V_a3 = 2.070300;V_a4 = 2.057100;
V_b1 = 6.865700;V_b2 = 0.438500;V_b3 = 26.893799;V_b4 = 102.477997;
V_ff(x,B) = \
  +V_a1*(4*pi/(V_b1+B))**1.5*safexp(-4*pi**2/(V_b1+B)*x*x) \
  +V_a2*(4*pi/(V_b2+B))**1.5*safexp(-4*pi**2/(V_b2+B)*x*x) \
  +V_a3*(4*pi/(V_b3+B))**1.5*safexp(-4*pi**2/(V_b3+B)*x*x) \
  +V_a4*(4*pi/(V_b4+B))**1.5*safexp(-4*pi**2/(V_b4+B)*x*x) \
  +V_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=21
V_plus_2_c = 1.229800;
V_plus_2_a1 = 10.106000;V_plus_2_a2 = 7.354100;V_plus_2_a3 = 2.288400;V_plus_2_a4 = 0.022300;
V_plus_2_b1 = 6.881800;V_plus_2_b2 = 0.440900;V_plus_2_b3 = 20.300400;V_plus_2_b4 = 115.122002;
V_plus_2_ff(x,B) = \
  +V_plus_2_a1*(4*pi/(V_plus_2_b1+B))**1.5*safexp(-4*pi**2/(V_plus_2_b1+B)*x*x) \
  +V_plus_2_a2*(4*pi/(V_plus_2_b2+B))**1.5*safexp(-4*pi**2/(V_plus_2_b2+B)*x*x) \
  +V_plus_2_a3*(4*pi/(V_plus_2_b3+B))**1.5*safexp(-4*pi**2/(V_plus_2_b3+B)*x*x) \
  +V_plus_2_a4*(4*pi/(V_plus_2_b4+B))**1.5*safexp(-4*pi**2/(V_plus_2_b4+B)*x*x) \
  +V_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=20
V_plus_3_c = 0.656565;
V_plus_3_a1 = 9.431410;V_plus_3_a2 = 7.741900;V_plus_3_a3 = 2.153430;V_plus_3_a4 = 0.016865;
V_plus_3_b1 = 6.395350;V_plus_3_b2 = 0.383349;V_plus_3_b3 = 15.190800;V_plus_3_b4 = 63.969002;
V_plus_3_ff(x,B) = \
  +V_plus_3_a1*(4*pi/(V_plus_3_b1+B))**1.5*safexp(-4*pi**2/(V_plus_3_b1+B)*x*x) \
  +V_plus_3_a2*(4*pi/(V_plus_3_b2+B))**1.5*safexp(-4*pi**2/(V_plus_3_b2+B)*x*x) \
  +V_plus_3_a3*(4*pi/(V_plus_3_b3+B))**1.5*safexp(-4*pi**2/(V_plus_3_b3+B)*x*x) \
  +V_plus_3_a4*(4*pi/(V_plus_3_b4+B))**1.5*safexp(-4*pi**2/(V_plus_3_b4+B)*x*x) \
  +V_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=18
V_plus_5_c = 1.714300;
V_plus_5_a1 = 15.688700;V_plus_5_a2 = 8.142080;V_plus_5_a3 = 2.030810;V_plus_5_a4 = -9.576000;
V_plus_5_b1 = 0.679003;V_plus_5_b2 = 5.401350;V_plus_5_b3 = 9.972780;V_plus_5_b4 = 0.940464;
V_plus_5_ff(x,B) = \
  +V_plus_5_a1*(4*pi/(V_plus_5_b1+B))**1.5*safexp(-4*pi**2/(V_plus_5_b1+B)*x*x) \
  +V_plus_5_a2*(4*pi/(V_plus_5_b2+B))**1.5*safexp(-4*pi**2/(V_plus_5_b2+B)*x*x) \
  +V_plus_5_a3*(4*pi/(V_plus_5_b3+B))**1.5*safexp(-4*pi**2/(V_plus_5_b3+B)*x*x) \
  +V_plus_5_a4*(4*pi/(V_plus_5_b4+B))**1.5*safexp(-4*pi**2/(V_plus_5_b4+B)*x*x) \
  +V_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=24
Cr_c = 1.183200;
Cr_a1 = 10.640600;Cr_a2 = 7.353700;Cr_a3 = 3.324000;Cr_a4 = 1.492200;
Cr_b1 = 6.103800;Cr_b2 = 0.392000;Cr_b3 = 20.262600;Cr_b4 = 98.739899;
Cr_ff(x,B) = \
  +Cr_a1*(4*pi/(Cr_b1+B))**1.5*safexp(-4*pi**2/(Cr_b1+B)*x*x) \
  +Cr_a2*(4*pi/(Cr_b2+B))**1.5*safexp(-4*pi**2/(Cr_b2+B)*x*x) \
  +Cr_a3*(4*pi/(Cr_b3+B))**1.5*safexp(-4*pi**2/(Cr_b3+B)*x*x) \
  +Cr_a4*(4*pi/(Cr_b4+B))**1.5*safexp(-4*pi**2/(Cr_b4+B)*x*x) \
  +Cr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=22
Cr_plus_2_c = 0.616898;
Cr_plus_2_a1 = 9.540340;Cr_plus_2_a2 = 7.750900;Cr_plus_2_a3 = 3.582740;Cr_plus_2_a4 = 0.509107;
Cr_plus_2_b1 = 5.660780;Cr_plus_2_b2 = 0.344261;Cr_plus_2_b3 = 13.307500;Cr_plus_2_b4 = 32.422401;
Cr_plus_2_ff(x,B) = \
  +Cr_plus_2_a1*(4*pi/(Cr_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Cr_plus_2_b1+B)*x*x) \
  +Cr_plus_2_a2*(4*pi/(Cr_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Cr_plus_2_b2+B)*x*x) \
  +Cr_plus_2_a3*(4*pi/(Cr_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Cr_plus_2_b3+B)*x*x) \
  +Cr_plus_2_a4*(4*pi/(Cr_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Cr_plus_2_b4+B)*x*x) \
  +Cr_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=21
Cr_plus_3_c = 0.518275;
Cr_plus_3_a1 = 9.680900;Cr_plus_3_a2 = 7.811360;Cr_plus_3_a3 = 2.876030;Cr_plus_3_a4 = 0.113575;
Cr_plus_3_b1 = 5.594630;Cr_plus_3_b2 = 0.334393;Cr_plus_3_b3 = 12.828800;Cr_plus_3_b4 = 32.876099;
Cr_plus_3_ff(x,B) = \
  +Cr_plus_3_a1*(4*pi/(Cr_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Cr_plus_3_b1+B)*x*x) \
  +Cr_plus_3_a2*(4*pi/(Cr_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Cr_plus_3_b2+B)*x*x) \
  +Cr_plus_3_a3*(4*pi/(Cr_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Cr_plus_3_b3+B)*x*x) \
  +Cr_plus_3_a4*(4*pi/(Cr_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Cr_plus_3_b4+B)*x*x) \
  +Cr_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=25
Mn_c = 1.089600;
Mn_a1 = 11.281900;Mn_a2 = 7.357300;Mn_a3 = 3.019300;Mn_a4 = 2.244100;
Mn_b1 = 5.340900;Mn_b2 = 0.343200;Mn_b3 = 17.867399;Mn_b4 = 83.754303;
Mn_ff(x,B) = \
  +Mn_a1*(4*pi/(Mn_b1+B))**1.5*safexp(-4*pi**2/(Mn_b1+B)*x*x) \
  +Mn_a2*(4*pi/(Mn_b2+B))**1.5*safexp(-4*pi**2/(Mn_b2+B)*x*x) \
  +Mn_a3*(4*pi/(Mn_b3+B))**1.5*safexp(-4*pi**2/(Mn_b3+B)*x*x) \
  +Mn_a4*(4*pi/(Mn_b4+B))**1.5*safexp(-4*pi**2/(Mn_b4+B)*x*x) \
  +Mn_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=23
Mn_plus_2_c = 1.087400;
Mn_plus_2_a1 = 10.806100;Mn_plus_2_a2 = 7.362000;Mn_plus_2_a3 = 3.526800;Mn_plus_2_a4 = 0.218400;
Mn_plus_2_b1 = 5.279600;Mn_plus_2_b2 = 0.343500;Mn_plus_2_b3 = 14.343000;Mn_plus_2_b4 = 41.323502;
Mn_plus_2_ff(x,B) = \
  +Mn_plus_2_a1*(4*pi/(Mn_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Mn_plus_2_b1+B)*x*x) \
  +Mn_plus_2_a2*(4*pi/(Mn_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Mn_plus_2_b2+B)*x*x) \
  +Mn_plus_2_a3*(4*pi/(Mn_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Mn_plus_2_b3+B)*x*x) \
  +Mn_plus_2_a4*(4*pi/(Mn_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Mn_plus_2_b4+B)*x*x) \
  +Mn_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=22
Mn_plus_3_c = 0.393974;
Mn_plus_3_a1 = 9.845210;Mn_plus_3_a2 = 7.871940;Mn_plus_3_a3 = 3.565310;Mn_plus_3_a4 = 0.323613;
Mn_plus_3_b1 = 4.917970;Mn_plus_3_b2 = 0.294393;Mn_plus_3_b3 = 10.817100;Mn_plus_3_b4 = 24.128099;
Mn_plus_3_ff(x,B) = \
  +Mn_plus_3_a1*(4*pi/(Mn_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Mn_plus_3_b1+B)*x*x) \
  +Mn_plus_3_a2*(4*pi/(Mn_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Mn_plus_3_b2+B)*x*x) \
  +Mn_plus_3_a3*(4*pi/(Mn_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Mn_plus_3_b3+B)*x*x) \
  +Mn_plus_3_a4*(4*pi/(Mn_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Mn_plus_3_b4+B)*x*x) \
  +Mn_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=21
Mn_plus_4_c = 0.251877;
Mn_plus_4_a1 = 9.962530;Mn_plus_4_a2 = 7.970570;Mn_plus_4_a3 = 2.760670;Mn_plus_4_a4 = 0.054447;
Mn_plus_4_b1 = 4.848500;Mn_plus_4_b2 = 0.283303;Mn_plus_4_b3 = 10.485200;Mn_plus_4_b4 = 27.573000;
Mn_plus_4_ff(x,B) = \
  +Mn_plus_4_a1*(4*pi/(Mn_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Mn_plus_4_b1+B)*x*x) \
  +Mn_plus_4_a2*(4*pi/(Mn_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Mn_plus_4_b2+B)*x*x) \
  +Mn_plus_4_a3*(4*pi/(Mn_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Mn_plus_4_b3+B)*x*x) \
  +Mn_plus_4_a4*(4*pi/(Mn_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Mn_plus_4_b4+B)*x*x) \
  +Mn_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=26
Fe_c = 1.036900;
Fe_a1 = 11.769500;Fe_a2 = 7.357300;Fe_a3 = 3.522200;Fe_a4 = 2.304500;
Fe_b1 = 4.761100;Fe_b2 = 0.307200;Fe_b3 = 15.353500;Fe_b4 = 76.880501;
Fe_ff(x,B) = \
  +Fe_a1*(4*pi/(Fe_b1+B))**1.5*safexp(-4*pi**2/(Fe_b1+B)*x*x) \
  +Fe_a2*(4*pi/(Fe_b2+B))**1.5*safexp(-4*pi**2/(Fe_b2+B)*x*x) \
  +Fe_a3*(4*pi/(Fe_b3+B))**1.5*safexp(-4*pi**2/(Fe_b3+B)*x*x) \
  +Fe_a4*(4*pi/(Fe_b4+B))**1.5*safexp(-4*pi**2/(Fe_b4+B)*x*x) \
  +Fe_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=24
Fe_plus_2_c = 1.009700;
Fe_plus_2_a1 = 11.042400;Fe_plus_2_a2 = 7.374000;Fe_plus_2_a3 = 4.134600;Fe_plus_2_a4 = 0.439900;
Fe_plus_2_b1 = 4.653800;Fe_plus_2_b2 = 0.305300;Fe_plus_2_b3 = 12.054600;Fe_plus_2_b4 = 31.280899;
Fe_plus_2_ff(x,B) = \
  +Fe_plus_2_a1*(4*pi/(Fe_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Fe_plus_2_b1+B)*x*x) \
  +Fe_plus_2_a2*(4*pi/(Fe_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Fe_plus_2_b2+B)*x*x) \
  +Fe_plus_2_a3*(4*pi/(Fe_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Fe_plus_2_b3+B)*x*x) \
  +Fe_plus_2_a4*(4*pi/(Fe_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Fe_plus_2_b4+B)*x*x) \
  +Fe_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=23
Fe_plus_3_c = 0.970700;
Fe_plus_3_a1 = 11.176400;Fe_plus_3_a2 = 7.386300;Fe_plus_3_a3 = 3.394800;Fe_plus_3_a4 = 0.072400;
Fe_plus_3_b1 = 4.614700;Fe_plus_3_b2 = 0.300500;Fe_plus_3_b3 = 11.672900;Fe_plus_3_b4 = 38.556599;
Fe_plus_3_ff(x,B) = \
  +Fe_plus_3_a1*(4*pi/(Fe_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Fe_plus_3_b1+B)*x*x) \
  +Fe_plus_3_a2*(4*pi/(Fe_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Fe_plus_3_b2+B)*x*x) \
  +Fe_plus_3_a3*(4*pi/(Fe_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Fe_plus_3_b3+B)*x*x) \
  +Fe_plus_3_a4*(4*pi/(Fe_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Fe_plus_3_b4+B)*x*x) \
  +Fe_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=27
Co_c = 1.011800;
Co_a1 = 12.284100;Co_a2 = 7.340900;Co_a3 = 4.003400;Co_a4 = 2.348800;
Co_b1 = 4.279100;Co_b2 = 0.278400;Co_b3 = 13.535900;Co_b4 = 71.169197;
Co_ff(x,B) = \
  +Co_a1*(4*pi/(Co_b1+B))**1.5*safexp(-4*pi**2/(Co_b1+B)*x*x) \
  +Co_a2*(4*pi/(Co_b2+B))**1.5*safexp(-4*pi**2/(Co_b2+B)*x*x) \
  +Co_a3*(4*pi/(Co_b3+B))**1.5*safexp(-4*pi**2/(Co_b3+B)*x*x) \
  +Co_a4*(4*pi/(Co_b4+B))**1.5*safexp(-4*pi**2/(Co_b4+B)*x*x) \
  +Co_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=25
Co_plus_2_c = 0.932400;
Co_plus_2_a1 = 11.229600;Co_plus_2_a2 = 7.388300;Co_plus_2_a3 = 4.739300;Co_plus_2_a4 = 0.710800;
Co_plus_2_b1 = 4.123100;Co_plus_2_b2 = 0.272600;Co_plus_2_b3 = 10.244300;Co_plus_2_b4 = 25.646599;
Co_plus_2_ff(x,B) = \
  +Co_plus_2_a1*(4*pi/(Co_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Co_plus_2_b1+B)*x*x) \
  +Co_plus_2_a2*(4*pi/(Co_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Co_plus_2_b2+B)*x*x) \
  +Co_plus_2_a3*(4*pi/(Co_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Co_plus_2_b3+B)*x*x) \
  +Co_plus_2_a4*(4*pi/(Co_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Co_plus_2_b4+B)*x*x) \
  +Co_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=24
Co_plus_3_c = 0.286667;
Co_plus_3_a1 = 10.337999;Co_plus_3_a2 = 7.881730;Co_plus_3_a3 = 4.767950;Co_plus_3_a4 = 0.725591;
Co_plus_3_b1 = 3.909690;Co_plus_3_b2 = 0.238668;Co_plus_3_b3 = 8.355830;Co_plus_3_b4 = 18.349100;
Co_plus_3_ff(x,B) = \
  +Co_plus_3_a1*(4*pi/(Co_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Co_plus_3_b1+B)*x*x) \
  +Co_plus_3_a2*(4*pi/(Co_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Co_plus_3_b2+B)*x*x) \
  +Co_plus_3_a3*(4*pi/(Co_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Co_plus_3_b3+B)*x*x) \
  +Co_plus_3_a4*(4*pi/(Co_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Co_plus_3_b4+B)*x*x) \
  +Co_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=28
Ni_c = 1.034100;
Ni_a1 = 12.837600;Ni_a2 = 7.292000;Ni_a3 = 4.443800;Ni_a4 = 2.380000;
Ni_b1 = 3.878500;Ni_b2 = 0.256500;Ni_b3 = 12.176300;Ni_b4 = 66.342102;
Ni_ff(x,B) = \
  +Ni_a1*(4*pi/(Ni_b1+B))**1.5*safexp(-4*pi**2/(Ni_b1+B)*x*x) \
  +Ni_a2*(4*pi/(Ni_b2+B))**1.5*safexp(-4*pi**2/(Ni_b2+B)*x*x) \
  +Ni_a3*(4*pi/(Ni_b3+B))**1.5*safexp(-4*pi**2/(Ni_b3+B)*x*x) \
  +Ni_a4*(4*pi/(Ni_b4+B))**1.5*safexp(-4*pi**2/(Ni_b4+B)*x*x) \
  +Ni_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=26
Ni_plus_2_c = 0.861400;
Ni_plus_2_a1 = 11.416600;Ni_plus_2_a2 = 7.400500;Ni_plus_2_a3 = 5.344200;Ni_plus_2_a4 = 0.977300;
Ni_plus_2_b1 = 3.676600;Ni_plus_2_b2 = 0.244900;Ni_plus_2_b3 = 8.873000;Ni_plus_2_b4 = 22.162600;
Ni_plus_2_ff(x,B) = \
  +Ni_plus_2_a1*(4*pi/(Ni_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ni_plus_2_b1+B)*x*x) \
  +Ni_plus_2_a2*(4*pi/(Ni_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ni_plus_2_b2+B)*x*x) \
  +Ni_plus_2_a3*(4*pi/(Ni_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ni_plus_2_b3+B)*x*x) \
  +Ni_plus_2_a4*(4*pi/(Ni_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ni_plus_2_b4+B)*x*x) \
  +Ni_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=25
Ni_plus_3_c = 0.386044;
Ni_plus_3_a1 = 10.780600;Ni_plus_3_a2 = 7.758680;Ni_plus_3_a3 = 5.227460;Ni_plus_3_a4 = 0.847114;
Ni_plus_3_b1 = 3.547700;Ni_plus_3_b2 = 0.223140;Ni_plus_3_b3 = 7.644680;Ni_plus_3_b4 = 16.967300;
Ni_plus_3_ff(x,B) = \
  +Ni_plus_3_a1*(4*pi/(Ni_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ni_plus_3_b1+B)*x*x) \
  +Ni_plus_3_a2*(4*pi/(Ni_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ni_plus_3_b2+B)*x*x) \
  +Ni_plus_3_a3*(4*pi/(Ni_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ni_plus_3_b3+B)*x*x) \
  +Ni_plus_3_a4*(4*pi/(Ni_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ni_plus_3_b4+B)*x*x) \
  +Ni_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=29
Cu_c = 1.191000;
Cu_a1 = 13.337999;Cu_a2 = 7.167600;Cu_a3 = 5.615800;Cu_a4 = 1.673500;
Cu_b1 = 3.582800;Cu_b2 = 0.247000;Cu_b3 = 11.396600;Cu_b4 = 64.812599;
Cu_ff(x,B) = \
  +Cu_a1*(4*pi/(Cu_b1+B))**1.5*safexp(-4*pi**2/(Cu_b1+B)*x*x) \
  +Cu_a2*(4*pi/(Cu_b2+B))**1.5*safexp(-4*pi**2/(Cu_b2+B)*x*x) \
  +Cu_a3*(4*pi/(Cu_b3+B))**1.5*safexp(-4*pi**2/(Cu_b3+B)*x*x) \
  +Cu_a4*(4*pi/(Cu_b4+B))**1.5*safexp(-4*pi**2/(Cu_b4+B)*x*x) \
  +Cu_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=28
Cu_plus_1_c = 0.890000;
Cu_plus_1_a1 = 11.947500;Cu_plus_1_a2 = 7.357300;Cu_plus_1_a3 = 6.245500;Cu_plus_1_a4 = 1.557800;
Cu_plus_1_b1 = 3.366900;Cu_plus_1_b2 = 0.227400;Cu_plus_1_b3 = 8.662500;Cu_plus_1_b4 = 25.848700;
Cu_plus_1_ff(x,B) = \
  +Cu_plus_1_a1*(4*pi/(Cu_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Cu_plus_1_b1+B)*x*x) \
  +Cu_plus_1_a2*(4*pi/(Cu_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Cu_plus_1_b2+B)*x*x) \
  +Cu_plus_1_a3*(4*pi/(Cu_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Cu_plus_1_b3+B)*x*x) \
  +Cu_plus_1_a4*(4*pi/(Cu_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Cu_plus_1_b4+B)*x*x) \
  +Cu_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=27
Cu_plus_2_c = 1.144310;
Cu_plus_2_a1 = 11.816800;Cu_plus_2_a2 = 7.111810;Cu_plus_2_a3 = 5.781350;Cu_plus_2_a4 = 1.145230;
Cu_plus_2_b1 = 3.374840;Cu_plus_2_b2 = 0.244078;Cu_plus_2_b3 = 7.987600;Cu_plus_2_b4 = 19.896999;
Cu_plus_2_ff(x,B) = \
  +Cu_plus_2_a1*(4*pi/(Cu_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Cu_plus_2_b1+B)*x*x) \
  +Cu_plus_2_a2*(4*pi/(Cu_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Cu_plus_2_b2+B)*x*x) \
  +Cu_plus_2_a3*(4*pi/(Cu_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Cu_plus_2_b3+B)*x*x) \
  +Cu_plus_2_a4*(4*pi/(Cu_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Cu_plus_2_b4+B)*x*x) \
  +Cu_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=30
Zn_c = 1.304100;
Zn_a1 = 14.074300;Zn_a2 = 7.031800;Zn_a3 = 5.162500;Zn_a4 = 2.410000;
Zn_b1 = 3.265500;Zn_b2 = 0.233300;Zn_b3 = 10.316299;Zn_b4 = 58.709702;
Zn_ff(x,B) = \
  +Zn_a1*(4*pi/(Zn_b1+B))**1.5*safexp(-4*pi**2/(Zn_b1+B)*x*x) \
  +Zn_a2*(4*pi/(Zn_b2+B))**1.5*safexp(-4*pi**2/(Zn_b2+B)*x*x) \
  +Zn_a3*(4*pi/(Zn_b3+B))**1.5*safexp(-4*pi**2/(Zn_b3+B)*x*x) \
  +Zn_a4*(4*pi/(Zn_b4+B))**1.5*safexp(-4*pi**2/(Zn_b4+B)*x*x) \
  +Zn_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=28
Zn_plus_2_c = 0.780700;
Zn_plus_2_a1 = 11.971900;Zn_plus_2_a2 = 7.386200;Zn_plus_2_a3 = 6.466800;Zn_plus_2_a4 = 1.394000;
Zn_plus_2_b1 = 2.994600;Zn_plus_2_b2 = 0.203100;Zn_plus_2_b3 = 7.082600;Zn_plus_2_b4 = 18.099499;
Zn_plus_2_ff(x,B) = \
  +Zn_plus_2_a1*(4*pi/(Zn_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Zn_plus_2_b1+B)*x*x) \
  +Zn_plus_2_a2*(4*pi/(Zn_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Zn_plus_2_b2+B)*x*x) \
  +Zn_plus_2_a3*(4*pi/(Zn_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Zn_plus_2_b3+B)*x*x) \
  +Zn_plus_2_a4*(4*pi/(Zn_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Zn_plus_2_b4+B)*x*x) \
  +Zn_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=31
Ga_c = 1.718900;
Ga_a1 = 15.235400;Ga_a2 = 6.700600;Ga_a3 = 4.359100;Ga_a4 = 2.962300;
Ga_b1 = 3.066900;Ga_b2 = 0.241200;Ga_b3 = 10.780500;Ga_b4 = 61.413498;
Ga_ff(x,B) = \
  +Ga_a1*(4*pi/(Ga_b1+B))**1.5*safexp(-4*pi**2/(Ga_b1+B)*x*x) \
  +Ga_a2*(4*pi/(Ga_b2+B))**1.5*safexp(-4*pi**2/(Ga_b2+B)*x*x) \
  +Ga_a3*(4*pi/(Ga_b3+B))**1.5*safexp(-4*pi**2/(Ga_b3+B)*x*x) \
  +Ga_a4*(4*pi/(Ga_b4+B))**1.5*safexp(-4*pi**2/(Ga_b4+B)*x*x) \
  +Ga_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=28
Ga_plus_3_c = 1.535450;
Ga_plus_3_a1 = 12.691999;Ga_plus_3_a2 = 6.698830;Ga_plus_3_a3 = 6.066920;Ga_plus_3_a4 = 1.006600;
Ga_plus_3_b1 = 2.812620;Ga_plus_3_b2 = 0.227890;Ga_plus_3_b3 = 6.364410;Ga_plus_3_b4 = 14.412200;
Ga_plus_3_ff(x,B) = \
  +Ga_plus_3_a1*(4*pi/(Ga_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ga_plus_3_b1+B)*x*x) \
  +Ga_plus_3_a2*(4*pi/(Ga_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ga_plus_3_b2+B)*x*x) \
  +Ga_plus_3_a3*(4*pi/(Ga_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ga_plus_3_b3+B)*x*x) \
  +Ga_plus_3_a4*(4*pi/(Ga_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ga_plus_3_b4+B)*x*x) \
  +Ga_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=32
Ge_c = 2.131300;
Ge_a1 = 16.081600;Ge_a2 = 6.374700;Ge_a3 = 3.706800;Ge_a4 = 3.683000;
Ge_b1 = 2.850900;Ge_b2 = 0.251600;Ge_b3 = 11.446800;Ge_b4 = 54.762501;
Ge_ff(x,B) = \
  +Ge_a1*(4*pi/(Ge_b1+B))**1.5*safexp(-4*pi**2/(Ge_b1+B)*x*x) \
  +Ge_a2*(4*pi/(Ge_b2+B))**1.5*safexp(-4*pi**2/(Ge_b2+B)*x*x) \
  +Ge_a3*(4*pi/(Ge_b3+B))**1.5*safexp(-4*pi**2/(Ge_b3+B)*x*x) \
  +Ge_a4*(4*pi/(Ge_b4+B))**1.5*safexp(-4*pi**2/(Ge_b4+B)*x*x) \
  +Ge_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=28
Ge_plus_4_c = 1.455720;
Ge_plus_4_a1 = 12.917200;Ge_plus_4_a2 = 6.700030;Ge_plus_4_a3 = 6.067910;Ge_plus_4_a4 = 0.859041;
Ge_plus_4_b1 = 2.537180;Ge_plus_4_b2 = 0.205855;Ge_plus_4_b3 = 5.479130;Ge_plus_4_b4 = 11.603000;
Ge_plus_4_ff(x,B) = \
  +Ge_plus_4_a1*(4*pi/(Ge_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Ge_plus_4_b1+B)*x*x) \
  +Ge_plus_4_a2*(4*pi/(Ge_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Ge_plus_4_b2+B)*x*x) \
  +Ge_plus_4_a3*(4*pi/(Ge_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Ge_plus_4_b3+B)*x*x) \
  +Ge_plus_4_a4*(4*pi/(Ge_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Ge_plus_4_b4+B)*x*x) \
  +Ge_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=33
As_c = 2.531000;
As_a1 = 16.672300;As_a2 = 6.070100;As_a3 = 3.431300;As_a4 = 4.277900;
As_b1 = 2.634500;As_b2 = 0.264700;As_b3 = 12.947900;As_b4 = 47.797199;
As_ff(x,B) = \
  +As_a1*(4*pi/(As_b1+B))**1.5*safexp(-4*pi**2/(As_b1+B)*x*x) \
  +As_a2*(4*pi/(As_b2+B))**1.5*safexp(-4*pi**2/(As_b2+B)*x*x) \
  +As_a3*(4*pi/(As_b3+B))**1.5*safexp(-4*pi**2/(As_b3+B)*x*x) \
  +As_a4*(4*pi/(As_b4+B))**1.5*safexp(-4*pi**2/(As_b4+B)*x*x) \
  +As_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=34
Se_c = 2.840900;
Se_a1 = 17.000599;Se_a2 = 5.819600;Se_a3 = 3.973100;Se_a4 = 4.354300;
Se_b1 = 2.409800;Se_b2 = 0.272600;Se_b3 = 15.237200;Se_b4 = 43.816299;
Se_ff(x,B) = \
  +Se_a1*(4*pi/(Se_b1+B))**1.5*safexp(-4*pi**2/(Se_b1+B)*x*x) \
  +Se_a2*(4*pi/(Se_b2+B))**1.5*safexp(-4*pi**2/(Se_b2+B)*x*x) \
  +Se_a3*(4*pi/(Se_b3+B))**1.5*safexp(-4*pi**2/(Se_b3+B)*x*x) \
  +Se_a4*(4*pi/(Se_b4+B))**1.5*safexp(-4*pi**2/(Se_b4+B)*x*x) \
  +Se_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=35
Br_c = 2.955700;
Br_a1 = 17.178900;Br_a2 = 5.235800;Br_a3 = 5.637700;Br_a4 = 3.985100;
Br_b1 = 2.172300;Br_b2 = 16.579599;Br_b3 = 0.260900;Br_b4 = 41.432800;
Br_ff(x,B) = \
  +Br_a1*(4*pi/(Br_b1+B))**1.5*safexp(-4*pi**2/(Br_b1+B)*x*x) \
  +Br_a2*(4*pi/(Br_b2+B))**1.5*safexp(-4*pi**2/(Br_b2+B)*x*x) \
  +Br_a3*(4*pi/(Br_b3+B))**1.5*safexp(-4*pi**2/(Br_b3+B)*x*x) \
  +Br_a4*(4*pi/(Br_b4+B))**1.5*safexp(-4*pi**2/(Br_b4+B)*x*x) \
  +Br_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Br_minus_1_c = 3.177600;
Br_minus_1_a1 = 17.171799;Br_minus_1_a2 = 6.333800;Br_minus_1_a3 = 5.575400;Br_minus_1_a4 = 3.727200;
Br_minus_1_b1 = 2.205900;Br_minus_1_b2 = 19.334499;Br_minus_1_b3 = 0.287100;Br_minus_1_b4 = 58.153500;
Br_minus_1_ff(x,B) = \
  +Br_minus_1_a1*(4*pi/(Br_minus_1_b1+B))**1.5*safexp(-4*pi**2/(Br_minus_1_b1+B)*x*x) \
  +Br_minus_1_a2*(4*pi/(Br_minus_1_b2+B))**1.5*safexp(-4*pi**2/(Br_minus_1_b2+B)*x*x) \
  +Br_minus_1_a3*(4*pi/(Br_minus_1_b3+B))**1.5*safexp(-4*pi**2/(Br_minus_1_b3+B)*x*x) \
  +Br_minus_1_a4*(4*pi/(Br_minus_1_b4+B))**1.5*safexp(-4*pi**2/(Br_minus_1_b4+B)*x*x) \
  +Br_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Kr_c = 2.825000;
Kr_a1 = 17.355499;Kr_a2 = 6.728600;Kr_a3 = 5.549300;Kr_a4 = 3.537500;
Kr_b1 = 1.938400;Kr_b2 = 16.562300;Kr_b3 = 0.226100;Kr_b4 = 39.397202;
Kr_ff(x,B) = \
  +Kr_a1*(4*pi/(Kr_b1+B))**1.5*safexp(-4*pi**2/(Kr_b1+B)*x*x) \
  +Kr_a2*(4*pi/(Kr_b2+B))**1.5*safexp(-4*pi**2/(Kr_b2+B)*x*x) \
  +Kr_a3*(4*pi/(Kr_b3+B))**1.5*safexp(-4*pi**2/(Kr_b3+B)*x*x) \
  +Kr_a4*(4*pi/(Kr_b4+B))**1.5*safexp(-4*pi**2/(Kr_b4+B)*x*x) \
  +Kr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=37
Rb_c = 3.487300;
Rb_a1 = 17.178400;Rb_a2 = 9.643499;Rb_a3 = 5.139900;Rb_a4 = 1.529200;
Rb_b1 = 1.788800;Rb_b2 = 17.315100;Rb_b3 = 0.274800;Rb_b4 = 164.933990;
Rb_ff(x,B) = \
  +Rb_a1*(4*pi/(Rb_b1+B))**1.5*safexp(-4*pi**2/(Rb_b1+B)*x*x) \
  +Rb_a2*(4*pi/(Rb_b2+B))**1.5*safexp(-4*pi**2/(Rb_b2+B)*x*x) \
  +Rb_a3*(4*pi/(Rb_b3+B))**1.5*safexp(-4*pi**2/(Rb_b3+B)*x*x) \
  +Rb_a4*(4*pi/(Rb_b4+B))**1.5*safexp(-4*pi**2/(Rb_b4+B)*x*x) \
  +Rb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Rb_plus_1_c = 2.078200;
Rb_plus_1_a1 = 17.581600;Rb_plus_1_a2 = 7.659800;Rb_plus_1_a3 = 5.898100;Rb_plus_1_a4 = 2.781700;
Rb_plus_1_b1 = 1.713900;Rb_plus_1_b2 = 14.795700;Rb_plus_1_b3 = 0.160300;Rb_plus_1_b4 = 31.208700;
Rb_plus_1_ff(x,B) = \
  +Rb_plus_1_a1*(4*pi/(Rb_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Rb_plus_1_b1+B)*x*x) \
  +Rb_plus_1_a2*(4*pi/(Rb_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Rb_plus_1_b2+B)*x*x) \
  +Rb_plus_1_a3*(4*pi/(Rb_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Rb_plus_1_b3+B)*x*x) \
  +Rb_plus_1_a4*(4*pi/(Rb_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Rb_plus_1_b4+B)*x*x) \
  +Rb_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=38
Sr_c = 2.506400;
Sr_a1 = 17.566299;Sr_a2 = 9.818399;Sr_a3 = 5.422000;Sr_a4 = 2.669400;
Sr_b1 = 1.556400;Sr_b2 = 14.098800;Sr_b3 = 0.166400;Sr_b4 = 132.376007;
Sr_ff(x,B) = \
  +Sr_a1*(4*pi/(Sr_b1+B))**1.5*safexp(-4*pi**2/(Sr_b1+B)*x*x) \
  +Sr_a2*(4*pi/(Sr_b2+B))**1.5*safexp(-4*pi**2/(Sr_b2+B)*x*x) \
  +Sr_a3*(4*pi/(Sr_b3+B))**1.5*safexp(-4*pi**2/(Sr_b3+B)*x*x) \
  +Sr_a4*(4*pi/(Sr_b4+B))**1.5*safexp(-4*pi**2/(Sr_b4+B)*x*x) \
  +Sr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Sr_plus_2_c = 41.402500;
Sr_plus_2_a1 = 18.087400;Sr_plus_2_a2 = 8.137300;Sr_plus_2_a3 = 2.565400;Sr_plus_2_a4 = -34.193001;
Sr_plus_2_b1 = 1.490700;Sr_plus_2_b2 = 12.696300;Sr_plus_2_b3 = 24.565100;Sr_plus_2_b4 = -0.013800;
Sr_plus_2_ff(x,B) = \
  +Sr_plus_2_a1*(4*pi/(Sr_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Sr_plus_2_b1+B)*x*x) \
  +Sr_plus_2_a2*(4*pi/(Sr_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Sr_plus_2_b2+B)*x*x) \
  +Sr_plus_2_a3*(4*pi/(Sr_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Sr_plus_2_b3+B)*x*x) \
  +Sr_plus_2_a4*(4*pi/(Sr_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Sr_plus_2_b4+B)*x*x) \
  +Sr_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=39
Y_c = 1.912130;
Y_a1 = 17.775999;Y_a2 = 10.294600;Y_a3 = 5.726290;Y_a4 = 3.265880;
Y_b1 = 1.402900;Y_b2 = 12.800600;Y_b3 = 0.125599;Y_b4 = 104.353996;
Y_ff(x,B) = \
  +Y_a1*(4*pi/(Y_b1+B))**1.5*safexp(-4*pi**2/(Y_b1+B)*x*x) \
  +Y_a2*(4*pi/(Y_b2+B))**1.5*safexp(-4*pi**2/(Y_b2+B)*x*x) \
  +Y_a3*(4*pi/(Y_b3+B))**1.5*safexp(-4*pi**2/(Y_b3+B)*x*x) \
  +Y_a4*(4*pi/(Y_b4+B))**1.5*safexp(-4*pi**2/(Y_b4+B)*x*x) \
  +Y_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Y_plus_3_c = 40.260201;
Y_plus_3_a1 = 17.926800;Y_plus_3_a2 = 9.153100;Y_plus_3_a3 = 1.767950;Y_plus_3_a4 = -33.108002;
Y_plus_3_b1 = 1.354170;Y_plus_3_b2 = 11.214500;Y_plus_3_b3 = 22.659901;Y_plus_3_b4 = -0.013190;
Y_plus_3_ff(x,B) = \
  +Y_plus_3_a1*(4*pi/(Y_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Y_plus_3_b1+B)*x*x) \
  +Y_plus_3_a2*(4*pi/(Y_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Y_plus_3_b2+B)*x*x) \
  +Y_plus_3_a3*(4*pi/(Y_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Y_plus_3_b3+B)*x*x) \
  +Y_plus_3_a4*(4*pi/(Y_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Y_plus_3_b4+B)*x*x) \
  +Y_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=40
Zr_c = 2.069290;
Zr_a1 = 17.876499;Zr_a2 = 10.948000;Zr_a3 = 5.417320;Zr_a4 = 3.657210;
Zr_b1 = 1.276180;Zr_b2 = 11.916000;Zr_b3 = 0.117622;Zr_b4 = 87.662697;
Zr_ff(x,B) = \
  +Zr_a1*(4*pi/(Zr_b1+B))**1.5*safexp(-4*pi**2/(Zr_b1+B)*x*x) \
  +Zr_a2*(4*pi/(Zr_b2+B))**1.5*safexp(-4*pi**2/(Zr_b2+B)*x*x) \
  +Zr_a3*(4*pi/(Zr_b3+B))**1.5*safexp(-4*pi**2/(Zr_b3+B)*x*x) \
  +Zr_a4*(4*pi/(Zr_b4+B))**1.5*safexp(-4*pi**2/(Zr_b4+B)*x*x) \
  +Zr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Zr_plus_4_c = 9.414539;
Zr_plus_4_a1 = 18.166800;Zr_plus_4_a2 = 10.056200;Zr_plus_4_a3 = 1.011180;Zr_plus_4_a4 = -2.647900;
Zr_plus_4_b1 = 1.214800;Zr_plus_4_b2 = 10.148300;Zr_plus_4_b3 = 21.605400;Zr_plus_4_b4 = -0.102760;
Zr_plus_4_ff(x,B) = \
  +Zr_plus_4_a1*(4*pi/(Zr_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Zr_plus_4_b1+B)*x*x) \
  +Zr_plus_4_a2*(4*pi/(Zr_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Zr_plus_4_b2+B)*x*x) \
  +Zr_plus_4_a3*(4*pi/(Zr_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Zr_plus_4_b3+B)*x*x) \
  +Zr_plus_4_a4*(4*pi/(Zr_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Zr_plus_4_b4+B)*x*x) \
  +Zr_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=41
Nb_c = 3.755910;
Nb_a1 = 17.614201;Nb_a2 = 12.014400;Nb_a3 = 4.041830;Nb_a4 = 3.533460;
Nb_b1 = 1.188650;Nb_b2 = 11.766000;Nb_b3 = 0.204785;Nb_b4 = 69.795700;
Nb_ff(x,B) = \
  +Nb_a1*(4*pi/(Nb_b1+B))**1.5*safexp(-4*pi**2/(Nb_b1+B)*x*x) \
  +Nb_a2*(4*pi/(Nb_b2+B))**1.5*safexp(-4*pi**2/(Nb_b2+B)*x*x) \
  +Nb_a3*(4*pi/(Nb_b3+B))**1.5*safexp(-4*pi**2/(Nb_b3+B)*x*x) \
  +Nb_a4*(4*pi/(Nb_b4+B))**1.5*safexp(-4*pi**2/(Nb_b4+B)*x*x) \
  +Nb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=38
Nb_plus_3_c = -12.912000;
Nb_plus_3_a1 = 19.881199;Nb_plus_3_a2 = 18.065300;Nb_plus_3_a3 = 11.017700;Nb_plus_3_a4 = 1.947150;
Nb_plus_3_b1 = 0.019175;Nb_plus_3_b2 = 1.133050;Nb_plus_3_b3 = 10.162100;Nb_plus_3_b4 = 28.338900;
Nb_plus_3_ff(x,B) = \
  +Nb_plus_3_a1*(4*pi/(Nb_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Nb_plus_3_b1+B)*x*x) \
  +Nb_plus_3_a2*(4*pi/(Nb_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Nb_plus_3_b2+B)*x*x) \
  +Nb_plus_3_a3*(4*pi/(Nb_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Nb_plus_3_b3+B)*x*x) \
  +Nb_plus_3_a4*(4*pi/(Nb_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Nb_plus_3_b4+B)*x*x) \
  +Nb_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Nb_plus_5_c = -6.393400;
Nb_plus_5_a1 = 17.916300;Nb_plus_5_a2 = 13.341700;Nb_plus_5_a3 = 10.799000;Nb_plus_5_a4 = 0.337905;
Nb_plus_5_b1 = 1.124460;Nb_plus_5_b2 = 0.028781;Nb_plus_5_b3 = 9.282060;Nb_plus_5_b4 = 25.722799;
Nb_plus_5_ff(x,B) = \
  +Nb_plus_5_a1*(4*pi/(Nb_plus_5_b1+B))**1.5*safexp(-4*pi**2/(Nb_plus_5_b1+B)*x*x) \
  +Nb_plus_5_a2*(4*pi/(Nb_plus_5_b2+B))**1.5*safexp(-4*pi**2/(Nb_plus_5_b2+B)*x*x) \
  +Nb_plus_5_a3*(4*pi/(Nb_plus_5_b3+B))**1.5*safexp(-4*pi**2/(Nb_plus_5_b3+B)*x*x) \
  +Nb_plus_5_a4*(4*pi/(Nb_plus_5_b4+B))**1.5*safexp(-4*pi**2/(Nb_plus_5_b4+B)*x*x) \
  +Nb_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=42
Mo_c = 4.387500;
Mo_a1 = 3.702500;Mo_a2 = 17.235600;Mo_a3 = 12.887600;Mo_a4 = 3.742900;
Mo_b1 = 0.277200;Mo_b2 = 1.095800;Mo_b3 = 11.004000;Mo_b4 = 61.658401;
Mo_ff(x,B) = \
  +Mo_a1*(4*pi/(Mo_b1+B))**1.5*safexp(-4*pi**2/(Mo_b1+B)*x*x) \
  +Mo_a2*(4*pi/(Mo_b2+B))**1.5*safexp(-4*pi**2/(Mo_b2+B)*x*x) \
  +Mo_a3*(4*pi/(Mo_b3+B))**1.5*safexp(-4*pi**2/(Mo_b3+B)*x*x) \
  +Mo_a4*(4*pi/(Mo_b4+B))**1.5*safexp(-4*pi**2/(Mo_b4+B)*x*x) \
  +Mo_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=39
Mo_plus_3_c = -14.421000;
Mo_plus_3_a1 = 21.166401;Mo_plus_3_a2 = 18.201700;Mo_plus_3_a3 = 11.742300;Mo_plus_3_a4 = 2.309510;
Mo_plus_3_b1 = 0.014734;Mo_plus_3_b2 = 1.030310;Mo_plus_3_b3 = 9.536590;Mo_plus_3_b4 = 26.630699;
Mo_plus_3_ff(x,B) = \
  +Mo_plus_3_a1*(4*pi/(Mo_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Mo_plus_3_b1+B)*x*x) \
  +Mo_plus_3_a2*(4*pi/(Mo_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Mo_plus_3_b2+B)*x*x) \
  +Mo_plus_3_a3*(4*pi/(Mo_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Mo_plus_3_b3+B)*x*x) \
  +Mo_plus_3_a4*(4*pi/(Mo_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Mo_plus_3_b4+B)*x*x) \
  +Mo_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=37
Mo_plus_5_c = -14.316000;
Mo_plus_5_a1 = 21.014900;Mo_plus_5_a2 = 18.099199;Mo_plus_5_a3 = 11.463200;Mo_plus_5_a4 = 0.740625;
Mo_plus_5_b1 = 0.014345;Mo_plus_5_b2 = 1.022380;Mo_plus_5_b3 = 8.788090;Mo_plus_5_b4 = 23.345200;
Mo_plus_5_ff(x,B) = \
  +Mo_plus_5_a1*(4*pi/(Mo_plus_5_b1+B))**1.5*safexp(-4*pi**2/(Mo_plus_5_b1+B)*x*x) \
  +Mo_plus_5_a2*(4*pi/(Mo_plus_5_b2+B))**1.5*safexp(-4*pi**2/(Mo_plus_5_b2+B)*x*x) \
  +Mo_plus_5_a3*(4*pi/(Mo_plus_5_b3+B))**1.5*safexp(-4*pi**2/(Mo_plus_5_b3+B)*x*x) \
  +Mo_plus_5_a4*(4*pi/(Mo_plus_5_b4+B))**1.5*safexp(-4*pi**2/(Mo_plus_5_b4+B)*x*x) \
  +Mo_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=36
Mo_plus_6_c = 0.344941;
Mo_plus_6_a1 = 17.887100;Mo_plus_6_a2 = 11.175000;Mo_plus_6_a3 = 6.578910;Mo_plus_6_a4 = 0.000000;
Mo_plus_6_b1 = 1.036490;Mo_plus_6_b2 = 8.480610;Mo_plus_6_b3 = 0.058881;Mo_plus_6_b4 = 0.000000;
Mo_plus_6_ff(x,B) = \
  +Mo_plus_6_a1*(4*pi/(Mo_plus_6_b1+B))**1.5*safexp(-4*pi**2/(Mo_plus_6_b1+B)*x*x) \
  +Mo_plus_6_a2*(4*pi/(Mo_plus_6_b2+B))**1.5*safexp(-4*pi**2/(Mo_plus_6_b2+B)*x*x) \
  +Mo_plus_6_a3*(4*pi/(Mo_plus_6_b3+B))**1.5*safexp(-4*pi**2/(Mo_plus_6_b3+B)*x*x) \
  +Mo_plus_6_a4*(4*pi/(Mo_plus_6_b4+B))**1.5*safexp(-4*pi**2/(Mo_plus_6_b4+B)*x*x) \
  +Mo_plus_6_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=43
Tc_c = 5.404280;
Tc_a1 = 19.130100;Tc_a2 = 11.094800;Tc_a3 = 4.649010;Tc_a4 = 2.712630;
Tc_b1 = 0.864132;Tc_b2 = 8.144870;Tc_b3 = 21.570700;Tc_b4 = 86.847198;
Tc_ff(x,B) = \
  +Tc_a1*(4*pi/(Tc_b1+B))**1.5*safexp(-4*pi**2/(Tc_b1+B)*x*x) \
  +Tc_a2*(4*pi/(Tc_b2+B))**1.5*safexp(-4*pi**2/(Tc_b2+B)*x*x) \
  +Tc_a3*(4*pi/(Tc_b3+B))**1.5*safexp(-4*pi**2/(Tc_b3+B)*x*x) \
  +Tc_a4*(4*pi/(Tc_b4+B))**1.5*safexp(-4*pi**2/(Tc_b4+B)*x*x) \
  +Tc_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=44
Ru_c = 5.378740;
Ru_a1 = 19.267399;Ru_a2 = 12.918200;Ru_a3 = 4.863370;Ru_a4 = 1.567560;
Ru_b1 = 0.808520;Ru_b2 = 8.434669;Ru_b3 = 24.799700;Ru_b4 = 94.292801;
Ru_ff(x,B) = \
  +Ru_a1*(4*pi/(Ru_b1+B))**1.5*safexp(-4*pi**2/(Ru_b1+B)*x*x) \
  +Ru_a2*(4*pi/(Ru_b2+B))**1.5*safexp(-4*pi**2/(Ru_b2+B)*x*x) \
  +Ru_a3*(4*pi/(Ru_b3+B))**1.5*safexp(-4*pi**2/(Ru_b3+B)*x*x) \
  +Ru_a4*(4*pi/(Ru_b4+B))**1.5*safexp(-4*pi**2/(Ru_b4+B)*x*x) \
  +Ru_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=41
Ru_plus_3_c = -3.189200;
Ru_plus_3_a1 = 18.563801;Ru_plus_3_a2 = 13.288500;Ru_plus_3_a3 = 9.326019;Ru_plus_3_a4 = 3.009640;
Ru_plus_3_b1 = 0.847329;Ru_plus_3_b2 = 8.371640;Ru_plus_3_b3 = 0.017662;Ru_plus_3_b4 = 22.886999;
Ru_plus_3_ff(x,B) = \
  +Ru_plus_3_a1*(4*pi/(Ru_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ru_plus_3_b1+B)*x*x) \
  +Ru_plus_3_a2*(4*pi/(Ru_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ru_plus_3_b2+B)*x*x) \
  +Ru_plus_3_a3*(4*pi/(Ru_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ru_plus_3_b3+B)*x*x) \
  +Ru_plus_3_a4*(4*pi/(Ru_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ru_plus_3_b4+B)*x*x) \
  +Ru_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=40
Ru_plus_4_c = 1.423570;
Ru_plus_4_a1 = 18.500299;Ru_plus_4_a2 = 13.178699;Ru_plus_4_a3 = 4.713040;Ru_plus_4_a4 = 2.185350;
Ru_plus_4_b1 = 0.844582;Ru_plus_4_b2 = 8.125340;Ru_plus_4_b3 = 0.364950;Ru_plus_4_b4 = 20.850399;
Ru_plus_4_ff(x,B) = \
  +Ru_plus_4_a1*(4*pi/(Ru_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Ru_plus_4_b1+B)*x*x) \
  +Ru_plus_4_a2*(4*pi/(Ru_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Ru_plus_4_b2+B)*x*x) \
  +Ru_plus_4_a3*(4*pi/(Ru_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Ru_plus_4_b3+B)*x*x) \
  +Ru_plus_4_a4*(4*pi/(Ru_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Ru_plus_4_b4+B)*x*x) \
  +Ru_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=45
Rh_c = 5.328000;
Rh_a1 = 19.295700;Rh_a2 = 14.350100;Rh_a3 = 4.734250;Rh_a4 = 1.289180;
Rh_b1 = 0.751536;Rh_b2 = 8.217580;Rh_b3 = 25.874901;Rh_b4 = 98.606201;
Rh_ff(x,B) = \
  +Rh_a1*(4*pi/(Rh_b1+B))**1.5*safexp(-4*pi**2/(Rh_b1+B)*x*x) \
  +Rh_a2*(4*pi/(Rh_b2+B))**1.5*safexp(-4*pi**2/(Rh_b2+B)*x*x) \
  +Rh_a3*(4*pi/(Rh_b3+B))**1.5*safexp(-4*pi**2/(Rh_b3+B)*x*x) \
  +Rh_a4*(4*pi/(Rh_b4+B))**1.5*safexp(-4*pi**2/(Rh_b4+B)*x*x) \
  +Rh_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=42
Rh_plus_3_c = 11.867800;
Rh_plus_3_a1 = 18.878500;Rh_plus_3_a2 = 14.125900;Rh_plus_3_a3 = 3.325150;Rh_plus_3_a4 = -6.198900;
Rh_plus_3_b1 = 0.764252;Rh_plus_3_b2 = 7.844380;Rh_plus_3_b3 = 21.248699;Rh_plus_3_b4 = -0.010360;
Rh_plus_3_ff(x,B) = \
  +Rh_plus_3_a1*(4*pi/(Rh_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Rh_plus_3_b1+B)*x*x) \
  +Rh_plus_3_a2*(4*pi/(Rh_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Rh_plus_3_b2+B)*x*x) \
  +Rh_plus_3_a3*(4*pi/(Rh_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Rh_plus_3_b3+B)*x*x) \
  +Rh_plus_3_a4*(4*pi/(Rh_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Rh_plus_3_b4+B)*x*x) \
  +Rh_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=41
Rh_plus_4_c = 11.283500;
Rh_plus_4_a1 = 18.854500;Rh_plus_4_a2 = 13.980600;Rh_plus_4_a3 = 2.534640;Rh_plus_4_a4 = -5.652600;
Rh_plus_4_b1 = 0.760825;Rh_plus_4_b2 = 7.624360;Rh_plus_4_b3 = 19.331699;Rh_plus_4_b4 = -0.010200;
Rh_plus_4_ff(x,B) = \
  +Rh_plus_4_a1*(4*pi/(Rh_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Rh_plus_4_b1+B)*x*x) \
  +Rh_plus_4_a2*(4*pi/(Rh_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Rh_plus_4_b2+B)*x*x) \
  +Rh_plus_4_a3*(4*pi/(Rh_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Rh_plus_4_b3+B)*x*x) \
  +Rh_plus_4_a4*(4*pi/(Rh_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Rh_plus_4_b4+B)*x*x) \
  +Rh_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
Pd_c = 5.265930;
Pd_a1 = 19.331900;Pd_a2 = 15.501699;Pd_a3 = 5.295370;Pd_a4 = 0.605844;
Pd_b1 = 0.698655;Pd_b2 = 7.989290;Pd_b3 = 25.205200;Pd_b4 = 76.898598;
Pd_ff(x,B) = \
  +Pd_a1*(4*pi/(Pd_b1+B))**1.5*safexp(-4*pi**2/(Pd_b1+B)*x*x) \
  +Pd_a2*(4*pi/(Pd_b2+B))**1.5*safexp(-4*pi**2/(Pd_b2+B)*x*x) \
  +Pd_a3*(4*pi/(Pd_b3+B))**1.5*safexp(-4*pi**2/(Pd_b3+B)*x*x) \
  +Pd_a4*(4*pi/(Pd_b4+B))**1.5*safexp(-4*pi**2/(Pd_b4+B)*x*x) \
  +Pd_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=44
Pd_plus_2_c = 5.291600;
Pd_plus_2_a1 = 19.170099;Pd_plus_2_a2 = 15.209600;Pd_plus_2_a3 = 4.322340;Pd_plus_2_a4 = 0.000000;
Pd_plus_2_b1 = 0.696219;Pd_plus_2_b2 = 7.555730;Pd_plus_2_b3 = 22.505699;Pd_plus_2_b4 = 0.000000;
Pd_plus_2_ff(x,B) = \
  +Pd_plus_2_a1*(4*pi/(Pd_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Pd_plus_2_b1+B)*x*x) \
  +Pd_plus_2_a2*(4*pi/(Pd_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Pd_plus_2_b2+B)*x*x) \
  +Pd_plus_2_a3*(4*pi/(Pd_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Pd_plus_2_b3+B)*x*x) \
  +Pd_plus_2_a4*(4*pi/(Pd_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Pd_plus_2_b4+B)*x*x) \
  +Pd_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=42
Pd_plus_4_c = 13.017400;
Pd_plus_4_a1 = 19.249300;Pd_plus_4_a2 = 14.790000;Pd_plus_4_a3 = 2.892890;Pd_plus_4_a4 = -7.949200;
Pd_plus_4_b1 = 0.683839;Pd_plus_4_b2 = 7.148330;Pd_plus_4_b3 = 17.914400;Pd_plus_4_b4 = 0.005127;
Pd_plus_4_ff(x,B) = \
  +Pd_plus_4_a1*(4*pi/(Pd_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Pd_plus_4_b1+B)*x*x) \
  +Pd_plus_4_a2*(4*pi/(Pd_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Pd_plus_4_b2+B)*x*x) \
  +Pd_plus_4_a3*(4*pi/(Pd_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Pd_plus_4_b3+B)*x*x) \
  +Pd_plus_4_a4*(4*pi/(Pd_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Pd_plus_4_b4+B)*x*x) \
  +Pd_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=47
Ag_c = 5.179000;
Ag_a1 = 19.280800;Ag_a2 = 16.688499;Ag_a3 = 4.804500;Ag_a4 = 1.046300;
Ag_b1 = 0.644600;Ag_b2 = 7.472600;Ag_b3 = 24.660500;Ag_b4 = 99.815598;
Ag_ff(x,B) = \
  +Ag_a1*(4*pi/(Ag_b1+B))**1.5*safexp(-4*pi**2/(Ag_b1+B)*x*x) \
  +Ag_a2*(4*pi/(Ag_b2+B))**1.5*safexp(-4*pi**2/(Ag_b2+B)*x*x) \
  +Ag_a3*(4*pi/(Ag_b3+B))**1.5*safexp(-4*pi**2/(Ag_b3+B)*x*x) \
  +Ag_a4*(4*pi/(Ag_b4+B))**1.5*safexp(-4*pi**2/(Ag_b4+B)*x*x) \
  +Ag_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
Ag_plus_1_c = 5.215720;
Ag_plus_1_a1 = 19.181200;Ag_plus_1_a2 = 15.971900;Ag_plus_1_a3 = 5.274750;Ag_plus_1_a4 = 0.357534;
Ag_plus_1_b1 = 0.646179;Ag_plus_1_b2 = 7.191230;Ag_plus_1_b3 = 21.732599;Ag_plus_1_b4 = 66.114700;
Ag_plus_1_ff(x,B) = \
  +Ag_plus_1_a1*(4*pi/(Ag_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Ag_plus_1_b1+B)*x*x) \
  +Ag_plus_1_a2*(4*pi/(Ag_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Ag_plus_1_b2+B)*x*x) \
  +Ag_plus_1_a3*(4*pi/(Ag_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Ag_plus_1_b3+B)*x*x) \
  +Ag_plus_1_a4*(4*pi/(Ag_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Ag_plus_1_b4+B)*x*x) \
  +Ag_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=45
Ag_plus_2_c = 5.214040;
Ag_plus_2_a1 = 19.164299;Ag_plus_2_a2 = 16.245600;Ag_plus_2_a3 = 4.370900;Ag_plus_2_a4 = 0.000000;
Ag_plus_2_b1 = 0.645643;Ag_plus_2_b2 = 7.185440;Ag_plus_2_b3 = 21.407200;Ag_plus_2_b4 = 0.000000;
Ag_plus_2_ff(x,B) = \
  +Ag_plus_2_a1*(4*pi/(Ag_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ag_plus_2_b1+B)*x*x) \
  +Ag_plus_2_a2*(4*pi/(Ag_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ag_plus_2_b2+B)*x*x) \
  +Ag_plus_2_a3*(4*pi/(Ag_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ag_plus_2_b3+B)*x*x) \
  +Ag_plus_2_a4*(4*pi/(Ag_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ag_plus_2_b4+B)*x*x) \
  +Ag_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=48
Cd_c = 5.069400;
Cd_a1 = 19.221399;Cd_a2 = 17.644400;Cd_a3 = 4.461000;Cd_a4 = 1.602900;
Cd_b1 = 0.594600;Cd_b2 = 6.908900;Cd_b3 = 24.700800;Cd_b4 = 87.482498;
Cd_ff(x,B) = \
  +Cd_a1*(4*pi/(Cd_b1+B))**1.5*safexp(-4*pi**2/(Cd_b1+B)*x*x) \
  +Cd_a2*(4*pi/(Cd_b2+B))**1.5*safexp(-4*pi**2/(Cd_b2+B)*x*x) \
  +Cd_a3*(4*pi/(Cd_b3+B))**1.5*safexp(-4*pi**2/(Cd_b3+B)*x*x) \
  +Cd_a4*(4*pi/(Cd_b4+B))**1.5*safexp(-4*pi**2/(Cd_b4+B)*x*x) \
  +Cd_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
Cd_plus_2_c = 5.119370;
Cd_plus_2_a1 = 19.151400;Cd_plus_2_a2 = 17.253500;Cd_plus_2_a3 = 4.471280;Cd_plus_2_a4 = 0.000000;
Cd_plus_2_b1 = 0.597922;Cd_plus_2_b2 = 6.806390;Cd_plus_2_b3 = 20.252100;Cd_plus_2_b4 = 0.000000;
Cd_plus_2_ff(x,B) = \
  +Cd_plus_2_a1*(4*pi/(Cd_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Cd_plus_2_b1+B)*x*x) \
  +Cd_plus_2_a2*(4*pi/(Cd_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Cd_plus_2_b2+B)*x*x) \
  +Cd_plus_2_a3*(4*pi/(Cd_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Cd_plus_2_b3+B)*x*x) \
  +Cd_plus_2_a4*(4*pi/(Cd_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Cd_plus_2_b4+B)*x*x) \
  +Cd_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=49
In_c = 4.939100;
In_a1 = 19.162399;In_a2 = 18.559601;In_a3 = 4.294800;In_a4 = 2.039600;
In_b1 = 0.547600;In_b2 = 6.377600;In_b3 = 25.849899;In_b4 = 92.802902;
In_ff(x,B) = \
  +In_a1*(4*pi/(In_b1+B))**1.5*safexp(-4*pi**2/(In_b1+B)*x*x) \
  +In_a2*(4*pi/(In_b2+B))**1.5*safexp(-4*pi**2/(In_b2+B)*x*x) \
  +In_a3*(4*pi/(In_b3+B))**1.5*safexp(-4*pi**2/(In_b3+B)*x*x) \
  +In_a4*(4*pi/(In_b4+B))**1.5*safexp(-4*pi**2/(In_b4+B)*x*x) \
  +In_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
In_plus_3_c = 4.996350;
In_plus_3_a1 = 19.104500;In_plus_3_a2 = 18.110800;In_plus_3_a3 = 3.788970;In_plus_3_a4 = 0.000000;
In_plus_3_b1 = 0.551522;In_plus_3_b2 = 6.324700;In_plus_3_b3 = 17.359501;In_plus_3_b4 = 0.000000;
In_plus_3_ff(x,B) = \
  +In_plus_3_a1*(4*pi/(In_plus_3_b1+B))**1.5*safexp(-4*pi**2/(In_plus_3_b1+B)*x*x) \
  +In_plus_3_a2*(4*pi/(In_plus_3_b2+B))**1.5*safexp(-4*pi**2/(In_plus_3_b2+B)*x*x) \
  +In_plus_3_a3*(4*pi/(In_plus_3_b3+B))**1.5*safexp(-4*pi**2/(In_plus_3_b3+B)*x*x) \
  +In_plus_3_a4*(4*pi/(In_plus_3_b4+B))**1.5*safexp(-4*pi**2/(In_plus_3_b4+B)*x*x) \
  +In_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=50
Sn_c = 4.782100;
Sn_a1 = 19.188900;Sn_a2 = 19.100500;Sn_a3 = 4.458500;Sn_a4 = 2.466300;
Sn_b1 = 5.830300;Sn_b2 = 0.503100;Sn_b3 = 26.890900;Sn_b4 = 83.957100;
Sn_ff(x,B) = \
  +Sn_a1*(4*pi/(Sn_b1+B))**1.5*safexp(-4*pi**2/(Sn_b1+B)*x*x) \
  +Sn_a2*(4*pi/(Sn_b2+B))**1.5*safexp(-4*pi**2/(Sn_b2+B)*x*x) \
  +Sn_a3*(4*pi/(Sn_b3+B))**1.5*safexp(-4*pi**2/(Sn_b3+B)*x*x) \
  +Sn_a4*(4*pi/(Sn_b4+B))**1.5*safexp(-4*pi**2/(Sn_b4+B)*x*x) \
  +Sn_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=48
Sn_plus_2_c = 4.786100;
Sn_plus_2_a1 = 19.109400;Sn_plus_2_a2 = 19.054800;Sn_plus_2_a3 = 4.564800;Sn_plus_2_a4 = 0.487000;
Sn_plus_2_b1 = 0.503600;Sn_plus_2_b2 = 5.837800;Sn_plus_2_b3 = 23.375200;Sn_plus_2_b4 = 62.206100;
Sn_plus_2_ff(x,B) = \
  +Sn_plus_2_a1*(4*pi/(Sn_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Sn_plus_2_b1+B)*x*x) \
  +Sn_plus_2_a2*(4*pi/(Sn_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Sn_plus_2_b2+B)*x*x) \
  +Sn_plus_2_a3*(4*pi/(Sn_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Sn_plus_2_b3+B)*x*x) \
  +Sn_plus_2_a4*(4*pi/(Sn_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Sn_plus_2_b4+B)*x*x) \
  +Sn_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
Sn_plus_4_c = 3.918200;
Sn_plus_4_a1 = 18.933300;Sn_plus_4_a2 = 19.713100;Sn_plus_4_a3 = 3.418200;Sn_plus_4_a4 = 0.019300;
Sn_plus_4_b1 = 5.764000;Sn_plus_4_b2 = 0.465500;Sn_plus_4_b3 = 14.004900;Sn_plus_4_b4 = -0.758300;
Sn_plus_4_ff(x,B) = \
  +Sn_plus_4_a1*(4*pi/(Sn_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Sn_plus_4_b1+B)*x*x) \
  +Sn_plus_4_a2*(4*pi/(Sn_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Sn_plus_4_b2+B)*x*x) \
  +Sn_plus_4_a3*(4*pi/(Sn_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Sn_plus_4_b3+B)*x*x) \
  +Sn_plus_4_a4*(4*pi/(Sn_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Sn_plus_4_b4+B)*x*x) \
  +Sn_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=51
Sb_c = 4.590900;
Sb_a1 = 19.641800;Sb_a2 = 19.045500;Sb_a3 = 5.037100;Sb_a4 = 2.682700;
Sb_b1 = 5.303400;Sb_b2 = 0.460700;Sb_b3 = 27.907400;Sb_b4 = 75.282501;
Sb_ff(x,B) = \
  +Sb_a1*(4*pi/(Sb_b1+B))**1.5*safexp(-4*pi**2/(Sb_b1+B)*x*x) \
  +Sb_a2*(4*pi/(Sb_b2+B))**1.5*safexp(-4*pi**2/(Sb_b2+B)*x*x) \
  +Sb_a3*(4*pi/(Sb_b3+B))**1.5*safexp(-4*pi**2/(Sb_b3+B)*x*x) \
  +Sb_a4*(4*pi/(Sb_b4+B))**1.5*safexp(-4*pi**2/(Sb_b4+B)*x*x) \
  +Sb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=48
Sb_plus_3_c = 4.696260;
Sb_plus_3_a1 = 18.975500;Sb_plus_3_a2 = 18.932999;Sb_plus_3_a3 = 5.107890;Sb_plus_3_a4 = 0.288753;
Sb_plus_3_b1 = 0.467196;Sb_plus_3_b2 = 5.221260;Sb_plus_3_b3 = 19.590200;Sb_plus_3_b4 = 55.511299;
Sb_plus_3_ff(x,B) = \
  +Sb_plus_3_a1*(4*pi/(Sb_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Sb_plus_3_b1+B)*x*x) \
  +Sb_plus_3_a2*(4*pi/(Sb_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Sb_plus_3_b2+B)*x*x) \
  +Sb_plus_3_a3*(4*pi/(Sb_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Sb_plus_3_b3+B)*x*x) \
  +Sb_plus_3_a4*(4*pi/(Sb_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Sb_plus_3_b4+B)*x*x) \
  +Sb_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=46
Sb_plus_5_c = 4.692630;
Sb_plus_5_a1 = 19.868500;Sb_plus_5_a2 = 19.030199;Sb_plus_5_a3 = 2.412530;Sb_plus_5_a4 = 0.000000;
Sb_plus_5_b1 = 5.448530;Sb_plus_5_b2 = 0.467973;Sb_plus_5_b3 = 14.125900;Sb_plus_5_b4 = 0.000000;
Sb_plus_5_ff(x,B) = \
  +Sb_plus_5_a1*(4*pi/(Sb_plus_5_b1+B))**1.5*safexp(-4*pi**2/(Sb_plus_5_b1+B)*x*x) \
  +Sb_plus_5_a2*(4*pi/(Sb_plus_5_b2+B))**1.5*safexp(-4*pi**2/(Sb_plus_5_b2+B)*x*x) \
  +Sb_plus_5_a3*(4*pi/(Sb_plus_5_b3+B))**1.5*safexp(-4*pi**2/(Sb_plus_5_b3+B)*x*x) \
  +Sb_plus_5_a4*(4*pi/(Sb_plus_5_b4+B))**1.5*safexp(-4*pi**2/(Sb_plus_5_b4+B)*x*x) \
  +Sb_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=52
Te_c = 4.352000;
Te_a1 = 19.964399;Te_a2 = 19.013800;Te_a3 = 6.144870;Te_a4 = 2.523900;
Te_b1 = 4.817420;Te_b2 = 0.420885;Te_b3 = 28.528400;Te_b4 = 70.840302;
Te_ff(x,B) = \
  +Te_a1*(4*pi/(Te_b1+B))**1.5*safexp(-4*pi**2/(Te_b1+B)*x*x) \
  +Te_a2*(4*pi/(Te_b2+B))**1.5*safexp(-4*pi**2/(Te_b2+B)*x*x) \
  +Te_a3*(4*pi/(Te_b3+B))**1.5*safexp(-4*pi**2/(Te_b3+B)*x*x) \
  +Te_a4*(4*pi/(Te_b4+B))**1.5*safexp(-4*pi**2/(Te_b4+B)*x*x) \
  +Te_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=53
I_c = 4.071200;
I_a1 = 20.147200;I_a2 = 18.994900;I_a3 = 7.513800;I_a4 = 2.273500;
I_b1 = 4.347000;I_b2 = 0.381400;I_b3 = 27.765999;I_b4 = 66.877602;
I_ff(x,B) = \
  +I_a1*(4*pi/(I_b1+B))**1.5*safexp(-4*pi**2/(I_b1+B)*x*x) \
  +I_a2*(4*pi/(I_b2+B))**1.5*safexp(-4*pi**2/(I_b2+B)*x*x) \
  +I_a3*(4*pi/(I_b3+B))**1.5*safexp(-4*pi**2/(I_b3+B)*x*x) \
  +I_a4*(4*pi/(I_b4+B))**1.5*safexp(-4*pi**2/(I_b4+B)*x*x) \
  +I_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
I_minus_1_c = 4.071400;
I_minus_1_a1 = 20.233200;I_minus_1_a2 = 18.997000;I_minus_1_a3 = 7.806900;I_minus_1_a4 = 2.886800;
I_minus_1_b1 = 4.357900;I_minus_1_b2 = 0.381500;I_minus_1_b3 = 29.525900;I_minus_1_b4 = 84.930397;
I_minus_1_ff(x,B) = \
  +I_minus_1_a1*(4*pi/(I_minus_1_b1+B))**1.5*safexp(-4*pi**2/(I_minus_1_b1+B)*x*x) \
  +I_minus_1_a2*(4*pi/(I_minus_1_b2+B))**1.5*safexp(-4*pi**2/(I_minus_1_b2+B)*x*x) \
  +I_minus_1_a3*(4*pi/(I_minus_1_b3+B))**1.5*safexp(-4*pi**2/(I_minus_1_b3+B)*x*x) \
  +I_minus_1_a4*(4*pi/(I_minus_1_b4+B))**1.5*safexp(-4*pi**2/(I_minus_1_b4+B)*x*x) \
  +I_minus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
Xe_c = 3.711800;
Xe_a1 = 20.293301;Xe_a2 = 19.029800;Xe_a3 = 8.976700;Xe_a4 = 1.990000;
Xe_b1 = 3.928200;Xe_b2 = 0.344000;Xe_b3 = 26.465900;Xe_b4 = 64.265800;
Xe_ff(x,B) = \
  +Xe_a1*(4*pi/(Xe_b1+B))**1.5*safexp(-4*pi**2/(Xe_b1+B)*x*x) \
  +Xe_a2*(4*pi/(Xe_b2+B))**1.5*safexp(-4*pi**2/(Xe_b2+B)*x*x) \
  +Xe_a3*(4*pi/(Xe_b3+B))**1.5*safexp(-4*pi**2/(Xe_b3+B)*x*x) \
  +Xe_a4*(4*pi/(Xe_b4+B))**1.5*safexp(-4*pi**2/(Xe_b4+B)*x*x) \
  +Xe_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=55
Cs_c = 3.335200;
Cs_a1 = 20.389200;Cs_a2 = 19.106199;Cs_a3 = 10.662000;Cs_a4 = 1.495300;
Cs_b1 = 3.569000;Cs_b2 = 0.310700;Cs_b3 = 24.387899;Cs_b4 = 213.903992;
Cs_ff(x,B) = \
  +Cs_a1*(4*pi/(Cs_b1+B))**1.5*safexp(-4*pi**2/(Cs_b1+B)*x*x) \
  +Cs_a2*(4*pi/(Cs_b2+B))**1.5*safexp(-4*pi**2/(Cs_b2+B)*x*x) \
  +Cs_a3*(4*pi/(Cs_b3+B))**1.5*safexp(-4*pi**2/(Cs_b3+B)*x*x) \
  +Cs_a4*(4*pi/(Cs_b4+B))**1.5*safexp(-4*pi**2/(Cs_b4+B)*x*x) \
  +Cs_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
Cs_plus_1_c = 3.279100;
Cs_plus_1_a1 = 20.352400;Cs_plus_1_a2 = 19.127800;Cs_plus_1_a3 = 10.282100;Cs_plus_1_a4 = 0.961500;
Cs_plus_1_b1 = 3.552000;Cs_plus_1_b2 = 0.308600;Cs_plus_1_b3 = 23.712799;Cs_plus_1_b4 = 59.456497;
Cs_plus_1_ff(x,B) = \
  +Cs_plus_1_a1*(4*pi/(Cs_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Cs_plus_1_b1+B)*x*x) \
  +Cs_plus_1_a2*(4*pi/(Cs_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Cs_plus_1_b2+B)*x*x) \
  +Cs_plus_1_a3*(4*pi/(Cs_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Cs_plus_1_b3+B)*x*x) \
  +Cs_plus_1_a4*(4*pi/(Cs_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Cs_plus_1_b4+B)*x*x) \
  +Cs_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=56
Ba_c = 2.773100;
Ba_a1 = 20.336100;Ba_a2 = 19.297001;Ba_a3 = 10.888000;Ba_a4 = 2.695900;
Ba_b1 = 3.216000;Ba_b2 = 0.275600;Ba_b3 = 20.207300;Ba_b4 = 167.201996;
Ba_ff(x,B) = \
  +Ba_a1*(4*pi/(Ba_b1+B))**1.5*safexp(-4*pi**2/(Ba_b1+B)*x*x) \
  +Ba_a2*(4*pi/(Ba_b2+B))**1.5*safexp(-4*pi**2/(Ba_b2+B)*x*x) \
  +Ba_a3*(4*pi/(Ba_b3+B))**1.5*safexp(-4*pi**2/(Ba_b3+B)*x*x) \
  +Ba_a4*(4*pi/(Ba_b4+B))**1.5*safexp(-4*pi**2/(Ba_b4+B)*x*x) \
  +Ba_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
Ba_plus_2_c = 3.029020;
Ba_plus_2_a1 = 20.180700;Ba_plus_2_a2 = 19.113600;Ba_plus_2_a3 = 10.905399;Ba_plus_2_a4 = 0.776340;
Ba_plus_2_b1 = 3.213670;Ba_plus_2_b2 = 0.283310;Ba_plus_2_b3 = 20.055799;Ba_plus_2_b4 = 51.745998;
Ba_plus_2_ff(x,B) = \
  +Ba_plus_2_a1*(4*pi/(Ba_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ba_plus_2_b1+B)*x*x) \
  +Ba_plus_2_a2*(4*pi/(Ba_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ba_plus_2_b2+B)*x*x) \
  +Ba_plus_2_a3*(4*pi/(Ba_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ba_plus_2_b3+B)*x*x) \
  +Ba_plus_2_a4*(4*pi/(Ba_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ba_plus_2_b4+B)*x*x) \
  +Ba_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=57
La_c = 2.146780;
La_a1 = 20.577999;La_a2 = 19.598999;La_a3 = 11.372700;La_a4 = 3.287190;
La_b1 = 2.948170;La_b2 = 0.244475;La_b3 = 18.772600;La_b4 = 133.123993;
La_ff(x,B) = \
  +La_a1*(4*pi/(La_b1+B))**1.5*safexp(-4*pi**2/(La_b1+B)*x*x) \
  +La_a2*(4*pi/(La_b2+B))**1.5*safexp(-4*pi**2/(La_b2+B)*x*x) \
  +La_a3*(4*pi/(La_b3+B))**1.5*safexp(-4*pi**2/(La_b3+B)*x*x) \
  +La_a4*(4*pi/(La_b4+B))**1.5*safexp(-4*pi**2/(La_b4+B)*x*x) \
  +La_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
La_plus_3_c = 2.408600;
La_plus_3_a1 = 20.248899;La_plus_3_a2 = 19.376301;La_plus_3_a3 = 11.632299;La_plus_3_a4 = 0.336048;
La_plus_3_b1 = 2.920700;La_plus_3_b2 = 0.250698;La_plus_3_b3 = 17.821100;La_plus_3_b4 = 54.945297;
La_plus_3_ff(x,B) = \
  +La_plus_3_a1*(4*pi/(La_plus_3_b1+B))**1.5*safexp(-4*pi**2/(La_plus_3_b1+B)*x*x) \
  +La_plus_3_a2*(4*pi/(La_plus_3_b2+B))**1.5*safexp(-4*pi**2/(La_plus_3_b2+B)*x*x) \
  +La_plus_3_a3*(4*pi/(La_plus_3_b3+B))**1.5*safexp(-4*pi**2/(La_plus_3_b3+B)*x*x) \
  +La_plus_3_a4*(4*pi/(La_plus_3_b4+B))**1.5*safexp(-4*pi**2/(La_plus_3_b4+B)*x*x) \
  +La_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=58
Ce_c = 1.862640;
Ce_a1 = 21.167099;Ce_a2 = 19.769501;Ce_a3 = 11.851299;Ce_a4 = 3.330490;
Ce_b1 = 2.812190;Ce_b2 = 0.226836;Ce_b3 = 17.608299;Ce_b4 = 127.112999;
Ce_ff(x,B) = \
  +Ce_a1*(4*pi/(Ce_b1+B))**1.5*safexp(-4*pi**2/(Ce_b1+B)*x*x) \
  +Ce_a2*(4*pi/(Ce_b2+B))**1.5*safexp(-4*pi**2/(Ce_b2+B)*x*x) \
  +Ce_a3*(4*pi/(Ce_b3+B))**1.5*safexp(-4*pi**2/(Ce_b3+B)*x*x) \
  +Ce_a4*(4*pi/(Ce_b4+B))**1.5*safexp(-4*pi**2/(Ce_b4+B)*x*x) \
  +Ce_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=55
Ce_plus_3_c = 2.090130;
Ce_plus_3_a1 = 20.803600;Ce_plus_3_a2 = 19.559000;Ce_plus_3_a3 = 11.936900;Ce_plus_3_a4 = 0.612376;
Ce_plus_3_b1 = 2.776910;Ce_plus_3_b2 = 0.231540;Ce_plus_3_b3 = 16.540800;Ce_plus_3_b4 = 43.169201;
Ce_plus_3_ff(x,B) = \
  +Ce_plus_3_a1*(4*pi/(Ce_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ce_plus_3_b1+B)*x*x) \
  +Ce_plus_3_a2*(4*pi/(Ce_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ce_plus_3_b2+B)*x*x) \
  +Ce_plus_3_a3*(4*pi/(Ce_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ce_plus_3_b3+B)*x*x) \
  +Ce_plus_3_a4*(4*pi/(Ce_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ce_plus_3_b4+B)*x*x) \
  +Ce_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=54
Ce_plus_4_c = 1.591800;
Ce_plus_4_a1 = 20.323500;Ce_plus_4_a2 = 19.818600;Ce_plus_4_a3 = 12.123300;Ce_plus_4_a4 = 0.144583;
Ce_plus_4_b1 = 2.659410;Ce_plus_4_b2 = 0.218850;Ce_plus_4_b3 = 15.799200;Ce_plus_4_b4 = 62.235500;
Ce_plus_4_ff(x,B) = \
  +Ce_plus_4_a1*(4*pi/(Ce_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Ce_plus_4_b1+B)*x*x) \
  +Ce_plus_4_a2*(4*pi/(Ce_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Ce_plus_4_b2+B)*x*x) \
  +Ce_plus_4_a3*(4*pi/(Ce_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Ce_plus_4_b3+B)*x*x) \
  +Ce_plus_4_a4*(4*pi/(Ce_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Ce_plus_4_b4+B)*x*x) \
  +Ce_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=59
Pr_c = 2.058300;
Pr_a1 = 22.043999;Pr_a2 = 19.669701;Pr_a3 = 12.385600;Pr_a4 = 2.824280;
Pr_b1 = 2.773930;Pr_b2 = 0.222087;Pr_b3 = 16.766899;Pr_b4 = 143.643997;
Pr_ff(x,B) = \
  +Pr_a1*(4*pi/(Pr_b1+B))**1.5*safexp(-4*pi**2/(Pr_b1+B)*x*x) \
  +Pr_a2*(4*pi/(Pr_b2+B))**1.5*safexp(-4*pi**2/(Pr_b2+B)*x*x) \
  +Pr_a3*(4*pi/(Pr_b3+B))**1.5*safexp(-4*pi**2/(Pr_b3+B)*x*x) \
  +Pr_a4*(4*pi/(Pr_b4+B))**1.5*safexp(-4*pi**2/(Pr_b4+B)*x*x) \
  +Pr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=56
Pr_plus_3_c = 1.771320;
Pr_plus_3_a1 = 21.372700;Pr_plus_3_a2 = 19.749100;Pr_plus_3_a3 = 12.132900;Pr_plus_3_a4 = 0.975180;
Pr_plus_3_b1 = 2.645200;Pr_plus_3_b2 = 0.214299;Pr_plus_3_b3 = 15.323000;Pr_plus_3_b4 = 36.406502;
Pr_plus_3_ff(x,B) = \
  +Pr_plus_3_a1*(4*pi/(Pr_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Pr_plus_3_b1+B)*x*x) \
  +Pr_plus_3_a2*(4*pi/(Pr_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Pr_plus_3_b2+B)*x*x) \
  +Pr_plus_3_a3*(4*pi/(Pr_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Pr_plus_3_b3+B)*x*x) \
  +Pr_plus_3_a4*(4*pi/(Pr_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Pr_plus_3_b4+B)*x*x) \
  +Pr_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=55
Pr_plus_4_c = 1.242850;
Pr_plus_4_a1 = 20.941299;Pr_plus_4_a2 = 20.053900;Pr_plus_4_a3 = 12.466800;Pr_plus_4_a4 = 0.296689;
Pr_plus_4_b1 = 2.544670;Pr_plus_4_b2 = 0.202481;Pr_plus_4_b3 = 14.813700;Pr_plus_4_b4 = 45.464298;
Pr_plus_4_ff(x,B) = \
  +Pr_plus_4_a1*(4*pi/(Pr_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Pr_plus_4_b1+B)*x*x) \
  +Pr_plus_4_a2*(4*pi/(Pr_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Pr_plus_4_b2+B)*x*x) \
  +Pr_plus_4_a3*(4*pi/(Pr_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Pr_plus_4_b3+B)*x*x) \
  +Pr_plus_4_a4*(4*pi/(Pr_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Pr_plus_4_b4+B)*x*x) \
  +Pr_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=60
Nd_c = 1.984860;
Nd_a1 = 22.684500;Nd_a2 = 19.684700;Nd_a3 = 12.774000;Nd_a4 = 2.851370;
Nd_b1 = 2.662480;Nd_b2 = 0.210628;Nd_b3 = 15.885000;Nd_b4 = 137.903000;
Nd_ff(x,B) = \
  +Nd_a1*(4*pi/(Nd_b1+B))**1.5*safexp(-4*pi**2/(Nd_b1+B)*x*x) \
  +Nd_a2*(4*pi/(Nd_b2+B))**1.5*safexp(-4*pi**2/(Nd_b2+B)*x*x) \
  +Nd_a3*(4*pi/(Nd_b3+B))**1.5*safexp(-4*pi**2/(Nd_b3+B)*x*x) \
  +Nd_a4*(4*pi/(Nd_b4+B))**1.5*safexp(-4*pi**2/(Nd_b4+B)*x*x) \
  +Nd_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=57
Nd_plus_3_c = 1.475880;
Nd_plus_3_a1 = 21.961000;Nd_plus_3_a2 = 19.933899;Nd_plus_3_a3 = 12.120000;Nd_plus_3_a4 = 1.510310;
Nd_plus_3_b1 = 2.527220;Nd_plus_3_b2 = 0.199237;Nd_plus_3_b3 = 14.178300;Nd_plus_3_b4 = 30.871700;
Nd_plus_3_ff(x,B) = \
  +Nd_plus_3_a1*(4*pi/(Nd_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Nd_plus_3_b1+B)*x*x) \
  +Nd_plus_3_a2*(4*pi/(Nd_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Nd_plus_3_b2+B)*x*x) \
  +Nd_plus_3_a3*(4*pi/(Nd_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Nd_plus_3_b3+B)*x*x) \
  +Nd_plus_3_a4*(4*pi/(Nd_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Nd_plus_3_b4+B)*x*x) \
  +Nd_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=61
Pm_c = 2.028760;
Pm_a1 = 23.340500;Pm_a2 = 19.609501;Pm_a3 = 13.123500;Pm_a4 = 2.875160;
Pm_b1 = 2.562700;Pm_b2 = 0.202088;Pm_b3 = 15.100900;Pm_b4 = 132.720993;
Pm_ff(x,B) = \
  +Pm_a1*(4*pi/(Pm_b1+B))**1.5*safexp(-4*pi**2/(Pm_b1+B)*x*x) \
  +Pm_a2*(4*pi/(Pm_b2+B))**1.5*safexp(-4*pi**2/(Pm_b2+B)*x*x) \
  +Pm_a3*(4*pi/(Pm_b3+B))**1.5*safexp(-4*pi**2/(Pm_b3+B)*x*x) \
  +Pm_a4*(4*pi/(Pm_b4+B))**1.5*safexp(-4*pi**2/(Pm_b4+B)*x*x) \
  +Pm_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=58
Pm_plus_3_c = 1.194990;
Pm_plus_3_a1 = 22.552700;Pm_plus_3_a2 = 20.110800;Pm_plus_3_a3 = 12.067100;Pm_plus_3_a4 = 2.074920;
Pm_plus_3_b1 = 2.417400;Pm_plus_3_b2 = 0.185769;Pm_plus_3_b3 = 13.127500;Pm_plus_3_b4 = 27.449100;
Pm_plus_3_ff(x,B) = \
  +Pm_plus_3_a1*(4*pi/(Pm_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Pm_plus_3_b1+B)*x*x) \
  +Pm_plus_3_a2*(4*pi/(Pm_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Pm_plus_3_b2+B)*x*x) \
  +Pm_plus_3_a3*(4*pi/(Pm_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Pm_plus_3_b3+B)*x*x) \
  +Pm_plus_3_a4*(4*pi/(Pm_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Pm_plus_3_b4+B)*x*x) \
  +Pm_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=62
Sm_c = 2.209630;
Sm_a1 = 24.004200;Sm_a2 = 19.425800;Sm_a3 = 13.439600;Sm_a4 = 2.896040;
Sm_b1 = 2.472740;Sm_b2 = 0.196451;Sm_b3 = 14.399600;Sm_b4 = 128.007004;
Sm_ff(x,B) = \
  +Sm_a1*(4*pi/(Sm_b1+B))**1.5*safexp(-4*pi**2/(Sm_b1+B)*x*x) \
  +Sm_a2*(4*pi/(Sm_b2+B))**1.5*safexp(-4*pi**2/(Sm_b2+B)*x*x) \
  +Sm_a3*(4*pi/(Sm_b3+B))**1.5*safexp(-4*pi**2/(Sm_b3+B)*x*x) \
  +Sm_a4*(4*pi/(Sm_b4+B))**1.5*safexp(-4*pi**2/(Sm_b4+B)*x*x) \
  +Sm_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=59
Sm_plus_3_c = 0.954586;
Sm_plus_3_a1 = 23.150400;Sm_plus_3_a2 = 20.259899;Sm_plus_3_a3 = 11.920200;Sm_plus_3_a4 = 2.714880;
Sm_plus_3_b1 = 2.316410;Sm_plus_3_b2 = 0.174081;Sm_plus_3_b3 = 12.157100;Sm_plus_3_b4 = 24.824200;
Sm_plus_3_ff(x,B) = \
  +Sm_plus_3_a1*(4*pi/(Sm_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Sm_plus_3_b1+B)*x*x) \
  +Sm_plus_3_a2*(4*pi/(Sm_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Sm_plus_3_b2+B)*x*x) \
  +Sm_plus_3_a3*(4*pi/(Sm_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Sm_plus_3_b3+B)*x*x) \
  +Sm_plus_3_a4*(4*pi/(Sm_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Sm_plus_3_b4+B)*x*x) \
  +Sm_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=63
Eu_c = 2.574500;
Eu_a1 = 24.627399;Eu_a2 = 19.088600;Eu_a3 = 13.760300;Eu_a4 = 2.922700;
Eu_b1 = 2.387900;Eu_b2 = 0.194200;Eu_b3 = 13.754600;Eu_b4 = 123.173996;
Eu_ff(x,B) = \
  +Eu_a1*(4*pi/(Eu_b1+B))**1.5*safexp(-4*pi**2/(Eu_b1+B)*x*x) \
  +Eu_a2*(4*pi/(Eu_b2+B))**1.5*safexp(-4*pi**2/(Eu_b2+B)*x*x) \
  +Eu_a3*(4*pi/(Eu_b3+B))**1.5*safexp(-4*pi**2/(Eu_b3+B)*x*x) \
  +Eu_a4*(4*pi/(Eu_b4+B))**1.5*safexp(-4*pi**2/(Eu_b4+B)*x*x) \
  +Eu_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=61
Eu_plus_2_c = 1.363890;
Eu_plus_2_a1 = 24.006300;Eu_plus_2_a2 = 19.950399;Eu_plus_2_a3 = 11.803400;Eu_plus_2_a4 = 3.872430;
Eu_plus_2_b1 = 2.277830;Eu_plus_2_b2 = 0.173530;Eu_plus_2_b3 = 11.609600;Eu_plus_2_b4 = 26.515600;
Eu_plus_2_ff(x,B) = \
  +Eu_plus_2_a1*(4*pi/(Eu_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Eu_plus_2_b1+B)*x*x) \
  +Eu_plus_2_a2*(4*pi/(Eu_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Eu_plus_2_b2+B)*x*x) \
  +Eu_plus_2_a3*(4*pi/(Eu_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Eu_plus_2_b3+B)*x*x) \
  +Eu_plus_2_a4*(4*pi/(Eu_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Eu_plus_2_b4+B)*x*x) \
  +Eu_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=60
Eu_plus_3_c = 0.759344;
Eu_plus_3_a1 = 23.749699;Eu_plus_3_a2 = 20.374500;Eu_plus_3_a3 = 11.850900;Eu_plus_3_a4 = 3.265030;
Eu_plus_3_b1 = 2.222580;Eu_plus_3_b2 = 0.163940;Eu_plus_3_b3 = 11.311000;Eu_plus_3_b4 = 22.996599;
Eu_plus_3_ff(x,B) = \
  +Eu_plus_3_a1*(4*pi/(Eu_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Eu_plus_3_b1+B)*x*x) \
  +Eu_plus_3_a2*(4*pi/(Eu_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Eu_plus_3_b2+B)*x*x) \
  +Eu_plus_3_a3*(4*pi/(Eu_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Eu_plus_3_b3+B)*x*x) \
  +Eu_plus_3_a4*(4*pi/(Eu_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Eu_plus_3_b4+B)*x*x) \
  +Eu_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=64
Gd_c = 2.419600;
Gd_a1 = 25.070900;Gd_a2 = 19.079800;Gd_a3 = 13.851800;Gd_a4 = 3.545450;
Gd_b1 = 2.253410;Gd_b2 = 0.181951;Gd_b3 = 12.933100;Gd_b4 = 101.397995;
Gd_ff(x,B) = \
  +Gd_a1*(4*pi/(Gd_b1+B))**1.5*safexp(-4*pi**2/(Gd_b1+B)*x*x) \
  +Gd_a2*(4*pi/(Gd_b2+B))**1.5*safexp(-4*pi**2/(Gd_b2+B)*x*x) \
  +Gd_a3*(4*pi/(Gd_b3+B))**1.5*safexp(-4*pi**2/(Gd_b3+B)*x*x) \
  +Gd_a4*(4*pi/(Gd_b4+B))**1.5*safexp(-4*pi**2/(Gd_b4+B)*x*x) \
  +Gd_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=61
Gd_plus_3_c = 0.645089;
Gd_plus_3_a1 = 24.346600;Gd_plus_3_a2 = 20.420799;Gd_plus_3_a3 = 11.870800;Gd_plus_3_a4 = 3.714900;
Gd_plus_3_b1 = 2.135530;Gd_plus_3_b2 = 0.155525;Gd_plus_3_b3 = 10.578199;Gd_plus_3_b4 = 21.702900;
Gd_plus_3_ff(x,B) = \
  +Gd_plus_3_a1*(4*pi/(Gd_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Gd_plus_3_b1+B)*x*x) \
  +Gd_plus_3_a2*(4*pi/(Gd_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Gd_plus_3_b2+B)*x*x) \
  +Gd_plus_3_a3*(4*pi/(Gd_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Gd_plus_3_b3+B)*x*x) \
  +Gd_plus_3_a4*(4*pi/(Gd_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Gd_plus_3_b4+B)*x*x) \
  +Gd_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=65
Tb_c = 3.582240;
Tb_a1 = 25.897600;Tb_a2 = 18.218500;Tb_a3 = 14.316700;Tb_a4 = 2.953540;
Tb_b1 = 2.242560;Tb_b2 = 0.196143;Tb_b3 = 12.664800;Tb_b4 = 115.362000;
Tb_ff(x,B) = \
  +Tb_a1*(4*pi/(Tb_b1+B))**1.5*safexp(-4*pi**2/(Tb_b1+B)*x*x) \
  +Tb_a2*(4*pi/(Tb_b2+B))**1.5*safexp(-4*pi**2/(Tb_b2+B)*x*x) \
  +Tb_a3*(4*pi/(Tb_b3+B))**1.5*safexp(-4*pi**2/(Tb_b3+B)*x*x) \
  +Tb_a4*(4*pi/(Tb_b4+B))**1.5*safexp(-4*pi**2/(Tb_b4+B)*x*x) \
  +Tb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=62
Tb_plus_3_c = 0.691967;
Tb_plus_3_a1 = 24.955900;Tb_plus_3_a2 = 20.327099;Tb_plus_3_a3 = 12.247100;Tb_plus_3_a4 = 3.773000;
Tb_plus_3_b1 = 2.056010;Tb_plus_3_b2 = 0.149525;Tb_plus_3_b3 = 10.049900;Tb_plus_3_b4 = 21.277300;
Tb_plus_3_ff(x,B) = \
  +Tb_plus_3_a1*(4*pi/(Tb_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Tb_plus_3_b1+B)*x*x) \
  +Tb_plus_3_a2*(4*pi/(Tb_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Tb_plus_3_b2+B)*x*x) \
  +Tb_plus_3_a3*(4*pi/(Tb_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Tb_plus_3_b3+B)*x*x) \
  +Tb_plus_3_a4*(4*pi/(Tb_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Tb_plus_3_b4+B)*x*x) \
  +Tb_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=66
Dy_c = 4.297280;
Dy_a1 = 26.507000;Dy_a2 = 17.638300;Dy_a3 = 14.559600;Dy_a4 = 2.965770;
Dy_b1 = 2.180200;Dy_b2 = 0.202172;Dy_b3 = 12.189899;Dy_b4 = 111.874001;
Dy_ff(x,B) = \
  +Dy_a1*(4*pi/(Dy_b1+B))**1.5*safexp(-4*pi**2/(Dy_b1+B)*x*x) \
  +Dy_a2*(4*pi/(Dy_b2+B))**1.5*safexp(-4*pi**2/(Dy_b2+B)*x*x) \
  +Dy_a3*(4*pi/(Dy_b3+B))**1.5*safexp(-4*pi**2/(Dy_b3+B)*x*x) \
  +Dy_a4*(4*pi/(Dy_b4+B))**1.5*safexp(-4*pi**2/(Dy_b4+B)*x*x) \
  +Dy_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=63
Dy_plus_3_c = 0.689690;
Dy_plus_3_a1 = 25.539499;Dy_plus_3_a2 = 20.286100;Dy_plus_3_a3 = 11.981200;Dy_plus_3_a4 = 4.500730;
Dy_plus_3_b1 = 1.980400;Dy_plus_3_b2 = 0.143384;Dy_plus_3_b3 = 9.349720;Dy_plus_3_b4 = 19.580999;
Dy_plus_3_ff(x,B) = \
  +Dy_plus_3_a1*(4*pi/(Dy_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Dy_plus_3_b1+B)*x*x) \
  +Dy_plus_3_a2*(4*pi/(Dy_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Dy_plus_3_b2+B)*x*x) \
  +Dy_plus_3_a3*(4*pi/(Dy_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Dy_plus_3_b3+B)*x*x) \
  +Dy_plus_3_a4*(4*pi/(Dy_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Dy_plus_3_b4+B)*x*x) \
  +Dy_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=67
Ho_c = 4.567960;
Ho_a1 = 26.904900;Ho_a2 = 17.293999;Ho_a3 = 14.558300;Ho_a4 = 3.638370;
Ho_b1 = 2.070510;Ho_b2 = 0.197940;Ho_b3 = 11.440700;Ho_b4 = 92.656601;
Ho_ff(x,B) = \
  +Ho_a1*(4*pi/(Ho_b1+B))**1.5*safexp(-4*pi**2/(Ho_b1+B)*x*x) \
  +Ho_a2*(4*pi/(Ho_b2+B))**1.5*safexp(-4*pi**2/(Ho_b2+B)*x*x) \
  +Ho_a3*(4*pi/(Ho_b3+B))**1.5*safexp(-4*pi**2/(Ho_b3+B)*x*x) \
  +Ho_a4*(4*pi/(Ho_b4+B))**1.5*safexp(-4*pi**2/(Ho_b4+B)*x*x) \
  +Ho_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=64
Ho_plus_3_c = 0.852795;
Ho_plus_3_a1 = 26.129601;Ho_plus_3_a2 = 20.099400;Ho_plus_3_a3 = 11.978800;Ho_plus_3_a4 = 4.936760;
Ho_plus_3_b1 = 1.910720;Ho_plus_3_b2 = 0.139358;Ho_plus_3_b3 = 8.800180;Ho_plus_3_b4 = 18.590799;
Ho_plus_3_ff(x,B) = \
  +Ho_plus_3_a1*(4*pi/(Ho_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ho_plus_3_b1+B)*x*x) \
  +Ho_plus_3_a2*(4*pi/(Ho_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ho_plus_3_b2+B)*x*x) \
  +Ho_plus_3_a3*(4*pi/(Ho_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ho_plus_3_b3+B)*x*x) \
  +Ho_plus_3_a4*(4*pi/(Ho_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ho_plus_3_b4+B)*x*x) \
  +Ho_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
Er_c = 5.920460;
Er_a1 = 27.656300;Er_a2 = 16.428499;Er_a3 = 14.977900;Er_a4 = 2.982330;
Er_b1 = 2.073560;Er_b2 = 0.223545;Er_b3 = 11.360400;Er_b4 = 105.703003;
Er_ff(x,B) = \
  +Er_a1*(4*pi/(Er_b1+B))**1.5*safexp(-4*pi**2/(Er_b1+B)*x*x) \
  +Er_a2*(4*pi/(Er_b2+B))**1.5*safexp(-4*pi**2/(Er_b2+B)*x*x) \
  +Er_a3*(4*pi/(Er_b3+B))**1.5*safexp(-4*pi**2/(Er_b3+B)*x*x) \
  +Er_a4*(4*pi/(Er_b4+B))**1.5*safexp(-4*pi**2/(Er_b4+B)*x*x) \
  +Er_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=65
Er_plus_3_c = 1.176130;
Er_plus_3_a1 = 26.722000;Er_plus_3_a2 = 19.774799;Er_plus_3_a3 = 12.150600;Er_plus_3_a4 = 5.173790;
Er_plus_3_b1 = 1.846590;Er_plus_3_b2 = 0.137290;Er_plus_3_b3 = 8.362249;Er_plus_3_b4 = 17.897400;
Er_plus_3_ff(x,B) = \
  +Er_plus_3_a1*(4*pi/(Er_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Er_plus_3_b1+B)*x*x) \
  +Er_plus_3_a2*(4*pi/(Er_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Er_plus_3_b2+B)*x*x) \
  +Er_plus_3_a3*(4*pi/(Er_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Er_plus_3_b3+B)*x*x) \
  +Er_plus_3_a4*(4*pi/(Er_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Er_plus_3_b4+B)*x*x) \
  +Er_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=69
Tm_c = 6.756210;
Tm_a1 = 28.181900;Tm_a2 = 15.885099;Tm_a3 = 15.154200;Tm_a4 = 2.987060;
Tm_b1 = 2.028590;Tm_b2 = 0.238849;Tm_b3 = 10.997499;Tm_b4 = 102.960999;
Tm_ff(x,B) = \
  +Tm_a1*(4*pi/(Tm_b1+B))**1.5*safexp(-4*pi**2/(Tm_b1+B)*x*x) \
  +Tm_a2*(4*pi/(Tm_b2+B))**1.5*safexp(-4*pi**2/(Tm_b2+B)*x*x) \
  +Tm_a3*(4*pi/(Tm_b3+B))**1.5*safexp(-4*pi**2/(Tm_b3+B)*x*x) \
  +Tm_a4*(4*pi/(Tm_b4+B))**1.5*safexp(-4*pi**2/(Tm_b4+B)*x*x) \
  +Tm_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=66
Tm_plus_3_c = 1.639290;
Tm_plus_3_a1 = 27.308300;Tm_plus_3_a2 = 19.332001;Tm_plus_3_a3 = 12.333900;Tm_plus_3_a4 = 5.383480;
Tm_plus_3_b1 = 1.787110;Tm_plus_3_b2 = 0.136974;Tm_plus_3_b3 = 7.967780;Tm_plus_3_b4 = 17.292200;
Tm_plus_3_ff(x,B) = \
  +Tm_plus_3_a1*(4*pi/(Tm_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Tm_plus_3_b1+B)*x*x) \
  +Tm_plus_3_a2*(4*pi/(Tm_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Tm_plus_3_b2+B)*x*x) \
  +Tm_plus_3_a3*(4*pi/(Tm_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Tm_plus_3_b3+B)*x*x) \
  +Tm_plus_3_a4*(4*pi/(Tm_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Tm_plus_3_b4+B)*x*x) \
  +Tm_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=70
Yb_c = 7.566720;
Yb_a1 = 28.664101;Yb_a2 = 15.434500;Yb_a3 = 15.308700;Yb_a4 = 2.989630;
Yb_b1 = 1.988900;Yb_b2 = 0.257119;Yb_b3 = 10.664700;Yb_b4 = 100.417000;
Yb_ff(x,B) = \
  +Yb_a1*(4*pi/(Yb_b1+B))**1.5*safexp(-4*pi**2/(Yb_b1+B)*x*x) \
  +Yb_a2*(4*pi/(Yb_b2+B))**1.5*safexp(-4*pi**2/(Yb_b2+B)*x*x) \
  +Yb_a3*(4*pi/(Yb_b3+B))**1.5*safexp(-4*pi**2/(Yb_b3+B)*x*x) \
  +Yb_a4*(4*pi/(Yb_b4+B))**1.5*safexp(-4*pi**2/(Yb_b4+B)*x*x) \
  +Yb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
Yb_plus_2_c = 3.709830;
Yb_plus_2_a1 = 28.120899;Yb_plus_2_a2 = 17.681700;Yb_plus_2_a3 = 13.333500;Yb_plus_2_a4 = 5.146570;
Yb_plus_2_b1 = 1.785030;Yb_plus_2_b2 = 0.159970;Yb_plus_2_b3 = 8.183040;Yb_plus_2_b4 = 20.389999;
Yb_plus_2_ff(x,B) = \
  +Yb_plus_2_a1*(4*pi/(Yb_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Yb_plus_2_b1+B)*x*x) \
  +Yb_plus_2_a2*(4*pi/(Yb_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Yb_plus_2_b2+B)*x*x) \
  +Yb_plus_2_a3*(4*pi/(Yb_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Yb_plus_2_b3+B)*x*x) \
  +Yb_plus_2_a4*(4*pi/(Yb_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Yb_plus_2_b4+B)*x*x) \
  +Yb_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=67
Yb_plus_3_c = 2.260010;
Yb_plus_3_a1 = 27.891700;Yb_plus_3_a2 = 18.761400;Yb_plus_3_a3 = 12.607200;Yb_plus_3_a4 = 5.476470;
Yb_plus_3_b1 = 1.732720;Yb_plus_3_b2 = 0.138790;Yb_plus_3_b3 = 7.644120;Yb_plus_3_b4 = 16.815300;
Yb_plus_3_ff(x,B) = \
  +Yb_plus_3_a1*(4*pi/(Yb_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Yb_plus_3_b1+B)*x*x) \
  +Yb_plus_3_a2*(4*pi/(Yb_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Yb_plus_3_b2+B)*x*x) \
  +Yb_plus_3_a3*(4*pi/(Yb_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Yb_plus_3_b3+B)*x*x) \
  +Yb_plus_3_a4*(4*pi/(Yb_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Yb_plus_3_b4+B)*x*x) \
  +Yb_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=71
Lu_c = 7.976280;
Lu_a1 = 28.947599;Lu_a2 = 15.220800;Lu_a3 = 15.100000;Lu_a4 = 3.716010;
Lu_b1 = 1.901820;Lu_b2 = 9.985189;Lu_b3 = 0.261033;Lu_b4 = 84.329803;
Lu_ff(x,B) = \
  +Lu_a1*(4*pi/(Lu_b1+B))**1.5*safexp(-4*pi**2/(Lu_b1+B)*x*x) \
  +Lu_a2*(4*pi/(Lu_b2+B))**1.5*safexp(-4*pi**2/(Lu_b2+B)*x*x) \
  +Lu_a3*(4*pi/(Lu_b3+B))**1.5*safexp(-4*pi**2/(Lu_b3+B)*x*x) \
  +Lu_a4*(4*pi/(Lu_b4+B))**1.5*safexp(-4*pi**2/(Lu_b4+B)*x*x) \
  +Lu_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
Lu_plus_3_c = 2.975730;
Lu_plus_3_a1 = 28.462799;Lu_plus_3_a2 = 18.121000;Lu_plus_3_a3 = 12.842899;Lu_plus_3_a4 = 5.594150;
Lu_plus_3_b1 = 1.682160;Lu_plus_3_b2 = 0.142292;Lu_plus_3_b3 = 7.337270;Lu_plus_3_b4 = 16.353500;
Lu_plus_3_ff(x,B) = \
  +Lu_plus_3_a1*(4*pi/(Lu_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Lu_plus_3_b1+B)*x*x) \
  +Lu_plus_3_a2*(4*pi/(Lu_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Lu_plus_3_b2+B)*x*x) \
  +Lu_plus_3_a3*(4*pi/(Lu_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Lu_plus_3_b3+B)*x*x) \
  +Lu_plus_3_a4*(4*pi/(Lu_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Lu_plus_3_b4+B)*x*x) \
  +Lu_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=72
Hf_c = 8.581540;
Hf_a1 = 29.143999;Hf_a2 = 15.172600;Hf_a3 = 14.758600;Hf_a4 = 4.300130;
Hf_b1 = 1.832620;Hf_b2 = 9.599899;Hf_b3 = 0.275116;Hf_b4 = 72.028999;
Hf_ff(x,B) = \
  +Hf_a1*(4*pi/(Hf_b1+B))**1.5*safexp(-4*pi**2/(Hf_b1+B)*x*x) \
  +Hf_a2*(4*pi/(Hf_b2+B))**1.5*safexp(-4*pi**2/(Hf_b2+B)*x*x) \
  +Hf_a3*(4*pi/(Hf_b3+B))**1.5*safexp(-4*pi**2/(Hf_b3+B)*x*x) \
  +Hf_a4*(4*pi/(Hf_b4+B))**1.5*safexp(-4*pi**2/(Hf_b4+B)*x*x) \
  +Hf_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
Hf_plus_4_c = 2.396990;
Hf_plus_4_a1 = 28.813099;Hf_plus_4_a2 = 18.460100;Hf_plus_4_a3 = 12.728500;Hf_plus_4_a4 = 5.599270;
Hf_plus_4_b1 = 1.591360;Hf_plus_4_b2 = 0.128903;Hf_plus_4_b3 = 6.762320;Hf_plus_4_b4 = 14.036600;
Hf_plus_4_ff(x,B) = \
  +Hf_plus_4_a1*(4*pi/(Hf_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Hf_plus_4_b1+B)*x*x) \
  +Hf_plus_4_a2*(4*pi/(Hf_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Hf_plus_4_b2+B)*x*x) \
  +Hf_plus_4_a3*(4*pi/(Hf_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Hf_plus_4_b3+B)*x*x) \
  +Hf_plus_4_a4*(4*pi/(Hf_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Hf_plus_4_b4+B)*x*x) \
  +Hf_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=73
Ta_c = 9.243540;
Ta_a1 = 29.202400;Ta_a2 = 15.229300;Ta_a3 = 14.513500;Ta_a4 = 4.764920;
Ta_b1 = 1.773330;Ta_b2 = 9.370460;Ta_b3 = 0.295977;Ta_b4 = 63.364399;
Ta_ff(x,B) = \
  +Ta_a1*(4*pi/(Ta_b1+B))**1.5*safexp(-4*pi**2/(Ta_b1+B)*x*x) \
  +Ta_a2*(4*pi/(Ta_b2+B))**1.5*safexp(-4*pi**2/(Ta_b2+B)*x*x) \
  +Ta_a3*(4*pi/(Ta_b3+B))**1.5*safexp(-4*pi**2/(Ta_b3+B)*x*x) \
  +Ta_a4*(4*pi/(Ta_b4+B))**1.5*safexp(-4*pi**2/(Ta_b4+B)*x*x) \
  +Ta_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
Ta_plus_5_c = 1.785550;
Ta_plus_5_a1 = 29.158699;Ta_plus_5_a2 = 18.840700;Ta_plus_5_a3 = 12.826799;Ta_plus_5_a4 = 5.386950;
Ta_plus_5_b1 = 1.507110;Ta_plus_5_b2 = 0.116741;Ta_plus_5_b3 = 6.315240;Ta_plus_5_b4 = 12.424400;
Ta_plus_5_ff(x,B) = \
  +Ta_plus_5_a1*(4*pi/(Ta_plus_5_b1+B))**1.5*safexp(-4*pi**2/(Ta_plus_5_b1+B)*x*x) \
  +Ta_plus_5_a2*(4*pi/(Ta_plus_5_b2+B))**1.5*safexp(-4*pi**2/(Ta_plus_5_b2+B)*x*x) \
  +Ta_plus_5_a3*(4*pi/(Ta_plus_5_b3+B))**1.5*safexp(-4*pi**2/(Ta_plus_5_b3+B)*x*x) \
  +Ta_plus_5_a4*(4*pi/(Ta_plus_5_b4+B))**1.5*safexp(-4*pi**2/(Ta_plus_5_b4+B)*x*x) \
  +Ta_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=74
W_c = 9.887500;
W_a1 = 29.081800;W_a2 = 15.430000;W_a3 = 14.432700;W_a4 = 5.119820;
W_b1 = 1.720290;W_b2 = 9.225900;W_b3 = 0.321703;W_b4 = 57.056000;
W_ff(x,B) = \
  +W_a1*(4*pi/(W_b1+B))**1.5*safexp(-4*pi**2/(W_b1+B)*x*x) \
  +W_a2*(4*pi/(W_b2+B))**1.5*safexp(-4*pi**2/(W_b2+B)*x*x) \
  +W_a3*(4*pi/(W_b3+B))**1.5*safexp(-4*pi**2/(W_b3+B)*x*x) \
  +W_a4*(4*pi/(W_b4+B))**1.5*safexp(-4*pi**2/(W_b4+B)*x*x) \
  +W_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=68
W_plus_6_c = 1.010740;
W_plus_6_a1 = 29.493599;W_plus_6_a2 = 19.376301;W_plus_6_a3 = 13.054399;W_plus_6_a4 = 5.064120;
W_plus_6_b1 = 1.427550;W_plus_6_b2 = 0.104621;W_plus_6_b3 = 5.936670;W_plus_6_b4 = 11.197200;
W_plus_6_ff(x,B) = \
  +W_plus_6_a1*(4*pi/(W_plus_6_b1+B))**1.5*safexp(-4*pi**2/(W_plus_6_b1+B)*x*x) \
  +W_plus_6_a2*(4*pi/(W_plus_6_b2+B))**1.5*safexp(-4*pi**2/(W_plus_6_b2+B)*x*x) \
  +W_plus_6_a3*(4*pi/(W_plus_6_b3+B))**1.5*safexp(-4*pi**2/(W_plus_6_b3+B)*x*x) \
  +W_plus_6_a4*(4*pi/(W_plus_6_b4+B))**1.5*safexp(-4*pi**2/(W_plus_6_b4+B)*x*x) \
  +W_plus_6_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=75
Re_c = 10.472000;
Re_a1 = 28.762100;Re_a2 = 15.718900;Re_a3 = 14.556400;Re_a4 = 5.441740;
Re_b1 = 1.671910;Re_b2 = 9.092270;Re_b3 = 0.350500;Re_b4 = 52.086098;
Re_ff(x,B) = \
  +Re_a1*(4*pi/(Re_b1+B))**1.5*safexp(-4*pi**2/(Re_b1+B)*x*x) \
  +Re_a2*(4*pi/(Re_b2+B))**1.5*safexp(-4*pi**2/(Re_b2+B)*x*x) \
  +Re_a3*(4*pi/(Re_b3+B))**1.5*safexp(-4*pi**2/(Re_b3+B)*x*x) \
  +Re_a4*(4*pi/(Re_b4+B))**1.5*safexp(-4*pi**2/(Re_b4+B)*x*x) \
  +Re_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=76
Os_c = 11.000500;
Os_a1 = 28.189400;Os_a2 = 16.154999;Os_a3 = 14.930500;Os_a4 = 5.675890;
Os_b1 = 1.629030;Os_b2 = 8.979480;Os_b3 = 0.382661;Os_b4 = 48.164700;
Os_ff(x,B) = \
  +Os_a1*(4*pi/(Os_b1+B))**1.5*safexp(-4*pi**2/(Os_b1+B)*x*x) \
  +Os_a2*(4*pi/(Os_b2+B))**1.5*safexp(-4*pi**2/(Os_b2+B)*x*x) \
  +Os_a3*(4*pi/(Os_b3+B))**1.5*safexp(-4*pi**2/(Os_b3+B)*x*x) \
  +Os_a4*(4*pi/(Os_b4+B))**1.5*safexp(-4*pi**2/(Os_b4+B)*x*x) \
  +Os_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=72
Os_plus_4_c = 6.498040;
Os_plus_4_a1 = 30.418999;Os_plus_4_a2 = 15.263700;Os_plus_4_a3 = 14.745800;Os_plus_4_a4 = 5.067950;
Os_plus_4_b1 = 1.371130;Os_plus_4_b2 = 6.847060;Os_plus_4_b3 = 0.165191;Os_plus_4_b4 = 18.003000;
Os_plus_4_ff(x,B) = \
  +Os_plus_4_a1*(4*pi/(Os_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Os_plus_4_b1+B)*x*x) \
  +Os_plus_4_a2*(4*pi/(Os_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Os_plus_4_b2+B)*x*x) \
  +Os_plus_4_a3*(4*pi/(Os_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Os_plus_4_b3+B)*x*x) \
  +Os_plus_4_a4*(4*pi/(Os_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Os_plus_4_b4+B)*x*x) \
  +Os_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=77
Ir_c = 11.472200;
Ir_a1 = 27.304899;Ir_a2 = 16.729599;Ir_a3 = 15.611500;Ir_a4 = 5.833770;
Ir_b1 = 1.592790;Ir_b2 = 8.865530;Ir_b3 = 0.417916;Ir_b4 = 45.001099;
Ir_ff(x,B) = \
  +Ir_a1*(4*pi/(Ir_b1+B))**1.5*safexp(-4*pi**2/(Ir_b1+B)*x*x) \
  +Ir_a2*(4*pi/(Ir_b2+B))**1.5*safexp(-4*pi**2/(Ir_b2+B)*x*x) \
  +Ir_a3*(4*pi/(Ir_b3+B))**1.5*safexp(-4*pi**2/(Ir_b3+B)*x*x) \
  +Ir_a4*(4*pi/(Ir_b4+B))**1.5*safexp(-4*pi**2/(Ir_b4+B)*x*x) \
  +Ir_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=74
Ir_plus_3_c = 8.279030;
Ir_plus_3_a1 = 30.415600;Ir_plus_3_a2 = 15.862000;Ir_plus_3_a3 = 13.614500;Ir_plus_3_a4 = 5.820080;
Ir_plus_3_b1 = 1.343230;Ir_plus_3_b2 = 7.109090;Ir_plus_3_b3 = 0.204633;Ir_plus_3_b4 = 20.325399;
Ir_plus_3_ff(x,B) = \
  +Ir_plus_3_a1*(4*pi/(Ir_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ir_plus_3_b1+B)*x*x) \
  +Ir_plus_3_a2*(4*pi/(Ir_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ir_plus_3_b2+B)*x*x) \
  +Ir_plus_3_a3*(4*pi/(Ir_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ir_plus_3_b3+B)*x*x) \
  +Ir_plus_3_a4*(4*pi/(Ir_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ir_plus_3_b4+B)*x*x) \
  +Ir_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=73
Ir_plus_4_c = 6.968240;
Ir_plus_4_a1 = 30.705799;Ir_plus_4_a2 = 15.551200;Ir_plus_4_a3 = 14.232600;Ir_plus_4_a4 = 5.536720;
Ir_plus_4_b1 = 1.309230;Ir_plus_4_b2 = 6.719830;Ir_plus_4_b3 = 0.167252;Ir_plus_4_b4 = 17.491100;
Ir_plus_4_ff(x,B) = \
  +Ir_plus_4_a1*(4*pi/(Ir_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Ir_plus_4_b1+B)*x*x) \
  +Ir_plus_4_a2*(4*pi/(Ir_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Ir_plus_4_b2+B)*x*x) \
  +Ir_plus_4_a3*(4*pi/(Ir_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Ir_plus_4_b3+B)*x*x) \
  +Ir_plus_4_a4*(4*pi/(Ir_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Ir_plus_4_b4+B)*x*x) \
  +Ir_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Pt_c = 11.688300;
Pt_a1 = 27.005899;Pt_a2 = 17.763901;Pt_a3 = 15.713100;Pt_a4 = 5.783700;
Pt_b1 = 1.512930;Pt_b2 = 8.811740;Pt_b3 = 0.424593;Pt_b4 = 38.610298;
Pt_ff(x,B) = \
  +Pt_a1*(4*pi/(Pt_b1+B))**1.5*safexp(-4*pi**2/(Pt_b1+B)*x*x) \
  +Pt_a2*(4*pi/(Pt_b2+B))**1.5*safexp(-4*pi**2/(Pt_b2+B)*x*x) \
  +Pt_a3*(4*pi/(Pt_b3+B))**1.5*safexp(-4*pi**2/(Pt_b3+B)*x*x) \
  +Pt_a4*(4*pi/(Pt_b4+B))**1.5*safexp(-4*pi**2/(Pt_b4+B)*x*x) \
  +Pt_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=76
Pt_plus_2_c = 9.853290;
Pt_plus_2_a1 = 29.842899;Pt_plus_2_a2 = 16.722401;Pt_plus_2_a3 = 13.215300;Pt_plus_2_a4 = 6.352340;
Pt_plus_2_b1 = 1.329270;Pt_plus_2_b2 = 7.389790;Pt_plus_2_b3 = 0.263297;Pt_plus_2_b4 = 22.942600;
Pt_plus_2_ff(x,B) = \
  +Pt_plus_2_a1*(4*pi/(Pt_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Pt_plus_2_b1+B)*x*x) \
  +Pt_plus_2_a2*(4*pi/(Pt_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Pt_plus_2_b2+B)*x*x) \
  +Pt_plus_2_a3*(4*pi/(Pt_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Pt_plus_2_b3+B)*x*x) \
  +Pt_plus_2_a4*(4*pi/(Pt_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Pt_plus_2_b4+B)*x*x) \
  +Pt_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=74
Pt_plus_4_c = 7.395340;
Pt_plus_4_a1 = 30.961201;Pt_plus_4_a2 = 15.982900;Pt_plus_4_a3 = 13.734800;Pt_plus_4_a4 = 5.920340;
Pt_plus_4_b1 = 1.248130;Pt_plus_4_b2 = 6.608340;Pt_plus_4_b3 = 0.168640;Pt_plus_4_b4 = 16.939199;
Pt_plus_4_ff(x,B) = \
  +Pt_plus_4_a1*(4*pi/(Pt_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Pt_plus_4_b1+B)*x*x) \
  +Pt_plus_4_a2*(4*pi/(Pt_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Pt_plus_4_b2+B)*x*x) \
  +Pt_plus_4_a3*(4*pi/(Pt_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Pt_plus_4_b3+B)*x*x) \
  +Pt_plus_4_a4*(4*pi/(Pt_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Pt_plus_4_b4+B)*x*x) \
  +Pt_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=79
Au_c = 12.065800;
Au_a1 = 16.881901;Au_a2 = 18.591299;Au_a3 = 25.558201;Au_a4 = 5.860000;
Au_b1 = 0.461100;Au_b2 = 8.621600;Au_b3 = 1.482600;Au_b4 = 36.395599;
Au_ff(x,B) = \
  +Au_a1*(4*pi/(Au_b1+B))**1.5*safexp(-4*pi**2/(Au_b1+B)*x*x) \
  +Au_a2*(4*pi/(Au_b2+B))**1.5*safexp(-4*pi**2/(Au_b2+B)*x*x) \
  +Au_a3*(4*pi/(Au_b3+B))**1.5*safexp(-4*pi**2/(Au_b3+B)*x*x) \
  +Au_a4*(4*pi/(Au_b4+B))**1.5*safexp(-4*pi**2/(Au_b4+B)*x*x) \
  +Au_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Au_plus_1_c = 11.229900;
Au_plus_1_a1 = 28.010899;Au_plus_1_a2 = 17.820400;Au_plus_1_a3 = 14.335899;Au_plus_1_a4 = 6.580770;
Au_plus_1_b1 = 1.353210;Au_plus_1_b2 = 7.739500;Au_plus_1_b3 = 0.356752;Au_plus_1_b4 = 26.404301;
Au_plus_1_ff(x,B) = \
  +Au_plus_1_a1*(4*pi/(Au_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Au_plus_1_b1+B)*x*x) \
  +Au_plus_1_a2*(4*pi/(Au_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Au_plus_1_b2+B)*x*x) \
  +Au_plus_1_a3*(4*pi/(Au_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Au_plus_1_b3+B)*x*x) \
  +Au_plus_1_a4*(4*pi/(Au_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Au_plus_1_b4+B)*x*x) \
  +Au_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=76
Au_plus_3_c = 9.096800;
Au_plus_3_a1 = 30.688599;Au_plus_3_a2 = 16.902901;Au_plus_3_a3 = 12.780100;Au_plus_3_a4 = 6.523540;
Au_plus_3_b1 = 1.219900;Au_plus_3_b2 = 6.828720;Au_plus_3_b3 = 0.212867;Au_plus_3_b4 = 18.659000;
Au_plus_3_ff(x,B) = \
  +Au_plus_3_a1*(4*pi/(Au_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Au_plus_3_b1+B)*x*x) \
  +Au_plus_3_a2*(4*pi/(Au_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Au_plus_3_b2+B)*x*x) \
  +Au_plus_3_a3*(4*pi/(Au_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Au_plus_3_b3+B)*x*x) \
  +Au_plus_3_a4*(4*pi/(Au_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Au_plus_3_b4+B)*x*x) \
  +Au_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=80
Hg_c = 12.608900;
Hg_a1 = 20.680901;Hg_a2 = 19.041700;Hg_a3 = 21.657499;Hg_a4 = 5.967600;
Hg_b1 = 0.545000;Hg_b2 = 8.448400;Hg_b3 = 1.572900;Hg_b4 = 38.324600;
Hg_ff(x,B) = \
  +Hg_a1*(4*pi/(Hg_b1+B))**1.5*safexp(-4*pi**2/(Hg_b1+B)*x*x) \
  +Hg_a2*(4*pi/(Hg_b2+B))**1.5*safexp(-4*pi**2/(Hg_b2+B)*x*x) \
  +Hg_a3*(4*pi/(Hg_b3+B))**1.5*safexp(-4*pi**2/(Hg_b3+B)*x*x) \
  +Hg_a4*(4*pi/(Hg_b4+B))**1.5*safexp(-4*pi**2/(Hg_b4+B)*x*x) \
  +Hg_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=79
Hg_plus_1_c = 12.020500;
Hg_plus_1_a1 = 25.085300;Hg_plus_1_a2 = 18.497299;Hg_plus_1_a3 = 16.888300;Hg_plus_1_a4 = 6.482160;
Hg_plus_1_b1 = 1.395070;Hg_plus_1_b2 = 7.651050;Hg_plus_1_b3 = 0.443378;Hg_plus_1_b4 = 28.226200;
Hg_plus_1_ff(x,B) = \
  +Hg_plus_1_a1*(4*pi/(Hg_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Hg_plus_1_b1+B)*x*x) \
  +Hg_plus_1_a2*(4*pi/(Hg_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Hg_plus_1_b2+B)*x*x) \
  +Hg_plus_1_a3*(4*pi/(Hg_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Hg_plus_1_b3+B)*x*x) \
  +Hg_plus_1_a4*(4*pi/(Hg_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Hg_plus_1_b4+B)*x*x) \
  +Hg_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Hg_plus_2_c = 10.626800;
Hg_plus_2_a1 = 29.564100;Hg_plus_2_a2 = 18.059999;Hg_plus_2_a3 = 12.837400;Hg_plus_2_a4 = 6.899120;
Hg_plus_2_b1 = 1.211520;Hg_plus_2_b2 = 7.056390;Hg_plus_2_b3 = 0.284738;Hg_plus_2_b4 = 20.748199;
Hg_plus_2_ff(x,B) = \
  +Hg_plus_2_a1*(4*pi/(Hg_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Hg_plus_2_b1+B)*x*x) \
  +Hg_plus_2_a2*(4*pi/(Hg_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Hg_plus_2_b2+B)*x*x) \
  +Hg_plus_2_a3*(4*pi/(Hg_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Hg_plus_2_b3+B)*x*x) \
  +Hg_plus_2_a4*(4*pi/(Hg_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Hg_plus_2_b4+B)*x*x) \
  +Hg_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=81
Tl_c = 13.174600;
Tl_a1 = 27.544600;Tl_a2 = 19.158400;Tl_a3 = 15.538000;Tl_a4 = 5.525930;
Tl_b1 = 0.655150;Tl_b2 = 8.707510;Tl_b3 = 1.963470;Tl_b4 = 45.814899;
Tl_ff(x,B) = \
  +Tl_a1*(4*pi/(Tl_b1+B))**1.5*safexp(-4*pi**2/(Tl_b1+B)*x*x) \
  +Tl_a2*(4*pi/(Tl_b2+B))**1.5*safexp(-4*pi**2/(Tl_b2+B)*x*x) \
  +Tl_a3*(4*pi/(Tl_b3+B))**1.5*safexp(-4*pi**2/(Tl_b3+B)*x*x) \
  +Tl_a4*(4*pi/(Tl_b4+B))**1.5*safexp(-4*pi**2/(Tl_b4+B)*x*x) \
  +Tl_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=80
Tl_plus_1_c = 12.525800;
Tl_plus_1_a1 = 21.398500;Tl_plus_1_a2 = 20.472300;Tl_plus_1_a3 = 18.747799;Tl_plus_1_a4 = 6.828470;
Tl_plus_1_b1 = 1.471100;Tl_plus_1_b2 = 0.517394;Tl_plus_1_b3 = 7.434630;Tl_plus_1_b4 = 28.848200;
Tl_plus_1_ff(x,B) = \
  +Tl_plus_1_a1*(4*pi/(Tl_plus_1_b1+B))**1.5*safexp(-4*pi**2/(Tl_plus_1_b1+B)*x*x) \
  +Tl_plus_1_a2*(4*pi/(Tl_plus_1_b2+B))**1.5*safexp(-4*pi**2/(Tl_plus_1_b2+B)*x*x) \
  +Tl_plus_1_a3*(4*pi/(Tl_plus_1_b3+B))**1.5*safexp(-4*pi**2/(Tl_plus_1_b3+B)*x*x) \
  +Tl_plus_1_a4*(4*pi/(Tl_plus_1_b4+B))**1.5*safexp(-4*pi**2/(Tl_plus_1_b4+B)*x*x) \
  +Tl_plus_1_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Tl_plus_3_c = 9.802700;
Tl_plus_3_a1 = 30.869499;Tl_plus_3_a2 = 18.384100;Tl_plus_3_a3 = 11.932800;Tl_plus_3_a4 = 7.005740;
Tl_plus_3_b1 = 1.100800;Tl_plus_3_b2 = 6.538520;Tl_plus_3_b3 = 0.219074;Tl_plus_3_b4 = 17.211399;
Tl_plus_3_ff(x,B) = \
  +Tl_plus_3_a1*(4*pi/(Tl_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Tl_plus_3_b1+B)*x*x) \
  +Tl_plus_3_a2*(4*pi/(Tl_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Tl_plus_3_b2+B)*x*x) \
  +Tl_plus_3_a3*(4*pi/(Tl_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Tl_plus_3_b3+B)*x*x) \
  +Tl_plus_3_a4*(4*pi/(Tl_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Tl_plus_3_b4+B)*x*x) \
  +Tl_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=82
Pb_c = 13.411800;
Pb_a1 = 31.061699;Pb_a2 = 13.063700;Pb_a3 = 18.441999;Pb_a4 = 5.969600;
Pb_b1 = 0.690200;Pb_b2 = 2.357600;Pb_b3 = 8.618000;Pb_b4 = 47.257900;
Pb_ff(x,B) = \
  +Pb_a1*(4*pi/(Pb_b1+B))**1.5*safexp(-4*pi**2/(Pb_b1+B)*x*x) \
  +Pb_a2*(4*pi/(Pb_b2+B))**1.5*safexp(-4*pi**2/(Pb_b2+B)*x*x) \
  +Pb_a3*(4*pi/(Pb_b3+B))**1.5*safexp(-4*pi**2/(Pb_b3+B)*x*x) \
  +Pb_a4*(4*pi/(Pb_b4+B))**1.5*safexp(-4*pi**2/(Pb_b4+B)*x*x) \
  +Pb_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=80
Pb_plus_2_c = 12.473400;
Pb_plus_2_a1 = 21.788601;Pb_plus_2_a2 = 19.568199;Pb_plus_2_a3 = 19.140600;Pb_plus_2_a4 = 7.011070;
Pb_plus_2_b1 = 1.336600;Pb_plus_2_b2 = 0.488383;Pb_plus_2_b3 = 6.772700;Pb_plus_2_b4 = 23.813200;
Pb_plus_2_ff(x,B) = \
  +Pb_plus_2_a1*(4*pi/(Pb_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Pb_plus_2_b1+B)*x*x) \
  +Pb_plus_2_a2*(4*pi/(Pb_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Pb_plus_2_b2+B)*x*x) \
  +Pb_plus_2_a3*(4*pi/(Pb_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Pb_plus_2_b3+B)*x*x) \
  +Pb_plus_2_a4*(4*pi/(Pb_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Pb_plus_2_b4+B)*x*x) \
  +Pb_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Pb_plus_4_c = 8.084280;
Pb_plus_4_a1 = 32.124397;Pb_plus_4_a2 = 18.800301;Pb_plus_4_a3 = 12.017500;Pb_plus_4_a4 = 6.968860;
Pb_plus_4_b1 = 1.005660;Pb_plus_4_b2 = 6.109260;Pb_plus_4_b3 = 0.147041;Pb_plus_4_b4 = 14.714000;
Pb_plus_4_ff(x,B) = \
  +Pb_plus_4_a1*(4*pi/(Pb_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Pb_plus_4_b1+B)*x*x) \
  +Pb_plus_4_a2*(4*pi/(Pb_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Pb_plus_4_b2+B)*x*x) \
  +Pb_plus_4_a3*(4*pi/(Pb_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Pb_plus_4_b3+B)*x*x) \
  +Pb_plus_4_a4*(4*pi/(Pb_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Pb_plus_4_b4+B)*x*x) \
  +Pb_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=83
Bi_c = 13.578199;
Bi_a1 = 33.368900;Bi_a2 = 12.951000;Bi_a3 = 16.587700;Bi_a4 = 6.469200;
Bi_b1 = 0.704000;Bi_b2 = 2.923800;Bi_b3 = 8.793700;Bi_b4 = 48.009300;
Bi_ff(x,B) = \
  +Bi_a1*(4*pi/(Bi_b1+B))**1.5*safexp(-4*pi**2/(Bi_b1+B)*x*x) \
  +Bi_a2*(4*pi/(Bi_b2+B))**1.5*safexp(-4*pi**2/(Bi_b2+B)*x*x) \
  +Bi_a3*(4*pi/(Bi_b3+B))**1.5*safexp(-4*pi**2/(Bi_b3+B)*x*x) \
  +Bi_a4*(4*pi/(Bi_b4+B))**1.5*safexp(-4*pi**2/(Bi_b4+B)*x*x) \
  +Bi_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=80
Bi_plus_3_c = 12.471100;
Bi_plus_3_a1 = 21.805300;Bi_plus_3_a2 = 19.502600;Bi_plus_3_a3 = 19.105301;Bi_plus_3_a4 = 7.102950;
Bi_plus_3_b1 = 1.235600;Bi_plus_3_b2 = 6.241490;Bi_plus_3_b3 = 0.469999;Bi_plus_3_b4 = 20.318501;
Bi_plus_3_ff(x,B) = \
  +Bi_plus_3_a1*(4*pi/(Bi_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Bi_plus_3_b1+B)*x*x) \
  +Bi_plus_3_a2*(4*pi/(Bi_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Bi_plus_3_b2+B)*x*x) \
  +Bi_plus_3_a3*(4*pi/(Bi_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Bi_plus_3_b3+B)*x*x) \
  +Bi_plus_3_a4*(4*pi/(Bi_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Bi_plus_3_b4+B)*x*x) \
  +Bi_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=78
Bi_plus_5_c = -6.799400;
Bi_plus_5_a1 = 33.536400;Bi_plus_5_a2 = 25.094601;Bi_plus_5_a3 = 19.249699;Bi_plus_5_a4 = 6.915550;
Bi_plus_5_b1 = 0.916540;Bi_plus_5_b2 = 0.390420;Bi_plus_5_b3 = 5.714140;Bi_plus_5_b4 = 12.828500;
Bi_plus_5_ff(x,B) = \
  +Bi_plus_5_a1*(4*pi/(Bi_plus_5_b1+B))**1.5*safexp(-4*pi**2/(Bi_plus_5_b1+B)*x*x) \
  +Bi_plus_5_a2*(4*pi/(Bi_plus_5_b2+B))**1.5*safexp(-4*pi**2/(Bi_plus_5_b2+B)*x*x) \
  +Bi_plus_5_a3*(4*pi/(Bi_plus_5_b3+B))**1.5*safexp(-4*pi**2/(Bi_plus_5_b3+B)*x*x) \
  +Bi_plus_5_a4*(4*pi/(Bi_plus_5_b4+B))**1.5*safexp(-4*pi**2/(Bi_plus_5_b4+B)*x*x) \
  +Bi_plus_5_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=84
Po_c = 13.677000;
Po_a1 = 34.672600;Po_a2 = 15.473300;Po_a3 = 13.113800;Po_a4 = 7.025880;
Po_b1 = 0.700999;Po_b2 = 3.550780;Po_b3 = 9.556419;Po_b4 = 47.004501;
Po_ff(x,B) = \
  +Po_a1*(4*pi/(Po_b1+B))**1.5*safexp(-4*pi**2/(Po_b1+B)*x*x) \
  +Po_a2*(4*pi/(Po_b2+B))**1.5*safexp(-4*pi**2/(Po_b2+B)*x*x) \
  +Po_a3*(4*pi/(Po_b3+B))**1.5*safexp(-4*pi**2/(Po_b3+B)*x*x) \
  +Po_a4*(4*pi/(Po_b4+B))**1.5*safexp(-4*pi**2/(Po_b4+B)*x*x) \
  +Po_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=85
At_c = 13.710800;
At_a1 = 35.316299;At_a2 = 19.021099;At_a3 = 9.498870;At_a4 = 7.425180;
At_b1 = 0.685870;At_b2 = 3.974580;At_b3 = 11.382400;At_b4 = 45.471500;
At_ff(x,B) = \
  +At_a1*(4*pi/(At_b1+B))**1.5*safexp(-4*pi**2/(At_b1+B)*x*x) \
  +At_a2*(4*pi/(At_b2+B))**1.5*safexp(-4*pi**2/(At_b2+B)*x*x) \
  +At_a3*(4*pi/(At_b3+B))**1.5*safexp(-4*pi**2/(At_b3+B)*x*x) \
  +At_a4*(4*pi/(At_b4+B))**1.5*safexp(-4*pi**2/(At_b4+B)*x*x) \
  +At_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=86
Rn_c = 13.690500;
Rn_a1 = 35.563099;Rn_a2 = 21.281601;Rn_a3 = 8.003700;Rn_a4 = 7.443300;
Rn_b1 = 0.663100;Rn_b2 = 4.069100;Rn_b3 = 14.042200;Rn_b4 = 44.247299;
Rn_ff(x,B) = \
  +Rn_a1*(4*pi/(Rn_b1+B))**1.5*safexp(-4*pi**2/(Rn_b1+B)*x*x) \
  +Rn_a2*(4*pi/(Rn_b2+B))**1.5*safexp(-4*pi**2/(Rn_b2+B)*x*x) \
  +Rn_a3*(4*pi/(Rn_b3+B))**1.5*safexp(-4*pi**2/(Rn_b3+B)*x*x) \
  +Rn_a4*(4*pi/(Rn_b4+B))**1.5*safexp(-4*pi**2/(Rn_b4+B)*x*x) \
  +Rn_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=87
Fr_c = 13.724700;
Fr_a1 = 35.929901;Fr_a2 = 23.054699;Fr_a3 = 12.143900;Fr_a4 = 2.112530;
Fr_b1 = 0.646453;Fr_b2 = 4.176190;Fr_b3 = 23.105200;Fr_b4 = 150.644989;
Fr_ff(x,B) = \
  +Fr_a1*(4*pi/(Fr_b1+B))**1.5*safexp(-4*pi**2/(Fr_b1+B)*x*x) \
  +Fr_a2*(4*pi/(Fr_b2+B))**1.5*safexp(-4*pi**2/(Fr_b2+B)*x*x) \
  +Fr_a3*(4*pi/(Fr_b3+B))**1.5*safexp(-4*pi**2/(Fr_b3+B)*x*x) \
  +Fr_a4*(4*pi/(Fr_b4+B))**1.5*safexp(-4*pi**2/(Fr_b4+B)*x*x) \
  +Fr_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=88
Ra_c = 13.621099;
Ra_a1 = 35.763000;Ra_a2 = 22.906399;Ra_a3 = 12.473900;Ra_a4 = 3.210970;
Ra_b1 = 0.616341;Ra_b2 = 3.871350;Ra_b3 = 19.988701;Ra_b4 = 142.324997;
Ra_ff(x,B) = \
  +Ra_a1*(4*pi/(Ra_b1+B))**1.5*safexp(-4*pi**2/(Ra_b1+B)*x*x) \
  +Ra_a2*(4*pi/(Ra_b2+B))**1.5*safexp(-4*pi**2/(Ra_b2+B)*x*x) \
  +Ra_a3*(4*pi/(Ra_b3+B))**1.5*safexp(-4*pi**2/(Ra_b3+B)*x*x) \
  +Ra_a4*(4*pi/(Ra_b4+B))**1.5*safexp(-4*pi**2/(Ra_b4+B)*x*x) \
  +Ra_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=86
Ra_plus_2_c = 13.543100;
Ra_plus_2_a1 = 35.215000;Ra_plus_2_a2 = 21.670000;Ra_plus_2_a3 = 7.913420;Ra_plus_2_a4 = 7.650780;
Ra_plus_2_b1 = 0.604909;Ra_plus_2_b2 = 3.576700;Ra_plus_2_b3 = 12.601000;Ra_plus_2_b4 = 29.843599;
Ra_plus_2_ff(x,B) = \
  +Ra_plus_2_a1*(4*pi/(Ra_plus_2_b1+B))**1.5*safexp(-4*pi**2/(Ra_plus_2_b1+B)*x*x) \
  +Ra_plus_2_a2*(4*pi/(Ra_plus_2_b2+B))**1.5*safexp(-4*pi**2/(Ra_plus_2_b2+B)*x*x) \
  +Ra_plus_2_a3*(4*pi/(Ra_plus_2_b3+B))**1.5*safexp(-4*pi**2/(Ra_plus_2_b3+B)*x*x) \
  +Ra_plus_2_a4*(4*pi/(Ra_plus_2_b4+B))**1.5*safexp(-4*pi**2/(Ra_plus_2_b4+B)*x*x) \
  +Ra_plus_2_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=89
Ac_c = 13.526600;
Ac_a1 = 35.659698;Ac_a2 = 23.103199;Ac_a3 = 12.597700;Ac_a4 = 4.086550;
Ac_b1 = 0.589092;Ac_b2 = 3.651550;Ac_b3 = 18.598999;Ac_b4 = 117.019997;
Ac_ff(x,B) = \
  +Ac_a1*(4*pi/(Ac_b1+B))**1.5*safexp(-4*pi**2/(Ac_b1+B)*x*x) \
  +Ac_a2*(4*pi/(Ac_b2+B))**1.5*safexp(-4*pi**2/(Ac_b2+B)*x*x) \
  +Ac_a3*(4*pi/(Ac_b3+B))**1.5*safexp(-4*pi**2/(Ac_b3+B)*x*x) \
  +Ac_a4*(4*pi/(Ac_b4+B))**1.5*safexp(-4*pi**2/(Ac_b4+B)*x*x) \
  +Ac_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=86
Ac_plus_3_c = 13.463699;
Ac_plus_3_a1 = 35.173599;Ac_plus_3_a2 = 22.111200;Ac_plus_3_a3 = 8.192160;Ac_plus_3_a4 = 7.055450;
Ac_plus_3_b1 = 0.579689;Ac_plus_3_b2 = 3.414370;Ac_plus_3_b3 = 12.918700;Ac_plus_3_b4 = 25.944300;
Ac_plus_3_ff(x,B) = \
  +Ac_plus_3_a1*(4*pi/(Ac_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Ac_plus_3_b1+B)*x*x) \
  +Ac_plus_3_a2*(4*pi/(Ac_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Ac_plus_3_b2+B)*x*x) \
  +Ac_plus_3_a3*(4*pi/(Ac_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Ac_plus_3_b3+B)*x*x) \
  +Ac_plus_3_a4*(4*pi/(Ac_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Ac_plus_3_b4+B)*x*x) \
  +Ac_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=90
Th_c = 13.431400;
Th_a1 = 35.564499;Th_a2 = 23.421900;Th_a3 = 12.747300;Th_a4 = 4.807030;
Th_b1 = 0.563359;Th_b2 = 3.462040;Th_b3 = 17.830900;Th_b4 = 99.172195;
Th_ff(x,B) = \
  +Th_a1*(4*pi/(Th_b1+B))**1.5*safexp(-4*pi**2/(Th_b1+B)*x*x) \
  +Th_a2*(4*pi/(Th_b2+B))**1.5*safexp(-4*pi**2/(Th_b2+B)*x*x) \
  +Th_a3*(4*pi/(Th_b3+B))**1.5*safexp(-4*pi**2/(Th_b3+B)*x*x) \
  +Th_a4*(4*pi/(Th_b4+B))**1.5*safexp(-4*pi**2/(Th_b4+B)*x*x) \
  +Th_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=86
Th_plus_4_c = 13.375999;
Th_plus_4_a1 = 35.100700;Th_plus_4_a2 = 22.441799;Th_plus_4_a3 = 9.785540;Th_plus_4_a4 = 5.294440;
Th_plus_4_b1 = 0.555054;Th_plus_4_b2 = 3.244980;Th_plus_4_b3 = 13.466100;Th_plus_4_b4 = 23.953300;
Th_plus_4_ff(x,B) = \
  +Th_plus_4_a1*(4*pi/(Th_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Th_plus_4_b1+B)*x*x) \
  +Th_plus_4_a2*(4*pi/(Th_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Th_plus_4_b2+B)*x*x) \
  +Th_plus_4_a3*(4*pi/(Th_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Th_plus_4_b3+B)*x*x) \
  +Th_plus_4_a4*(4*pi/(Th_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Th_plus_4_b4+B)*x*x) \
  +Th_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=91
Pa_c = 13.428699;
Pa_a1 = 35.884701;Pa_a2 = 23.294800;Pa_a3 = 14.189100;Pa_a4 = 4.172870;
Pa_b1 = 0.547751;Pa_b2 = 3.415190;Pa_b3 = 16.923500;Pa_b4 = 105.250999;
Pa_ff(x,B) = \
  +Pa_a1*(4*pi/(Pa_b1+B))**1.5*safexp(-4*pi**2/(Pa_b1+B)*x*x) \
  +Pa_a2*(4*pi/(Pa_b2+B))**1.5*safexp(-4*pi**2/(Pa_b2+B)*x*x) \
  +Pa_a3*(4*pi/(Pa_b3+B))**1.5*safexp(-4*pi**2/(Pa_b3+B)*x*x) \
  +Pa_a4*(4*pi/(Pa_b4+B))**1.5*safexp(-4*pi**2/(Pa_b4+B)*x*x) \
  +Pa_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=92
U_c = 13.396600;
U_a1 = 36.022800;U_a2 = 23.412800;U_a3 = 14.949100;U_a4 = 4.188000;
U_b1 = 0.529300;U_b2 = 3.325300;U_b3 = 16.092699;U_b4 = 100.612999;
U_ff(x,B) = \
  +U_a1*(4*pi/(U_b1+B))**1.5*safexp(-4*pi**2/(U_b1+B)*x*x) \
  +U_a2*(4*pi/(U_b2+B))**1.5*safexp(-4*pi**2/(U_b2+B)*x*x) \
  +U_a3*(4*pi/(U_b3+B))**1.5*safexp(-4*pi**2/(U_b3+B)*x*x) \
  +U_a4*(4*pi/(U_b4+B))**1.5*safexp(-4*pi**2/(U_b4+B)*x*x) \
  +U_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=89
U_plus_3_c = 13.309200;
U_plus_3_a1 = 35.574699;U_plus_3_a2 = 22.525900;U_plus_3_a3 = 12.216499;U_plus_3_a4 = 5.370730;
U_plus_3_b1 = 0.520480;U_plus_3_b2 = 3.122930;U_plus_3_b3 = 12.714800;U_plus_3_b4 = 26.339399;
U_plus_3_ff(x,B) = \
  +U_plus_3_a1*(4*pi/(U_plus_3_b1+B))**1.5*safexp(-4*pi**2/(U_plus_3_b1+B)*x*x) \
  +U_plus_3_a2*(4*pi/(U_plus_3_b2+B))**1.5*safexp(-4*pi**2/(U_plus_3_b2+B)*x*x) \
  +U_plus_3_a3*(4*pi/(U_plus_3_b3+B))**1.5*safexp(-4*pi**2/(U_plus_3_b3+B)*x*x) \
  +U_plus_3_a4*(4*pi/(U_plus_3_b4+B))**1.5*safexp(-4*pi**2/(U_plus_3_b4+B)*x*x) \
  +U_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=88
U_plus_4_c = 13.267099;
U_plus_4_a1 = 35.371498;U_plus_4_a2 = 22.532600;U_plus_4_a3 = 12.029100;U_plus_4_a4 = 4.798400;
U_plus_4_b1 = 0.516598;U_plus_4_b2 = 3.050530;U_plus_4_b3 = 12.572300;U_plus_4_b4 = 23.458200;
U_plus_4_ff(x,B) = \
  +U_plus_4_a1*(4*pi/(U_plus_4_b1+B))**1.5*safexp(-4*pi**2/(U_plus_4_b1+B)*x*x) \
  +U_plus_4_a2*(4*pi/(U_plus_4_b2+B))**1.5*safexp(-4*pi**2/(U_plus_4_b2+B)*x*x) \
  +U_plus_4_a3*(4*pi/(U_plus_4_b3+B))**1.5*safexp(-4*pi**2/(U_plus_4_b3+B)*x*x) \
  +U_plus_4_a4*(4*pi/(U_plus_4_b4+B))**1.5*safexp(-4*pi**2/(U_plus_4_b4+B)*x*x) \
  +U_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=86
U_plus_6_c = 13.166500;
U_plus_6_a1 = 34.850899;U_plus_6_a2 = 22.758400;U_plus_6_a3 = 14.009900;U_plus_6_a4 = 1.214570;
U_plus_6_b1 = 0.507079;U_plus_6_b2 = 2.890300;U_plus_6_b3 = 13.176700;U_plus_6_b4 = 25.201700;
U_plus_6_ff(x,B) = \
  +U_plus_6_a1*(4*pi/(U_plus_6_b1+B))**1.5*safexp(-4*pi**2/(U_plus_6_b1+B)*x*x) \
  +U_plus_6_a2*(4*pi/(U_plus_6_b2+B))**1.5*safexp(-4*pi**2/(U_plus_6_b2+B)*x*x) \
  +U_plus_6_a3*(4*pi/(U_plus_6_b3+B))**1.5*safexp(-4*pi**2/(U_plus_6_b3+B)*x*x) \
  +U_plus_6_a4*(4*pi/(U_plus_6_b4+B))**1.5*safexp(-4*pi**2/(U_plus_6_b4+B)*x*x) \
  +U_plus_6_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=93
Np_c = 13.357300;
Np_a1 = 36.187401;Np_a2 = 23.596399;Np_a3 = 15.640200;Np_a4 = 4.185500;
Np_b1 = 0.511929;Np_b2 = 3.253960;Np_b3 = 15.362200;Np_b4 = 97.490799;
Np_ff(x,B) = \
  +Np_a1*(4*pi/(Np_b1+B))**1.5*safexp(-4*pi**2/(Np_b1+B)*x*x) \
  +Np_a2*(4*pi/(Np_b2+B))**1.5*safexp(-4*pi**2/(Np_b2+B)*x*x) \
  +Np_a3*(4*pi/(Np_b3+B))**1.5*safexp(-4*pi**2/(Np_b3+B)*x*x) \
  +Np_a4*(4*pi/(Np_b4+B))**1.5*safexp(-4*pi**2/(Np_b4+B)*x*x) \
  +Np_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=90
Np_plus_3_c = 13.254400;
Np_plus_3_a1 = 35.707397;Np_plus_3_a2 = 22.612999;Np_plus_3_a3 = 12.989799;Np_plus_3_a4 = 5.432270;
Np_plus_3_b1 = 0.502322;Np_plus_3_b2 = 3.038070;Np_plus_3_b3 = 12.144899;Np_plus_3_b4 = 25.492800;
Np_plus_3_ff(x,B) = \
  +Np_plus_3_a1*(4*pi/(Np_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Np_plus_3_b1+B)*x*x) \
  +Np_plus_3_a2*(4*pi/(Np_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Np_plus_3_b2+B)*x*x) \
  +Np_plus_3_a3*(4*pi/(Np_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Np_plus_3_b3+B)*x*x) \
  +Np_plus_3_a4*(4*pi/(Np_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Np_plus_3_b4+B)*x*x) \
  +Np_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=89
Np_plus_4_c = 13.211599;
Np_plus_4_a1 = 35.510300;Np_plus_4_a2 = 22.578699;Np_plus_4_a3 = 12.776600;Np_plus_4_a4 = 4.921590;
Np_plus_4_b1 = 0.498626;Np_plus_4_b2 = 2.966270;Np_plus_4_b3 = 11.948400;Np_plus_4_b4 = 22.750200;
Np_plus_4_ff(x,B) = \
  +Np_plus_4_a1*(4*pi/(Np_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Np_plus_4_b1+B)*x*x) \
  +Np_plus_4_a2*(4*pi/(Np_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Np_plus_4_b2+B)*x*x) \
  +Np_plus_4_a3*(4*pi/(Np_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Np_plus_4_b3+B)*x*x) \
  +Np_plus_4_a4*(4*pi/(Np_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Np_plus_4_b4+B)*x*x) \
  +Np_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=87
Np_plus_6_c = 13.113000;
Np_plus_6_a1 = 35.013599;Np_plus_6_a2 = 22.728600;Np_plus_6_a3 = 14.388400;Np_plus_6_a4 = 1.756690;
Np_plus_6_b1 = 0.489810;Np_plus_6_b2 = 2.810990;Np_plus_6_b3 = 12.330000;Np_plus_6_b4 = 22.658100;
Np_plus_6_ff(x,B) = \
  +Np_plus_6_a1*(4*pi/(Np_plus_6_b1+B))**1.5*safexp(-4*pi**2/(Np_plus_6_b1+B)*x*x) \
  +Np_plus_6_a2*(4*pi/(Np_plus_6_b2+B))**1.5*safexp(-4*pi**2/(Np_plus_6_b2+B)*x*x) \
  +Np_plus_6_a3*(4*pi/(Np_plus_6_b3+B))**1.5*safexp(-4*pi**2/(Np_plus_6_b3+B)*x*x) \
  +Np_plus_6_a4*(4*pi/(Np_plus_6_b4+B))**1.5*safexp(-4*pi**2/(Np_plus_6_b4+B)*x*x) \
  +Np_plus_6_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=94
Pu_c = 13.381200;
Pu_a1 = 36.525398;Pu_a2 = 23.808300;Pu_a3 = 16.770700;Pu_a4 = 3.479470;
Pu_b1 = 0.499384;Pu_b2 = 3.263710;Pu_b3 = 14.945499;Pu_b4 = 105.979996;
Pu_ff(x,B) = \
  +Pu_a1*(4*pi/(Pu_b1+B))**1.5*safexp(-4*pi**2/(Pu_b1+B)*x*x) \
  +Pu_a2*(4*pi/(Pu_b2+B))**1.5*safexp(-4*pi**2/(Pu_b2+B)*x*x) \
  +Pu_a3*(4*pi/(Pu_b3+B))**1.5*safexp(-4*pi**2/(Pu_b3+B)*x*x) \
  +Pu_a4*(4*pi/(Pu_b4+B))**1.5*safexp(-4*pi**2/(Pu_b4+B)*x*x) \
  +Pu_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=91
Pu_plus_3_c = 13.199100;
Pu_plus_3_a1 = 35.840000;Pu_plus_3_a2 = 22.716900;Pu_plus_3_a3 = 13.580700;Pu_plus_3_a4 = 5.660160;
Pu_plus_3_b1 = 0.484936;Pu_plus_3_b2 = 2.961180;Pu_plus_3_b3 = 11.533100;Pu_plus_3_b4 = 24.399200;
Pu_plus_3_ff(x,B) = \
  +Pu_plus_3_a1*(4*pi/(Pu_plus_3_b1+B))**1.5*safexp(-4*pi**2/(Pu_plus_3_b1+B)*x*x) \
  +Pu_plus_3_a2*(4*pi/(Pu_plus_3_b2+B))**1.5*safexp(-4*pi**2/(Pu_plus_3_b2+B)*x*x) \
  +Pu_plus_3_a3*(4*pi/(Pu_plus_3_b3+B))**1.5*safexp(-4*pi**2/(Pu_plus_3_b3+B)*x*x) \
  +Pu_plus_3_a4*(4*pi/(Pu_plus_3_b4+B))**1.5*safexp(-4*pi**2/(Pu_plus_3_b4+B)*x*x) \
  +Pu_plus_3_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=90
Pu_plus_4_c = 13.155500;
Pu_plus_4_a1 = 35.649300;Pu_plus_4_a2 = 22.646000;Pu_plus_4_a3 = 13.359500;Pu_plus_4_a4 = 5.188310;
Pu_plus_4_b1 = 0.481422;Pu_plus_4_b2 = 2.890200;Pu_plus_4_b3 = 11.316000;Pu_plus_4_b4 = 21.830099;
Pu_plus_4_ff(x,B) = \
  +Pu_plus_4_a1*(4*pi/(Pu_plus_4_b1+B))**1.5*safexp(-4*pi**2/(Pu_plus_4_b1+B)*x*x) \
  +Pu_plus_4_a2*(4*pi/(Pu_plus_4_b2+B))**1.5*safexp(-4*pi**2/(Pu_plus_4_b2+B)*x*x) \
  +Pu_plus_4_a3*(4*pi/(Pu_plus_4_b3+B))**1.5*safexp(-4*pi**2/(Pu_plus_4_b3+B)*x*x) \
  +Pu_plus_4_a4*(4*pi/(Pu_plus_4_b4+B))**1.5*safexp(-4*pi**2/(Pu_plus_4_b4+B)*x*x) \
  +Pu_plus_4_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=88
Pu_plus_6_c = 13.058200;
Pu_plus_6_a1 = 35.173599;Pu_plus_6_a2 = 22.718100;Pu_plus_6_a3 = 14.763500;Pu_plus_6_a4 = 2.286780;
Pu_plus_6_b1 = 0.473204;Pu_plus_6_b2 = 2.738480;Pu_plus_6_b3 = 11.552999;Pu_plus_6_b4 = 20.930300;
Pu_plus_6_ff(x,B) = \
  +Pu_plus_6_a1*(4*pi/(Pu_plus_6_b1+B))**1.5*safexp(-4*pi**2/(Pu_plus_6_b1+B)*x*x) \
  +Pu_plus_6_a2*(4*pi/(Pu_plus_6_b2+B))**1.5*safexp(-4*pi**2/(Pu_plus_6_b2+B)*x*x) \
  +Pu_plus_6_a3*(4*pi/(Pu_plus_6_b3+B))**1.5*safexp(-4*pi**2/(Pu_plus_6_b3+B)*x*x) \
  +Pu_plus_6_a4*(4*pi/(Pu_plus_6_b4+B))**1.5*safexp(-4*pi**2/(Pu_plus_6_b4+B)*x*x) \
  +Pu_plus_6_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=95
Am_c = 13.359200;
Am_a1 = 36.670601;Am_a2 = 24.099199;Am_a3 = 17.341499;Am_a4 = 3.493310;
Am_b1 = 0.483629;Am_b2 = 3.206470;Am_b3 = 14.313600;Am_b4 = 102.272995;
Am_ff(x,B) = \
  +Am_a1*(4*pi/(Am_b1+B))**1.5*safexp(-4*pi**2/(Am_b1+B)*x*x) \
  +Am_a2*(4*pi/(Am_b2+B))**1.5*safexp(-4*pi**2/(Am_b2+B)*x*x) \
  +Am_a3*(4*pi/(Am_b3+B))**1.5*safexp(-4*pi**2/(Am_b3+B)*x*x) \
  +Am_a4*(4*pi/(Am_b4+B))**1.5*safexp(-4*pi**2/(Am_b4+B)*x*x) \
  +Am_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=96
Cm_c = 13.288700;
Cm_a1 = 36.648800;Cm_a2 = 24.409599;Cm_a3 = 17.399000;Cm_a4 = 4.216650;
Cm_b1 = 0.465154;Cm_b2 = 3.089970;Cm_b3 = 13.434600;Cm_b4 = 88.483398;
Cm_ff(x,B) = \
  +Cm_a1*(4*pi/(Cm_b1+B))**1.5*safexp(-4*pi**2/(Cm_b1+B)*x*x) \
  +Cm_a2*(4*pi/(Cm_b2+B))**1.5*safexp(-4*pi**2/(Cm_b2+B)*x*x) \
  +Cm_a3*(4*pi/(Cm_b3+B))**1.5*safexp(-4*pi**2/(Cm_b3+B)*x*x) \
  +Cm_a4*(4*pi/(Cm_b4+B))**1.5*safexp(-4*pi**2/(Cm_b4+B)*x*x) \
  +Cm_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=97
Bk_c = 13.275400;
Bk_a1 = 36.788101;Bk_a2 = 24.773600;Bk_a3 = 17.891899;Bk_a4 = 4.232840;
Bk_b1 = 0.451018;Bk_b2 = 3.046190;Bk_b3 = 12.894600;Bk_b4 = 86.002998;
Bk_ff(x,B) = \
  +Bk_a1*(4*pi/(Bk_b1+B))**1.5*safexp(-4*pi**2/(Bk_b1+B)*x*x) \
  +Bk_a2*(4*pi/(Bk_b2+B))**1.5*safexp(-4*pi**2/(Bk_b2+B)*x*x) \
  +Bk_a3*(4*pi/(Bk_b3+B))**1.5*safexp(-4*pi**2/(Bk_b3+B)*x*x) \
  +Bk_a4*(4*pi/(Bk_b4+B))**1.5*safexp(-4*pi**2/(Bk_b4+B)*x*x) \
  +Bk_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);
Z=98
Cf_c = 13.267400;
Cf_a1 = 36.918499;Cf_a2 = 25.199499;Cf_a3 = 18.331699;Cf_a4 = 4.243910;
Cf_b1 = 0.437533;Cf_b2 = 3.007750;Cf_b3 = 12.404400;Cf_b4 = 83.788101;
Cf_ff(x,B) = \
  +Cf_a1*(4*pi/(Cf_b1+B))**1.5*safexp(-4*pi**2/(Cf_b1+B)*x*x) \
  +Cf_a2*(4*pi/(Cf_b2+B))**1.5*safexp(-4*pi**2/(Cf_b2+B)*x*x) \
  +Cf_a3*(4*pi/(Cf_b3+B))**1.5*safexp(-4*pi**2/(Cf_b3+B)*x*x) \
  +Cf_a4*(4*pi/(Cf_b4+B))**1.5*safexp(-4*pi**2/(Cf_b4+B)*x*x) \
  +Cf_c *(4*pi/(1e-99+B))**1.5*safexp(-4*pi**2/(1e-99+B)*x*x);