function GARF_distribution = GARF_distgen(protein_atomtype_index_number,ligand_atomtype_index_number) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % This octave code is for generating a GARF potential-driven distribution % from 1 to 20 angstroms given a selected atom pair, using the index number % "protein_atomtype_index_number" and "ligand_atomtype_index_number". The % parameter set for all the 480 nonbonding atom pairwise interactions is % included in this file and assigned to the variable "Paraset". %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % protein atom type list, the index number from 1 to 20 shows all the % defined protein atom types in GARF. C=1 ; Cs=2 ; CA=3 ; CB=4 ; CC=5 ; CN=6 ; CR=7 ; CT=8 ; CV=9 ; CW=10; N=11; N2=12; N3=13; NA=14; NB=15; O=16 ; O2=17; OH=18; S=19; SH=20; % ligand atom type list, the index number from 1 to 24 shows all the % defined ligand atom types in GARF. C_3=1; C_2=2; C_1=3; C_ar=4; O_3=5; O_3p=6; O_2=7; O_co2=8; O_2v=9; N_2=10; N_am=11; N_pl3=12; N_4=13; P=14; F=15; Cl=16; Br=17; I=18; C_cat=19; S_3=20; S_o=21; C_3X=22; C_2X=23; C_arX=24; % Paraset_index_num is the row number locating the parameter set for % generating the GARF potential-driven distribution regarding a selected % protein atom type (using "protein_atomtype_index_number", from 1 to 20) % and a ligand atom type (using "ligand_atomtype_index_number", from 1 to % 24). Paraset_index_num=(protein_atomtype_index_number-1)*24+ligand_atomtype_index_number; Paraset=[1.846 3.5 9 -1.591 3.5 6 2.4723 3.3 12 -1.1634 3.3 6 51.496 2.3 8 -11.254 2.3 4.9 25.133 2.6 8 -9.1649 2.6 5.1 1.0565 3 12 -0.88233 3 6 0.85047 3.3 8 -0.53819 3.3 4 33.314 2.4 12 -4.9627 2.4 6 0.19654 3.6 11 -0.32283 3.6 6 6.4713 2.9 12 -2.439 2.9 6 1.625 3.2 10 -1.1666 3.2 6 0.62148 3.5 12 -0.63944 3.5 6 35.699 2.4 12 -5.318 2.4 6 0.64491 3.1 12 -0.74101 3.1 6 28.996 2.9 12 -4.1542 2.9 6 69.658 2.1 12 -6.5 2.1 6 189.73 2.1 12 -12.684 2.1 6 2.1619 3.5 12 -1.4481 3.5 6 2.1619 3.5 12 -1.4481 3.5 6 2.4955 3.4 12 -1.9693 3.4 6 0.070199 3.9 12 -0.22695 3.9 6 228.47 2.6 12 -11.541 2.6 6 3.5515 3.3 12 -2.2735 3.3 6 74.481 2.3 8 -19.507 2.3 5.3 1.9191 3.7 8 -1.9858 3.7 6 1.1846 3.4 11 -1.0937 3.4 6 447.12 2.4 12 -18.038 2.4 6 103.18 2.4 12 -8.3024 2.4 6 1.055 3.6 12 -0.8735 3.6 6 17.823 2.9 9 -7.6023 2.9 6 29.857 3 12 -5.2425 3 6 4.1091 3.5 12 -1.6803 3.5 6 0.78165 3.7 12 -0.97924 3.7 6 19312 2 12 -158.2 2 6 341.88 2 12 -20.103 2 6 0.42074 3.7 12 -0.45451 3.7 6 30.099 2.8 12 -6.107 2.8 6 50.732 2.6 12 -7.6522 2.6 6 9.273 3.3 12 -3.789 3.3 6 712.38 2.5 12 -32.27 2.5 6 68.511 2.4 12 -6.9998 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 104.16 3 12 -14.088 3 6 2.7914 3.4 12 -2.451 3.4 6 2.7914 3.4 12 -2.451 3.4 6 2.9516 3.4 10 -2.1126 3.4 6 28.912 3 12 -4.0119 3 6 3.5358 3 12 -1.7122 3 6 5.7767 3.4 8 -5.5649 3.4 6 14.51 3 8 -7.8838 3 5 103.18 2.4 12 -8.3024 2.4 6 5.5477 3 12 -3.4207 3 6 7.4897 2.3 8 -0.64565 2.3 2.3 39.764 2.4 12 -5.9235 2.4 6 78.713 2.1 12 -7.3448 2.1 6 0.20735 3.6 12 -0.37554 3.6 6 1.4512 3.3 8 -1.3854 3.3 6 3.7065 3.4 8 -3.5756 3.4 5.8 4.5682 3.1 8 -4.0536 3.1 6 33.495 2 8 -0.70109 2 2 142.41 2.1 12 -9.5207 2.1 6 9.273 3.3 12 -3.789 3.3 6 0.84836 2.9 12 -0.78147 2.9 6 68.511 2.4 12 -6.9998 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 2.7161 3 12 -1.8193 3 6 0.82267 3.4 12 -0.84417 3.4 6 369.15 2.4 12 -18.169 2.4 6 0.62125 3.6 12 -1.0043 3.6 6 4.4767 3.2 12 -3.5697 3.2 6 100.13 2.4 12 -16.185 2.4 6 0.43282 3.7 12 -0.79867 3.7 6 2.0185 3.4 12 -1.5225 3.4 6 103.18 2.4 12 -8.3024 2.4 6 6.2232 2.7 8 -1.0162 2.7 2.9 2.0697 3.5 12 -1.2038 3.5 6 1.0495 3.5 12 -0.94198 3.5 6 1.6068 3.4 12 -1.4109 3.4 6 107.41 3 12 -10.022 3 6 146.51 3 12 -13.671 3 6 36.806 2.4 12 -5.4829 2.4 6 88.552 2.4 12 -8.263 2.4 6 64.776 2.4 12 -9.6494 2.4 6 18.765 3 12 -3.7704 3 6 9.273 3.3 12 -3.789 3.3 6 2185 2.5 12 -60.632 2.5 6 66.197 3 12 -8.9538 3 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 2.2121 3.5 12 -1.9856 3.5 6 35.114 3 12 -7.0555 3 6 35.114 3 12 -7.0555 3 6 356.51 2.1 11 -26.068 2.1 6 8.7516 2.6 8 -2.022 2.6 3.8 244.73 2.1 12 -16.361 2.1 6 2.4345 3.1 12 -1.2281 3.1 6 3.9752 3.4 8 -3.6449 3.4 6 103.18 2.4 12 -8.3024 2.4 6 3.0864 3 8 -1.2388 3 4.1 1.998 3.3 12 -1.6096 3.3 6 92.478 2.1 12 -10.269 2.1 6 20.503 2.7 12 -3.2959 2.7 6 37.238 2.4 12 -5.5472 2.4 6 4.3177 3.2 12 -3.0795 3.2 6 0.73961 3.1 12 -0.71416 3.1 6 12.903 2.6 8 -3.87 2.6 4.7 57.214 2.4 12 -10.046 2.4 6 6572.7 2 12 -101.32 2 6 43.425 3 12 -8.7253 3 6 492.63 2.6 12 -24.885 2.6 6 68.511 2.4 12 -6.9998 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 3.076 3.2 12 -2.2645 3.2 6 34.942 3 12 -4.979 3 5.7 231.44 2.9 12 -22.511 2.9 6 304.5 2 12 -11.077 2 6 104.36 2.1 12 -6.9769 2.1 6 3.4124 3 8 -1.3277 3 3.9 0.079074 3.9 12 -0.25564 3.9 6 1.3776 3.5 12 -1.0662 3.5 6 47.182 3 12 -6.3819 3 6 5.3622 3 12 -1.5085 3 6 0.24492 3.7 12 -0.41742 3.7 6 124.91 2 8 -3.9646 2 2.9 5.7784 3 8 -3.1909 3 5.1 3.5294 3.2 12 -1.7693 3.2 6 1857.1 2.5 12 -51.532 2.5 6 40.884 2.4 8 -12.526 2.4 5.4 68.713 2.3 8 -13.751 2.3 4.7 0.44991 3.6 12 -0.67146 3.6 6 17.976 3 12 -4.7747 3 6 9.273 3.3 12 -3.789 3.3 6 116.08 2.6 12 -10.298 2.6 6 2.1141 3.6 8 -1.5905 3.6 4.4 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 122.87 3 12 -13.595 3 6 56.729 3 12 -8.73 3 6 56.729 3 12 -8.73 3 6 4.6984 3.1 12 -2.7265 3.1 6 311.51 2.1 12 -16.408 2.1 6 2.8012 3 12 -1.5918 3 6 5.6034 2.8 8 -1.8378 2.8 4 34.589 2.7 12 -9.1875 2.7 6 26.742 2.7 8 -5.0236 2.7 3.8 10.209 2.7 12 -3.0831 2.7 6 1.588 3 8 -1.3855 3 3.8 0.54189 3.7 12 -0.67888 3.7 6 7.946 2.4 12 -2.6253 2.4 6 0.71067 3.2 12 -1.0055 3.2 6 150.41 2.1 12 -16.701 2.1 6 132.08 2.2 12 -10.77 2.2 6 1.4036 3.5 8 -1.369 3.5 5.7 35.054 2.6 12 -6.3477 2.6 6 65.393 2.4 8 -10.783 2.4 4.3 28.045 3 12 -5.7913 3 6 18.922 3.3 12 -4.2019 3.3 6 103.97 2.4 12 -13.198 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 55.111 2.4 12 -8.2097 2.4 6 5.5071 3.2 12 -2.4897 3.2 6 1.6076 3.7 12 -2.0141 3.7 6 0.754 3.6 8 -0.83752 3.6 5.5 19.908 2.6 8 -7.9971 2.6 5.2 1.1988 3.2 12 -0.98729 3.2 6 0.64009 3.8 8 -0.86529 3.8 5.1 1.6028 3.4 12 -1.4073 3.4 6 18.859 2.7 12 -3.6764 2.7 6 12.02 2.6 12 -2.1029 2.6 3.8 93.006 2.1 12 -7.9705 2.1 6 0.2096 3.7 10 -0.33108 3.7 6 4.9479 2.5 12 -1.5639 2.5 6 0.6534 3.1 12 -0.63091 3.1 6 0.49862 3.1 12 -0.48146 3.1 6 1.6805 3.3 12 -1.6918 3.3 6 1.8635 3.2 8 -1.4891 3.2 5.3 46.835 2.2 10 -7.9041 2.2 6 1.5194 3.1 12 -0.93932 3.1 6 100.35 2.4 10 -14.805 2.4 6 0.37074 3.1 12 -0.35798 3.1 6 0.84039 3.4 12 -0.86235 3.4 6 1.6063 3.3 12 -0.98307 3.3 6 1.6063 3.3 12 -0.98307 3.3 6 86.555 2.4 12 -9.4016 2.4 6 1.8824 3.3 12 -1.4032 3.3 6 22.873 2.9 12 -3.7499 2.9 6 1.0715 3.7 8 -1.2155 3.7 5.4 1.0647 3.7 8 -1.2393 3.7 4.3 0.53005 3.6 12 -0.91213 3.6 5.9 0.29441 3.6 12 -0.42569 3.6 6 260.89 2.1 12 -17.442 2.1 6 103.18 2.4 12 -8.3024 2.4 6 2.0991 3 12 -1.406 3 6 0.6145 3.4 12 -0.98124 3.4 6 0.17678 3.6 12 -0.35357 3.6 6 18.413 2.4 12 -4.6369 2.4 6 2.9091 2.9 12 -2.0884 2.9 6 976.68 2.6 12 -38.619 2.6 6 35.32 2.4 12 -5.2614 2.4 6 65.099 2.1 12 -6.0745 2.1 6 1.749 3.3 12 -1.5003 3.3 6 7.6682 3.8 12 -3.7755 3.8 6 20.041 3.3 12 -4.8022 3.3 6 1.6935 2.9 12 -1.56 2.9 6 98.732 2.4 12 -13.57 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 343.87 2.4 12 -20.791 2.4 6 1049.6 2.5 12 -37.029 2.5 6 3639.4 2.2 12 -87.106 2.2 6 0.67646 3.7 9 -0.75064 3.7 5.9 86.766 2.3 10 -11.307 2.3 5.9 1.5509 3.4 8 -1.5716 3.4 6 73.591 2.3 12 -6.192 2.3 6 428.39 2.2 8 -95.317 2.2 5.7 103.18 2.4 12 -8.3024 2.4 6 1.0121 3.5 12 -0.88167 3.5 6 892.85 2.1 12 -30.417 2.1 6 230.24 2.5 12 -13.539 2.5 6 24.48 2.3 8 -5.5185 2.3 4.1 153.17 2.1 12 -14.293 2.1 6 153.17 2.1 12 -14.293 2.1 6 0.41016 3.7 12 -0.69902 3.7 6 91.271 2.3 8 -30.214 2.3 5.8 1.1565 3.5 12 -1.19 3.5 6 341.66 2 8 -73.566 2 5.6 9.273 3.3 12 -3.789 3.3 6 0.84836 2.9 12 -0.78147 2.9 6 68.511 2.4 12 -6.9998 2.4 6 1707 2 12 -36.559 2 6 1707 2 12 -36.559 2 6 2205.2 2 12 -47.229 2 6 783.12 2.3 12 -38.794 2.3 6 27.574 3 12 -6.7267 3 6 60.419 2.4 12 -6.5627 2.4 6 6.063 3 9 -3.008 3 5.9 65.655 2.1 8 -8.7666 2.1 4.5 0.61257 3.5 12 -0.9068 3.5 6 2.281 3 8 -1.0399 3 3.6 0.81601 3.4 12 -0.71649 3.4 6 3.6367 3 8 -1.6016 3 2.6 1.1462 2.7 12 -1.3204 2.7 6 1.0035 2.8 12 -1.0897 2.8 6 7.8056 2.3 12 -5.0065 2.3 6 1.0261 2.7 12 -1.7223 2.7 6 0.70608 2.8 12 -0.93347 2.8 6 0.027183 3.9 12 -0.29624 3.9 6 6.6943 2.4 12 -2.8829 2.4 6 0.42022 3.7 8 -0.52682 3.7 3.7 0.42454 3.7 12 -0.72354 3.7 6 5.2842 3 12 -2.3292 3 6 0.65952 2.9 12 -0.60752 2.9 6 3.5377 3 12 -1.7131 3 6 3.5238 3.4 12 -1.9834 3.4 6 3.5238 3.4 12 -1.9834 3.4 6 1.5779 3.5 12 -1.2213 3.5 6 0.52165 3.6 12 -0.85697 3.6 6 0.79725 3.7 12 -0.9988 3.7 6 0.70288 3.6 8 -0.8056 3.6 5.9 3.1038 3.1 12 -1.8298 3.1 6 3.0416 3 11 -1.374 3 6 7.3139 2.4 8 -2.2574 2.4 4.3 67.12 2.1 12 -5.2891 2.1 6 33.007 2.7 8 -6.5818 2.7 3.7 4.0041 2.8 8 -1.4296 2.8 3.2 0.81691 2.7 12 -0.86828 2.7 6 10.194 2.4 12 -4.3903 2.4 6 0.66206 2.7 12 -0.91795 2.7 6 1.1821 2.7 12 -1.9843 2.7 6 95.241 2.1 12 -10.575 2.1 6 20.299 2 8 -3.011 2 3.6 0.51123 3.3 11 -0.52895 3.3 6 1.0659 2.9 12 -0.82087 2.9 6 7.4455 3.7 11 -3.7148 3.7 6 15.678 2.8 12 -4.4287 2.8 6 495.16 2.1 15 -30.723 2.1 7.9 25.108 2.7 8 -5.3066 2.7 3.9 212.53 2 12 -7.7311 2 6 212.53 2 12 -7.7311 2 6 1.4669 3.7 8 -1.5933 3.7 5.1 64.376 2.4 12 -9.5898 2.4 6 1.5099 3.6 12 -1.3638 3.6 6 147.35 2.1 12 -9.8514 2.1 6 2.4098 3 12 -1.3694 3 6 98 2.1 12 -6.5518 2.1 6 0.84216 3.1 8 -0.62787 3.1 3 13.387 2.5 8 -2.7362 2.5 3.6 776.33 2 12 -24.26 2 6 181.31 2.1 12 -12.122 2.1 6 0.53165 2.8 12 -0.93628 2.8 6 0.6749 2.8 12 -0.84055 2.8 6 1.3855 2.4 12 -0.89023 2.4 6 2.2549 2.4 12 -1.867 2.4 6 10.786 2.3 8 -3.997 2.3 4.5 1.3475 2.7 12 -1.4322 2.7 6 5.8855 2.4 12 -2.5346 2.4 6 1.7873 3.5 12 -1.6043 3.5 6 1.278 3.3 12 -1.7354 3.3 6 2.4654 3.3 12 -2.183 3.3 6 52.631 3 12 -9.8818 3 6 2259.4 2.4 12 -58.702 2.4 6 2259.4 2.4 12 -58.702 2.4 6 2259.4 2.4 12 -58.702 2.4 6 1863.8 2.5 12 -51.718 2.5 6 54.107 3 12 -7.3186 3 6 213.02 2.6 12 -18.138 2.6 6 147.35 2.1 12 -9.8514 2.1 6 197.17 2.1 12 -15.029 2.1 6 0.96703 3.3 12 -1.0223 3.3 6 14.01 2.3 8 -2.0038 2.3 3.3 6.8572 2.7 8 -1.1197 2.7 2.9 198.97 2.6 9 -41.541 2.6 5.7 1.7816 3.1 12 -1.2733 3.1 6 13.701 2.1 12 -5.0914 2.1 6 1.8956 2.9 12 -1.4598 2.9 6 1.0617 2.7 12 -1.2739 2.7 6 0.10702 2.8 12 -0.26622 2.8 6 92.013 2.1 12 -10.217 2.1 6 4.8266 2.8 8 -3.8502 2.8 6 26.458 2.3 8 -5.3033 2.3 4.3 34.228 2.3 12 -7.8467 2.3 6 1.2981 3.4 12 -1.1398 3.4 6 4.4558 3.1 12 -3.133 3.1 6 24.406 3 12 -6.3909 3 6 1.3596 3.8 8 -1.0994 3.8 3.3 1.3596 3.8 8 -1.0994 3.8 3.3 1.3596 3.8 8 -1.0994 3.8 3.3 18.851 2.7 12 -5.693 2.7 6 1634.2 2.5 12 -45.346 2.5 6 7.4944 3.3 12 -4.5866 3.3 6 124.96 2.3 12 -12.483 2.3 6 2.0907 3.1 12 -1.2325 3.1 6 0.59859 3.4 12 -0.61423 3.4 6 0.6524 3.8 8 -0.86987 3.8 6 0.78413 3.7 12 -0.98236 3.7 6 4528.6 2.2 12 -108.39 2.2 6 0.45392 3.6 12 -0.64607 3.6 6 0.14868 3 12 -0.44984 3 6 0.65898 3 12 -0.74395 3 6 0.62491 2.7 12 -0.81403 2.7 6 0.42495 2.7 12 -0.68329 2.7 6 1.3477 3.8 8 -1.2142 3.8 2.4 0.79986 2.8 12 -1.0575 2.8 6 0.54543 3.2 8 -0.6317 3.2 3.7 86.104 2.1 12 -11.436 2.1 6 257.77 2.5 12 -19.911 2.5 6 346.03 2.1 12 -19.713 2.1 6 4.4164 3.1 8 -3.8922 3.1 5.8 83.488 2.7 12 -11.785 2.7 6 83.488 2.7 12 -11.785 2.7 6 83.488 2.7 12 -11.785 2.7 6 1.1531 3.5 12 -1.035 3.5 6 310.86 2.1 12 -20.783 2.1 6 662.34 2.2 12 -42.033 2.2 6 0.42305 3.6 9 -0.52698 3.6 6 23.318 2.4 11 -4.5902 2.4 6 4.4931 2.8 12 -1.8617 2.8 6 2.1818 2.9 12 -1.6802 2.9 6 0.02744 3.9 12 -0.10505 3.9 6 1.4473 3.4 8 -1.0633 3.4 3 7.2536 2.5 12 -2.6448 2.5 6 0.14888 3 12 -0.5729 3 6 0.82943 2.8 12 -0.61898 2.8 6 1.5013 2.5 8 -0.74436 2.5 3 0.22857 3 8 -0.24627 3 2.2 0.49513 3.2 8 -0.5352 3.2 3.3 1.0233 2.7 12 -1.2278 2.7 6 12.11 2.2 12 -6.4896 2.2 6 0.5502 2.8 12 -1.1004 2.8 6 0.012027 3.8 12 -0.1503 3.8 6 62.833 2.4 12 -7.976 2.4 6 1.0465 2.8 12 -0.6529 2.8 6 112.87 2.1 12 -10.532 2.1 6 14.114 3 8 -6.5609 3 5.2 14.114 3 8 -6.5609 3 5.2 113.79 2.1 12 -15.113 2.1 6 36.507 2.1 11 -5.2045 2.1 6 73.666 2.1 12 -5.8049 2.1 6 25.39 2.3 8 -13.818 2.3 6 53.385 2.1 10 -11.531 2.1 6 1.1048 3 8 -0.94661 3 5.3 0.11688 3.7 12 -0.40224 3.7 6 0.98697 3.1 12 -0.95301 3.1 6 0.44837 3.7 12 -0.76415 3.7 6 8.8969 2.6 12 -3.3756 2.6 6 0.458 2.7 12 -1.3857 2.7 6 59.487 2.1 12 -5.5508 2.1 6 0.12279 3.3 8 -0.33058 3.3 2.1 4.7732 2.2 8 -2.62 2.2 4.7 60.199 2 8 -2.5201 2 2.5 5.3778 2.2 12 -2.5306 2.2 6 3.1221 2.4 12 -2.5851 2.4 6 0.065401 3.4 12 -0.42841 3.4 6 18.704 2.1 12 -6.5169 2.1 6 1.4341 3.5 12 -1.11 3.5 6 0.19009 3.5 8 -0.26124 3.5 2.4 3.1602 3.8 8 -3.1427 3.8 6 25.83 3 12 -6.861 3 6 25.83 3 12 -6.861 3 6 3.6405 3 12 -2.282 3 6 3.2056 2.7 12 -1.8548 2.7 6 182.24 2.4 11 -20.356 2.4 6 0.094078 3.6 12 -0.27961 3.6 6 11.172 2.3 8 -5.3017 2.3 5.5 11.368 2.4 12 -1.9961 2.4 6 8.7154 2.2 8 -3.1748 2.2 4.2 0.01802 3.9 12 -0.098197 3.9 6 134.6 2.1 12 -14.945 2.1 6 0.051493 3.9 12 -0.19713 3.9 6 0.41405 2.7 12 -0.79223 2.7 6 0.57448 3 12 -0.68509 3 6 0.56161 2.6 12 -1.0728 2.6 6 5.2246 2.2 12 -3.8351 2.2 6 14.927 2.5 11 -5.722 2.5 5.8 8.1415 2.1 12 -2.4473 2.1 6 22.172 2.1 12 -6.6649 2.1 6 7.3779 2.2 12 -4.0413 2.2 6 0.4081 2.8 12 -0.81621 2.8 6 3.4998 3 12 -2.1229 3 6 1.358 2.5 8 -0.78079 2.5 2.9 0.994 3.3 12 -0.95334 3.3 6 1.3618 3 12 -0.9121 3 6 1.3618 3 12 -0.9121 3 6 2.6336 3 8 -1.0348 3 3 0.58175 3.2 8 -0.69187 3.2 3.5 537.66 2.2 12 -38.026 2.2 6 0.2193 2.7 15 -0.64248 2.7 7.9 55.368 2.1 10 -13.43 2.1 6 0.35097 3.4 10 -0.44528 3.4 6 0.26097 3.7 12 -0.32694 3.7 6 19.795 2.8 12 -4.0162 2.8 6 366.47 2.4 12 -22.157 2.4 6 1.1956 3.4 12 -0.94352 3.4 6 6.7482 3 8 -3.302 3 4.8 8.8525 3.1 12 -2.8627 3.1 6 0.060924 3.9 12 -0.19697 3.9 6 20.761 3.4 12 -4.5223 3.4 6 292.35 2.4 12 -17.676 2.4 6 1.3998 3.3 12 -1.2008 3.3 6 103.79 2.4 12 -11.273 2.4 6 13.085 3 12 -4.0164 3 6 848.21 2 12 -26.506 2 6 1752.5 2.3 12 -48.136 2.3 6 273.25 2.3 11 -27.831 2.3 6 1.6793 3.4 12 -1.7231 3.4 6 1.6793 3.4 12 -1.7231 3.4 6 1.6793 3.4 12 -1.7231 3.4 6 62.622 2 8 -2.3423 2 3.2 45.758 3 12 -8.0343 3 6 64.16 3 12 -8.6783 3 6 70.583 2.4 12 -7.6667 2.4 6 1.9166 3.3 12 -1.2086 3.3 6 30.53 2.6 12 -4.6051 2.6 6 2.935 3.2 8 -2.9271 3.2 6 17.047 2.7 12 -5.1482 2.7 6 17.047 2.7 12 -5.1482 2.7 6 20.271 2.1 8 -4.1216 2.1 4.3 34.557 2.2 8 -1.7947 2.2 2.3 24.232 3 12 -6.4364 3 6 0.32163 3.4 8 -0.44072 3.4 2.9 140.69 2.1 12 -13.128 2.1 6 12.728 3.1 8 -8.1702 3.1 5.5 150.48 2.1 12 -16.709 2.1 6 2.633 2.8 12 -2.5008 2.8 6 16.039 3 12 -4.6444 3 6 848.21 2 12 -26.506 2 6 42.931 3 12 -7.538 3 6 357.17 2.6 11 -26.74 2.6 5.7 1.6793 3.4 12 -1.7231 3.4 6 1.6793 3.4 12 -1.7231 3.4 6 1.6793 3.4 12 -1.7231 3.4 6 2122 2.4 12 -64.677 2.4 6 820.75 2.4 12 -43.231 2.4 6 7062.2 2 11 -109.81 2 5.6 120.86 2.2 12 -8.5478 2.2 6 1.0125 3 8 -1.1533 3 5.1 3.489 3.4 9 -2.7764 3.4 6 ]; R_int=1; R_end=20; rd = 0.005; Ruler=(R_int:rd:R_end)'; GARF_distribution(:,i) = exp((Paraset(Paraset_index_num,1).*(Paraset(Paraset_index_num,2)./Ruler).^Paraset(Paraset_index_num,3)+Paraset(Paraset_index_num,4).*(Paraset(Paraset_index_num,5)./Ruler).^Paraset(Paraset_index_num,6))./-0.5918);