# here "s" is sin(theta)/lambda where theta is the Bragg angle: half the angle between
#the incident and scattered rays, and lambda is the wavelength in Angstrom
# B is the Debye-Waller-Ott "temperature factor" (same units as a Protein Data Bank file)
# retun values of the Ee_sf(s,B) are in "electrons" (ratio to the scattering from a point electron)

# single electron (becomes a point electron when B=0)
Ano_sf(s,B) = exp(-B*s*s);
Z=1
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_sf(s)=H_a1*exp(-H_b1*s*s)+H_a2*exp(-H_b2*s*s)+H_a3*exp(-H_b3*s*s)+H_a4*exp(-H_b4*s*s)+H_c
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_sf(s)=H_minus_1_a1*exp(-H_minus_1_b1*s*s)+H_minus_1_a2*exp(-H_minus_1_b2*s*s)+H_minus_1_a3*exp(-H_minus_1_b3*s*s)+H_minus_1_a4*exp(-H_minus_1_b4*s*s)+H_minus_1_c
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_sf(s)=He_a1*exp(-He_b1*s*s)+He_a2*exp(-He_b2*s*s)+He_a3*exp(-He_b3*s*s)+He_a4*exp(-He_b4*s*s)+He_c
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_sf(s)=Li_a1*exp(-Li_b1*s*s)+Li_a2*exp(-Li_b2*s*s)+Li_a3*exp(-Li_b3*s*s)+Li_a4*exp(-Li_b4*s*s)+Li_c
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_sf(s)=Li_plus_1_a1*exp(-Li_plus_1_b1*s*s)+Li_plus_1_a2*exp(-Li_plus_1_b2*s*s)+Li_plus_1_a3*exp(-Li_plus_1_b3*s*s)+Li_plus_1_a4*exp(-Li_plus_1_b4*s*s)+Li_plus_1_c
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_sf(s)=Be_a1*exp(-Be_b1*s*s)+Be_a2*exp(-Be_b2*s*s)+Be_a3*exp(-Be_b3*s*s)+Be_a4*exp(-Be_b4*s*s)+Be_c
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_sf(s)=Be_plus_2_a1*exp(-Be_plus_2_b1*s*s)+Be_plus_2_a2*exp(-Be_plus_2_b2*s*s)+Be_plus_2_a3*exp(-Be_plus_2_b3*s*s)+Be_plus_2_a4*exp(-Be_plus_2_b4*s*s)+Be_plus_2_c
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_sf(s)=B_a1*exp(-B_b1*s*s)+B_a2*exp(-B_b2*s*s)+B_a3*exp(-B_b3*s*s)+B_a4*exp(-B_b4*s*s)+B_c
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_sf(s)=C_a1*exp(-C_b1*s*s)+C_a2*exp(-C_b2*s*s)+C_a3*exp(-C_b3*s*s)+C_a4*exp(-C_b4*s*s)+C_c
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_sf(s)=Cv_a1*exp(-Cv_b1*s*s)+Cv_a2*exp(-Cv_b2*s*s)+Cv_a3*exp(-Cv_b3*s*s)+Cv_a4*exp(-Cv_b4*s*s)+Cv_c
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_sf(s)=N_a1*exp(-N_b1*s*s)+N_a2*exp(-N_b2*s*s)+N_a3*exp(-N_b3*s*s)+N_a4*exp(-N_b4*s*s)+N_c
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_sf(s)=O_a1*exp(-O_b1*s*s)+O_a2*exp(-O_b2*s*s)+O_a3*exp(-O_b3*s*s)+O_a4*exp(-O_b4*s*s)+O_c
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_sf(s)=O_minus_1_a1*exp(-O_minus_1_b1*s*s)+O_minus_1_a2*exp(-O_minus_1_b2*s*s)+O_minus_1_a3*exp(-O_minus_1_b3*s*s)+O_minus_1_a4*exp(-O_minus_1_b4*s*s)+O_minus_1_c
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_sf(s)=F_a1*exp(-F_b1*s*s)+F_a2*exp(-F_b2*s*s)+F_a3*exp(-F_b3*s*s)+F_a4*exp(-F_b4*s*s)+F_c
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_sf(s)=F_minus_1_a1*exp(-F_minus_1_b1*s*s)+F_minus_1_a2*exp(-F_minus_1_b2*s*s)+F_minus_1_a3*exp(-F_minus_1_b3*s*s)+F_minus_1_a4*exp(-F_minus_1_b4*s*s)+F_minus_1_c
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_sf(s)=Ne_a1*exp(-Ne_b1*s*s)+Ne_a2*exp(-Ne_b2*s*s)+Ne_a3*exp(-Ne_b3*s*s)+Ne_a4*exp(-Ne_b4*s*s)+Ne_c
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_sf(s)=Na_a1*exp(-Na_b1*s*s)+Na_a2*exp(-Na_b2*s*s)+Na_a3*exp(-Na_b3*s*s)+Na_a4*exp(-Na_b4*s*s)+Na_c
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_sf(s)=Na_plus_1_a1*exp(-Na_plus_1_b1*s*s)+Na_plus_1_a2*exp(-Na_plus_1_b2*s*s)+Na_plus_1_a3*exp(-Na_plus_1_b3*s*s)+Na_plus_1_a4*exp(-Na_plus_1_b4*s*s)+Na_plus_1_c
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_sf(s)=Mg_a1*exp(-Mg_b1*s*s)+Mg_a2*exp(-Mg_b2*s*s)+Mg_a3*exp(-Mg_b3*s*s)+Mg_a4*exp(-Mg_b4*s*s)+Mg_c
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_sf(s)=Mg_plus_2_a1*exp(-Mg_plus_2_b1*s*s)+Mg_plus_2_a2*exp(-Mg_plus_2_b2*s*s)+Mg_plus_2_a3*exp(-Mg_plus_2_b3*s*s)+Mg_plus_2_a4*exp(-Mg_plus_2_b4*s*s)+Mg_plus_2_c
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_sf(s)=Al_a1*exp(-Al_b1*s*s)+Al_a2*exp(-Al_b2*s*s)+Al_a3*exp(-Al_b3*s*s)+Al_a4*exp(-Al_b4*s*s)+Al_c
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_sf(s)=Al_plus_3_a1*exp(-Al_plus_3_b1*s*s)+Al_plus_3_a2*exp(-Al_plus_3_b2*s*s)+Al_plus_3_a3*exp(-Al_plus_3_b3*s*s)+Al_plus_3_a4*exp(-Al_plus_3_b4*s*s)+Al_plus_3_c
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_sf(s)=Si_a1*exp(-Si_b1*s*s)+Si_a2*exp(-Si_b2*s*s)+Si_a3*exp(-Si_b3*s*s)+Si_a4*exp(-Si_b4*s*s)+Si_c
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_sf(s)=Siv_a1*exp(-Siv_b1*s*s)+Siv_a2*exp(-Siv_b2*s*s)+Siv_a3*exp(-Siv_b3*s*s)+Siv_a4*exp(-Siv_b4*s*s)+Siv_c
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_sf(s)=Si_plus_4_a1*exp(-Si_plus_4_b1*s*s)+Si_plus_4_a2*exp(-Si_plus_4_b2*s*s)+Si_plus_4_a3*exp(-Si_plus_4_b3*s*s)+Si_plus_4_a4*exp(-Si_plus_4_b4*s*s)+Si_plus_4_c
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_sf(s)=P_a1*exp(-P_b1*s*s)+P_a2*exp(-P_b2*s*s)+P_a3*exp(-P_b3*s*s)+P_a4*exp(-P_b4*s*s)+P_c
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_sf(s)=S_a1*exp(-S_b1*s*s)+S_a2*exp(-S_b2*s*s)+S_a3*exp(-S_b3*s*s)+S_a4*exp(-S_b4*s*s)+S_c
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_sf(s)=Cl_a1*exp(-Cl_b1*s*s)+Cl_a2*exp(-Cl_b2*s*s)+Cl_a3*exp(-Cl_b3*s*s)+Cl_a4*exp(-Cl_b4*s*s)+Cl_c
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_sf(s)=Cl_minus_1_a1*exp(-Cl_minus_1_b1*s*s)+Cl_minus_1_a2*exp(-Cl_minus_1_b2*s*s)+Cl_minus_1_a3*exp(-Cl_minus_1_b3*s*s)+Cl_minus_1_a4*exp(-Cl_minus_1_b4*s*s)+Cl_minus_1_c
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_sf(s)=Ar_a1*exp(-Ar_b1*s*s)+Ar_a2*exp(-Ar_b2*s*s)+Ar_a3*exp(-Ar_b3*s*s)+Ar_a4*exp(-Ar_b4*s*s)+Ar_c
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_sf(s)=K_a1*exp(-K_b1*s*s)+K_a2*exp(-K_b2*s*s)+K_a3*exp(-K_b3*s*s)+K_a4*exp(-K_b4*s*s)+K_c
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_sf(s)=K_plus_1_a1*exp(-K_plus_1_b1*s*s)+K_plus_1_a2*exp(-K_plus_1_b2*s*s)+K_plus_1_a3*exp(-K_plus_1_b3*s*s)+K_plus_1_a4*exp(-K_plus_1_b4*s*s)+K_plus_1_c
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_sf(s)=Ca_a1*exp(-Ca_b1*s*s)+Ca_a2*exp(-Ca_b2*s*s)+Ca_a3*exp(-Ca_b3*s*s)+Ca_a4*exp(-Ca_b4*s*s)+Ca_c
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_sf(s)=Ca_plus_2_a1*exp(-Ca_plus_2_b1*s*s)+Ca_plus_2_a2*exp(-Ca_plus_2_b2*s*s)+Ca_plus_2_a3*exp(-Ca_plus_2_b3*s*s)+Ca_plus_2_a4*exp(-Ca_plus_2_b4*s*s)+Ca_plus_2_c
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_sf(s)=Sc_a1*exp(-Sc_b1*s*s)+Sc_a2*exp(-Sc_b2*s*s)+Sc_a3*exp(-Sc_b3*s*s)+Sc_a4*exp(-Sc_b4*s*s)+Sc_c
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_sf(s)=Sc_plus_3_a1*exp(-Sc_plus_3_b1*s*s)+Sc_plus_3_a2*exp(-Sc_plus_3_b2*s*s)+Sc_plus_3_a3*exp(-Sc_plus_3_b3*s*s)+Sc_plus_3_a4*exp(-Sc_plus_3_b4*s*s)+Sc_plus_3_c
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_sf(s)=Ti_a1*exp(-Ti_b1*s*s)+Ti_a2*exp(-Ti_b2*s*s)+Ti_a3*exp(-Ti_b3*s*s)+Ti_a4*exp(-Ti_b4*s*s)+Ti_c
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_sf(s)=Ti_plus_2_a1*exp(-Ti_plus_2_b1*s*s)+Ti_plus_2_a2*exp(-Ti_plus_2_b2*s*s)+Ti_plus_2_a3*exp(-Ti_plus_2_b3*s*s)+Ti_plus_2_a4*exp(-Ti_plus_2_b4*s*s)+Ti_plus_2_c
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_sf(s)=Ti_plus_3_a1*exp(-Ti_plus_3_b1*s*s)+Ti_plus_3_a2*exp(-Ti_plus_3_b2*s*s)+Ti_plus_3_a3*exp(-Ti_plus_3_b3*s*s)+Ti_plus_3_a4*exp(-Ti_plus_3_b4*s*s)+Ti_plus_3_c
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_sf(s)=Ti_plus_4_a1*exp(-Ti_plus_4_b1*s*s)+Ti_plus_4_a2*exp(-Ti_plus_4_b2*s*s)+Ti_plus_4_a3*exp(-Ti_plus_4_b3*s*s)+Ti_plus_4_a4*exp(-Ti_plus_4_b4*s*s)+Ti_plus_4_c
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_sf(s)=V_a1*exp(-V_b1*s*s)+V_a2*exp(-V_b2*s*s)+V_a3*exp(-V_b3*s*s)+V_a4*exp(-V_b4*s*s)+V_c
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_sf(s)=V_plus_2_a1*exp(-V_plus_2_b1*s*s)+V_plus_2_a2*exp(-V_plus_2_b2*s*s)+V_plus_2_a3*exp(-V_plus_2_b3*s*s)+V_plus_2_a4*exp(-V_plus_2_b4*s*s)+V_plus_2_c
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_sf(s)=V_plus_3_a1*exp(-V_plus_3_b1*s*s)+V_plus_3_a2*exp(-V_plus_3_b2*s*s)+V_plus_3_a3*exp(-V_plus_3_b3*s*s)+V_plus_3_a4*exp(-V_plus_3_b4*s*s)+V_plus_3_c
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_sf(s)=V_plus_5_a1*exp(-V_plus_5_b1*s*s)+V_plus_5_a2*exp(-V_plus_5_b2*s*s)+V_plus_5_a3*exp(-V_plus_5_b3*s*s)+V_plus_5_a4*exp(-V_plus_5_b4*s*s)+V_plus_5_c
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_sf(s)=Cr_a1*exp(-Cr_b1*s*s)+Cr_a2*exp(-Cr_b2*s*s)+Cr_a3*exp(-Cr_b3*s*s)+Cr_a4*exp(-Cr_b4*s*s)+Cr_c
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_sf(s)=Cr_plus_2_a1*exp(-Cr_plus_2_b1*s*s)+Cr_plus_2_a2*exp(-Cr_plus_2_b2*s*s)+Cr_plus_2_a3*exp(-Cr_plus_2_b3*s*s)+Cr_plus_2_a4*exp(-Cr_plus_2_b4*s*s)+Cr_plus_2_c
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_sf(s)=Cr_plus_3_a1*exp(-Cr_plus_3_b1*s*s)+Cr_plus_3_a2*exp(-Cr_plus_3_b2*s*s)+Cr_plus_3_a3*exp(-Cr_plus_3_b3*s*s)+Cr_plus_3_a4*exp(-Cr_plus_3_b4*s*s)+Cr_plus_3_c
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_sf(s)=Mn_a1*exp(-Mn_b1*s*s)+Mn_a2*exp(-Mn_b2*s*s)+Mn_a3*exp(-Mn_b3*s*s)+Mn_a4*exp(-Mn_b4*s*s)+Mn_c
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_sf(s)=Mn_plus_2_a1*exp(-Mn_plus_2_b1*s*s)+Mn_plus_2_a2*exp(-Mn_plus_2_b2*s*s)+Mn_plus_2_a3*exp(-Mn_plus_2_b3*s*s)+Mn_plus_2_a4*exp(-Mn_plus_2_b4*s*s)+Mn_plus_2_c
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_sf(s)=Mn_plus_3_a1*exp(-Mn_plus_3_b1*s*s)+Mn_plus_3_a2*exp(-Mn_plus_3_b2*s*s)+Mn_plus_3_a3*exp(-Mn_plus_3_b3*s*s)+Mn_plus_3_a4*exp(-Mn_plus_3_b4*s*s)+Mn_plus_3_c
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_sf(s)=Mn_plus_4_a1*exp(-Mn_plus_4_b1*s*s)+Mn_plus_4_a2*exp(-Mn_plus_4_b2*s*s)+Mn_plus_4_a3*exp(-Mn_plus_4_b3*s*s)+Mn_plus_4_a4*exp(-Mn_plus_4_b4*s*s)+Mn_plus_4_c
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_sf(s)=Fe_a1*exp(-Fe_b1*s*s)+Fe_a2*exp(-Fe_b2*s*s)+Fe_a3*exp(-Fe_b3*s*s)+Fe_a4*exp(-Fe_b4*s*s)+Fe_c
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_sf(s)=Fe_plus_2_a1*exp(-Fe_plus_2_b1*s*s)+Fe_plus_2_a2*exp(-Fe_plus_2_b2*s*s)+Fe_plus_2_a3*exp(-Fe_plus_2_b3*s*s)+Fe_plus_2_a4*exp(-Fe_plus_2_b4*s*s)+Fe_plus_2_c
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_sf(s)=Fe_plus_3_a1*exp(-Fe_plus_3_b1*s*s)+Fe_plus_3_a2*exp(-Fe_plus_3_b2*s*s)+Fe_plus_3_a3*exp(-Fe_plus_3_b3*s*s)+Fe_plus_3_a4*exp(-Fe_plus_3_b4*s*s)+Fe_plus_3_c
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_sf(s)=Co_a1*exp(-Co_b1*s*s)+Co_a2*exp(-Co_b2*s*s)+Co_a3*exp(-Co_b3*s*s)+Co_a4*exp(-Co_b4*s*s)+Co_c
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_sf(s)=Co_plus_2_a1*exp(-Co_plus_2_b1*s*s)+Co_plus_2_a2*exp(-Co_plus_2_b2*s*s)+Co_plus_2_a3*exp(-Co_plus_2_b3*s*s)+Co_plus_2_a4*exp(-Co_plus_2_b4*s*s)+Co_plus_2_c
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_sf(s)=Co_plus_3_a1*exp(-Co_plus_3_b1*s*s)+Co_plus_3_a2*exp(-Co_plus_3_b2*s*s)+Co_plus_3_a3*exp(-Co_plus_3_b3*s*s)+Co_plus_3_a4*exp(-Co_plus_3_b4*s*s)+Co_plus_3_c
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_sf(s)=Ni_a1*exp(-Ni_b1*s*s)+Ni_a2*exp(-Ni_b2*s*s)+Ni_a3*exp(-Ni_b3*s*s)+Ni_a4*exp(-Ni_b4*s*s)+Ni_c
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_sf(s)=Ni_plus_2_a1*exp(-Ni_plus_2_b1*s*s)+Ni_plus_2_a2*exp(-Ni_plus_2_b2*s*s)+Ni_plus_2_a3*exp(-Ni_plus_2_b3*s*s)+Ni_plus_2_a4*exp(-Ni_plus_2_b4*s*s)+Ni_plus_2_c
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_sf(s)=Ni_plus_3_a1*exp(-Ni_plus_3_b1*s*s)+Ni_plus_3_a2*exp(-Ni_plus_3_b2*s*s)+Ni_plus_3_a3*exp(-Ni_plus_3_b3*s*s)+Ni_plus_3_a4*exp(-Ni_plus_3_b4*s*s)+Ni_plus_3_c
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_sf(s)=Cu_a1*exp(-Cu_b1*s*s)+Cu_a2*exp(-Cu_b2*s*s)+Cu_a3*exp(-Cu_b3*s*s)+Cu_a4*exp(-Cu_b4*s*s)+Cu_c
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_sf(s)=Cu_plus_1_a1*exp(-Cu_plus_1_b1*s*s)+Cu_plus_1_a2*exp(-Cu_plus_1_b2*s*s)+Cu_plus_1_a3*exp(-Cu_plus_1_b3*s*s)+Cu_plus_1_a4*exp(-Cu_plus_1_b4*s*s)+Cu_plus_1_c
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_sf(s)=Cu_plus_2_a1*exp(-Cu_plus_2_b1*s*s)+Cu_plus_2_a2*exp(-Cu_plus_2_b2*s*s)+Cu_plus_2_a3*exp(-Cu_plus_2_b3*s*s)+Cu_plus_2_a4*exp(-Cu_plus_2_b4*s*s)+Cu_plus_2_c
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_sf(s)=Zn_a1*exp(-Zn_b1*s*s)+Zn_a2*exp(-Zn_b2*s*s)+Zn_a3*exp(-Zn_b3*s*s)+Zn_a4*exp(-Zn_b4*s*s)+Zn_c
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_sf(s)=Zn_plus_2_a1*exp(-Zn_plus_2_b1*s*s)+Zn_plus_2_a2*exp(-Zn_plus_2_b2*s*s)+Zn_plus_2_a3*exp(-Zn_plus_2_b3*s*s)+Zn_plus_2_a4*exp(-Zn_plus_2_b4*s*s)+Zn_plus_2_c
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_sf(s)=Ga_a1*exp(-Ga_b1*s*s)+Ga_a2*exp(-Ga_b2*s*s)+Ga_a3*exp(-Ga_b3*s*s)+Ga_a4*exp(-Ga_b4*s*s)+Ga_c
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_sf(s)=Ga_plus_3_a1*exp(-Ga_plus_3_b1*s*s)+Ga_plus_3_a2*exp(-Ga_plus_3_b2*s*s)+Ga_plus_3_a3*exp(-Ga_plus_3_b3*s*s)+Ga_plus_3_a4*exp(-Ga_plus_3_b4*s*s)+Ga_plus_3_c
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_sf(s)=Ge_a1*exp(-Ge_b1*s*s)+Ge_a2*exp(-Ge_b2*s*s)+Ge_a3*exp(-Ge_b3*s*s)+Ge_a4*exp(-Ge_b4*s*s)+Ge_c
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_sf(s)=Ge_plus_4_a1*exp(-Ge_plus_4_b1*s*s)+Ge_plus_4_a2*exp(-Ge_plus_4_b2*s*s)+Ge_plus_4_a3*exp(-Ge_plus_4_b3*s*s)+Ge_plus_4_a4*exp(-Ge_plus_4_b4*s*s)+Ge_plus_4_c
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_sf(s)=As_a1*exp(-As_b1*s*s)+As_a2*exp(-As_b2*s*s)+As_a3*exp(-As_b3*s*s)+As_a4*exp(-As_b4*s*s)+As_c
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_sf(s)=Se_a1*exp(-Se_b1*s*s)+Se_a2*exp(-Se_b2*s*s)+Se_a3*exp(-Se_b3*s*s)+Se_a4*exp(-Se_b4*s*s)+Se_c
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_sf(s)=Br_a1*exp(-Br_b1*s*s)+Br_a2*exp(-Br_b2*s*s)+Br_a3*exp(-Br_b3*s*s)+Br_a4*exp(-Br_b4*s*s)+Br_c
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_sf(s)=Br_minus_1_a1*exp(-Br_minus_1_b1*s*s)+Br_minus_1_a2*exp(-Br_minus_1_b2*s*s)+Br_minus_1_a3*exp(-Br_minus_1_b3*s*s)+Br_minus_1_a4*exp(-Br_minus_1_b4*s*s)+Br_minus_1_c
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_sf(s)=Kr_a1*exp(-Kr_b1*s*s)+Kr_a2*exp(-Kr_b2*s*s)+Kr_a3*exp(-Kr_b3*s*s)+Kr_a4*exp(-Kr_b4*s*s)+Kr_c
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_sf(s)=Rb_a1*exp(-Rb_b1*s*s)+Rb_a2*exp(-Rb_b2*s*s)+Rb_a3*exp(-Rb_b3*s*s)+Rb_a4*exp(-Rb_b4*s*s)+Rb_c
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_sf(s)=Rb_plus_1_a1*exp(-Rb_plus_1_b1*s*s)+Rb_plus_1_a2*exp(-Rb_plus_1_b2*s*s)+Rb_plus_1_a3*exp(-Rb_plus_1_b3*s*s)+Rb_plus_1_a4*exp(-Rb_plus_1_b4*s*s)+Rb_plus_1_c
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_sf(s)=Sr_a1*exp(-Sr_b1*s*s)+Sr_a2*exp(-Sr_b2*s*s)+Sr_a3*exp(-Sr_b3*s*s)+Sr_a4*exp(-Sr_b4*s*s)+Sr_c
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_sf(s)=Sr_plus_2_a1*exp(-Sr_plus_2_b1*s*s)+Sr_plus_2_a2*exp(-Sr_plus_2_b2*s*s)+Sr_plus_2_a3*exp(-Sr_plus_2_b3*s*s)+Sr_plus_2_a4*exp(-Sr_plus_2_b4*s*s)+Sr_plus_2_c
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_sf(s)=Y_a1*exp(-Y_b1*s*s)+Y_a2*exp(-Y_b2*s*s)+Y_a3*exp(-Y_b3*s*s)+Y_a4*exp(-Y_b4*s*s)+Y_c
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_sf(s)=Y_plus_3_a1*exp(-Y_plus_3_b1*s*s)+Y_plus_3_a2*exp(-Y_plus_3_b2*s*s)+Y_plus_3_a3*exp(-Y_plus_3_b3*s*s)+Y_plus_3_a4*exp(-Y_plus_3_b4*s*s)+Y_plus_3_c
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_sf(s)=Zr_a1*exp(-Zr_b1*s*s)+Zr_a2*exp(-Zr_b2*s*s)+Zr_a3*exp(-Zr_b3*s*s)+Zr_a4*exp(-Zr_b4*s*s)+Zr_c
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_sf(s)=Zr_plus_4_a1*exp(-Zr_plus_4_b1*s*s)+Zr_plus_4_a2*exp(-Zr_plus_4_b2*s*s)+Zr_plus_4_a3*exp(-Zr_plus_4_b3*s*s)+Zr_plus_4_a4*exp(-Zr_plus_4_b4*s*s)+Zr_plus_4_c
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_sf(s)=Nb_a1*exp(-Nb_b1*s*s)+Nb_a2*exp(-Nb_b2*s*s)+Nb_a3*exp(-Nb_b3*s*s)+Nb_a4*exp(-Nb_b4*s*s)+Nb_c
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_sf(s)=Nb_plus_3_a1*exp(-Nb_plus_3_b1*s*s)+Nb_plus_3_a2*exp(-Nb_plus_3_b2*s*s)+Nb_plus_3_a3*exp(-Nb_plus_3_b3*s*s)+Nb_plus_3_a4*exp(-Nb_plus_3_b4*s*s)+Nb_plus_3_c
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_sf(s)=Nb_plus_5_a1*exp(-Nb_plus_5_b1*s*s)+Nb_plus_5_a2*exp(-Nb_plus_5_b2*s*s)+Nb_plus_5_a3*exp(-Nb_plus_5_b3*s*s)+Nb_plus_5_a4*exp(-Nb_plus_5_b4*s*s)+Nb_plus_5_c
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_sf(s)=Mo_a1*exp(-Mo_b1*s*s)+Mo_a2*exp(-Mo_b2*s*s)+Mo_a3*exp(-Mo_b3*s*s)+Mo_a4*exp(-Mo_b4*s*s)+Mo_c
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_sf(s)=Mo_plus_3_a1*exp(-Mo_plus_3_b1*s*s)+Mo_plus_3_a2*exp(-Mo_plus_3_b2*s*s)+Mo_plus_3_a3*exp(-Mo_plus_3_b3*s*s)+Mo_plus_3_a4*exp(-Mo_plus_3_b4*s*s)+Mo_plus_3_c
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_sf(s)=Mo_plus_5_a1*exp(-Mo_plus_5_b1*s*s)+Mo_plus_5_a2*exp(-Mo_plus_5_b2*s*s)+Mo_plus_5_a3*exp(-Mo_plus_5_b3*s*s)+Mo_plus_5_a4*exp(-Mo_plus_5_b4*s*s)+Mo_plus_5_c
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_sf(s)=Mo_plus_6_a1*exp(-Mo_plus_6_b1*s*s)+Mo_plus_6_a2*exp(-Mo_plus_6_b2*s*s)+Mo_plus_6_a3*exp(-Mo_plus_6_b3*s*s)+Mo_plus_6_a4*exp(-Mo_plus_6_b4*s*s)+Mo_plus_6_c
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_sf(s)=Tc_a1*exp(-Tc_b1*s*s)+Tc_a2*exp(-Tc_b2*s*s)+Tc_a3*exp(-Tc_b3*s*s)+Tc_a4*exp(-Tc_b4*s*s)+Tc_c
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_sf(s)=Ru_a1*exp(-Ru_b1*s*s)+Ru_a2*exp(-Ru_b2*s*s)+Ru_a3*exp(-Ru_b3*s*s)+Ru_a4*exp(-Ru_b4*s*s)+Ru_c
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_sf(s)=Ru_plus_3_a1*exp(-Ru_plus_3_b1*s*s)+Ru_plus_3_a2*exp(-Ru_plus_3_b2*s*s)+Ru_plus_3_a3*exp(-Ru_plus_3_b3*s*s)+Ru_plus_3_a4*exp(-Ru_plus_3_b4*s*s)+Ru_plus_3_c
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_sf(s)=Ru_plus_4_a1*exp(-Ru_plus_4_b1*s*s)+Ru_plus_4_a2*exp(-Ru_plus_4_b2*s*s)+Ru_plus_4_a3*exp(-Ru_plus_4_b3*s*s)+Ru_plus_4_a4*exp(-Ru_plus_4_b4*s*s)+Ru_plus_4_c
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_sf(s)=Rh_a1*exp(-Rh_b1*s*s)+Rh_a2*exp(-Rh_b2*s*s)+Rh_a3*exp(-Rh_b3*s*s)+Rh_a4*exp(-Rh_b4*s*s)+Rh_c
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_sf(s)=Rh_plus_3_a1*exp(-Rh_plus_3_b1*s*s)+Rh_plus_3_a2*exp(-Rh_plus_3_b2*s*s)+Rh_plus_3_a3*exp(-Rh_plus_3_b3*s*s)+Rh_plus_3_a4*exp(-Rh_plus_3_b4*s*s)+Rh_plus_3_c
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_sf(s)=Rh_plus_4_a1*exp(-Rh_plus_4_b1*s*s)+Rh_plus_4_a2*exp(-Rh_plus_4_b2*s*s)+Rh_plus_4_a3*exp(-Rh_plus_4_b3*s*s)+Rh_plus_4_a4*exp(-Rh_plus_4_b4*s*s)+Rh_plus_4_c
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_sf(s)=Pd_a1*exp(-Pd_b1*s*s)+Pd_a2*exp(-Pd_b2*s*s)+Pd_a3*exp(-Pd_b3*s*s)+Pd_a4*exp(-Pd_b4*s*s)+Pd_c
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_sf(s)=Pd_plus_2_a1*exp(-Pd_plus_2_b1*s*s)+Pd_plus_2_a2*exp(-Pd_plus_2_b2*s*s)+Pd_plus_2_a3*exp(-Pd_plus_2_b3*s*s)+Pd_plus_2_a4*exp(-Pd_plus_2_b4*s*s)+Pd_plus_2_c
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_sf(s)=Pd_plus_4_a1*exp(-Pd_plus_4_b1*s*s)+Pd_plus_4_a2*exp(-Pd_plus_4_b2*s*s)+Pd_plus_4_a3*exp(-Pd_plus_4_b3*s*s)+Pd_plus_4_a4*exp(-Pd_plus_4_b4*s*s)+Pd_plus_4_c
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_sf(s)=Ag_a1*exp(-Ag_b1*s*s)+Ag_a2*exp(-Ag_b2*s*s)+Ag_a3*exp(-Ag_b3*s*s)+Ag_a4*exp(-Ag_b4*s*s)+Ag_c
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_sf(s)=Ag_plus_1_a1*exp(-Ag_plus_1_b1*s*s)+Ag_plus_1_a2*exp(-Ag_plus_1_b2*s*s)+Ag_plus_1_a3*exp(-Ag_plus_1_b3*s*s)+Ag_plus_1_a4*exp(-Ag_plus_1_b4*s*s)+Ag_plus_1_c
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_sf(s)=Ag_plus_2_a1*exp(-Ag_plus_2_b1*s*s)+Ag_plus_2_a2*exp(-Ag_plus_2_b2*s*s)+Ag_plus_2_a3*exp(-Ag_plus_2_b3*s*s)+Ag_plus_2_a4*exp(-Ag_plus_2_b4*s*s)+Ag_plus_2_c
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_sf(s)=Cd_a1*exp(-Cd_b1*s*s)+Cd_a2*exp(-Cd_b2*s*s)+Cd_a3*exp(-Cd_b3*s*s)+Cd_a4*exp(-Cd_b4*s*s)+Cd_c
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_sf(s)=Cd_plus_2_a1*exp(-Cd_plus_2_b1*s*s)+Cd_plus_2_a2*exp(-Cd_plus_2_b2*s*s)+Cd_plus_2_a3*exp(-Cd_plus_2_b3*s*s)+Cd_plus_2_a4*exp(-Cd_plus_2_b4*s*s)+Cd_plus_2_c
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_sf(s)=In_a1*exp(-In_b1*s*s)+In_a2*exp(-In_b2*s*s)+In_a3*exp(-In_b3*s*s)+In_a4*exp(-In_b4*s*s)+In_c
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_sf(s)=In_plus_3_a1*exp(-In_plus_3_b1*s*s)+In_plus_3_a2*exp(-In_plus_3_b2*s*s)+In_plus_3_a3*exp(-In_plus_3_b3*s*s)+In_plus_3_a4*exp(-In_plus_3_b4*s*s)+In_plus_3_c
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_sf(s)=Sn_a1*exp(-Sn_b1*s*s)+Sn_a2*exp(-Sn_b2*s*s)+Sn_a3*exp(-Sn_b3*s*s)+Sn_a4*exp(-Sn_b4*s*s)+Sn_c
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_sf(s)=Sn_plus_2_a1*exp(-Sn_plus_2_b1*s*s)+Sn_plus_2_a2*exp(-Sn_plus_2_b2*s*s)+Sn_plus_2_a3*exp(-Sn_plus_2_b3*s*s)+Sn_plus_2_a4*exp(-Sn_plus_2_b4*s*s)+Sn_plus_2_c
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_sf(s)=Sn_plus_4_a1*exp(-Sn_plus_4_b1*s*s)+Sn_plus_4_a2*exp(-Sn_plus_4_b2*s*s)+Sn_plus_4_a3*exp(-Sn_plus_4_b3*s*s)+Sn_plus_4_a4*exp(-Sn_plus_4_b4*s*s)+Sn_plus_4_c
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_sf(s)=Sb_a1*exp(-Sb_b1*s*s)+Sb_a2*exp(-Sb_b2*s*s)+Sb_a3*exp(-Sb_b3*s*s)+Sb_a4*exp(-Sb_b4*s*s)+Sb_c
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_sf(s)=Sb_plus_3_a1*exp(-Sb_plus_3_b1*s*s)+Sb_plus_3_a2*exp(-Sb_plus_3_b2*s*s)+Sb_plus_3_a3*exp(-Sb_plus_3_b3*s*s)+Sb_plus_3_a4*exp(-Sb_plus_3_b4*s*s)+Sb_plus_3_c
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_sf(s)=Sb_plus_5_a1*exp(-Sb_plus_5_b1*s*s)+Sb_plus_5_a2*exp(-Sb_plus_5_b2*s*s)+Sb_plus_5_a3*exp(-Sb_plus_5_b3*s*s)+Sb_plus_5_a4*exp(-Sb_plus_5_b4*s*s)+Sb_plus_5_c
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_sf(s)=Te_a1*exp(-Te_b1*s*s)+Te_a2*exp(-Te_b2*s*s)+Te_a3*exp(-Te_b3*s*s)+Te_a4*exp(-Te_b4*s*s)+Te_c
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_sf(s)=I_a1*exp(-I_b1*s*s)+I_a2*exp(-I_b2*s*s)+I_a3*exp(-I_b3*s*s)+I_a4*exp(-I_b4*s*s)+I_c
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_sf(s)=I_minus_1_a1*exp(-I_minus_1_b1*s*s)+I_minus_1_a2*exp(-I_minus_1_b2*s*s)+I_minus_1_a3*exp(-I_minus_1_b3*s*s)+I_minus_1_a4*exp(-I_minus_1_b4*s*s)+I_minus_1_c
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_sf(s)=Xe_a1*exp(-Xe_b1*s*s)+Xe_a2*exp(-Xe_b2*s*s)+Xe_a3*exp(-Xe_b3*s*s)+Xe_a4*exp(-Xe_b4*s*s)+Xe_c
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_sf(s)=Cs_a1*exp(-Cs_b1*s*s)+Cs_a2*exp(-Cs_b2*s*s)+Cs_a3*exp(-Cs_b3*s*s)+Cs_a4*exp(-Cs_b4*s*s)+Cs_c
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_sf(s)=Cs_plus_1_a1*exp(-Cs_plus_1_b1*s*s)+Cs_plus_1_a2*exp(-Cs_plus_1_b2*s*s)+Cs_plus_1_a3*exp(-Cs_plus_1_b3*s*s)+Cs_plus_1_a4*exp(-Cs_plus_1_b4*s*s)+Cs_plus_1_c
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_sf(s)=Ba_a1*exp(-Ba_b1*s*s)+Ba_a2*exp(-Ba_b2*s*s)+Ba_a3*exp(-Ba_b3*s*s)+Ba_a4*exp(-Ba_b4*s*s)+Ba_c
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_sf(s)=Ba_plus_2_a1*exp(-Ba_plus_2_b1*s*s)+Ba_plus_2_a2*exp(-Ba_plus_2_b2*s*s)+Ba_plus_2_a3*exp(-Ba_plus_2_b3*s*s)+Ba_plus_2_a4*exp(-Ba_plus_2_b4*s*s)+Ba_plus_2_c
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_sf(s)=La_a1*exp(-La_b1*s*s)+La_a2*exp(-La_b2*s*s)+La_a3*exp(-La_b3*s*s)+La_a4*exp(-La_b4*s*s)+La_c
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_sf(s)=La_plus_3_a1*exp(-La_plus_3_b1*s*s)+La_plus_3_a2*exp(-La_plus_3_b2*s*s)+La_plus_3_a3*exp(-La_plus_3_b3*s*s)+La_plus_3_a4*exp(-La_plus_3_b4*s*s)+La_plus_3_c
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_sf(s)=Ce_a1*exp(-Ce_b1*s*s)+Ce_a2*exp(-Ce_b2*s*s)+Ce_a3*exp(-Ce_b3*s*s)+Ce_a4*exp(-Ce_b4*s*s)+Ce_c
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_sf(s)=Ce_plus_3_a1*exp(-Ce_plus_3_b1*s*s)+Ce_plus_3_a2*exp(-Ce_plus_3_b2*s*s)+Ce_plus_3_a3*exp(-Ce_plus_3_b3*s*s)+Ce_plus_3_a4*exp(-Ce_plus_3_b4*s*s)+Ce_plus_3_c
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_sf(s)=Ce_plus_4_a1*exp(-Ce_plus_4_b1*s*s)+Ce_plus_4_a2*exp(-Ce_plus_4_b2*s*s)+Ce_plus_4_a3*exp(-Ce_plus_4_b3*s*s)+Ce_plus_4_a4*exp(-Ce_plus_4_b4*s*s)+Ce_plus_4_c
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_sf(s)=Pr_a1*exp(-Pr_b1*s*s)+Pr_a2*exp(-Pr_b2*s*s)+Pr_a3*exp(-Pr_b3*s*s)+Pr_a4*exp(-Pr_b4*s*s)+Pr_c
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_sf(s)=Pr_plus_3_a1*exp(-Pr_plus_3_b1*s*s)+Pr_plus_3_a2*exp(-Pr_plus_3_b2*s*s)+Pr_plus_3_a3*exp(-Pr_plus_3_b3*s*s)+Pr_plus_3_a4*exp(-Pr_plus_3_b4*s*s)+Pr_plus_3_c
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_sf(s)=Pr_plus_4_a1*exp(-Pr_plus_4_b1*s*s)+Pr_plus_4_a2*exp(-Pr_plus_4_b2*s*s)+Pr_plus_4_a3*exp(-Pr_plus_4_b3*s*s)+Pr_plus_4_a4*exp(-Pr_plus_4_b4*s*s)+Pr_plus_4_c
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_sf(s)=Nd_a1*exp(-Nd_b1*s*s)+Nd_a2*exp(-Nd_b2*s*s)+Nd_a3*exp(-Nd_b3*s*s)+Nd_a4*exp(-Nd_b4*s*s)+Nd_c
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_sf(s)=Nd_plus_3_a1*exp(-Nd_plus_3_b1*s*s)+Nd_plus_3_a2*exp(-Nd_plus_3_b2*s*s)+Nd_plus_3_a3*exp(-Nd_plus_3_b3*s*s)+Nd_plus_3_a4*exp(-Nd_plus_3_b4*s*s)+Nd_plus_3_c
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_sf(s)=Pm_a1*exp(-Pm_b1*s*s)+Pm_a2*exp(-Pm_b2*s*s)+Pm_a3*exp(-Pm_b3*s*s)+Pm_a4*exp(-Pm_b4*s*s)+Pm_c
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_sf(s)=Pm_plus_3_a1*exp(-Pm_plus_3_b1*s*s)+Pm_plus_3_a2*exp(-Pm_plus_3_b2*s*s)+Pm_plus_3_a3*exp(-Pm_plus_3_b3*s*s)+Pm_plus_3_a4*exp(-Pm_plus_3_b4*s*s)+Pm_plus_3_c
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_sf(s)=Sm_a1*exp(-Sm_b1*s*s)+Sm_a2*exp(-Sm_b2*s*s)+Sm_a3*exp(-Sm_b3*s*s)+Sm_a4*exp(-Sm_b4*s*s)+Sm_c
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_sf(s)=Sm_plus_3_a1*exp(-Sm_plus_3_b1*s*s)+Sm_plus_3_a2*exp(-Sm_plus_3_b2*s*s)+Sm_plus_3_a3*exp(-Sm_plus_3_b3*s*s)+Sm_plus_3_a4*exp(-Sm_plus_3_b4*s*s)+Sm_plus_3_c
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_sf(s)=Eu_a1*exp(-Eu_b1*s*s)+Eu_a2*exp(-Eu_b2*s*s)+Eu_a3*exp(-Eu_b3*s*s)+Eu_a4*exp(-Eu_b4*s*s)+Eu_c
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_sf(s)=Eu_plus_2_a1*exp(-Eu_plus_2_b1*s*s)+Eu_plus_2_a2*exp(-Eu_plus_2_b2*s*s)+Eu_plus_2_a3*exp(-Eu_plus_2_b3*s*s)+Eu_plus_2_a4*exp(-Eu_plus_2_b4*s*s)+Eu_plus_2_c
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_sf(s)=Eu_plus_3_a1*exp(-Eu_plus_3_b1*s*s)+Eu_plus_3_a2*exp(-Eu_plus_3_b2*s*s)+Eu_plus_3_a3*exp(-Eu_plus_3_b3*s*s)+Eu_plus_3_a4*exp(-Eu_plus_3_b4*s*s)+Eu_plus_3_c
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_sf(s)=Gd_a1*exp(-Gd_b1*s*s)+Gd_a2*exp(-Gd_b2*s*s)+Gd_a3*exp(-Gd_b3*s*s)+Gd_a4*exp(-Gd_b4*s*s)+Gd_c
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_sf(s)=Gd_plus_3_a1*exp(-Gd_plus_3_b1*s*s)+Gd_plus_3_a2*exp(-Gd_plus_3_b2*s*s)+Gd_plus_3_a3*exp(-Gd_plus_3_b3*s*s)+Gd_plus_3_a4*exp(-Gd_plus_3_b4*s*s)+Gd_plus_3_c
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_sf(s)=Tb_a1*exp(-Tb_b1*s*s)+Tb_a2*exp(-Tb_b2*s*s)+Tb_a3*exp(-Tb_b3*s*s)+Tb_a4*exp(-Tb_b4*s*s)+Tb_c
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_sf(s)=Tb_plus_3_a1*exp(-Tb_plus_3_b1*s*s)+Tb_plus_3_a2*exp(-Tb_plus_3_b2*s*s)+Tb_plus_3_a3*exp(-Tb_plus_3_b3*s*s)+Tb_plus_3_a4*exp(-Tb_plus_3_b4*s*s)+Tb_plus_3_c
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_sf(s)=Dy_a1*exp(-Dy_b1*s*s)+Dy_a2*exp(-Dy_b2*s*s)+Dy_a3*exp(-Dy_b3*s*s)+Dy_a4*exp(-Dy_b4*s*s)+Dy_c
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_sf(s)=Dy_plus_3_a1*exp(-Dy_plus_3_b1*s*s)+Dy_plus_3_a2*exp(-Dy_plus_3_b2*s*s)+Dy_plus_3_a3*exp(-Dy_plus_3_b3*s*s)+Dy_plus_3_a4*exp(-Dy_plus_3_b4*s*s)+Dy_plus_3_c
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_sf(s)=Ho_a1*exp(-Ho_b1*s*s)+Ho_a2*exp(-Ho_b2*s*s)+Ho_a3*exp(-Ho_b3*s*s)+Ho_a4*exp(-Ho_b4*s*s)+Ho_c
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_sf(s)=Ho_plus_3_a1*exp(-Ho_plus_3_b1*s*s)+Ho_plus_3_a2*exp(-Ho_plus_3_b2*s*s)+Ho_plus_3_a3*exp(-Ho_plus_3_b3*s*s)+Ho_plus_3_a4*exp(-Ho_plus_3_b4*s*s)+Ho_plus_3_c
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_sf(s)=Er_a1*exp(-Er_b1*s*s)+Er_a2*exp(-Er_b2*s*s)+Er_a3*exp(-Er_b3*s*s)+Er_a4*exp(-Er_b4*s*s)+Er_c
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_sf(s)=Er_plus_3_a1*exp(-Er_plus_3_b1*s*s)+Er_plus_3_a2*exp(-Er_plus_3_b2*s*s)+Er_plus_3_a3*exp(-Er_plus_3_b3*s*s)+Er_plus_3_a4*exp(-Er_plus_3_b4*s*s)+Er_plus_3_c
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_sf(s)=Tm_a1*exp(-Tm_b1*s*s)+Tm_a2*exp(-Tm_b2*s*s)+Tm_a3*exp(-Tm_b3*s*s)+Tm_a4*exp(-Tm_b4*s*s)+Tm_c
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_sf(s)=Tm_plus_3_a1*exp(-Tm_plus_3_b1*s*s)+Tm_plus_3_a2*exp(-Tm_plus_3_b2*s*s)+Tm_plus_3_a3*exp(-Tm_plus_3_b3*s*s)+Tm_plus_3_a4*exp(-Tm_plus_3_b4*s*s)+Tm_plus_3_c
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_sf(s)=Yb_a1*exp(-Yb_b1*s*s)+Yb_a2*exp(-Yb_b2*s*s)+Yb_a3*exp(-Yb_b3*s*s)+Yb_a4*exp(-Yb_b4*s*s)+Yb_c
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_sf(s)=Yb_plus_2_a1*exp(-Yb_plus_2_b1*s*s)+Yb_plus_2_a2*exp(-Yb_plus_2_b2*s*s)+Yb_plus_2_a3*exp(-Yb_plus_2_b3*s*s)+Yb_plus_2_a4*exp(-Yb_plus_2_b4*s*s)+Yb_plus_2_c
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_sf(s)=Yb_plus_3_a1*exp(-Yb_plus_3_b1*s*s)+Yb_plus_3_a2*exp(-Yb_plus_3_b2*s*s)+Yb_plus_3_a3*exp(-Yb_plus_3_b3*s*s)+Yb_plus_3_a4*exp(-Yb_plus_3_b4*s*s)+Yb_plus_3_c
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_sf(s)=Lu_a1*exp(-Lu_b1*s*s)+Lu_a2*exp(-Lu_b2*s*s)+Lu_a3*exp(-Lu_b3*s*s)+Lu_a4*exp(-Lu_b4*s*s)+Lu_c
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_sf(s)=Lu_plus_3_a1*exp(-Lu_plus_3_b1*s*s)+Lu_plus_3_a2*exp(-Lu_plus_3_b2*s*s)+Lu_plus_3_a3*exp(-Lu_plus_3_b3*s*s)+Lu_plus_3_a4*exp(-Lu_plus_3_b4*s*s)+Lu_plus_3_c
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_sf(s)=Hf_a1*exp(-Hf_b1*s*s)+Hf_a2*exp(-Hf_b2*s*s)+Hf_a3*exp(-Hf_b3*s*s)+Hf_a4*exp(-Hf_b4*s*s)+Hf_c
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_sf(s)=Hf_plus_4_a1*exp(-Hf_plus_4_b1*s*s)+Hf_plus_4_a2*exp(-Hf_plus_4_b2*s*s)+Hf_plus_4_a3*exp(-Hf_plus_4_b3*s*s)+Hf_plus_4_a4*exp(-Hf_plus_4_b4*s*s)+Hf_plus_4_c
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_sf(s)=Ta_a1*exp(-Ta_b1*s*s)+Ta_a2*exp(-Ta_b2*s*s)+Ta_a3*exp(-Ta_b3*s*s)+Ta_a4*exp(-Ta_b4*s*s)+Ta_c
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_sf(s)=Ta_plus_5_a1*exp(-Ta_plus_5_b1*s*s)+Ta_plus_5_a2*exp(-Ta_plus_5_b2*s*s)+Ta_plus_5_a3*exp(-Ta_plus_5_b3*s*s)+Ta_plus_5_a4*exp(-Ta_plus_5_b4*s*s)+Ta_plus_5_c
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_sf(s)=W_a1*exp(-W_b1*s*s)+W_a2*exp(-W_b2*s*s)+W_a3*exp(-W_b3*s*s)+W_a4*exp(-W_b4*s*s)+W_c
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_sf(s)=W_plus_6_a1*exp(-W_plus_6_b1*s*s)+W_plus_6_a2*exp(-W_plus_6_b2*s*s)+W_plus_6_a3*exp(-W_plus_6_b3*s*s)+W_plus_6_a4*exp(-W_plus_6_b4*s*s)+W_plus_6_c
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_sf(s)=Re_a1*exp(-Re_b1*s*s)+Re_a2*exp(-Re_b2*s*s)+Re_a3*exp(-Re_b3*s*s)+Re_a4*exp(-Re_b4*s*s)+Re_c
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_sf(s)=Os_a1*exp(-Os_b1*s*s)+Os_a2*exp(-Os_b2*s*s)+Os_a3*exp(-Os_b3*s*s)+Os_a4*exp(-Os_b4*s*s)+Os_c
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_sf(s)=Os_plus_4_a1*exp(-Os_plus_4_b1*s*s)+Os_plus_4_a2*exp(-Os_plus_4_b2*s*s)+Os_plus_4_a3*exp(-Os_plus_4_b3*s*s)+Os_plus_4_a4*exp(-Os_plus_4_b4*s*s)+Os_plus_4_c
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_sf(s)=Ir_a1*exp(-Ir_b1*s*s)+Ir_a2*exp(-Ir_b2*s*s)+Ir_a3*exp(-Ir_b3*s*s)+Ir_a4*exp(-Ir_b4*s*s)+Ir_c
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_sf(s)=Ir_plus_3_a1*exp(-Ir_plus_3_b1*s*s)+Ir_plus_3_a2*exp(-Ir_plus_3_b2*s*s)+Ir_plus_3_a3*exp(-Ir_plus_3_b3*s*s)+Ir_plus_3_a4*exp(-Ir_plus_3_b4*s*s)+Ir_plus_3_c
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_sf(s)=Ir_plus_4_a1*exp(-Ir_plus_4_b1*s*s)+Ir_plus_4_a2*exp(-Ir_plus_4_b2*s*s)+Ir_plus_4_a3*exp(-Ir_plus_4_b3*s*s)+Ir_plus_4_a4*exp(-Ir_plus_4_b4*s*s)+Ir_plus_4_c
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_sf(s)=Pt_a1*exp(-Pt_b1*s*s)+Pt_a2*exp(-Pt_b2*s*s)+Pt_a3*exp(-Pt_b3*s*s)+Pt_a4*exp(-Pt_b4*s*s)+Pt_c
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_sf(s)=Pt_plus_2_a1*exp(-Pt_plus_2_b1*s*s)+Pt_plus_2_a2*exp(-Pt_plus_2_b2*s*s)+Pt_plus_2_a3*exp(-Pt_plus_2_b3*s*s)+Pt_plus_2_a4*exp(-Pt_plus_2_b4*s*s)+Pt_plus_2_c
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_sf(s)=Pt_plus_4_a1*exp(-Pt_plus_4_b1*s*s)+Pt_plus_4_a2*exp(-Pt_plus_4_b2*s*s)+Pt_plus_4_a3*exp(-Pt_plus_4_b3*s*s)+Pt_plus_4_a4*exp(-Pt_plus_4_b4*s*s)+Pt_plus_4_c
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_sf(s)=Au_a1*exp(-Au_b1*s*s)+Au_a2*exp(-Au_b2*s*s)+Au_a3*exp(-Au_b3*s*s)+Au_a4*exp(-Au_b4*s*s)+Au_c
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_sf(s)=Au_plus_1_a1*exp(-Au_plus_1_b1*s*s)+Au_plus_1_a2*exp(-Au_plus_1_b2*s*s)+Au_plus_1_a3*exp(-Au_plus_1_b3*s*s)+Au_plus_1_a4*exp(-Au_plus_1_b4*s*s)+Au_plus_1_c
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_sf(s)=Au_plus_3_a1*exp(-Au_plus_3_b1*s*s)+Au_plus_3_a2*exp(-Au_plus_3_b2*s*s)+Au_plus_3_a3*exp(-Au_plus_3_b3*s*s)+Au_plus_3_a4*exp(-Au_plus_3_b4*s*s)+Au_plus_3_c
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_sf(s)=Hg_a1*exp(-Hg_b1*s*s)+Hg_a2*exp(-Hg_b2*s*s)+Hg_a3*exp(-Hg_b3*s*s)+Hg_a4*exp(-Hg_b4*s*s)+Hg_c
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_sf(s)=Hg_plus_1_a1*exp(-Hg_plus_1_b1*s*s)+Hg_plus_1_a2*exp(-Hg_plus_1_b2*s*s)+Hg_plus_1_a3*exp(-Hg_plus_1_b3*s*s)+Hg_plus_1_a4*exp(-Hg_plus_1_b4*s*s)+Hg_plus_1_c
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_sf(s)=Hg_plus_2_a1*exp(-Hg_plus_2_b1*s*s)+Hg_plus_2_a2*exp(-Hg_plus_2_b2*s*s)+Hg_plus_2_a3*exp(-Hg_plus_2_b3*s*s)+Hg_plus_2_a4*exp(-Hg_plus_2_b4*s*s)+Hg_plus_2_c
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_sf(s)=Tl_a1*exp(-Tl_b1*s*s)+Tl_a2*exp(-Tl_b2*s*s)+Tl_a3*exp(-Tl_b3*s*s)+Tl_a4*exp(-Tl_b4*s*s)+Tl_c
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_sf(s)=Tl_plus_1_a1*exp(-Tl_plus_1_b1*s*s)+Tl_plus_1_a2*exp(-Tl_plus_1_b2*s*s)+Tl_plus_1_a3*exp(-Tl_plus_1_b3*s*s)+Tl_plus_1_a4*exp(-Tl_plus_1_b4*s*s)+Tl_plus_1_c
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_sf(s)=Tl_plus_3_a1*exp(-Tl_plus_3_b1*s*s)+Tl_plus_3_a2*exp(-Tl_plus_3_b2*s*s)+Tl_plus_3_a3*exp(-Tl_plus_3_b3*s*s)+Tl_plus_3_a4*exp(-Tl_plus_3_b4*s*s)+Tl_plus_3_c
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_sf(s)=Pb_a1*exp(-Pb_b1*s*s)+Pb_a2*exp(-Pb_b2*s*s)+Pb_a3*exp(-Pb_b3*s*s)+Pb_a4*exp(-Pb_b4*s*s)+Pb_c
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_sf(s)=Pb_plus_2_a1*exp(-Pb_plus_2_b1*s*s)+Pb_plus_2_a2*exp(-Pb_plus_2_b2*s*s)+Pb_plus_2_a3*exp(-Pb_plus_2_b3*s*s)+Pb_plus_2_a4*exp(-Pb_plus_2_b4*s*s)+Pb_plus_2_c
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_sf(s)=Pb_plus_4_a1*exp(-Pb_plus_4_b1*s*s)+Pb_plus_4_a2*exp(-Pb_plus_4_b2*s*s)+Pb_plus_4_a3*exp(-Pb_plus_4_b3*s*s)+Pb_plus_4_a4*exp(-Pb_plus_4_b4*s*s)+Pb_plus_4_c
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_sf(s)=Bi_a1*exp(-Bi_b1*s*s)+Bi_a2*exp(-Bi_b2*s*s)+Bi_a3*exp(-Bi_b3*s*s)+Bi_a4*exp(-Bi_b4*s*s)+Bi_c
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_sf(s)=Bi_plus_3_a1*exp(-Bi_plus_3_b1*s*s)+Bi_plus_3_a2*exp(-Bi_plus_3_b2*s*s)+Bi_plus_3_a3*exp(-Bi_plus_3_b3*s*s)+Bi_plus_3_a4*exp(-Bi_plus_3_b4*s*s)+Bi_plus_3_c
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_sf(s)=Bi_plus_5_a1*exp(-Bi_plus_5_b1*s*s)+Bi_plus_5_a2*exp(-Bi_plus_5_b2*s*s)+Bi_plus_5_a3*exp(-Bi_plus_5_b3*s*s)+Bi_plus_5_a4*exp(-Bi_plus_5_b4*s*s)+Bi_plus_5_c
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_sf(s)=Po_a1*exp(-Po_b1*s*s)+Po_a2*exp(-Po_b2*s*s)+Po_a3*exp(-Po_b3*s*s)+Po_a4*exp(-Po_b4*s*s)+Po_c
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_sf(s)=At_a1*exp(-At_b1*s*s)+At_a2*exp(-At_b2*s*s)+At_a3*exp(-At_b3*s*s)+At_a4*exp(-At_b4*s*s)+At_c
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_sf(s)=Rn_a1*exp(-Rn_b1*s*s)+Rn_a2*exp(-Rn_b2*s*s)+Rn_a3*exp(-Rn_b3*s*s)+Rn_a4*exp(-Rn_b4*s*s)+Rn_c
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_sf(s)=Fr_a1*exp(-Fr_b1*s*s)+Fr_a2*exp(-Fr_b2*s*s)+Fr_a3*exp(-Fr_b3*s*s)+Fr_a4*exp(-Fr_b4*s*s)+Fr_c
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_sf(s)=Ra_a1*exp(-Ra_b1*s*s)+Ra_a2*exp(-Ra_b2*s*s)+Ra_a3*exp(-Ra_b3*s*s)+Ra_a4*exp(-Ra_b4*s*s)+Ra_c
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_sf(s)=Ra_plus_2_a1*exp(-Ra_plus_2_b1*s*s)+Ra_plus_2_a2*exp(-Ra_plus_2_b2*s*s)+Ra_plus_2_a3*exp(-Ra_plus_2_b3*s*s)+Ra_plus_2_a4*exp(-Ra_plus_2_b4*s*s)+Ra_plus_2_c
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_sf(s)=Ac_a1*exp(-Ac_b1*s*s)+Ac_a2*exp(-Ac_b2*s*s)+Ac_a3*exp(-Ac_b3*s*s)+Ac_a4*exp(-Ac_b4*s*s)+Ac_c
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_sf(s)=Ac_plus_3_a1*exp(-Ac_plus_3_b1*s*s)+Ac_plus_3_a2*exp(-Ac_plus_3_b2*s*s)+Ac_plus_3_a3*exp(-Ac_plus_3_b3*s*s)+Ac_plus_3_a4*exp(-Ac_plus_3_b4*s*s)+Ac_plus_3_c
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_sf(s)=Th_a1*exp(-Th_b1*s*s)+Th_a2*exp(-Th_b2*s*s)+Th_a3*exp(-Th_b3*s*s)+Th_a4*exp(-Th_b4*s*s)+Th_c
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_sf(s)=Th_plus_4_a1*exp(-Th_plus_4_b1*s*s)+Th_plus_4_a2*exp(-Th_plus_4_b2*s*s)+Th_plus_4_a3*exp(-Th_plus_4_b3*s*s)+Th_plus_4_a4*exp(-Th_plus_4_b4*s*s)+Th_plus_4_c
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_sf(s)=Pa_a1*exp(-Pa_b1*s*s)+Pa_a2*exp(-Pa_b2*s*s)+Pa_a3*exp(-Pa_b3*s*s)+Pa_a4*exp(-Pa_b4*s*s)+Pa_c
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_sf(s)=U_a1*exp(-U_b1*s*s)+U_a2*exp(-U_b2*s*s)+U_a3*exp(-U_b3*s*s)+U_a4*exp(-U_b4*s*s)+U_c
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_sf(s)=U_plus_3_a1*exp(-U_plus_3_b1*s*s)+U_plus_3_a2*exp(-U_plus_3_b2*s*s)+U_plus_3_a3*exp(-U_plus_3_b3*s*s)+U_plus_3_a4*exp(-U_plus_3_b4*s*s)+U_plus_3_c
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_sf(s)=U_plus_4_a1*exp(-U_plus_4_b1*s*s)+U_plus_4_a2*exp(-U_plus_4_b2*s*s)+U_plus_4_a3*exp(-U_plus_4_b3*s*s)+U_plus_4_a4*exp(-U_plus_4_b4*s*s)+U_plus_4_c
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_sf(s)=U_plus_6_a1*exp(-U_plus_6_b1*s*s)+U_plus_6_a2*exp(-U_plus_6_b2*s*s)+U_plus_6_a3*exp(-U_plus_6_b3*s*s)+U_plus_6_a4*exp(-U_plus_6_b4*s*s)+U_plus_6_c
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_sf(s)=Np_a1*exp(-Np_b1*s*s)+Np_a2*exp(-Np_b2*s*s)+Np_a3*exp(-Np_b3*s*s)+Np_a4*exp(-Np_b4*s*s)+Np_c
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_sf(s)=Np_plus_3_a1*exp(-Np_plus_3_b1*s*s)+Np_plus_3_a2*exp(-Np_plus_3_b2*s*s)+Np_plus_3_a3*exp(-Np_plus_3_b3*s*s)+Np_plus_3_a4*exp(-Np_plus_3_b4*s*s)+Np_plus_3_c
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_sf(s)=Np_plus_4_a1*exp(-Np_plus_4_b1*s*s)+Np_plus_4_a2*exp(-Np_plus_4_b2*s*s)+Np_plus_4_a3*exp(-Np_plus_4_b3*s*s)+Np_plus_4_a4*exp(-Np_plus_4_b4*s*s)+Np_plus_4_c
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_sf(s)=Np_plus_6_a1*exp(-Np_plus_6_b1*s*s)+Np_plus_6_a2*exp(-Np_plus_6_b2*s*s)+Np_plus_6_a3*exp(-Np_plus_6_b3*s*s)+Np_plus_6_a4*exp(-Np_plus_6_b4*s*s)+Np_plus_6_c
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_sf(s)=Pu_a1*exp(-Pu_b1*s*s)+Pu_a2*exp(-Pu_b2*s*s)+Pu_a3*exp(-Pu_b3*s*s)+Pu_a4*exp(-Pu_b4*s*s)+Pu_c
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_sf(s)=Pu_plus_3_a1*exp(-Pu_plus_3_b1*s*s)+Pu_plus_3_a2*exp(-Pu_plus_3_b2*s*s)+Pu_plus_3_a3*exp(-Pu_plus_3_b3*s*s)+Pu_plus_3_a4*exp(-Pu_plus_3_b4*s*s)+Pu_plus_3_c
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_sf(s)=Pu_plus_4_a1*exp(-Pu_plus_4_b1*s*s)+Pu_plus_4_a2*exp(-Pu_plus_4_b2*s*s)+Pu_plus_4_a3*exp(-Pu_plus_4_b3*s*s)+Pu_plus_4_a4*exp(-Pu_plus_4_b4*s*s)+Pu_plus_4_c
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_sf(s)=Pu_plus_6_a1*exp(-Pu_plus_6_b1*s*s)+Pu_plus_6_a2*exp(-Pu_plus_6_b2*s*s)+Pu_plus_6_a3*exp(-Pu_plus_6_b3*s*s)+Pu_plus_6_a4*exp(-Pu_plus_6_b4*s*s)+Pu_plus_6_c
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_sf(s)=Am_a1*exp(-Am_b1*s*s)+Am_a2*exp(-Am_b2*s*s)+Am_a3*exp(-Am_b3*s*s)+Am_a4*exp(-Am_b4*s*s)+Am_c
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_sf(s)=Cm_a1*exp(-Cm_b1*s*s)+Cm_a2*exp(-Cm_b2*s*s)+Cm_a3*exp(-Cm_b3*s*s)+Cm_a4*exp(-Cm_b4*s*s)+Cm_c
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_sf(s)=Bk_a1*exp(-Bk_b1*s*s)+Bk_a2*exp(-Bk_b2*s*s)+Bk_a3*exp(-Bk_b3*s*s)+Bk_a4*exp(-Bk_b4*s*s)+Bk_c
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_sf(s)=Cf_a1*exp(-Cf_b1*s*s)+Cf_a2*exp(-Cf_b2*s*s)+Cf_a3*exp(-Cf_b3*s*s)+Cf_a4*exp(-Cf_b4*s*s)+Cf_c