function sf = scattering_factor_parameters(a)
a=char(a);
if strcmp(a,'H')
sf=[0.493002 10.5109 0.322912 26.1257 0.140191 3.14236 0.040810 57.7997 0.003038];
elseif strcmp(a,'Li')%This is for Li(0)
sf=[1.1282 3.9546 0.7508 1.0524 0.6175 85.3905 0.4653 168.261 0.0377];
%elseif strcmp(a,'Li')%This is for Li+
%sf=[0.6968 4.6237 0.7888 1.9557 0.3414 0.6316 0.1563 10.0953 0.0167];
elseif strcmp(a,'C')
sf=[2.31 20.8439 1.02 10.2075 1.5886 0.5687 0.865 51.6512 0.2156];
elseif strcmp(a,'N')
sf=[12.2126 0.0057 3.1322 9.8933 2.0125 28.9975 1.1663 0.5826 -11.529];
elseif strcmp(a,'O')
sf=[3.04850 13.2771 2.2868 5.7011 1.5463 0.3239 0.867 32.9089 0.2508];
elseif strcmp(a,'B')
sf=[2.0545 23.2185 1.3326 1.021 1.0979 60.3498 0.7068 0.1403 -0.1932];
%elseif strcmp(a(b),'Na')%this is for Na+
% sf=[3.2565 2.6671 3.93620 6.1153 1.3998 0.2001 1.0032 14.0390 0.404];
elseif strcmp(a,'Na')
sf=[4.7626 3.285 3.1736 8.8422 1.2674 0.3136 1.1128 129.424 0.676];
elseif strcmp(a,'Si')
sf=[6.2915 2.4386 3.0353 32.3337 1.9891 0.6785 1.541 81.6937 1.1407];
elseif strcmp(a,'P')
sf=[6.4345 1.9067 4.1791 27.1570 1.7800 0.526 1.4908 68.1645 11.1149];
elseif strcmp(a,'S')
sf=[6.9053 1.46790 5.2034 22.2151 1.4379 0.2536 1.5863 56.172 0.8669];
%elseif strcmp(a,'Cl')%this is for Cl-
% sf=[18.2915 0.0066 7.2084 1.1717 6.5337 19.5424 2.3386 60.4486 -16.378];
elseif strcmp(a,'Cl')
sf=[11.4604 0.0104 7.1964 1.16620 6.25560 18.5194 1.6455 47.7784 -9.5574];
elseif strcmp(a,'Cr')
sf=[10.6406 6.1038 7.3537 0.392 3.324 20.2626 1.4922 98.7399 1.1832];
elseif strcmp(a,'Mn')
sf=[11.2819 5.3409 7.3573 0.3432 3.01930 17.8674 2.24410 83.7543 1.08960];
elseif strcmp(a,'Fe')
sf=[11.7695 4.7611 7.3573 0.3072 3.5222 15.3535 2.3045 76.8805 1.0369];
elseif strcmp(a,'Co')
sf=[12.2841 4.2791 7.3409 0.2784 4.0034 13.5359 2.2488 71.1692 1.0118];
elseif strcmp(a,'Ni')
sf=[12.8376 3.8785 7.292 0.2565 4.4438 12.1763 2.38 66.3421 1.0341];
elseif strcmp(a,'Cu')
sf=[13.338 3.5828 7.1676 0.247 5.6158 11.3966 1.6735 64.8126 1.191];
elseif strcmp(a,'Zn')
sf=[14.0743 3.2655 7.0318 0.2333 5.1652 10.3163 2.41 58.7097 1.3041];
elseif strcmp(a,'Se')
sf=[17.0006 2.4098 5.8196 0.2726 3.9731 15.2372 4.3543 43.8163 2.8409];
elseif strcmp(a,'Cs')
sf=[20.3892 3.569 19.1062 0.3107 10.662 24.3879 1.4953 213.904 3.3352];
%elseif strcmp(a,'Cs')%this is for Cs+
%sf=[20.3524 3.552 19.1278 0.3086 10.2821 23.7128 0.9615 59.4565 3.2791];
elseif strcmp(a,'Mo')
sf=[3.7025 0.2772 17.2356 1.0958 12.8876 11.004 3.7429 61.6584 4.3875];
end