├── README.md ├── examples ├── NiTe2 │ ├── EIGENVAL │ ├── KPOINTS │ ├── POSCAR │ ├── ebands3d.bxsf │ └── fs.png ├── SrVO3 │ ├── EIGENVAL │ ├── KPOINTS │ ├── POSCAR │ ├── ebands3d.bxsf │ └── fs.png └── copper │ ├── EIGENVAL │ ├── KPOINTS │ ├── OUTCAR │ ├── POSCAR │ ├── ebands3d.bxsf │ ├── fs.png │ ├── fs_matplotlib_screenshot.jpg │ ├── fs_matplotlib_screenshot.png │ ├── fs_mayavi_screenshot.jpg │ ├── fs_mayavi_screenshot.png │ ├── fs_uc_mayavi_screenshot.png │ └── fs_xcrysden_screenshot.png └── fs.py /README.md: -------------------------------------------------------------------------------- 1 | # Fermi Surface from VASP calculation 2 | 3 | ## Fermi Surface 4 | 5 | In condensed matter physics, the [Fermi Surface](https://en.wikipedia.org/wiki/Fermi_surface) is the surface of iso-energy in the reciprocal space, the energy of which is referred to as the *Fermi energy* or *Fermi level*. To visualize the Fermi surface, first you have to obtain the band energies in the reciprocal space. For example, one can sample the reciprocal space primitive cell using a uniform k-points grid and then use any *ab initio* code, e.g. *VASP*, to generate the band energies on these uniform grid points. With the band energies on the uniform grid, one can either plot the iso-surface with some coding or just use [Xcrysden](http://www.xcrysden.org/doc/fermi.html) to do the plotting work. 6 | 7 | ## Prerequisites 8 | 9 | The following python library must be installed to run the script 10 | 11 | * [ase](https://wiki.fysik.dtu.dk/ase/ase/io/io.html) 12 | * [spglib](https://atztogo.github.io/spglib/) 13 | * [scikit-image](https://scikit-image.org/) 14 | * [matplotlib [optional]](https://matplotlib.org/) 15 | * [mayavi [optional]](https://docs.enthought.com/mayavi/mayavi/) 16 | 17 | ## Steps to visualize the Fermi surface 18 | 19 | **1. Calculate the band energies using VASP.** 20 | 21 | Note that Gamma-centered `KPOINTS` is required. 22 | 23 | ``` 24 | Monkhorst-Pack 25 | 0 26 | Gamma 27 | 21 21 21 28 | 0 0 0 29 | ``` 30 | 31 | Moreover, there is *no need* to turn off the symmetry in the `INCAR` file. As is well-known, VASP utilize the symmetry to fold the k-points grid to symmetry-irreducible k-points in the BZ, which will significantly reduce the computation cost. The script use `spglib` to unfold the irreducible k-points in BZ to the uniform k-points grid. 32 | 33 | **2. Visualizing the Fermi surface** 34 | 35 | At this step, use the script `fs.py` to process the `VASP` output files. 36 | 37 | ``` 38 | usage: fs.py [-h] [-i EIGVENVAL] [--poscar POSCAR] [--kpoints KPOINTS] 39 | [--plot {xcrys,mpl,mayavi}] [--cell {uc,bz}] [--efermi EFERMI] 40 | [--symprec SYMPREC] [--kmesh KMESH KMESH KMESH] 41 | 42 | optional arguments: 43 | -h, --help show this help message and exit 44 | -i EIGVENVAL Location of the EIGENVAL file. By default, the 45 | KPOINTS/POSCAR are also read from the directory where 46 | 47 | 48 | --kpoints KPOINTS location of VASP KPOINTS 49 | --plot {xcrys,mpl,mayavi} 50 | Fermi surface plotting method 51 | 52 | --efermi EFERMI the Fermi energy 53 | --symprec SYMPREC the symmetry precision used in spglib 54 | --kmesh KMESH KMESH KMESH 55 | the kmesh in the KPOINTS 56 | --interp NEW_KMESH NEW_KMESH NEW_KMESH 57 | the new grid size 58 | ``` 59 | 60 | There are **three** different ways to visualize the Fermi surface: 61 | 62 | + Using Xcrysden 63 | 64 | First, you need to generate the `bxsf` file as input for `xcrysden` 65 | 66 | ```bash 67 | fs.py -i EIGENVAL 68 | ``` 69 | 70 | then, launch xcrysden 71 | 72 | ```bash 73 | xcrysden --bxsf ebands3d.bxsf 74 | ``` 75 | 76 | For example, below is the the Fermi surface of `fcc Cu` 77 | 78 |

79 | 80 |

81 | 82 | + Using Matplotlib 83 | 84 | ```bash 85 | fs.py --plot mpl 86 | ``` 87 | 88 |

89 | 90 |

91 | 92 | + Using Mayavi 93 | 94 | ```bash 95 | fs.py --plot mayavi # Fermi Surface in Brillouin Zone 96 | fs.py --plot mayavi --cell uc # Fermi Surface in primitive cell 97 | ``` 98 | 99 |

100 | 101 | 102 |

103 | 104 | You can also use `vaspkit` ([tutorial](https://vaspkit.com/tutorials.html#fermi-surface)) to generate `KPOINTS` with given k-spacing. Since the content in `KPOINTS` is quiet different, the script requires to input `kmesh` manually. 105 | 106 | ```bash 107 | fs.py --kmesh 42 42 24 -i EIGENVAL 108 | fs.py --kmesh 42 42 24 --interp 210 210 120 --plot mayavi --efermi 9.05973773 109 | ``` 110 | 111 |

112 | 113 |

114 | 115 | *** 116 | 117 | _This file was generated by Remarkable on Thursday, 27. February 2020 04:45PM_ 118 | -------------------------------------------------------------------------------- /examples/NiTe2/POSCAR: -------------------------------------------------------------------------------- 1 | generated by phonopy 2 | 1.0 3 | 3.7403137504547965 0.0000000000000000 0.0000000000000000 4 | -1.8701568752273983 3.2392067260181032 0.0000000000000000 5 | 0.0000000000000000 0.0000000000000000 5.2659797668000010 6 | Ni Te 7 | 1 2 8 | Direct 9 | 0.0000000000000000 0.0000000000000000 0.4999999999999999 10 | 0.6666666666666666 0.3333333333333333 0.7665059859999999 11 | 0.3333333333333334 0.6666666666666667 0.2334940139999999 12 | -------------------------------------------------------------------------------- /examples/NiTe2/fs.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/NiTe2/fs.png -------------------------------------------------------------------------------- /examples/SrVO3/KPOINTS: -------------------------------------------------------------------------------- 1 | Monkhorst-Pack 2 | 0 3 | Gamma 4 | 21 21 21 5 | 0 0 0 6 | -------------------------------------------------------------------------------- /examples/SrVO3/POSCAR: -------------------------------------------------------------------------------- 1 | Sr V O 2 | 1.0000000000000000 3 | 3.9008910000000001 0.0000000000000000 0.0000000000000000 4 | 0.0000000000000000 3.9008910000000001 0.0000000000000000 5 | 0.0000000000000000 0.0000000000000000 3.9008910000000001 6 | Sr V O 7 | 1 1 3 8 | Direct 9 | 0.5000000000000000 0.5000000000000000 0.5000000000000000 10 | 0.0000000000000000 0.0000000000000000 0.0000000000000000 11 | 0.0000000000000000 0.5000000000000000 0.0000000000000000 12 | 0.0000000000000000 0.0000000000000000 0.5000000000000000 13 | 0.5000000000000000 0.0000000000000000 0.0000000000000000 14 | -------------------------------------------------------------------------------- /examples/SrVO3/fs.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/SrVO3/fs.png -------------------------------------------------------------------------------- /examples/copper/EIGENVAL: -------------------------------------------------------------------------------- 1 | 1 1 1000 1 2 | 0.1197732E+02 0.2568176E-09 0.2568176E-09 0.2568176E-09 0.5000000E-15 3 | 1.000000000000000E-004 4 | CAR 5 | fcc copper 6 | 11 286 12 7 | 8 | 0.3742002E-15 0.3742002E-15 0.3742002E-15 0.1079797E-03 9 | 1 -1.835058 10 | 2 4.425966 11 | 3 4.425970 12 | 4 4.425970 13 | 5 5.235326 14 | 6 5.235326 15 | 7 29.895653 16 | 8 32.425637 17 | 9 32.425674 18 | 10 32.425685 19 | 11 34.907028 20 | 12 34.907063 21 | 22 | 0.4761905E-01 0.3712308E-15 0.3712308E-15 0.8638376E-03 23 | 1 -1.753154 24 | 2 4.408741 25 | 3 4.441552 26 | 4 4.441555 27 | 5 5.229121 28 | 6 5.229121 29 | 7 29.445080 30 | 8 31.280182 31 | 9 32.145006 32 | 10 32.145040 33 | 11 35.305804 34 | 12 35.305966 35 | 36 | 0.9523810E-01 0.3816392E-15 0.3816392E-15 0.8638376E-03 37 | 1 -1.509372 38 | 2 4.360421 39 | 3 4.483628 40 | 4 4.483630 41 | 5 5.214598 42 | 6 5.214600 43 | 7 27.566975 44 | 8 30.632301 45 | 9 31.587902 46 | 10 31.587934 47 | 11 36.194588 48 | 12 36.196697 49 | 50 | 0.1428571E+00 0.3747003E-15 0.3747003E-15 0.8638376E-03 51 | 1 -1.110126 52 | 2 4.291853 53 | 3 4.537692 54 | 4 4.537692 55 | 5 5.205372 56 | 6 5.205377 57 | 7 25.138539 58 | 8 30.411216 59 | 9 31.018839 60 | 10 31.019432 61 | 11 36.791820 62 | 12 37.135042 63 | 64 | 0.1904762E+00 0.3747003E-15 0.3747003E-15 0.8638376E-03 65 | 1 -0.568525 66 | 2 4.221556 67 | 3 4.578882 68 | 4 4.578884 69 | 5 5.224137 70 | 6 5.224146 71 | 7 22.606361 72 | 8 30.170305 73 | 9 30.516791 74 | 10 30.516819 75 | 11 36.744348 76 | 12 37.035444 77 | 78 | 0.2380952E+00 0.3191891E-15 0.3191891E-15 0.8638376E-03 79 | 1 0.090024 80 | 2 4.180326 81 | 3 4.584475 82 | 4 4.584478 83 | 5 5.290277 84 | 6 5.290289 85 | 7 19.516059 86 | 8 29.849184 87 | 9 30.102341 88 | 10 30.102347 89 | 11 35.958211 90 | 12 35.958640 91 | 92 | 0.2857143E+00 0.3330669E-15 0.3330669E-15 0.8638376E-03 93 | 1 0.814545 94 | 2 4.219293 95 | 3 4.554659 96 | 4 4.554663 97 | 5 5.398575 98 | 6 5.398590 99 | 7 16.683079 100 | 8 29.473900 101 | 9 29.779113 102 | 10 29.779121 103 | 11 34.854426 104 | 12 34.854505 105 | 106 | 0.3333333E+00 0.3608225E-15 0.3608225E-15 0.8638376E-03 107 | 1 1.504859 108 | 2 4.429809 109 | 3 4.506699 110 | 4 4.506703 111 | 5 5.524138 112 | 6 5.524155 113 | 7 15.810787 114 | 8 29.097676 115 | 9 29.543781 116 | 10 29.543805 117 | 11 33.896447 118 | 12 33.896506 119 | 120 | 0.3809524E+00 0.3053113E-15 0.3053113E-15 0.8638376E-03 121 | 1 2.018655 122 | 2 4.458060 123 | 3 4.458063 124 | 4 4.923506 125 | 5 5.642112 126 | 6 5.642129 127 | 7 13.990619 128 | 8 28.772010 129 | 9 29.388488 130 | 10 29.388531 131 | 11 33.133715 132 | 12 33.133793 133 | 134 | 0.4285714E+00 0.3330669E-15 0.3330669E-15 0.8638376E-03 135 | 1 2.294735 136 | 2 4.419817 137 | 3 4.419821 138 | 4 5.671521 139 | 5 5.732203 140 | 6 5.732222 141 | 7 12.342859 142 | 8 28.535094 143 | 9 29.300381 144 | 10 29.300441 145 | 11 32.599290 146 | 12 32.599392 147 | 148 | 0.4761905E+00 0.3885781E-15 0.3885781E-15 0.8638376E-03 149 | 1 2.400532 150 | 2 4.399172 151 | 3 4.399175 152 | 4 5.780973 153 | 5 5.780993 154 | 6 6.338733 155 | 7 11.255931 156 | 8 28.410922 157 | 9 29.262930 158 | 10 29.263001 159 | 11 32.322336 160 | 12 32.322442 161 | 162 | 0.4761905E-01 0.4761905E-01 0.3729628E-15 0.6478782E-03 163 | 1 -1.725877 164 | 2 4.393979 165 | 3 4.451612 166 | 4 4.451614 167 | 5 5.206392 168 | 6 5.247084 169 | 7 29.510504 170 | 8 31.324889 171 | 9 31.326021 172 | 10 31.973981 173 | 11 35.524491 174 | 12 36.063281 175 | 176 | 0.9523810E-01 0.4761905E-01 0.3712308E-15 0.2591513E-02 177 | 1 -1.536153 178 | 2 4.350749 179 | 3 4.487356 180 | 4 4.489848 181 | 5 5.173795 182 | 6 5.254250 183 | 7 28.328770 184 | 8 30.270869 185 | 9 30.528115 186 | 10 32.242543 187 | 11 35.513852 188 | 12 37.034229 189 | 190 | 0.1428571E+00 0.4761905E-01 0.3747003E-15 0.2591513E-02 191 | 1 -1.188757 192 | 2 4.286449 193 | 3 4.539321 194 | 4 4.543551 195 | 5 5.142541 196 | 6 5.260736 197 | 7 26.013346 198 | 8 29.823048 199 | 9 29.837819 200 | 10 32.675095 201 | 11 35.178116 202 | 12 37.522321 203 | 204 | 0.1904762E+00 0.4761905E-01 0.3885781E-15 0.2591513E-02 205 | 1 -0.694238 206 | 2 4.213534 207 | 3 4.585829 208 | 4 4.591088 209 | 5 5.138037 210 | 6 5.284588 211 | 7 23.530119 212 | 8 29.229828 213 | 9 29.503825 214 | 10 32.805543 215 | 11 34.645575 216 | 12 37.568843 217 | 218 | 0.2380952E+00 0.4761905E-01 0.3747003E-15 0.2591513E-02 219 | 1 -0.073039 220 | 2 4.155954 221 | 3 4.595197 222 | 4 4.609906 223 | 5 5.192615 224 | 6 5.344260 225 | 7 20.732685 226 | 8 28.756736 227 | 9 29.178442 228 | 10 32.391353 229 | 11 34.178714 230 | 12 37.410651 231 | 232 | 0.2857143E+00 0.4761905E-01 0.3191891E-15 0.2591513E-02 233 | 1 0.634416 234 | 2 4.153409 235 | 3 4.560977 236 | 4 4.593057 237 | 5 5.312119 238 | 6 5.440818 239 | 7 17.609459 240 | 8 28.409150 241 | 9 28.865688 242 | 10 31.591365 243 | 11 33.822821 244 | 12 36.541653 245 | 246 | 0.3333333E+00 0.4761905E-01 0.3330669E-15 0.2591513E-02 247 | 1 1.346799 248 | 2 4.274268 249 | 3 4.504726 250 | 4 4.552185 251 | 5 5.477826 252 | 6 5.555294 253 | 7 16.238126 254 | 8 28.192955 255 | 9 28.577748 256 | 10 30.708755 257 | 11 33.512584 258 | 12 35.523560 259 | 260 | 0.3809524E+00 0.4761905E-01 0.3608225E-15 0.2591513E-02 261 | 1 1.928226 262 | 2 4.446721 263 | 3 4.503668 264 | 4 4.607835 265 | 5 5.663696 266 | 6 5.687861 267 | 7 14.659246 268 | 8 28.115437 269 | 9 28.322371 270 | 10 29.896875 271 | 11 33.284847 272 | 12 34.630921 273 | 274 | 0.4285714E+00 0.4761905E-01 0.3885781E-15 0.2591513E-02 275 | 1 2.274677 276 | 2 4.400968 277 | 3 4.460340 278 | 4 5.096330 279 | 5 5.745501 280 | 6 6.019886 281 | 7 12.893465 282 | 8 28.088716 283 | 9 28.184420 284 | 10 29.247132 285 | 11 33.177626 286 | 12 33.884867 287 | 288 | 0.4761905E+00 0.4761905E-01 0.3608225E-15 0.2591513E-02 289 | 1 2.424048 290 | 2 4.374150 291 | 3 4.431156 292 | 4 5.414930 293 | 5 5.786301 294 | 6 6.561485 295 | 7 11.553904 296 | 8 27.833330 297 | 9 28.405954 298 | 10 28.852347 299 | 11 33.207903 300 | 12 33.294720 301 | 302 | -0.4761905E+00 0.4761905E-01 0.4440892E-15 0.2591513E-02 303 | 1 2.455218 304 | 2 4.369667 305 | 3 4.421125 306 | 4 5.495749 307 | 5 5.778657 308 | 6 6.807556 309 | 7 11.132304 310 | 8 27.547196 311 | 9 28.751493 312 | 10 28.780352 313 | 11 32.866921 314 | 12 33.376619 315 | 316 | -0.4285714E+00 0.4761905E-01 0.4163336E-15 0.2591513E-02 317 | 1 2.390809 318 | 2 4.387689 319 | 3 4.432209 320 | 4 5.373327 321 | 5 5.721820 322 | 6 6.395682 323 | 7 11.930654 324 | 8 27.307858 325 | 9 28.880223 326 | 10 29.295590 327 | 11 32.613558 328 | 12 33.673115 329 | 330 | -0.3809524E+00 0.4761905E-01 0.3608225E-15 0.2591513E-02 331 | 1 2.186042 332 | 2 4.425285 333 | 3 4.463101 334 | 4 4.958343 335 | 5 5.622706 336 | 6 5.959299 337 | 7 13.450219 338 | 8 27.174494 339 | 9 29.175336 340 | 10 29.913745 341 | 11 32.566153 342 | 12 34.072460 343 | 344 | -0.3333333E+00 0.4761905E-01 0.3608225E-15 0.2591513E-02 345 | 1 1.758682 346 | 2 4.470109 347 | 3 4.486726 348 | 4 4.509401 349 | 5 5.494245 350 | 6 5.733699 351 | 7 15.292365 352 | 8 27.165573 353 | 9 29.596376 354 | 10 30.550604 355 | 11 32.797463 356 | 12 34.524228 357 | 358 | -0.2857143E+00 0.4761905E-01 0.3608225E-15 0.2591513E-02 359 | 1 1.116769 360 | 2 4.221212 361 | 3 4.530598 362 | 4 4.561437 363 | 5 5.354888 364 | 6 5.572239 365 | 7 16.606175 366 | 8 27.281418 367 | 9 30.100301 368 | 10 31.094270 369 | 11 33.402560 370 | 12 34.941371 371 | 372 | -0.2380952E+00 0.4761905E-01 0.3330669E-15 0.2591513E-02 373 | 1 0.393105 374 | 2 4.150357 375 | 3 4.573983 376 | 4 4.602806 377 | 5 5.230664 378 | 6 5.436041 379 | 7 18.631181 380 | 8 27.516010 381 | 9 30.624536 382 | 10 31.528461 383 | 11 34.373309 384 | 12 35.234622 385 | 386 | -0.1904762E+00 0.4761905E-01 0.3469447E-15 0.2591513E-02 387 | 1 -0.292217 388 | 2 4.177102 389 | 3 4.590067 390 | 4 4.609677 391 | 5 5.152690 392 | 6 5.333040 393 | 7 21.723195 394 | 8 27.862044 395 | 9 31.079841 396 | 10 31.930416 397 | 11 35.413998 398 | 12 35.568059 399 | 400 | -0.1428571E+00 0.4761905E-01 0.3885781E-15 0.2591513E-02 401 | 1 -0.874471 402 | 2 4.243937 403 | 3 4.569639 404 | 4 4.571650 405 | 5 5.136909 406 | 6 5.273600 407 | 7 24.363694 408 | 8 28.315725 409 | 9 31.374646 410 | 10 32.342666 411 | 11 35.590829 412 | 12 36.354045 413 | 414 | -0.9523810E-01 0.4761905E-01 0.3747003E-15 0.2591513E-02 415 | 1 -1.321639 416 | 2 4.317147 417 | 3 4.512712 418 | 4 4.523216 419 | 5 5.162218 420 | 6 5.249933 421 | 7 26.805517 422 | 8 28.903367 423 | 9 31.486007 424 | 10 32.711082 425 | 11 35.560318 426 | 12 35.793027 427 | 428 | -0.4761905E-01 0.4761905E-01 0.3816392E-15 0.1295756E-02 429 | 1 -1.617311 430 | 2 4.376467 431 | 3 4.463696 432 | 4 4.473683 433 | 5 5.198718 434 | 6 5.240920 435 | 7 28.731150 436 | 8 29.988374 437 | 9 31.522028 438 | 10 32.735641 439 | 11 35.018388 440 | 12 35.874053 441 | 442 | 0.9523810E-01 0.9523810E-01 0.3761130E-15 0.6478782E-03 443 | 1 -1.401383 444 | 2 4.301351 445 | 3 4.527403 446 | 4 4.527405 447 | 5 5.123887 448 | 6 5.281124 449 | 7 28.694203 450 | 8 29.581543 451 | 9 29.581563 452 | 10 30.955438 453 | 11 37.034082 454 | 12 37.090803 455 | 456 | 0.1428571E+00 0.9523810E-01 0.3677614E-15 0.2591513E-02 457 | 1 -1.108576 458 | 2 4.233508 459 | 3 4.582885 460 | 4 4.584795 461 | 5 5.074029 462 | 6 5.302291 463 | 7 26.961889 464 | 8 28.762390 465 | 9 29.140132 466 | 10 31.232846 467 | 11 35.477593 468 | 12 37.273005 469 | 470 | 0.1904762E+00 0.9523810E-01 0.3747003E-15 0.2591513E-02 471 | 1 -0.666257 472 | 2 4.156154 473 | 3 4.629221 474 | 4 4.642240 475 | 5 5.055596 476 | 6 5.332353 477 | 7 24.528059 478 | 8 28.077303 479 | 9 28.821216 480 | 10 32.049980 481 | 11 34.100005 482 | 12 37.139859 483 | 484 | 0.2380952E+00 0.9523810E-01 0.3885781E-15 0.2591513E-02 485 | 1 -0.091235 486 | 2 4.084598 487 | 3 4.622809 488 | 4 4.676447 489 | 5 5.119786 490 | 6 5.388743 491 | 7 21.948274 492 | 8 27.531050 493 | 9 28.490368 494 | 10 31.855258 495 | 11 34.018962 496 | 12 36.939859 497 | 498 | 0.2857143E+00 0.9523810E-01 0.3747003E-15 0.2591513E-02 499 | 1 0.583862 500 | 2 4.045720 501 | 3 4.562623 502 | 4 4.674777 503 | 5 5.276342 504 | 6 5.477324 505 | 7 18.844339 506 | 8 27.129049 507 | 9 28.165159 508 | 10 30.732195 509 | 11 34.883470 510 | 12 36.909824 511 | 512 | 0.3333333E+00 0.9523810E-01 0.3191891E-15 0.2591513E-02 513 | 1 1.293746 514 | 2 4.083900 515 | 3 4.483983 516 | 4 4.643714 517 | 5 5.505648 518 | 6 5.583664 519 | 7 16.657272 520 | 8 26.878276 521 | 9 27.826960 522 | 10 29.659901 523 | 11 35.283865 524 | 12 37.371438 525 | 526 | 0.3809524E+00 0.9523810E-01 0.3885781E-15 0.2591513E-02 527 | 1 1.918670 528 | 2 4.253629 529 | 3 4.420545 530 | 4 4.598339 531 | 5 5.684956 532 | 6 5.821250 533 | 7 15.442462 534 | 8 26.787325 535 | 9 27.406744 536 | 10 28.848678 537 | 11 35.228189 538 | 12 36.372621 539 | 540 | 0.4285714E+00 0.9523810E-01 0.3608225E-15 0.2591513E-02 541 | 1 2.328900 542 | 2 4.326503 543 | 3 4.552476 544 | 4 4.617884 545 | 5 5.760296 546 | 6 6.289392 547 | 7 13.592527 548 | 8 26.837872 549 | 9 26.865514 550 | 10 28.409898 551 | 11 35.097202 552 | 12 35.558152 553 | 554 | 0.4761905E+00 0.9523810E-01 0.3608225E-15 0.2591513E-02 555 | 1 2.521993 556 | 2 4.294913 557 | 3 4.516163 558 | 4 4.925592 559 | 5 5.794477 560 | 6 6.943897 561 | 7 12.036922 562 | 8 26.216060 563 | 9 27.120833 564 | 10 28.278532 565 | 11 34.875486 566 | 12 35.008918 567 | 568 | -0.4761905E+00 0.9523810E-01 0.4996004E-15 0.2591513E-02 569 | 1 2.582706 570 | 2 4.285390 571 | 3 4.496335 572 | 4 5.078113 573 | 5 5.779957 574 | 6 7.445473 575 | 7 11.193830 576 | 8 25.666661 577 | 9 27.558544 578 | 10 28.353716 579 | 11 34.327866 580 | 12 34.983177 581 | 582 | -0.4285714E+00 0.9523810E-01 0.4440892E-15 0.2591513E-02 583 | 1 2.555540 584 | 2 4.302360 585 | 3 4.496611 586 | 4 5.039731 587 | 5 5.715765 588 | 6 7.191007 589 | 7 11.647325 590 | 8 25.244829 591 | 9 28.180047 592 | 10 28.597306 593 | 11 33.914173 594 | 12 35.002924 595 | 596 | -0.3809524E+00 0.9523810E-01 0.3885781E-15 0.2591513E-02 597 | 1 2.417988 598 | 2 4.343399 599 | 3 4.517933 600 | 4 4.791177 601 | 5 5.608102 602 | 6 6.582367 603 | 7 13.034104 604 | 8 24.969843 605 | 9 28.981584 606 | 10 28.998372 607 | 11 33.631853 608 | 12 35.026935 609 | 610 | -0.3333333E+00 0.9523810E-01 0.3053113E-15 0.2591513E-02 611 | 1 2.089094 612 | 2 4.401610 613 | 3 4.427148 614 | 4 4.557974 615 | 5 5.469121 616 | 6 6.134973 617 | 7 14.823355 618 | 8 24.846370 619 | 9 29.549851 620 | 10 29.954065 621 | 11 33.478461 622 | 12 34.996204 623 | 624 | -0.2857143E+00 0.9523810E-01 0.3608225E-15 0.2591513E-02 625 | 1 1.524045 626 | 2 4.170218 627 | 3 4.471170 628 | 4 4.610599 629 | 5 5.315460 630 | 6 5.853046 631 | 7 16.514609 632 | 8 24.871901 633 | 9 30.237177 634 | 10 31.081990 635 | 11 33.454206 636 | 12 34.850710 637 | 638 | -0.2380952E+00 0.9523810E-01 0.3608225E-15 0.2591513E-02 639 | 1 0.824003 640 | 2 4.084219 641 | 3 4.535996 642 | 4 4.662124 643 | 5 5.171376 644 | 6 5.647161 645 | 7 17.954781 646 | 8 25.040103 647 | 9 31.024288 648 | 10 32.327899 649 | 11 33.576556 650 | 12 34.570130 651 | 652 | -0.1904762E+00 0.9523810E-01 0.3608225E-15 0.2591513E-02 653 | 1 0.125813 654 | 2 4.102190 655 | 3 4.582002 656 | 4 4.682910 657 | 5 5.074686 658 | 6 5.485476 659 | 7 20.905537 660 | 8 25.346150 661 | 9 31.803509 662 | 10 33.402467 663 | 11 34.101513 664 | 12 34.245657 665 | 666 | -0.1428571E+00 0.9523810E-01 0.3469447E-15 0.2591513E-02 667 | 1 -0.488398 668 | 2 4.167720 669 | 3 4.594544 670 | 4 4.643322 671 | 5 5.061018 672 | 6 5.367695 673 | 7 23.682050 674 | 8 25.801113 675 | 9 32.223410 676 | 10 33.772651 677 | 11 34.249262 678 | 12 35.408591 679 | 680 | -0.9523810E-01 0.9523810E-01 0.3885781E-15 0.1295756E-02 681 | 1 -0.977765 682 | 2 4.246534 683 | 3 4.569766 684 | 4 4.572197 685 | 5 5.104317 686 | 6 5.293905 687 | 7 25.868799 688 | 8 26.678391 689 | 9 31.983358 690 | 10 33.640587 691 | 11 34.902743 692 | 12 35.359459 693 | 694 | 0.1428571E+00 0.1428571E+00 0.3761130E-15 0.6478782E-03 695 | 1 -0.871939 696 | 2 4.157477 697 | 3 4.650632 698 | 4 4.650634 699 | 5 5.002122 700 | 6 5.334634 701 | 7 27.824841 702 | 8 27.841636 703 | 9 27.841658 704 | 10 29.750883 705 | 11 34.074405 706 | 12 39.044686 707 | 708 | 0.1904762E+00 0.1428571E+00 0.3747003E-15 0.2591513E-02 709 | 1 -0.486090 710 | 2 4.071883 711 | 3 4.707115 712 | 4 4.722975 713 | 5 4.965442 714 | 6 5.369275 715 | 7 25.577945 716 | 8 27.065375 717 | 9 28.099084 718 | 10 29.806218 719 | 11 32.761380 720 | 12 39.241694 721 | 722 | 0.2380952E+00 0.1428571E+00 0.3747003E-15 0.2591513E-02 723 | 1 0.034599 724 | 2 3.988104 725 | 3 4.671522 726 | 4 4.777007 727 | 5 5.055522 728 | 6 5.422714 729 | 7 23.115183 730 | 8 26.438615 731 | 9 27.832530 732 | 10 29.894301 733 | 11 32.472769 734 | 12 38.885887 735 | 736 | 0.2857143E+00 0.1428571E+00 0.3747003E-15 0.2591513E-02 737 | 1 0.662730 738 | 2 3.923368 739 | 3 4.575097 740 | 4 4.794469 741 | 5 5.257209 742 | 6 5.505073 743 | 7 20.259500 744 | 8 25.967754 745 | 9 27.410466 746 | 10 29.114944 747 | 11 33.577665 748 | 12 38.435136 749 | 750 | 0.3333333E+00 0.1428571E+00 0.3885781E-15 0.2591513E-02 751 | 1 1.345178 752 | 2 3.904698 753 | 3 4.470744 754 | 4 4.776519 755 | 5 5.544002 756 | 6 5.606237 757 | 7 17.342204 758 | 8 25.661095 759 | 9 26.794566 760 | 10 28.373266 761 | 11 35.178546 762 | 12 37.976435 763 | 764 | 0.3809524E+00 0.1428571E+00 0.3191891E-15 0.2591513E-02 765 | 1 1.982458 766 | 2 3.968585 767 | 3 4.386206 768 | 4 4.737135 769 | 5 5.703428 770 | 6 5.940854 771 | 7 16.159098 772 | 8 25.528205 773 | 9 25.973989 774 | 10 27.961210 775 | 11 36.964067 776 | 12 37.571442 777 | 778 | 0.4285714E+00 0.1428571E+00 0.3330669E-15 0.2591513E-02 779 | 1 2.441274 780 | 2 4.100701 781 | 3 4.381036 782 | 4 4.690903 783 | 5 5.774711 784 | 6 6.507423 785 | 7 14.443682 786 | 8 25.112033 787 | 9 25.578066 788 | 10 27.786139 789 | 11 37.206156 790 | 12 37.391702 791 | 792 | 0.4761905E+00 0.1428571E+00 0.3330669E-15 0.2591513E-02 793 | 1 2.677180 794 | 2 4.143158 795 | 3 4.548766 796 | 4 4.649229 797 | 5 5.804645 798 | 6 7.277684 799 | 7 12.717880 800 | 8 24.321056 801 | 9 25.816511 802 | 10 27.769765 803 | 11 36.660939 804 | 12 36.918274 805 | 806 | -0.4761905E+00 0.1428571E+00 0.4440892E-15 0.2591513E-02 807 | 1 2.766741 808 | 2 4.143412 809 | 3 4.620004 810 | 4 4.710885 811 | 5 5.784874 812 | 6 8.048297 813 | 7 11.484503 814 | 8 23.642951 815 | 9 26.245143 816 | 10 27.892971 817 | 11 36.059559 818 | 12 36.678619 819 | 820 | -0.4285714E+00 0.1428571E+00 0.3885781E-15 0.2591513E-02 821 | 1 2.772594 822 | 2 4.163371 823 | 3 4.608651 824 | 4 4.736015 825 | 5 5.715042 826 | 6 8.072617 827 | 7 11.491467 828 | 8 23.096138 829 | 9 26.863440 830 | 10 28.157843 831 | 11 35.589906 832 | 12 36.460080 833 | 834 | -0.3809524E+00 0.1428571E+00 0.3885781E-15 0.2591513E-02 835 | 1 2.694284 836 | 2 4.211075 837 | 3 4.591760 838 | 4 4.617696 839 | 5 5.600723 840 | 6 7.370175 841 | 7 12.734989 842 | 8 22.696812 843 | 9 27.665804 844 | 10 28.571553 845 | 11 35.251629 846 | 12 36.219029 847 | 848 | -0.3333333E+00 0.1428571E+00 0.3608225E-15 0.2591513E-02 849 | 1 2.466504 850 | 2 4.281611 851 | 3 4.325257 852 | 4 4.647536 853 | 5 5.453283 854 | 6 6.704081 855 | 7 14.466314 856 | 8 22.465260 857 | 9 28.641707 858 | 10 29.140732 859 | 11 35.043734 860 | 12 35.899086 861 | 862 | -0.2857143E+00 0.1428571E+00 0.3885781E-15 0.2591513E-02 863 | 1 2.005953 864 | 2 4.095706 865 | 3 4.369606 866 | 4 4.695332 867 | 5 5.287438 868 | 6 6.246966 869 | 7 16.287920 870 | 8 22.411181 871 | 9 29.778008 872 | 10 29.868779 873 | 11 34.965027 874 | 12 35.440606 875 | 876 | -0.2380952E+00 0.1428571E+00 0.3608225E-15 0.2591513E-02 877 | 1 1.359109 878 | 2 4.007309 879 | 3 4.458503 880 | 4 4.753914 881 | 5 5.122210 882 | 6 5.931513 883 | 7 17.506184 884 | 8 22.535344 885 | 9 30.754607 886 | 10 31.059200 887 | 11 34.807305 888 | 12 35.013746 889 | 890 | -0.1904762E+00 0.1428571E+00 0.3330669E-15 0.2591513E-02 891 | 1 0.667896 892 | 2 4.021218 893 | 3 4.535196 894 | 4 4.794669 895 | 5 4.995473 896 | 6 5.692181 897 | 7 20.182832 898 | 8 22.839915 899 | 9 31.790913 900 | 10 32.465388 901 | 11 34.014964 902 | 12 35.178171 903 | 904 | -0.1428571E+00 0.1428571E+00 0.3330669E-15 0.1295756E-02 905 | 1 0.034532 906 | 2 4.085843 907 | 3 4.584375 908 | 4 4.739270 909 | 5 4.992583 910 | 6 5.505364 911 | 7 22.844548 912 | 8 23.561580 913 | 9 32.946899 914 | 10 33.145531 915 | 11 33.953900 916 | 12 35.132008 917 | 918 | 0.1904762E+00 0.1904762E+00 0.3761130E-15 0.6478782E-03 919 | 1 -0.160227 920 | 2 3.976694 921 | 3 4.816876 922 | 4 4.816878 923 | 5 4.871152 924 | 6 5.403187 925 | 7 26.200163 926 | 8 26.200151 927 | 9 27.040837 928 | 10 28.323761 929 | 11 31.351833 930 | 12 41.422209 931 | 932 | 0.2380952E+00 0.1904762E+00 0.3747003E-15 0.2591513E-02 933 | 1 0.297822 934 | 2 3.881921 935 | 3 4.743732 936 | 4 4.898411 937 | 5 4.999449 938 | 6 5.449476 939 | 7 24.230838 940 | 8 25.487682 941 | 9 27.121738 942 | 10 27.911443 943 | 11 30.669485 944 | 12 41.477309 945 | 946 | 0.2857143E+00 0.1904762E+00 0.3747003E-15 0.2591513E-02 947 | 1 0.863122 948 | 2 3.800056 949 | 3 4.623315 950 | 4 4.942029 951 | 5 5.235773 952 | 6 5.523645 953 | 7 21.674200 954 | 8 24.935623 955 | 9 26.396704 956 | 10 27.768354 957 | 11 31.317115 958 | 12 41.061593 959 | 960 | 0.3333333E+00 0.1904762E+00 0.3747003E-15 0.2591513E-02 961 | 1 1.489908 962 | 2 3.747516 963 | 3 4.505027 964 | 4 4.941857 965 | 5 5.563985 966 | 6 5.620731 967 | 7 18.542298 968 | 8 24.555575 969 | 9 25.367370 970 | 10 27.527699 971 | 11 32.771114 972 | 12 40.502689 973 | 974 | 0.3809524E+00 0.1904762E+00 0.3747003E-15 0.2591513E-02 975 | 1 2.097365 976 | 2 3.745998 977 | 3 4.405176 978 | 4 4.911517 979 | 5 5.716827 980 | 6 6.023666 981 | 7 16.494592 982 | 8 24.302055 983 | 9 24.363405 984 | 10 27.357874 985 | 11 34.544495 986 | 12 40.034614 987 | 988 | 0.4285714E+00 0.1904762E+00 0.3469447E-15 0.2591513E-02 989 | 1 2.571873 990 | 2 3.809101 991 | 3 4.351738 992 | 4 4.867463 993 | 5 5.787410 994 | 6 6.670942 995 | 7 15.424064 996 | 8 23.272500 997 | 9 24.374242 998 | 10 27.279844 999 | 11 36.470269 1000 | 12 39.338037 1001 | 1002 | 0.4761905E+00 0.1904762E+00 0.3330669E-15 0.2591513E-02 1003 | 1 2.845955 1004 | 2 3.892114 1005 | 3 4.382843 1006 | 4 4.822182 1007 | 5 5.815744 1008 | 6 7.539117 1009 | 7 13.597587 1010 | 8 22.296525 1011 | 9 24.582315 1012 | 10 27.299870 1013 | 11 38.168093 1014 | 12 38.566591 1015 | 1016 | -0.4761905E+00 0.1904762E+00 0.3885781E-15 0.2591513E-02 1017 | 1 2.969409 1018 | 2 3.937465 1019 | 3 4.468974 1020 | 4 4.784527 1021 | 5 5.793550 1022 | 6 8.531826 1023 | 7 12.061345 1024 | 8 21.410737 1025 | 9 24.982852 1026 | 10 27.427777 1027 | 11 37.921437 1028 | 12 38.382103 1029 | 1030 | -0.4285714E+00 0.1904762E+00 0.3885781E-15 0.2591513E-02 1031 | 1 3.010261 1032 | 2 3.974654 1033 | 3 4.502379 1034 | 4 4.760792 1035 | 5 5.720572 1036 | 6 9.003302 1037 | 7 11.488553 1038 | 8 20.580559 1039 | 9 25.575449 1040 | 10 27.675642 1041 | 11 37.407846 1042 | 12 38.084363 1043 | 1044 | -0.3809524E+00 0.1904762E+00 0.3885781E-15 0.2591513E-02 1045 | 1 2.984795 1046 | 2 4.034280 1047 | 3 4.417451 1048 | 4 4.754909 1049 | 5 5.602730 1050 | 6 8.297991 1051 | 7 12.542466 1052 | 8 19.910131 1053 | 9 26.354594 1054 | 10 28.055277 1055 | 11 37.028294 1056 | 12 37.691445 1057 | 1058 | -0.3333333E+00 0.1904762E+00 0.4440892E-15 0.2591513E-02 1059 | 1 2.850824 1060 | 2 4.119316 1061 | 3 4.224696 1062 | 4 4.768658 1063 | 5 5.451432 1064 | 6 7.431016 1065 | 7 14.224518 1066 | 8 19.497811 1067 | 9 27.305245 1068 | 10 28.576059 1069 | 11 36.784254 1070 | 12 37.189942 1071 | 1072 | -0.2857143E+00 0.1904762E+00 0.3885781E-15 0.2591513E-02 1073 | 1 2.513399 1074 | 2 4.025155 1075 | 3 4.232784 1076 | 4 4.802123 1077 | 5 5.280999 1078 | 6 6.768149 1079 | 7 16.071114 1080 | 8 19.340778 1081 | 9 28.411090 1082 | 10 29.242419 1083 | 11 36.530073 1084 | 12 36.677673 1085 | 1086 | -0.2380952E+00 0.1904762E+00 0.3330669E-15 0.2591513E-02 1087 | 1 1.954845 1088 | 2 3.943639 1089 | 3 4.346340 1090 | 4 4.853025 1091 | 5 5.106034 1092 | 6 6.301085 1093 | 7 17.139496 1094 | 8 19.406902 1095 | 9 29.654589 1096 | 10 30.048634 1097 | 11 35.697420 1098 | 12 36.711265 1099 | 1100 | -0.1904762E+00 0.1904762E+00 0.3608225E-15 0.1295756E-02 1101 | 1 1.299779 1102 | 2 3.957680 1103 | 3 4.452345 1104 | 4 4.912232 1105 | 5 4.946382 1106 | 6 5.958166 1107 | 7 19.338873 1108 | 8 20.068998 1109 | 9 30.959799 1110 | 10 31.015622 1111 | 11 34.769667 1112 | 12 36.821588 1113 | 1114 | 0.2380952E+00 0.2380952E+00 0.3635121E-15 0.6478782E-03 1115 | 1 0.684161 1116 | 2 3.776267 1117 | 3 4.789028 1118 | 4 5.017272 1119 | 5 5.017274 1120 | 6 5.480476 1121 | 7 24.687642 1122 | 8 24.682673 1123 | 9 26.238037 1124 | 10 26.379907 1125 | 11 29.549562 1126 | 12 43.728678 1127 | 1128 | 0.2857143E+00 0.2380952E+00 0.3642919E-15 0.2591513E-02 1129 | 1 1.167702 1130 | 2 3.681546 1131 | 3 4.730947 1132 | 4 5.100325 1133 | 5 5.211727 1134 | 6 5.536905 1135 | 7 22.899943 1136 | 8 24.028644 1137 | 9 25.038879 1138 | 10 26.682374 1139 | 11 29.478111 1140 | 12 43.677134 1141 | 1142 | 0.3333333E+00 0.2380952E+00 0.3608225E-15 0.2591513E-02 1143 | 1 1.704668 1144 | 2 3.609589 1145 | 3 4.636670 1146 | 4 5.126688 1147 | 5 5.554228 1148 | 6 5.626264 1149 | 7 20.061546 1150 | 8 23.545570 1151 | 9 23.801670 1152 | 10 26.878575 1153 | 11 30.479217 1154 | 12 43.133382 1155 | 1156 | 0.3809524E+00 0.2380952E+00 0.3885781E-15 0.2591513E-02 1157 | 1 2.228261 1158 | 2 3.573366 1159 | 3 4.542039 1160 | 4 5.109343 1161 | 5 5.723331 1162 | 6 6.064489 1163 | 7 17.207388 1164 | 8 22.552220 1165 | 9 23.287054 1166 | 10 26.867557 1167 | 11 32.080139 1168 | 12 42.318239 1169 | 1170 | 0.4285714E+00 0.2380952E+00 0.3747003E-15 0.2591513E-02 1171 | 1 2.658800 1172 | 2 3.583755 1173 | 3 4.466751 1174 | 4 5.070018 1175 | 5 5.796861 1176 | 6 6.784760 1177 | 7 16.262970 1178 | 8 21.297071 1179 | 9 23.257290 1180 | 10 26.848049 1181 | 11 33.967136 1182 | 12 41.436925 1183 | 1184 | 0.4761905E+00 0.2380952E+00 0.3608225E-15 0.2591513E-02 1185 | 1 2.946380 1186 | 2 3.635673 1187 | 3 4.425787 1188 | 4 5.023335 1189 | 5 5.827121 1190 | 6 7.734437 1191 | 7 14.659522 1192 | 8 20.023406 1193 | 9 23.445057 1194 | 10 26.884978 1195 | 11 36.025986 1196 | 12 40.564019 1197 | 1198 | -0.4761905E+00 0.2380952E+00 0.3885781E-15 0.2591513E-02 1199 | 1 3.111905 1200 | 2 3.698381 1201 | 3 4.413897 1202 | 4 4.978688 1203 | 5 5.805573 1204 | 6 8.869537 1205 | 7 12.937218 1206 | 8 18.790522 1207 | 9 23.811035 1208 | 10 27.004313 1209 | 11 38.147367 1210 | 12 39.870947 1211 | 1212 | -0.4285714E+00 0.2380952E+00 0.3885781E-15 0.2591513E-02 1213 | 1 3.203985 1214 | 2 3.757409 1215 | 3 4.391398 1216 | 4 4.942352 1217 | 5 5.732672 1218 | 6 9.859172 1219 | 7 11.748063 1220 | 8 17.701789 1221 | 9 24.368078 1222 | 10 27.222805 1223 | 11 39.303111 1224 | 12 39.579262 1225 | 1226 | -0.3809524E+00 0.2380952E+00 0.3885781E-15 0.2591513E-02 1227 | 1 3.238440 1228 | 2 3.832084 1229 | 3 4.309055 1230 | 4 4.918412 1231 | 5 5.616268 1232 | 6 9.363601 1233 | 7 12.436898 1234 | 8 17.129540 1235 | 9 25.120857 1236 | 10 27.553489 1237 | 11 38.865326 1238 | 12 39.362853 1239 | 1240 | -0.3333333E+00 0.2380952E+00 0.3885781E-15 0.2591513E-02 1241 | 1 3.188477 1242 | 2 3.937181 1243 | 3 4.153084 1244 | 4 4.908087 1245 | 5 5.469730 1246 | 6 8.322486 1247 | 7 14.070112 1248 | 8 16.895534 1249 | 9 26.045398 1250 | 10 28.006110 1251 | 11 38.551674 1252 | 12 38.721324 1253 | 1254 | -0.2857143E+00 0.2380952E+00 0.4440892E-15 0.2591513E-02 1255 | 1 2.973395 1256 | 2 3.987126 1257 | 3 4.073270 1258 | 4 4.908585 1259 | 5 5.311679 1260 | 6 7.446428 1261 | 7 15.857577 1262 | 8 16.892532 1263 | 9 27.121227 1264 | 10 28.584307 1265 | 11 37.849471 1266 | 12 38.335850 1267 | 1268 | -0.2380952E+00 0.2380952E+00 0.4163336E-15 0.1295756E-02 1269 | 1 2.533188 1270 | 2 3.923011 1271 | 3 4.212206 1272 | 4 4.904757 1273 | 5 5.172005 1274 | 6 6.787825 1275 | 7 16.738970 1276 | 8 17.166376 1277 | 9 28.331037 1278 | 10 29.278787 1279 | 11 36.794912 1280 | 12 38.007501 1281 | 1282 | 0.2857143E+00 0.2857143E+00 0.3635121E-15 0.6478782E-03 1283 | 1 1.546524 1284 | 2 3.574621 1285 | 3 4.878390 1286 | 4 5.239130 1287 | 5 5.239131 1288 | 6 5.560441 1289 | 7 23.252418 1290 | 8 23.236189 1291 | 9 23.575042 1292 | 10 25.847353 1293 | 11 28.602935 1294 | 12 45.591279 1295 | 1296 | 0.3333333E+00 0.2857143E+00 0.3642919E-15 0.2591513E-02 1297 | 1 1.955563 1298 | 2 3.488726 1299 | 3 4.916172 1300 | 4 5.312588 1301 | 5 5.511997 1302 | 6 5.625182 1303 | 7 21.450679 1304 | 8 22.201601 1305 | 9 22.615952 1306 | 10 26.107359 1307 | 11 28.763255 1308 | 12 45.420707 1309 | 1310 | 0.3809524E+00 0.2857143E+00 0.3955170E-15 0.2591513E-02 1311 | 1 2.343511 1312 | 2 3.433680 1313 | 3 4.852687 1314 | 4 5.315194 1315 | 5 5.721685 1316 | 6 6.068144 1317 | 7 18.648120 1318 | 8 20.700376 1319 | 9 22.251366 1320 | 10 26.373687 1321 | 11 29.875471 1322 | 12 44.539712 1323 | 1324 | 0.4285714E+00 0.2857143E+00 0.3885781E-15 0.2591513E-02 1325 | 1 2.669451 1326 | 2 3.417138 1327 | 3 4.763266 1328 | 4 5.283065 1329 | 5 5.801904 1330 | 6 6.864409 1331 | 7 16.813843 1332 | 8 19.267667 1333 | 9 22.213015 1334 | 10 26.458123 1335 | 11 31.559497 1336 | 12 43.504510 1337 | 1338 | 0.4761905E+00 0.2857143E+00 0.3608225E-15 0.2591513E-02 1339 | 1 2.919062 1340 | 2 3.440353 1341 | 3 4.668072 1342 | 4 5.237114 1343 | 5 5.838022 1344 | 6 7.883309 1345 | 7 15.840592 1346 | 8 18.026912 1347 | 9 22.410250 1348 | 10 26.524506 1349 | 11 33.536105 1350 | 12 42.634631 1351 | 1352 | -0.4761905E+00 0.2857143E+00 0.3885781E-15 0.2591513E-02 1353 | 1 3.105726 1354 | 2 3.492730 1355 | 3 4.570503 1356 | 4 5.187210 1357 | 5 5.820765 1358 | 6 9.093324 1359 | 7 14.064914 1360 | 8 17.113794 1361 | 9 22.733624 1362 | 10 26.639100 1363 | 11 35.702616 1364 | 12 41.874770 1365 | 1366 | -0.4285714E+00 0.2857143E+00 0.3885781E-15 0.2591513E-02 1367 | 1 3.253322 1368 | 2 3.558225 1369 | 3 4.462151 1370 | 4 5.138422 1371 | 5 5.752248 1372 | 6 10.398320 1373 | 7 12.483881 1374 | 8 16.631047 1375 | 9 23.241976 1376 | 10 26.829733 1377 | 11 38.007162 1378 | 12 41.229189 1379 | 1380 | -0.3809524E+00 0.2857143E+00 0.3885781E-15 0.2591513E-02 1381 | 1 3.365408 1382 | 2 3.637367 1383 | 3 4.328182 1384 | 4 5.092656 1385 | 5 5.643878 1386 | 6 10.562566 1387 | 7 12.369924 1388 | 8 16.406611 1389 | 9 23.979150 1390 | 10 27.112743 1391 | 11 40.302490 1392 | 12 40.684939 1393 | 1394 | -0.3333333E+00 0.2857143E+00 0.3885781E-15 0.2591513E-02 1395 | 1 3.401770 1396 | 2 3.752746 1397 | 3 4.166098 1398 | 4 5.047022 1399 | 5 5.514969 1400 | 6 9.386245 1401 | 7 13.886449 1402 | 8 16.343488 1403 | 9 24.889393 1404 | 10 27.498691 1405 | 11 40.185559 1406 | 12 40.355235 1407 | 1408 | -0.2857143E+00 0.2857143E+00 0.3885781E-15 0.1295756E-02 1409 | 1 3.289742 1410 | 2 3.905277 1411 | 3 4.028318 1412 | 4 4.990789 1413 | 5 5.393554 1414 | 6 8.315803 1415 | 7 15.313151 1416 | 8 16.507754 1417 | 9 25.941598 1418 | 10 27.991778 1419 | 11 39.192865 1420 | 12 39.544100 1421 | 1422 | 0.3333333E+00 0.3333333E+00 0.3635121E-15 0.6478782E-03 1423 | 1 2.207043 1424 | 2 3.389607 1425 | 3 5.358645 1426 | 4 5.464340 1427 | 5 5.464342 1428 | 6 5.635707 1429 | 7 20.497928 1430 | 8 21.769544 1431 | 9 21.743739 1432 | 10 25.437311 1433 | 11 28.056096 1434 | 12 47.594491 1435 | 1436 | 0.3809524E+00 0.3333333E+00 0.3642919E-15 0.2591513E-02 1437 | 1 2.435904 1438 | 2 3.320033 1439 | 3 5.324910 1440 | 4 5.510696 1441 | 5 5.711249 1442 | 6 6.071573 1443 | 7 18.938170 1444 | 8 20.092690 1445 | 9 21.272922 1446 | 10 25.701584 1447 | 11 28.308519 1448 | 12 46.716630 1449 | 1450 | 0.4285714E+00 0.3333333E+00 0.3538836E-15 0.2591513E-02 1451 | 1 2.636857 1452 | 2 3.287459 1453 | 3 5.198136 1454 | 4 5.486946 1455 | 5 5.802016 1456 | 6 6.950598 1457 | 7 17.791777 1458 | 8 17.905763 1459 | 9 21.221395 1460 | 10 26.032915 1461 | 11 29.463799 1462 | 12 45.709972 1463 | 1464 | 0.4761905E+00 0.3333333E+00 0.4024558E-15 0.2591513E-02 1465 | 1 2.815321 1466 | 2 3.295025 1467 | 3 5.049583 1468 | 4 5.443671 1469 | 5 5.848555 1470 | 6 8.015263 1471 | 7 16.740660 1472 | 8 16.948650 1473 | 9 21.460520 1474 | 10 26.193455 1475 | 11 31.204960 1476 | 12 44.736612 1477 | 1478 | -0.4761905E+00 0.3333333E+00 0.4440892E-15 0.2591513E-02 1479 | 1 2.980606 1480 | 2 3.339281 1481 | 3 4.882560 1482 | 4 5.390762 1483 | 5 5.839293 1484 | 6 9.239596 1485 | 7 15.369588 1486 | 8 16.476750 1487 | 9 21.715390 1488 | 10 26.327623 1489 | 11 33.256118 1490 | 12 43.888488 1491 | 1492 | -0.4285714E+00 0.3333333E+00 0.3885781E-15 0.2591513E-02 1493 | 1 3.141950 1494 | 2 3.409844 1495 | 3 4.698180 1496 | 4 5.330207 1497 | 5 5.780100 1498 | 6 10.584791 1499 | 7 13.664109 1500 | 8 16.031565 1501 | 9 22.150188 1502 | 10 26.503215 1503 | 11 35.509748 1504 | 12 43.124013 1505 | 1506 | -0.3809524E+00 0.3333333E+00 0.3885781E-15 0.2591513E-02 1507 | 1 3.301484 1508 | 2 3.494015 1509 | 3 4.498899 1510 | 4 5.258112 1511 | 5 5.688551 1512 | 6 11.781144 1513 | 7 12.301704 1514 | 8 15.690886 1515 | 9 22.891033 1516 | 10 26.748480 1517 | 11 37.920565 1518 | 12 42.355973 1519 | 1520 | -0.3333333E+00 0.3333333E+00 0.3885781E-15 0.1295756E-02 1521 | 1 3.424297 1522 | 2 3.596423 1523 | 3 4.302021 1524 | 4 5.166188 1525 | 5 5.591431 1526 | 6 10.619691 1527 | 7 13.284212 1528 | 8 15.754492 1529 | 9 23.838018 1530 | 10 27.078135 1531 | 11 40.454539 1532 | 12 41.457199 1533 | 1534 | 0.3809524E+00 0.3809524E+00 0.3635121E-15 0.6478782E-03 1535 | 1 2.518581 1536 | 2 3.236257 1537 | 3 5.668998 1538 | 4 5.668998 1539 | 5 5.699291 1540 | 6 6.353794 1541 | 7 17.687101 1542 | 8 20.436826 1543 | 9 20.414398 1544 | 10 25.140323 1545 | 11 27.723916 1546 | 12 48.503293 1547 | 1548 | 0.4285714E+00 0.3809524E+00 0.3642919E-15 0.2591513E-02 1549 | 1 2.606908 1550 | 2 3.189148 1551 | 3 5.622951 1552 | 4 5.658558 1553 | 5 5.799763 1554 | 6 7.154038 1555 | 7 16.888350 1556 | 8 19.311501 1557 | 9 20.369015 1558 | 10 25.426034 1559 | 11 28.029148 1560 | 12 47.744690 1561 | 1562 | 0.4761905E+00 0.3809524E+00 0.3608225E-15 0.2591513E-02 1563 | 1 2.710635 1564 | 2 3.183249 1565 | 3 5.461311 1566 | 4 5.618936 1567 | 5 5.860717 1568 | 6 8.173120 1569 | 7 16.211063 1570 | 8 17.918206 1571 | 9 20.627767 1572 | 10 25.817645 1573 | 11 29.211106 1574 | 12 46.774196 1575 | 1576 | -0.4761905E+00 0.3809524E+00 0.3885781E-15 0.2591513E-02 1577 | 1 2.834417 1578 | 2 3.218564 1579 | 3 5.250641 1580 | 4 5.566447 1581 | 5 5.862548 1582 | 6 9.311336 1583 | 7 15.535104 1584 | 8 16.747770 1585 | 9 20.735805 1586 | 10 26.043330 1587 | 11 31.001290 1588 | 12 45.828471 1589 | 1590 | -0.4285714E+00 0.3809524E+00 0.3885781E-15 0.2591513E-02 1591 | 1 2.979413 1592 | 2 3.289681 1593 | 3 5.007418 1594 | 4 5.495552 1595 | 5 5.817435 1596 | 6 10.483392 1597 | 7 14.892965 1598 | 8 15.300601 1599 | 9 21.015240 1600 | 10 26.234137 1601 | 11 33.119330 1602 | 12 44.858256 1603 | 1604 | -0.3809524E+00 0.3809524E+00 0.4440892E-15 0.1295756E-02 1605 | 1 3.142082 1606 | 2 3.385505 1607 | 3 4.747027 1608 | 4 5.395221 1609 | 5 5.750469 1610 | 6 11.498335 1611 | 7 13.461964 1612 | 8 15.074611 1613 | 9 21.826534 1614 | 10 26.459716 1615 | 11 35.445914 1616 | 12 43.744870 1617 | 1618 | 0.4285714E+00 0.4285714E+00 0.3635121E-15 0.6478782E-03 1619 | 1 2.604738 1620 | 2 3.127038 1621 | 3 5.745347 1622 | 4 5.826885 1623 | 5 5.826886 1624 | 6 7.645696 1625 | 7 16.115715 1626 | 8 19.820417 1627 | 9 19.819854 1628 | 10 24.947792 1629 | 11 27.533465 1630 | 12 48.896224 1631 | 1632 | 0.4761905E+00 0.4285714E+00 0.3608225E-15 0.2591513E-02 1633 | 1 2.640272 1634 | 2 3.106389 1635 | 3 5.732910 1636 | 4 5.782145 1637 | 5 5.882842 1638 | 6 8.408891 1639 | 7 15.316994 1640 | 8 19.353004 1641 | 9 20.074013 1642 | 10 25.265866 1643 | 11 27.879034 1644 | 12 48.466869 1645 | 1646 | -0.4761905E+00 0.4285714E+00 0.3885781E-15 0.2591513E-02 1647 | 1 2.716234 1648 | 2 3.129175 1649 | 3 5.591748 1650 | 4 5.690563 1651 | 5 5.892130 1652 | 6 9.268839 1653 | 7 14.754561 1654 | 8 18.174747 1655 | 9 20.031105 1656 | 10 25.713060 1657 | 11 29.091798 1658 | 12 47.518102 1659 | 1660 | -0.4285714E+00 0.4285714E+00 0.3885781E-15 0.1295756E-02 1661 | 1 2.831213 1662 | 2 3.192700 1663 | 3 5.319960 1664 | 4 5.616080 1665 | 5 5.860980 1666 | 6 10.028201 1667 | 7 14.608215 1668 | 8 16.700117 1669 | 9 20.153768 1670 | 10 25.994613 1671 | 11 30.935106 1672 | 12 46.326768 1673 | 1674 | 0.4761905E+00 0.4761905E+00 0.3887138E-15 0.6478782E-03 1675 | 1 2.616252 1676 | 2 3.070269 1677 | 3 5.769516 1678 | 4 5.913344 1679 | 5 5.913345 1680 | 6 8.735480 1681 | 7 14.662285 1682 | 8 20.163567 1683 | 9 20.163074 1684 | 10 24.853230 1685 | 11 27.447058 1686 | 12 49.180076 1687 | 1688 | -0.4761905E+00 0.4761905E+00 0.3885781E-15 0.1295756E-02 1689 | 1 2.641795 1690 | 2 3.078065 1691 | 3 5.753873 1692 | 4 5.834620 1693 | 5 5.917416 1694 | 6 9.064675 1695 | 7 14.446049 1696 | 8 19.644387 1697 | 9 20.073088 1698 | 10 25.213293 1699 | 11 27.831891 1700 | 12 48.732523 1701 | 1702 | 0.1428571E+00 0.9523810E-01 0.4761905E-01 0.2591513E-02 1703 | 1 -1.294658 1704 | 2 4.293767 1705 | 3 4.532047 1706 | 4 4.538366 1707 | 5 5.128728 1708 | 6 5.273690 1709 | 7 27.413334 1710 | 8 28.330522 1711 | 9 30.862747 1712 | 10 32.227121 1713 | 11 35.968415 1714 | 12 36.207297 1715 | 1716 | 0.1904762E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1717 | 1 -0.899337 1718 | 2 4.220446 1719 | 3 4.585803 1720 | 4 4.594467 1721 | 5 5.089157 1722 | 6 5.303383 1723 | 7 25.284327 1724 | 8 27.362817 1725 | 9 30.824265 1726 | 10 32.783888 1727 | 11 34.859710 1728 | 12 36.597687 1729 | 1730 | 0.2380952E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1731 | 1 -0.363557 1732 | 2 4.145745 1733 | 3 4.605219 1734 | 4 4.647504 1735 | 5 5.089655 1736 | 6 5.363156 1737 | 7 22.780350 1738 | 8 26.827494 1739 | 9 30.498282 1740 | 10 33.304229 1741 | 11 33.750180 1742 | 12 36.512266 1743 | 1744 | 0.2857143E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1745 | 1 0.286933 1746 | 2 4.094641 1747 | 3 4.580208 1748 | 4 4.658467 1749 | 5 5.164827 1750 | 6 5.465771 1751 | 7 19.742648 1752 | 8 26.433159 1753 | 9 30.083627 1754 | 10 32.197224 1755 | 11 34.134539 1756 | 12 35.957510 1757 | 1758 | 0.3333333E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1759 | 1 1.000610 1760 | 2 4.108720 1761 | 3 4.522075 1762 | 4 4.629342 1763 | 5 5.304017 1764 | 6 5.609430 1765 | 7 17.197794 1766 | 8 26.166175 1767 | 9 29.609024 1768 | 10 31.102247 1769 | 11 34.176573 1770 | 12 35.668156 1771 | 1772 | 0.3809524E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1773 | 1 1.676627 1774 | 2 4.254291 1775 | 3 4.453716 1776 | 4 4.582454 1777 | 5 5.469431 1778 | 6 5.799376 1779 | 7 15.991722 1780 | 8 26.030895 1781 | 9 29.038469 1782 | 10 30.246208 1783 | 11 33.909536 1784 | 12 35.501308 1785 | 1786 | 0.4285714E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1787 | 1 2.173099 1788 | 2 4.387386 1789 | 3 4.525867 1790 | 4 4.593580 1791 | 5 5.619843 1792 | 6 6.093890 1793 | 7 14.143349 1794 | 8 26.032995 1795 | 9 28.363966 1796 | 10 29.686636 1797 | 11 33.635353 1798 | 12 35.287122 1799 | 1800 | 0.4761905E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1801 | 1 2.433928 1802 | 2 4.341274 1803 | 3 4.494140 1804 | 4 4.960084 1805 | 5 5.721889 1806 | 6 6.598869 1807 | 7 12.465536 1808 | 8 26.170537 1809 | 9 27.697818 1810 | 10 29.361408 1811 | 11 33.451018 1812 | 12 35.074572 1813 | 1814 | -0.4761905E+00 0.9523810E-01 0.4761905E-01 0.5183026E-02 1815 | 1 2.531805 1816 | 2 4.316344 1817 | 3 4.473834 1818 | 4 5.165785 1819 | 5 5.770182 1820 | 6 7.170809 1821 | 7 11.327666 1822 | 8 26.394980 1823 | 9 27.189120 1824 | 10 29.208157 1825 | 11 33.363761 1826 | 12 34.945154 1827 | 1828 | 0.1904762E+00 0.1428571E+00 0.4761905E-01 0.2591513E-02 1829 | 1 -0.768506 1830 | 2 4.156130 1831 | 3 4.635327 1832 | 4 4.653041 1833 | 5 5.024722 1834 | 6 5.336133 1835 | 7 26.058835 1836 | 8 26.730985 1837 | 9 29.728661 1838 | 10 30.885488 1839 | 11 34.160544 1840 | 12 37.857720 1841 | 1842 | 0.2380952E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1843 | 1 -0.285578 1844 | 2 4.073949 1845 | 3 4.645417 1846 | 4 4.723056 1847 | 5 5.026827 1848 | 6 5.392621 1849 | 7 23.846950 1850 | 8 25.891106 1851 | 9 29.657897 1852 | 10 31.111811 1853 | 11 33.353888 1854 | 12 37.748240 1855 | 1856 | 0.2857143E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1857 | 1 0.318796 1858 | 2 4.003898 1859 | 3 4.592831 1860 | 4 4.755096 1861 | 5 5.125830 1862 | 6 5.490220 1863 | 7 21.127315 1864 | 8 25.423444 1865 | 9 29.167553 1866 | 10 30.616551 1867 | 11 33.984062 1868 | 12 37.147959 1869 | 1870 | 0.3333333E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1871 | 1 1.002945 1872 | 2 3.973248 1873 | 3 4.506274 1874 | 4 4.742124 1875 | 5 5.299905 1876 | 6 5.637383 1877 | 7 17.818265 1878 | 8 25.103029 1879 | 9 28.451955 1880 | 10 29.901194 1881 | 11 35.308264 1882 | 12 36.474438 1883 | 1884 | 0.3809524E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1885 | 1 1.684786 1886 | 2 4.024210 1887 | 3 4.417453 1888 | 4 4.703488 1889 | 5 5.495770 1890 | 6 5.860233 1891 | 7 16.465116 1892 | 8 24.926537 1893 | 9 27.555402 1894 | 10 29.438341 1895 | 11 35.815694 1896 | 12 36.641177 1897 | 1898 | 0.4285714E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1899 | 1 2.230521 1900 | 2 4.182883 1901 | 3 4.368008 1902 | 4 4.657438 1903 | 5 5.650496 1904 | 6 6.239694 1905 | 7 14.973384 1906 | 8 24.896494 1907 | 9 26.649089 1908 | 10 29.175609 1909 | 11 35.252536 1910 | 12 37.230293 1911 | 1912 | 0.4761905E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1913 | 1 2.544699 1914 | 2 4.229757 1915 | 3 4.540028 1916 | 4 4.625341 1917 | 5 5.740935 1918 | 6 6.847423 1919 | 7 13.168241 1920 | 8 24.989471 1921 | 9 25.864872 1922 | 10 29.059158 1923 | 11 34.777117 1924 | 12 37.157486 1925 | 1926 | -0.4761905E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1927 | 1 2.675658 1928 | 2 4.206053 1929 | 3 4.577481 1930 | 4 4.780460 1931 | 5 5.773991 1932 | 6 7.595033 1933 | 7 11.740915 1934 | 8 24.895860 1935 | 9 25.556265 1936 | 10 29.080546 1937 | 11 34.399699 1938 | 12 36.976647 1939 | 1940 | -0.4285714E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1941 | 1 2.703771 1942 | 2 4.204300 1943 | 3 4.567154 1944 | 4 4.852961 1945 | 5 5.756147 1946 | 6 7.917411 1947 | 7 11.261718 1948 | 8 24.424189 1949 | 9 25.931015 1950 | 10 29.245759 1951 | 11 34.123472 1952 | 12 36.821829 1953 | 1954 | -0.3809524E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1955 | 1 2.650917 1956 | 2 4.230991 1957 | 3 4.573295 1958 | 4 4.749661 1959 | 5 5.691794 1960 | 6 7.382636 1961 | 7 12.157671 1962 | 8 24.027238 1963 | 9 26.539678 1964 | 10 29.560429 1965 | 11 33.953499 1966 | 12 36.693593 1967 | 1968 | -0.3333333E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1969 | 1 2.472896 1970 | 2 4.281899 1971 | 3 4.473650 1972 | 4 4.615921 1973 | 5 5.581739 1974 | 6 6.696017 1975 | 7 13.750001 1976 | 8 23.776479 1977 | 9 27.306099 1978 | 10 30.021683 1979 | 11 33.900008 1980 | 12 36.542401 1981 | 1982 | -0.2857143E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1983 | 1 2.081071 1984 | 2 4.186752 1985 | 3 4.364379 1986 | 4 4.654522 1987 | 5 5.432108 1988 | 6 6.212109 1989 | 7 15.619965 1990 | 8 23.684145 1991 | 9 28.204940 1992 | 10 30.609842 1993 | 11 33.981170 1994 | 12 36.284417 1995 | 1996 | -0.2380952E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 1997 | 1 1.471282 1998 | 2 4.038520 1999 | 3 4.446114 2000 | 4 4.705744 2001 | 5 5.263870 2002 | 6 5.898126 2003 | 7 16.883523 2004 | 8 23.749957 2005 | 9 29.210461 2006 | 10 31.275596 2007 | 11 34.223646 2008 | 12 35.826141 2009 | 2010 | -0.1904762E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 2011 | 1 0.772709 2012 | 2 4.016002 2013 | 3 4.527245 2014 | 4 4.748719 2015 | 5 5.111251 2016 | 6 5.672238 2017 | 7 18.879630 2018 | 8 23.968972 2019 | 9 30.281014 2020 | 10 31.916252 2021 | 11 34.594486 2022 | 12 35.233793 2023 | 2024 | -0.1428571E+00 0.1428571E+00 0.4761905E-01 0.5183026E-02 2025 | 1 0.106303 2026 | 2 4.063366 2027 | 3 4.587799 2028 | 4 4.744022 2029 | 5 5.026158 2030 | 6 5.500455 2031 | 7 22.095348 2032 | 8 24.345080 2033 | 9 31.314214 2034 | 10 32.293630 2035 | 11 34.375290 2036 | 12 35.555477 2037 | 2038 | -0.9523810E-01 0.1428571E+00 0.4761905E-01 0.2591513E-02 2039 | 1 -0.462025 2040 | 2 4.139947 2041 | 3 4.608710 2042 | 4 4.673294 2043 | 5 5.036822 2044 | 6 5.378647 2045 | 7 24.411336 2046 | 8 25.118916 2047 | 9 31.673601 2048 | 10 32.405212 2049 | 11 34.465126 2050 | 12 36.466372 2051 | 2052 | 0.2380952E+00 0.1904762E+00 0.4761905E-01 0.2591513E-02 2053 | 1 -0.061907 2054 | 2 3.981404 2055 | 3 4.707411 2056 | 4 4.817819 2057 | 5 4.974568 2058 | 6 5.417870 2059 | 7 24.738076 2060 | 8 25.232534 2061 | 9 28.643686 2062 | 10 29.049935 2063 | 11 31.800729 2064 | 12 40.307354 2065 | 2066 | 0.2857143E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2067 | 1 0.485143 2068 | 2 3.897623 2069 | 3 4.622951 2070 | 4 4.884134 2071 | 5 5.100410 2072 | 6 5.505295 2073 | 7 22.436640 2074 | 8 24.508466 2075 | 9 28.024728 2076 | 10 29.093083 2077 | 11 31.975980 2078 | 12 39.538510 2079 | 2080 | 0.3333333E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2081 | 1 1.119603 2082 | 2 3.837908 2083 | 3 4.512262 2084 | 4 4.892518 2085 | 5 5.307211 2086 | 6 5.651334 2087 | 7 19.263731 2088 | 8 24.110938 2089 | 9 26.941191 2090 | 10 28.933690 2091 | 11 33.206138 2092 | 12 38.751480 2093 | 2094 | 0.3809524E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2095 | 1 1.774610 2096 | 2 3.827467 2097 | 3 4.405028 2098 | 4 4.864759 2099 | 5 5.525080 2100 | 6 5.905373 2101 | 7 16.599713 2102 | 8 23.873160 2103 | 9 25.845516 2104 | 10 28.724488 2105 | 11 34.825584 2106 | 12 38.047165 2107 | 2108 | 0.4285714E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2109 | 1 2.336927 2110 | 2 3.893181 2111 | 3 4.328684 2112 | 4 4.821034 2113 | 5 5.676520 2114 | 6 6.362524 2115 | 7 15.857978 2116 | 8 23.788476 2117 | 9 24.838343 2118 | 10 28.588142 2119 | 11 36.294788 2120 | 12 37.775583 2121 | 2122 | 0.4761905E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2123 | 1 2.695114 2124 | 2 3.995380 2125 | 3 4.344538 2126 | 4 4.774958 2127 | 5 5.758561 2128 | 6 7.062824 2129 | 7 14.041331 2130 | 8 23.662239 2131 | 9 24.152317 2132 | 10 28.553307 2133 | 11 36.318818 2134 | 12 39.174134 2135 | 2136 | -0.4761905E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2137 | 1 2.860889 2138 | 2 4.027788 2139 | 3 4.476473 2140 | 4 4.735067 2141 | 5 5.782791 2142 | 6 7.959950 2143 | 7 12.394325 2144 | 8 22.993820 2145 | 9 24.308418 2146 | 10 28.638019 2147 | 11 35.877178 2148 | 12 39.064126 2149 | 2150 | -0.4285714E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2151 | 1 2.917381 2152 | 2 4.038124 2153 | 3 4.573850 2154 | 4 4.706554 2155 | 5 5.754224 2156 | 6 8.666313 2157 | 7 11.407538 2158 | 8 22.295345 2159 | 9 24.750027 2160 | 10 28.859699 2161 | 11 35.464502 2162 | 12 38.738031 2163 | 2164 | -0.3809524E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2165 | 1 2.904174 2166 | 2 4.070901 2167 | 3 4.534586 2168 | 4 4.702866 2169 | 5 5.677955 2170 | 6 8.298379 2171 | 7 11.962334 2172 | 8 21.702738 2173 | 9 25.364320 2174 | 10 29.234374 2175 | 11 35.123131 2176 | 12 38.434406 2177 | 2178 | -0.3333333E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2179 | 1 2.799846 2180 | 2 4.131197 2181 | 3 4.355214 2182 | 4 4.723149 2183 | 5 5.559775 2184 | 6 7.455831 2185 | 7 13.475653 2186 | 8 21.273602 2187 | 9 26.132756 2188 | 10 29.774575 2189 | 11 34.849943 2190 | 12 38.107036 2191 | 2192 | -0.2857143E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2193 | 1 2.521702 2194 | 2 4.103262 2195 | 3 4.238309 2196 | 4 4.759832 2197 | 5 5.406257 2198 | 6 6.764127 2199 | 7 15.311388 2200 | 8 21.047566 2201 | 9 27.038085 2202 | 10 30.486365 2203 | 11 34.616556 2204 | 12 37.684912 2205 | 2206 | -0.2380952E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2207 | 1 2.009949 2208 | 2 3.967027 2209 | 3 4.333892 2210 | 4 4.810668 2211 | 5 5.231973 2212 | 6 6.286313 2213 | 7 16.631325 2214 | 8 21.038381 2215 | 9 28.065499 2216 | 10 31.365579 2217 | 11 34.345618 2218 | 12 37.109618 2219 | 2220 | -0.1904762E+00 0.1904762E+00 0.4761905E-01 0.5183026E-02 2221 | 1 1.355703 2222 | 2 3.942179 2223 | 3 4.437350 2224 | 4 4.868375 2225 | 5 5.059895 2226 | 6 5.949918 2227 | 7 18.197502 2228 | 8 21.272635 2229 | 9 29.199135 2230 | 10 32.385823 2231 | 11 33.876381 2232 | 12 36.468803 2233 | 2234 | -0.1428571E+00 0.1904762E+00 0.4761905E-01 0.2591513E-02 2235 | 1 0.693548 2236 | 2 3.987086 2237 | 3 4.527679 2238 | 4 4.863604 2239 | 5 4.971319 2240 | 6 5.695363 2241 | 7 21.221820 2242 | 8 22.004398 2243 | 9 30.404904 2244 | 10 32.946855 2245 | 11 33.620920 2246 | 12 35.902588 2247 | 2248 | 0.2857143E+00 0.2380952E+00 0.4761905E-01 0.2591513E-02 2249 | 1 0.775075 2250 | 2 3.787133 2251 | 3 4.678208 2252 | 4 5.017334 2253 | 5 5.104395 2254 | 6 5.514516 2255 | 7 23.446985 2256 | 8 23.824767 2257 | 9 26.554935 2258 | 10 27.729944 2259 | 11 30.311963 2260 | 12 42.023420 2261 | 2262 | 0.3333333E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2263 | 1 1.336501 2264 | 2 3.709777 2265 | 3 4.574093 2266 | 4 5.069986 2267 | 5 5.310613 2268 | 6 5.649574 2269 | 7 20.881286 2270 | 8 23.180288 2271 | 9 25.338757 2272 | 10 28.084143 2273 | 11 30.995379 2274 | 12 41.039838 2275 | 2276 | 0.3809524E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2277 | 1 1.924549 2278 | 2 3.665158 2279 | 3 4.466197 2280 | 4 5.057964 2281 | 5 5.543142 2282 | 6 5.926842 2283 | 7 17.721676 2284 | 8 22.846844 2285 | 9 24.143559 2286 | 10 28.084906 2287 | 11 32.450274 2288 | 12 40.175676 2289 | 2290 | 0.4285714E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2291 | 1 2.450941 2292 | 2 3.670185 2293 | 3 4.376578 2294 | 4 5.018792 2295 | 5 5.693787 2296 | 6 6.451209 2297 | 7 16.370245 2298 | 8 22.622986 2299 | 9 23.107027 2300 | 10 28.030785 2301 | 11 34.234309 2302 | 12 39.390854 2303 | 2304 | 0.4761905E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2305 | 1 2.823519 2306 | 2 3.726626 2307 | 3 4.333591 2308 | 4 4.969867 2309 | 5 5.773403 2310 | 6 7.232279 2311 | 7 15.070206 2312 | 8 21.785512 2313 | 9 22.915496 2314 | 10 28.033854 2315 | 11 36.115245 2316 | 12 38.799760 2317 | 2318 | -0.4761905E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2319 | 1 3.026926 2320 | 2 3.789770 2321 | 3 4.356157 2322 | 4 4.921191 2323 | 5 5.794783 2324 | 6 8.239237 2325 | 7 13.285517 2326 | 8 20.692579 2327 | 9 23.172611 2328 | 10 28.132695 2329 | 11 37.295045 2330 | 12 39.133878 2331 | 2332 | -0.4285714E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2333 | 1 3.123404 2334 | 2 3.829742 2335 | 3 4.398391 2336 | 4 4.880568 2337 | 5 5.761623 2338 | 6 9.283725 2339 | 7 11.872702 2340 | 8 19.605236 2341 | 9 23.634139 2342 | 10 28.350179 2343 | 11 37.092447 2344 | 12 40.519530 2345 | 2346 | -0.3809524E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2347 | 1 3.155271 2348 | 2 3.876009 2349 | 3 4.375578 2350 | 4 4.857959 2351 | 5 5.679731 2352 | 6 9.326739 2353 | 7 11.877137 2354 | 8 18.771783 2355 | 9 24.236106 2356 | 10 28.704766 2357 | 11 36.691438 2358 | 12 40.265550 2359 | 2360 | -0.3333333E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2361 | 1 3.116516 2362 | 2 3.951376 2363 | 3 4.247412 2364 | 4 4.859581 2365 | 5 5.557264 2366 | 6 8.368419 2367 | 7 13.306881 2368 | 8 18.183654 2369 | 9 24.995675 2370 | 10 29.211401 2371 | 11 36.299853 2372 | 12 39.829684 2373 | 2374 | -0.2857143E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2375 | 1 2.946308 2376 | 2 4.021301 2377 | 3 4.094954 2378 | 4 4.881258 2379 | 5 5.404050 2380 | 6 7.470981 2381 | 7 15.100316 2382 | 8 17.937153 2383 | 9 25.889810 2384 | 10 29.879024 2385 | 11 35.890868 2386 | 12 39.267942 2387 | 2388 | -0.2380952E+00 0.2380952E+00 0.4761905E-01 0.5183026E-02 2389 | 1 2.555929 2390 | 2 3.918518 2391 | 3 4.194583 2392 | 4 4.915912 2393 | 5 5.235202 2394 | 6 6.796405 2395 | 7 16.406728 2396 | 8 17.961584 2397 | 9 26.901352 2398 | 10 30.704941 2399 | 11 35.353345 2400 | 12 38.521768 2401 | 2402 | -0.1904762E+00 0.2380952E+00 0.4761905E-01 0.2591513E-02 2403 | 1 1.982367 2404 | 2 3.897900 2405 | 3 4.321527 2406 | 4 4.944673 2407 | 5 5.083183 2408 | 6 6.310966 2409 | 7 17.495154 2410 | 8 18.449822 2411 | 9 28.014794 2412 | 10 31.650097 2413 | 11 34.568443 2414 | 12 37.605394 2415 | 2416 | 0.3333333E+00 0.2857143E+00 0.4761905E-01 0.2591513E-02 2417 | 1 1.627751 2418 | 2 3.591356 2419 | 3 4.719129 2420 | 4 5.238643 2421 | 5 5.337916 2422 | 6 5.633935 2423 | 7 22.110636 2424 | 8 22.402742 2425 | 9 23.776164 2426 | 10 27.005340 2427 | 11 29.417037 2428 | 12 43.322989 2429 | 2430 | 0.3809524E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2431 | 1 2.102497 2432 | 2 3.527894 2433 | 3 4.654176 2434 | 4 5.271482 2435 | 5 5.548298 2436 | 6 5.923453 2437 | 7 19.331972 2438 | 8 21.794892 2439 | 9 22.492084 2440 | 10 27.409129 2441 | 11 30.284275 2442 | 12 42.407427 2443 | 2444 | 0.4285714E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2445 | 1 2.528614 2446 | 2 3.501875 2447 | 3 4.573737 2448 | 4 5.238351 2449 | 5 5.701253 2450 | 6 6.509977 2451 | 7 16.972533 2452 | 8 20.949159 2453 | 9 21.742958 2454 | 10 27.508409 2455 | 11 31.855019 2456 | 12 41.509440 2457 | 2458 | 0.4761905E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2459 | 1 2.856079 2460 | 2 3.518731 2461 | 3 4.500148 2462 | 4 5.186959 2463 | 5 5.784910 2464 | 6 7.364504 2465 | 7 16.134091 2466 | 8 19.600925 2467 | 9 21.782513 2468 | 10 27.558245 2469 | 11 33.755695 2470 | 12 40.690456 2471 | 2472 | -0.4761905E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2473 | 1 3.077207 2474 | 2 3.566828 2475 | 3 4.443559 2476 | 4 5.130221 2477 | 5 5.809090 2478 | 6 8.447830 2479 | 7 14.385575 2480 | 8 18.339139 2481 | 9 22.057096 2482 | 10 27.662575 2483 | 11 35.829969 2484 | 12 39.991830 2485 | 2486 | -0.4285714E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2487 | 1 3.225113 2488 | 2 3.620603 2489 | 3 4.395043 2490 | 4 5.076151 2491 | 5 5.777256 2492 | 6 9.682835 2493 | 7 12.724977 2494 | 8 17.381633 2495 | 9 22.512605 2496 | 10 27.859958 2497 | 11 37.859596 2498 | 12 39.574502 2499 | 2500 | -0.3809524E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2501 | 1 3.325614 2502 | 2 3.676637 2503 | 3 4.323090 2504 | 4 5.033127 2505 | 5 5.696714 2506 | 6 10.414300 2507 | 7 11.929386 2508 | 8 16.726048 2509 | 9 23.139831 2510 | 10 28.172396 2511 | 11 38.394777 2512 | 12 40.798287 2513 | 2514 | -0.3333333E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2515 | 1 3.362772 2516 | 2 3.762187 2517 | 3 4.196550 2518 | 4 5.007105 2519 | 5 5.578792 2520 | 6 9.437053 2521 | 7 13.202310 2522 | 8 16.473734 2523 | 9 23.919201 2524 | 10 28.615251 2525 | 11 37.955493 2526 | 12 41.617948 2527 | 2528 | -0.2857143E+00 0.2857143E+00 0.4761905E-01 0.5183026E-02 2529 | 1 3.286371 2530 | 2 3.888931 2531 | 3 4.034142 2532 | 4 4.994688 2533 | 5 5.440015 2534 | 6 8.354603 2535 | 7 14.910052 2536 | 8 16.574917 2537 | 9 24.810297 2538 | 10 29.196559 2539 | 11 37.348151 2540 | 12 40.888539 2541 | 2542 | -0.2380952E+00 0.2857143E+00 0.4761905E-01 0.2591513E-02 2543 | 1 3.016218 2544 | 2 3.926153 2545 | 3 4.041849 2546 | 4 4.976515 2547 | 5 5.309298 2548 | 6 7.472014 2549 | 7 16.057327 2550 | 8 16.750588 2551 | 9 25.808448 2552 | 10 29.909376 2553 | 11 36.483419 2554 | 12 39.831545 2555 | 2556 | 0.3809524E+00 0.3333333E+00 0.4761905E-01 0.2591513E-02 2557 | 1 2.277749 2558 | 2 3.411288 2559 | 3 4.985981 2560 | 4 5.463625 2561 | 5 5.576924 2562 | 6 5.912210 2563 | 7 20.517545 2564 | 8 20.724271 2565 | 9 21.082643 2566 | 10 26.457364 2567 | 11 28.838229 2568 | 12 44.714873 2569 | 2570 | 0.4285714E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2571 | 1 2.562466 2572 | 2 3.367961 2573 | 3 4.929317 2574 | 4 5.465575 2575 | 5 5.697987 2576 | 6 6.565836 2577 | 7 18.468872 2578 | 8 18.970628 2579 | 9 20.640437 2580 | 10 26.922588 2581 | 11 29.786332 2582 | 12 43.760163 2583 | 2584 | 0.4761905E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2585 | 1 2.798050 2586 | 2 3.364239 2587 | 3 4.832260 2588 | 4 5.408777 2589 | 5 5.792607 2590 | 6 7.485070 2591 | 7 16.834068 2592 | 8 17.802240 2593 | 9 20.717262 2594 | 10 27.112076 2595 | 11 31.429283 2596 | 12 42.847488 2597 | 2598 | -0.4761905E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2599 | 1 2.990417 2600 | 2 3.397871 2601 | 3 4.717500 2602 | 4 5.342786 2603 | 5 5.825856 2604 | 6 8.613197 2605 | 7 15.637478 2606 | 8 16.967459 2607 | 9 20.969928 2608 | 10 27.248306 2609 | 11 33.418507 2610 | 12 42.014095 2611 | 2612 | -0.4285714E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2613 | 1 3.155962 2614 | 2 3.456039 2615 | 3 4.588567 2616 | 4 5.274789 2617 | 5 5.801376 2618 | 6 9.907148 2619 | 7 13.890502 2620 | 8 16.566770 2621 | 9 21.393296 2622 | 10 27.433802 2623 | 11 35.612340 2624 | 12 41.254889 2625 | 2626 | -0.3809524E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2627 | 1 3.310226 2628 | 2 3.522054 2629 | 3 4.440724 2630 | 4 5.209389 2631 | 5 5.729568 2632 | 6 11.227269 2633 | 7 12.391294 2634 | 8 16.284990 2635 | 9 22.040596 2636 | 10 27.707128 2637 | 11 37.940403 2638 | 12 40.546225 2639 | 2640 | -0.3333333E+00 0.3333333E+00 0.4761905E-01 0.5183026E-02 2641 | 1 3.441401 2642 | 2 3.595915 2643 | 3 4.271531 2644 | 4 5.147691 2645 | 5 5.627590 2646 | 6 10.663858 2647 | 7 13.022511 2648 | 8 16.082116 2649 | 9 22.870900 2650 | 10 28.088051 2651 | 11 39.579859 2652 | 12 40.552758 2653 | 2654 | -0.2857143E+00 0.3333333E+00 0.4761905E-01 0.2591513E-02 2655 | 1 3.458867 2656 | 2 3.723647 2657 | 3 4.108883 2658 | 4 5.081571 2659 | 5 5.520763 2660 | 6 9.426664 2661 | 7 14.328782 2662 | 8 16.163969 2663 | 9 23.809323 2664 | 10 28.586185 2665 | 11 38.668140 2666 | 12 42.053148 2667 | 2668 | 0.4285714E+00 0.3809524E+00 0.4761905E-01 0.2591513E-02 2669 | 1 2.579596 2670 | 2 3.262065 2671 | 3 5.340527 2672 | 4 5.668502 2673 | 5 5.690942 2674 | 6 6.712987 2675 | 7 17.710623 2676 | 8 19.649727 2677 | 9 19.925590 2678 | 10 26.066854 2679 | 11 28.467688 2680 | 12 46.052184 2681 | 2682 | 0.4761905E+00 0.3809524E+00 0.4761905E-01 0.5183026E-02 2683 | 1 2.717300 2684 | 2 3.243601 2685 | 3 5.234579 2686 | 4 5.612353 2687 | 5 5.798359 2688 | 6 7.645584 2689 | 7 16.909657 2690 | 8 17.927687 2691 | 9 19.864094 2692 | 10 26.592159 2693 | 11 29.462984 2694 | 12 45.088251 2695 | 2696 | -0.4761905E+00 0.3809524E+00 0.4761905E-01 0.5183026E-02 2697 | 1 2.856621 2698 | 2 3.266243 2699 | 3 5.083341 2700 | 4 5.534621 2701 | 5 5.847292 2702 | 6 8.759621 2703 | 7 16.295868 2704 | 8 16.825240 2705 | 9 20.030563 2706 | 10 26.862608 2707 | 11 31.159561 2708 | 12 44.155491 2709 | 2710 | -0.4285714E+00 0.3809524E+00 0.4761905E-01 0.5183026E-02 2711 | 1 3.003229 2712 | 2 3.324920 2713 | 3 4.896876 2714 | 4 5.453993 2715 | 5 5.834937 2716 | 6 10.000721 2717 | 7 15.232010 2718 | 8 15.834181 2719 | 9 20.311487 2720 | 10 27.072502 2721 | 11 33.221006 2722 | 12 43.245649 2723 | 2724 | -0.3809524E+00 0.3809524E+00 0.4761905E-01 0.5183026E-02 2725 | 1 3.160655 2726 | 2 3.407701 2727 | 3 4.684093 2728 | 4 5.367608 2729 | 5 5.777778 2730 | 6 11.289538 2731 | 7 13.580851 2732 | 8 15.391298 2733 | 9 20.913126 2734 | 10 27.324853 2735 | 11 35.499977 2736 | 12 42.283565 2737 | 2738 | -0.3333333E+00 0.3809524E+00 0.4761905E-01 0.2591513E-02 2739 | 1 3.324866 2740 | 2 3.497133 2741 | 3 4.461213 2742 | 4 5.268030 2743 | 5 5.699890 2744 | 6 12.043545 2745 | 7 12.422666 2746 | 8 15.392604 2747 | 9 21.836445 2748 | 10 27.658010 2749 | 11 37.932392 2750 | 12 41.125973 2751 | 2752 | 0.4761905E+00 0.4285714E+00 0.4761905E-01 0.2591513E-02 2753 | 1 2.658850 2754 | 2 3.155635 2755 | 3 5.575255 2756 | 4 5.744563 2757 | 5 5.826546 2758 | 6 7.931821 2759 | 7 16.161069 2760 | 8 19.297414 2761 | 9 19.463400 2762 | 10 25.816668 2763 | 11 28.249242 2764 | 12 47.225237 2765 | 2766 | -0.4761905E+00 0.4285714E+00 0.4761905E-01 0.5183026E-02 2767 | 1 2.741123 2768 | 2 3.165004 2769 | 3 5.441990 2770 | 4 5.675536 2771 | 5 5.878953 2772 | 6 8.898746 2773 | 7 15.427982 2774 | 8 18.070318 2775 | 9 19.450066 2776 | 10 26.400348 2777 | 11 29.283134 2778 | 12 46.242634 2779 | 2780 | -0.4285714E+00 0.4285714E+00 0.4761905E-01 0.5183026E-02 2781 | 1 2.854649 2782 | 2 3.215842 2783 | 3 5.228940 2784 | 4 5.591746 2785 | 5 5.876773 2786 | 6 9.918545 2787 | 7 14.896202 2788 | 8 16.744154 2789 | 9 19.470965 2790 | 10 26.742010 2791 | 11 31.029757 2792 | 12 45.148579 2793 | 2794 | -0.3809524E+00 0.4285714E+00 0.4761905E-01 0.2591513E-02 2795 | 1 2.997280 2796 | 2 3.301227 2797 | 3 4.964412 2798 | 4 5.494257 2799 | 5 5.833865 2800 | 6 10.802682 2801 | 7 14.703239 2802 | 8 15.246322 2803 | 9 19.934252 2804 | 10 27.015519 2805 | 11 33.156217 2806 | 12 43.859553 2807 | 2808 | -0.4761905E+00 0.4761905E+00 0.4761905E-01 0.2591513E-02 2809 | 1 2.667032 2810 | 2 3.100278 2811 | 3 5.703932 2812 | 4 5.753452 2813 | 5 5.913279 2814 | 6 9.015640 2815 | 7 14.708027 2816 | 8 19.555443 2817 | 9 19.297711 2818 | 10 25.694636 2819 | 11 28.148661 2820 | 12 47.978589 2821 | 2822 | -0.4285714E+00 0.4761905E+00 0.4761905E-01 0.2591513E-02 2823 | 1 2.739531 2824 | 2 3.138119 2825 | 3 5.522073 2826 | 4 5.681092 2827 | 5 5.912472 2828 | 6 9.577998 2829 | 7 14.527850 2830 | 8 18.211770 2831 | 9 19.164616 2832 | 10 26.337431 2833 | 11 29.225666 2834 | 12 46.769302 2835 | 2836 | 0.2857143E+00 0.1904762E+00 0.9523810E-01 0.2591513E-02 2837 | 1 0.228172 2838 | 2 3.989763 2839 | 3 4.602337 2840 | 4 4.827359 2841 | 5 5.022019 2842 | 6 5.505787 2843 | 7 23.019519 2844 | 8 23.670442 2845 | 9 29.835249 2846 | 10 30.738805 2847 | 11 33.092326 2848 | 12 37.650761 2849 | 2850 | 0.3333333E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2851 | 1 0.847450 2852 | 2 3.921954 2853 | 3 4.516023 2854 | 4 4.880010 2855 | 5 5.096573 2856 | 6 5.681384 2857 | 7 20.344820 2858 | 8 22.968550 2859 | 9 28.646899 2860 | 10 31.047181 2861 | 11 33.681541 2862 | 12 36.909187 2863 | 2864 | 0.3809524E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2865 | 1 1.518867 2866 | 2 3.897009 2867 | 3 4.412681 2868 | 4 4.854115 2869 | 5 5.262442 2870 | 6 5.927308 2871 | 7 17.163082 2872 | 8 22.663302 2873 | 9 27.471291 2874 | 10 30.881040 2875 | 11 34.633804 2876 | 12 36.602238 2877 | 2878 | 0.4285714E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2879 | 1 2.142336 2880 | 2 3.952613 2881 | 3 4.315205 2882 | 4 4.806762 2883 | 5 5.444163 2884 | 6 6.283885 2885 | 7 16.212620 2886 | 8 22.539472 2887 | 9 26.386744 2888 | 10 30.601183 2889 | 11 34.785977 2890 | 12 37.470512 2891 | 2892 | 0.4761905E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2893 | 1 2.584027 2894 | 2 4.083002 2895 | 3 4.277467 2896 | 4 4.759323 2897 | 5 5.595566 2898 | 6 6.830213 2899 | 7 14.615292 2900 | 8 22.600328 2901 | 9 25.420589 2902 | 10 30.335504 2903 | 11 34.512136 2904 | 12 38.628069 2905 | 2906 | -0.4761905E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2907 | 1 2.803144 2908 | 2 4.095256 2909 | 3 4.434827 2910 | 4 4.720400 2911 | 5 5.695072 2912 | 6 7.607847 2913 | 7 12.865546 2914 | 8 22.845069 2915 | 9 24.589683 2916 | 10 30.134167 2917 | 11 34.276677 2918 | 12 39.013076 2919 | 2920 | -0.4285714E+00 0.1904762E+00 0.9523810E-01 0.5183026E-02 2921 | 1 2.881431 2922 | 2 4.069289 2923 | 3 4.583901 2924 | 4 4.691400 2925 | 5 5.742113 2926 | 6 8.414434 2927 | 7 11.578502 2928 | 8 23.213111 2929 | 9 23.933195 2930 | 10 30.026149 2931 | 11 34.153889 2932 | 12 39.075632 2933 | 2934 | 0.3333333E+00 0.2380952E+00 0.9523810E-01 0.2591513E-02 2935 | 1 1.041006 2936 | 2 3.816956 2937 | 3 4.526947 2938 | 4 5.020875 2939 | 5 5.132132 2940 | 6 5.673561 2941 | 7 21.583684 2942 | 8 22.231481 2943 | 9 26.959232 2944 | 10 29.816185 2945 | 11 32.003774 2946 | 12 38.921873 2947 | 2948 | 0.3809524E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 2949 | 1 1.656520 2950 | 2 3.763473 2951 | 3 4.413412 2952 | 4 5.044388 2953 | 5 5.275185 2954 | 6 5.937581 2955 | 7 18.591522 2956 | 8 21.560922 2957 | 9 25.730388 2958 | 10 30.150621 2959 | 11 32.928811 2960 | 12 38.124872 2961 | 2962 | 0.4285714E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 2963 | 1 2.251537 2964 | 2 3.759858 2965 | 3 4.303893 2966 | 4 5.002314 2967 | 5 5.464349 2968 | 6 6.339295 2969 | 7 16.468724 2970 | 8 21.356483 2971 | 9 24.583423 2972 | 10 30.113321 2973 | 11 34.342853 2974 | 12 37.569518 2975 | 2976 | 0.4761905E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 2977 | 1 2.712402 2978 | 2 3.821174 2979 | 3 4.230833 2980 | 4 4.949783 2981 | 5 5.615689 2982 | 6 6.952155 2983 | 7 15.589489 2984 | 8 21.382914 2985 | 9 23.538394 2986 | 10 30.036186 2987 | 11 35.335633 2988 | 12 37.841982 2989 | 2990 | -0.4761905E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 2991 | 1 2.971574 2992 | 2 3.885060 2993 | 3 4.258757 2994 | 4 4.897562 2995 | 5 5.708205 2996 | 6 7.812079 2997 | 7 13.766904 2998 | 8 21.589696 2999 | 9 22.601160 3000 | 10 30.019778 3001 | 11 35.333939 3002 | 12 39.416741 3003 | 3004 | -0.4285714E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3005 | 1 3.082918 3006 | 2 3.889274 3007 | 3 4.364225 3008 | 4 4.852926 3009 | 5 5.743279 3010 | 6 8.824117 3011 | 7 12.194118 3012 | 8 21.651960 3013 | 9 22.162202 3014 | 10 30.102304 3015 | 11 35.047661 3016 | 12 41.358651 3017 | 3018 | -0.3809524E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3019 | 1 3.114266 3020 | 2 3.894846 3021 | 3 4.403394 3022 | 4 4.830919 3023 | 5 5.727460 3024 | 6 9.313326 3025 | 7 11.579006 3026 | 8 20.924436 3027 | 9 22.716674 3028 | 10 30.304542 3029 | 11 34.772541 3030 | 12 41.169915 3031 | 3032 | -0.3333333E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3033 | 1 3.080863 3034 | 2 3.931834 3035 | 3 4.313852 3036 | 4 4.845804 3037 | 5 5.665375 3038 | 6 8.523243 3039 | 7 12.676451 3040 | 8 20.197576 3041 | 9 23.411535 3042 | 10 30.635695 3043 | 11 34.553055 3044 | 12 40.870509 3045 | 3046 | -0.2857143E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3047 | 1 2.944933 3048 | 2 3.970931 3049 | 3 4.164654 3050 | 4 4.884645 3051 | 5 5.557138 3052 | 6 7.574247 3053 | 7 14.376103 3054 | 8 19.769461 3055 | 9 24.228682 3056 | 10 31.090661 3057 | 11 34.372467 3058 | 12 39.738863 3059 | 3060 | -0.2380952E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3061 | 1 2.619440 3062 | 2 3.891364 3063 | 3 4.181007 3064 | 4 4.934460 3065 | 5 5.405056 3066 | 6 6.835192 3067 | 7 16.072989 3068 | 8 19.630065 3069 | 9 25.179245 3070 | 10 31.639299 3071 | 11 34.150201 3072 | 12 38.098082 3073 | 3074 | -0.1904762E+00 0.2380952E+00 0.9523810E-01 0.5183026E-02 3075 | 1 2.086875 3076 | 2 3.837057 3077 | 3 4.290809 3078 | 4 4.990331 3079 | 5 5.224449 3080 | 6 6.316928 3081 | 7 16.878945 3082 | 8 19.698093 3083 | 9 26.249487 3084 | 10 32.179309 3085 | 11 33.700819 3086 | 12 36.960230 3087 | 3088 | -0.1428571E+00 0.2380952E+00 0.9523810E-01 0.2591513E-02 3089 | 1 1.458974 3090 | 2 3.857967 3091 | 3 4.409732 3092 | 4 5.037757 3093 | 5 5.050828 3094 | 6 5.951238 3095 | 7 19.444723 3096 | 8 20.334123 3097 | 9 27.424859 3098 | 10 32.113273 3099 | 11 33.347866 3100 | 12 36.648911 3101 | 3102 | 0.3809524E+00 0.2857143E+00 0.9523810E-01 0.2591513E-02 3103 | 1 1.862016 3104 | 2 3.640146 3105 | 3 4.488092 3106 | 4 5.238801 3107 | 5 5.317635 3108 | 6 5.930432 3109 | 7 20.045470 3110 | 8 20.682112 3111 | 9 24.064539 3112 | 10 29.030342 3113 | 11 31.195767 3114 | 12 40.178356 3115 | 3116 | 0.4285714E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3117 | 1 2.374116 3118 | 2 3.605654 3119 | 3 4.389637 3120 | 4 5.231113 3121 | 5 5.477251 3122 | 6 6.380034 3123 | 7 17.405749 3124 | 8 20.145423 3125 | 9 22.810568 3126 | 10 29.440912 3127 | 11 32.309763 3128 | 12 39.321526 3129 | 3130 | 0.4761905E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3131 | 1 2.793647 3132 | 2 3.617973 3133 | 3 4.303873 3134 | 4 5.171603 3135 | 5 5.631808 3136 | 6 7.059061 3137 | 7 16.337793 3138 | 8 20.132604 3139 | 9 21.567761 3140 | 10 29.525617 3141 | 11 33.973571 3142 | 12 38.611515 3143 | 3144 | -0.4761905E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3145 | 1 3.069512 3146 | 2 3.664839 3147 | 3 4.260947 3148 | 4 5.105639 3149 | 5 5.725090 3150 | 6 7.988672 3151 | 7 14.841966 3152 | 8 20.052276 3153 | 9 20.613739 3154 | 10 29.597591 3155 | 11 35.643084 3156 | 12 38.272970 3157 | 3158 | -0.4285714E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3159 | 1 3.227044 3160 | 2 3.696941 3161 | 3 4.271743 3162 | 4 5.042199 3163 | 5 5.757267 3164 | 6 9.122282 3165 | 7 13.104127 3166 | 8 19.065422 3167 | 9 20.821294 3168 | 10 29.749506 3169 | 11 36.350038 3170 | 12 39.166561 3171 | 3172 | -0.3809524E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3173 | 1 3.315955 3174 | 2 3.710147 3175 | 3 4.280094 3176 | 4 4.995661 3177 | 5 5.734944 3178 | 6 10.141231 3179 | 7 11.864248 3180 | 8 18.074587 3181 | 9 21.431110 3182 | 10 30.016795 3183 | 11 36.138607 3184 | 12 41.155689 3185 | 3186 | -0.3333333E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3187 | 1 3.339729 3188 | 2 3.754001 3189 | 3 4.212198 3190 | 4 4.983509 3191 | 5 5.664900 3192 | 6 9.584259 3193 | 7 12.584221 3194 | 8 17.396596 3195 | 9 22.183290 3196 | 10 30.418724 3197 | 11 35.758078 3198 | 12 42.271488 3199 | 3200 | -0.2857143E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3201 | 1 3.272665 3202 | 2 3.843634 3203 | 3 4.066624 3204 | 4 5.005931 3205 | 5 5.552832 3206 | 6 8.482074 3207 | 7 14.234344 3208 | 8 17.241234 3209 | 9 23.091874 3210 | 10 30.964768 3211 | 11 35.308907 3212 | 12 40.295225 3213 | 3214 | -0.2380952E+00 0.2857143E+00 0.9523810E-01 0.5183026E-02 3215 | 1 3.053718 3216 | 2 3.863547 3217 | 3 4.024576 3218 | 4 5.044780 3219 | 5 5.407988 3220 | 6 7.548378 3221 | 7 15.936869 3222 | 8 17.165689 3223 | 9 24.048745 3224 | 10 31.650111 3225 | 11 34.694906 3226 | 12 38.658233 3227 | 3228 | -0.1904762E+00 0.2857143E+00 0.9523810E-01 0.2591513E-02 3229 | 1 2.633442 3230 | 2 3.817215 3231 | 3 4.153218 3232 | 4 5.068935 3233 | 5 5.263788 3234 | 6 6.838667 3235 | 7 16.490627 3236 | 8 17.477213 3237 | 9 25.099989 3238 | 10 32.403982 3239 | 11 33.809755 3240 | 12 37.557232 3241 | 3242 | 0.4285714E+00 0.3333333E+00 0.9523810E-01 0.2591513E-02 3243 | 1 2.479560 3244 | 2 3.476243 3245 | 3 4.610849 3246 | 4 5.462196 3247 | 5 5.492209 3248 | 6 6.428911 3249 | 7 18.721445 3250 | 8 19.400475 3251 | 9 20.936923 3252 | 10 28.407515 3253 | 11 30.614729 3254 | 12 41.467897 3255 | 3256 | 0.4761905E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3257 | 1 2.794274 3258 | 2 3.464178 3259 | 3 4.537522 3260 | 4 5.405425 3261 | 5 5.641278 3262 | 6 7.170409 3263 | 7 16.977259 3264 | 8 18.939882 3265 | 9 19.704901 3266 | 10 28.904214 3267 | 11 31.849487 3268 | 12 40.578277 3269 | 3270 | -0.4761905E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3271 | 1 3.030466 3272 | 2 3.491512 3273 | 3 4.458876 3274 | 4 5.322526 3275 | 5 5.744841 3276 | 6 8.154049 3277 | 7 16.003853 3278 | 8 18.181075 3279 | 9 19.221135 3280 | 10 29.100437 3281 | 11 33.662259 3282 | 12 39.791584 3283 | 3284 | -0.4285714E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3285 | 1 3.204835 3286 | 2 3.537912 3287 | 3 4.385532 3288 | 4 5.238190 3289 | 5 5.783627 3290 | 6 9.344848 3291 | 7 14.252851 3292 | 8 17.223318 3293 | 9 19.580342 3294 | 10 29.289971 3295 | 11 35.678242 3296 | 12 39.158216 3297 | 3298 | -0.3809524E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3299 | 1 3.352092 3300 | 2 3.573018 3301 | 3 4.308305 3302 | 4 5.164146 3303 | 5 5.764704 3304 | 6 10.650796 3305 | 7 12.650025 3306 | 8 16.683582 3307 | 9 20.142670 3308 | 10 29.561420 3309 | 11 37.358467 3310 | 12 39.081357 3311 | 3312 | -0.3333333E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3313 | 1 3.485789 3314 | 2 3.593100 3315 | 3 4.199722 3316 | 4 5.116234 3317 | 5 5.697976 3318 | 6 10.784106 3319 | 7 12.544561 3320 | 8 16.450306 3321 | 9 20.940067 3322 | 10 29.945352 3323 | 11 37.370994 3324 | 12 40.661644 3325 | 3326 | -0.2857143E+00 0.3333333E+00 0.9523810E-01 0.5183026E-02 3327 | 1 3.505570 3328 | 2 3.693437 3329 | 3 4.047292 3330 | 4 5.099198 3331 | 5 5.596282 3332 | 6 9.561676 3333 | 7 14.081794 3334 | 8 16.378224 3335 | 9 21.962673 3336 | 10 30.450654 3337 | 11 36.745830 3338 | 12 40.732450 3339 | 3340 | -0.2380952E+00 0.3333333E+00 0.9523810E-01 0.2591513E-02 3341 | 1 3.367294 3342 | 2 3.854137 3343 | 3 3.910625 3344 | 4 5.091552 3345 | 5 5.483175 3346 | 6 8.455181 3347 | 7 15.484384 3348 | 8 16.542028 3349 | 9 22.988311 3350 | 10 31.058614 3351 | 11 35.815398 3352 | 12 39.489303 3353 | 3354 | 0.4761905E+00 0.3809524E+00 0.9523810E-01 0.2591513E-02 3355 | 1 2.752977 3356 | 2 3.339766 3357 | 3 4.883542 3358 | 4 5.597393 3359 | 5 5.666975 3360 | 6 7.336604 3361 | 7 17.976402 3362 | 8 18.007929 3363 | 9 18.485781 3364 | 10 27.948937 3365 | 11 30.221243 3366 | 12 42.799868 3367 | 3368 | -0.4761905E+00 0.3809524E+00 0.9523810E-01 0.5183026E-02 3369 | 1 2.918632 3370 | 2 3.353278 3371 | 3 4.786823 3372 | 4 5.515491 3373 | 5 5.773283 3374 | 6 8.338128 3375 | 7 16.853831 3376 | 8 16.937763 3377 | 9 18.408598 3378 | 10 28.530262 3379 | 11 31.539914 3380 | 12 41.887360 3381 | 3382 | -0.4285714E+00 0.3809524E+00 0.9523810E-01 0.5183026E-02 3383 | 1 3.068276 3384 | 2 3.403149 3385 | 3 4.662331 3386 | 4 5.414688 3387 | 5 5.823419 3388 | 6 9.526102 3389 | 7 15.562244 3390 | 8 16.525347 3391 | 9 18.585696 3392 | 10 28.825066 3393 | 11 33.450493 3394 | 12 41.034293 3395 | 3396 | -0.3809524E+00 0.3809524E+00 0.9523810E-01 0.5183026E-02 3397 | 1 3.214368 3398 | 2 3.473272 3399 | 3 4.513893 3400 | 4 5.317582 3401 | 5 5.813856 3402 | 6 10.851468 3403 | 7 13.861354 3404 | 8 16.125031 3405 | 9 18.962952 3406 | 10 29.108822 3407 | 11 35.629797 3408 | 12 40.207004 3409 | 3410 | -0.3333333E+00 0.3809524E+00 0.9523810E-01 0.5183026E-02 3411 | 1 3.367487 3412 | 2 3.541732 3413 | 3 4.341730 3414 | 4 5.234699 3415 | 5 5.758144 3416 | 6 12.027208 3417 | 7 12.504056 3418 | 8 15.794286 3419 | 9 19.696680 3420 | 10 29.466296 3421 | 11 37.925358 3422 | 12 39.331957 3423 | 3424 | -0.2857143E+00 0.3809524E+00 0.9523810E-01 0.2591513E-02 3425 | 1 3.529644 3426 | 2 3.576935 3427 | 3 4.161884 3428 | 4 5.167134 3429 | 5 5.675863 3430 | 6 10.809951 3431 | 7 13.520910 3432 | 8 15.838614 3433 | 9 20.747068 3434 | 10 29.918819 3435 | 11 38.134543 3436 | 12 40.197178 3437 | 3438 | -0.4761905E+00 0.4285714E+00 0.9523810E-01 0.2591513E-02 3439 | 1 2.811862 3440 | 2 3.242157 3441 | 3 5.143024 3442 | 4 5.635804 3443 | 5 5.825894 3444 | 6 8.578279 3445 | 7 16.290502 3446 | 8 17.974049 3447 | 9 18.101594 3448 | 10 27.648796 3449 | 11 29.981351 3450 | 12 44.109186 3451 | 3452 | -0.4285714E+00 0.4285714E+00 0.9523810E-01 0.5183026E-02 3453 | 1 2.923783 3454 | 2 3.283266 3455 | 3 4.992610 3456 | 4 5.547316 3457 | 5 5.873759 3458 | 6 9.665156 3459 | 7 15.653118 3460 | 8 16.830716 3461 | 9 17.952547 3462 | 10 28.310008 3463 | 11 31.364021 3464 | 12 43.111682 3465 | 3466 | -0.3809524E+00 0.4285714E+00 0.9523810E-01 0.5183026E-02 3467 | 1 3.054473 3468 | 2 3.359905 3469 | 3 4.795318 3470 | 4 5.444118 3471 | 5 5.870826 3472 | 6 10.806483 3473 | 7 15.059471 3474 | 8 15.436690 3475 | 9 18.082875 3476 | 10 28.689926 3477 | 11 33.345265 3478 | 12 42.043197 3479 | 3480 | -0.3333333E+00 0.4285714E+00 0.9523810E-01 0.2591513E-02 3481 | 1 3.203913 3482 | 2 3.456321 3483 | 3 4.566716 3484 | 4 5.341047 3485 | 5 5.824553 3486 | 6 11.792794 3487 | 7 13.672003 3488 | 8 15.162877 3489 | 9 18.642987 3490 | 10 29.049104 3491 | 11 35.605131 3492 | 12 40.775965 3493 | 3494 | -0.4285714E+00 0.4761905E+00 0.9523810E-01 0.2591513E-02 3495 | 1 2.810186 3496 | 2 3.191341 3497 | 3 5.299934 3498 | 4 5.630823 3499 | 5 5.913330 3500 | 6 9.709582 3501 | 7 14.851209 3502 | 8 18.196480 3503 | 9 17.644025 3504 | 10 27.500627 3505 | 11 29.868662 3506 | 12 45.110945 3507 | 3508 | -0.3809524E+00 0.4761905E+00 0.9523810E-01 0.2591513E-02 3509 | 1 2.916092 3510 | 2 3.259390 3511 | 3 5.073427 3512 | 4 5.545804 3513 | 5 5.912042 3514 | 6 10.417021 3515 | 7 14.699861 3516 | 8 16.814190 3517 | 9 17.646550 3518 | 10 28.237300 3519 | 11 31.307239 3520 | 12 43.723483 3521 | 3522 | 0.4285714E+00 0.2857143E+00 0.1428571E+00 0.2591513E-02 3523 | 1 2.217889 3524 | 2 3.717821 3525 | 3 4.294896 3526 | 4 5.204598 3527 | 5 5.244364 3528 | 6 6.340876 3529 | 7 17.903509 3530 | 8 18.820528 3531 | 9 24.501874 3532 | 10 31.321423 3533 | 11 32.971404 3534 | 12 37.515302 3535 | 3536 | 0.4761905E+00 0.2857143E+00 0.1428571E+00 0.5183026E-02 3537 | 1 2.719563 3538 | 2 3.720653 3539 | 3 4.176253 3540 | 4 5.150423 3541 | 5 5.419842 3542 | 6 6.910232 3543 | 7 16.438479 3544 | 8 18.407249 3545 | 9 23.350963 3546 | 10 31.671231 3547 | 11 33.775116 3548 | 12 37.310077 3549 | 3550 | -0.4761905E+00 0.2857143E+00 0.1428571E+00 0.5183026E-02 3551 | 1 3.068809 3552 | 2 3.771387 3553 | 3 4.100962 3554 | 4 5.083890 3555 | 5 5.572762 3556 | 6 7.703986 3557 | 7 15.403506 3558 | 8 18.492801 3559 | 9 22.254239 3560 | 10 31.580293 3561 | 11 34.498214 3562 | 12 37.986711 3563 | 3564 | -0.4285714E+00 0.2857143E+00 0.1428571E+00 0.5183026E-02 3565 | 1 3.248562 3566 | 2 3.775852 3567 | 3 4.148590 3568 | 4 5.022043 3569 | 5 5.675557 3570 | 6 8.727241 3571 | 7 13.612896 3572 | 8 18.810931 3573 | 9 21.184763 3574 | 10 31.453073 3575 | 11 34.673954 3576 | 12 39.602014 3577 | 3578 | -0.3809524E+00 0.2857143E+00 0.1428571E+00 0.5183026E-02 3579 | 1 3.324116 3580 | 2 3.735774 3581 | 3 4.236780 3582 | 4 4.980197 3583 | 5 5.725652 3584 | 6 9.829601 3585 | 7 12.117387 3586 | 8 19.342699 3587 | 9 20.143015 3588 | 10 31.377026 3589 | 11 34.670446 3590 | 12 41.622124 3591 | 3592 | 0.4761905E+00 0.3333333E+00 0.1428571E+00 0.2591513E-02 3593 | 1 2.778622 3594 | 2 3.583790 3595 | 3 4.292769 3596 | 4 5.358277 3597 | 5 5.461905 3598 | 6 7.005840 3599 | 7 16.993951 3600 | 8 17.909998 3601 | 9 21.435897 3602 | 10 30.670625 3603 | 11 32.623125 3604 | 12 38.545027 3605 | 3606 | -0.4761905E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3607 | 1 3.084723 3608 | 2 3.607415 3609 | 3 4.202097 3610 | 4 5.292312 3611 | 5 5.601668 3612 | 6 7.852481 3613 | 7 16.268724 3614 | 8 17.688065 3615 | 9 20.031676 3616 | 10 31.185594 3617 | 11 33.800389 3618 | 12 37.991794 3619 | 3620 | -0.4285714E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3621 | 1 3.281100 3622 | 2 3.645971 3623 | 3 4.153090 3624 | 4 5.206784 3625 | 5 5.706088 3626 | 6 8.923857 3627 | 7 14.732229 3628 | 8 17.846805 3629 | 9 18.712620 3630 | 10 31.359098 3631 | 11 35.111267 3632 | 12 37.999407 3633 | 3634 | -0.3809524E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3635 | 1 3.417423 3636 | 2 3.634720 3637 | 3 4.154357 3638 | 4 5.133434 3639 | 5 5.754034 3640 | 6 10.166349 3641 | 7 13.044159 3642 | 8 17.757939 3643 | 9 18.047383 3644 | 10 31.532792 3645 | 11 35.542212 3646 | 12 39.189299 3647 | 3648 | -0.3333333E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3649 | 1 3.533554 3650 | 2 3.584416 3651 | 3 4.137440 3652 | 4 5.097525 3653 | 5 5.747474 3654 | 6 10.917237 3655 | 7 12.191715 3656 | 8 16.982796 3657 | 9 18.622388 3658 | 10 31.790306 3659 | 11 35.356770 3660 | 12 40.786453 3661 | 3662 | -0.2857143E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3663 | 1 3.484005 3664 | 2 3.679438 3665 | 3 4.048006 3666 | 4 5.113147 3667 | 5 5.690262 3668 | 6 9.803782 3669 | 7 13.547029 3670 | 8 16.773768 3671 | 9 19.744536 3672 | 10 32.154602 3673 | 11 34.959612 3674 | 12 39.797462 3675 | 3676 | -0.2380952E+00 0.3333333E+00 0.1428571E+00 0.5183026E-02 3677 | 1 3.357147 3678 | 2 3.768830 3679 | 3 3.955608 3680 | 4 5.163032 3681 | 5 5.588439 3682 | 6 8.636115 3683 | 7 15.257798 3684 | 8 16.721648 3685 | 9 20.939478 3686 | 10 32.625143 3687 | 11 34.287256 3688 | 12 38.314788 3689 | 3690 | -0.1904762E+00 0.3333333E+00 0.1428571E+00 0.2591513E-02 3691 | 1 3.109271 3692 | 2 3.730903 3693 | 3 4.039070 3694 | 4 5.205924 3695 | 5 5.466840 3696 | 6 7.666155 3697 | 7 16.260337 3698 | 8 16.962578 3699 | 9 22.150075 3700 | 10 33.147006 3701 | 11 33.188442 3702 | 12 37.486222 3703 | 3704 | -0.4761905E+00 0.3809524E+00 0.1428571E+00 0.2591513E-02 3705 | 1 3.008005 3706 | 2 3.470223 3707 | 3 4.462205 3708 | 4 5.452644 3709 | 5 5.665148 3710 | 6 8.048923 3711 | 7 16.724208 3712 | 8 17.636827 3713 | 9 17.914786 3714 | 10 30.181274 3715 | 11 32.357572 3716 | 12 39.718754 3717 | 3718 | -0.4285714E+00 0.3809524E+00 0.1428571E+00 0.5183026E-02 3719 | 1 3.169097 3720 | 2 3.516988 3721 | 3 4.373366 3722 | 4 5.368130 3723 | 5 5.758204 3724 | 6 9.136956 3725 | 7 15.930067 3726 | 8 17.007315 3727 | 9 17.335033 3728 | 10 30.819772 3729 | 11 33.763299 3730 | 12 38.972397 3731 | 3732 | -0.3809524E+00 0.3809524E+00 0.1428571E+00 0.5183026E-02 3733 | 1 3.297249 3734 | 2 3.578839 3735 | 3 4.281442 3736 | 4 5.270431 3737 | 5 5.808269 3738 | 6 10.404248 3739 | 7 14.243954 3740 | 8 16.667837 3741 | 9 17.266536 3742 | 10 31.158116 3743 | 11 35.588635 3744 | 12 38.414524 3745 | 3746 | -0.3333333E+00 0.3809524E+00 0.1428571E+00 0.5183026E-02 3747 | 1 3.412870 3748 | 2 3.630677 3749 | 3 4.177227 3750 | 4 5.196839 3751 | 5 5.803687 3752 | 6 11.730589 3753 | 7 12.714662 3754 | 8 16.404100 3755 | 9 17.577466 3756 | 10 31.501343 3757 | 11 36.719155 3758 | 12 38.694377 3759 | 3760 | -0.2857143E+00 0.3809524E+00 0.1428571E+00 0.5183026E-02 3761 | 1 3.500322 3762 | 2 3.689563 3763 | 3 4.040185 3764 | 4 5.167010 3765 | 5 5.751037 3766 | 6 11.039157 3767 | 7 13.446364 3768 | 8 16.200022 3769 | 9 18.328080 3770 | 10 31.923821 3771 | 11 36.293358 3772 | 12 39.632030 3773 | 3774 | -0.2380952E+00 0.3809524E+00 0.1428571E+00 0.2591513E-02 3775 | 1 3.501472 3776 | 2 3.792877 3777 | 3 3.886366 3778 | 4 5.171338 3779 | 5 5.666470 3780 | 6 9.754484 3781 | 7 14.756195 3782 | 8 16.293590 3783 | 9 19.614827 3784 | 10 32.407598 3785 | 11 35.295806 3786 | 12 38.952099 3787 | 3788 | -0.4285714E+00 0.4285714E+00 0.1428571E+00 0.2591513E-02 3789 | 1 3.033464 3790 | 2 3.388118 3791 | 3 4.684476 3792 | 4 5.480143 3793 | 5 5.824804 3794 | 6 9.385570 3795 | 7 16.481436 3796 | 8 16.755672 3797 | 9 16.910284 3798 | 10 29.856380 3799 | 11 32.185533 3800 | 12 40.955400 3801 | 3802 | -0.3809524E+00 0.4285714E+00 0.1428571E+00 0.5183026E-02 3803 | 1 3.146409 3804 | 2 3.461632 3805 | 3 4.546543 3806 | 4 5.381651 3807 | 5 5.870882 3808 | 6 10.581352 3809 | 7 15.573584 3810 | 8 16.029292 3811 | 9 16.771008 3812 | 10 30.596671 3813 | 11 33.719944 3814 | 12 40.034259 3815 | 3816 | -0.3333333E+00 0.4285714E+00 0.1428571E+00 0.5183026E-02 3817 | 1 3.263512 3818 | 2 3.556226 3819 | 3 4.378543 3820 | 4 5.285936 3821 | 5 5.868750 3822 | 6 11.842636 3823 | 7 13.958060 3824 | 8 15.596286 3825 | 9 16.878722 3826 | 10 31.047310 3827 | 11 35.782767 3828 | 12 39.027005 3829 | 3830 | -0.2857143E+00 0.4285714E+00 0.1428571E+00 0.2591513E-02 3831 | 1 3.385763 3832 | 2 3.649091 3833 | 3 4.190898 3834 | 4 5.215329 3835 | 5 5.822166 3836 | 6 12.424608 3837 | 7 12.940454 3838 | 8 15.579125 3839 | 9 17.340029 3840 | 10 31.479961 3841 | 11 37.778062 3842 | 12 37.985933 3843 | 3844 | -0.3809524E+00 0.4761905E+00 0.1428571E+00 0.2591513E-02 3845 | 1 3.015343 3846 | 2 3.345175 3847 | 3 4.830726 3848 | 4 5.471059 3849 | 5 5.913270 3850 | 6 10.631904 3851 | 7 15.106746 3852 | 8 16.756196 3853 | 9 16.469939 3854 | 10 29.694673 3855 | 11 32.102307 3856 | 12 41.959510 3857 | 3858 | -0.3333333E+00 0.4761905E+00 0.1428571E+00 0.2591513E-02 3859 | 1 3.129411 3860 | 2 3.443045 3861 | 3 4.616936 3862 | 4 5.376464 3863 | 5 5.912323 3864 | 6 11.454194 3865 | 7 14.955418 3866 | 8 15.508663 3867 | 9 16.523181 3868 | 10 30.521946 3869 | 11 33.702557 3870 | 12 40.599047 3871 | 3872 | -0.4285714E+00 0.3809524E+00 0.1904762E+00 0.2591513E-02 3873 | 1 3.287111 3874 | 2 3.653179 3875 | 3 4.097669 3876 | 4 5.289904 3877 | 5 5.664155 3878 | 6 8.846189 3879 | 7 16.000652 3880 | 8 16.683604 3881 | 9 18.657121 3882 | 10 32.640255 3883 | 11 33.874467 3884 | 12 37.687579 3885 | 3886 | -0.3809524E+00 0.3809524E+00 0.1904762E+00 0.5183026E-02 3887 | 1 3.399162 3888 | 2 3.716903 3889 | 3 4.031766 3890 | 4 5.224564 3891 | 5 5.751568 3892 | 6 10.032490 3893 | 7 14.659383 3894 | 8 16.529689 3895 | 9 17.587029 3896 | 10 33.124865 3897 | 11 34.604568 3898 | 12 37.999573 3899 | 3900 | -0.3333333E+00 0.3809524E+00 0.1904762E+00 0.5183026E-02 3901 | 1 3.456916 3902 | 2 3.724504 3903 | 3 4.030234 3904 | 4 5.173605 3905 | 5 5.802151 3906 | 6 11.357042 3907 | 7 13.071269 3908 | 8 16.575875 3909 | 9 16.741846 3910 | 10 33.157367 3911 | 11 35.088242 3912 | 12 39.059375 3913 | 3914 | -0.3809524E+00 0.4285714E+00 0.1904762E+00 0.2591513E-02 3915 | 1 3.265724 3916 | 2 3.595146 3917 | 3 4.267865 3918 | 4 5.311357 3919 | 5 5.823520 3920 | 6 10.295949 3921 | 7 15.581466 3922 | 8 16.433189 3923 | 9 16.624134 3924 | 10 32.301019 3925 | 11 34.318465 3926 | 12 38.262854 3927 | 3928 | -0.3333333E+00 0.4285714E+00 0.1904762E+00 0.5183026E-02 3929 | 1 3.334978 3930 | 2 3.697830 3931 | 3 4.151571 3932 | 4 5.236089 3933 | 5 5.869272 3934 | 6 11.599416 3935 | 7 14.268251 3936 | 8 16.317940 3937 | 9 16.331009 3938 | 10 33.081834 3939 | 11 35.547462 3940 | 12 37.744675 3941 | 3942 | -0.2857143E+00 0.4285714E+00 0.1904762E+00 0.5183026E-02 3943 | 1 3.386914 3944 | 2 3.802484 3945 | 3 4.020500 3946 | 4 5.194216 3947 | 5 5.868454 3948 | 6 12.652463 3949 | 7 13.085393 3950 | 8 16.054927 3951 | 9 16.337646 3952 | 10 33.529768 3953 | 11 35.947635 3954 | 12 38.022765 3955 | 3956 | -0.2380952E+00 0.4285714E+00 0.1904762E+00 0.2591513E-02 3957 | 1 3.415420 3958 | 2 3.863656 3959 | 3 3.916244 3960 | 4 5.199391 3961 | 5 5.822352 3962 | 6 11.303828 3963 | 7 14.214951 3964 | 8 16.072078 3965 | 9 16.760043 3966 | 10 33.914345 3967 | 11 34.946517 3968 | 12 38.207734 3969 | 3970 | -0.3333333E+00 0.4761905E+00 0.1904762E+00 0.2591513E-02 3971 | 1 3.224138 3972 | 2 3.564203 3973 | 3 4.389599 3974 | 4 5.302052 3975 | 5 5.913092 3976 | 6 11.695916 3977 | 7 15.325065 3978 | 8 15.485069 3979 | 9 16.286145 3980 | 10 32.131367 3981 | 11 34.506939 3982 | 12 38.932649 3983 | 3984 | -0.2857143E+00 0.4761905E+00 0.1904762E+00 0.2591513E-02 3985 | 1 3.304929 3986 | 2 3.690464 3987 | 3 4.203138 3988 | 4 5.233145 3989 | 5 5.912636 3990 | 6 12.633979 3991 | 7 14.113548 3992 | 8 15.420527 3993 | 9 16.278651 3994 | 10 33.034314 3995 | 11 36.080010 3996 | 12 37.599707 3997 | 3998 | -0.2857143E+00 0.4761905E+00 0.2380952E+00 0.2591513E-02 3999 | 1 3.346197 4000 | 2 3.850309 4001 | 3 4.006059 4002 | 4 5.199797 4003 | 5 5.912939 4004 | 6 12.863850 4005 | 7 13.930542 4006 | 8 15.786482 4007 | 9 16.135597 4008 | 10 34.771540 4009 | 11 35.500311 4010 | 12 37.054203 4011 | -------------------------------------------------------------------------------- /examples/copper/KPOINTS: -------------------------------------------------------------------------------- 1 | Monkhorst-Pack 2 | 0 3 | Gamma 4 | 21 21 21 5 | 0 0 0 6 | -------------------------------------------------------------------------------- /examples/copper/POSCAR: -------------------------------------------------------------------------------- 1 | Cu 2 | 1.00000000000000 3 | 0.0000000000000000 1.8159749723196192 1.8159749723196192 4 | 1.8159749723196192 0.0000000000000000 1.8159749723196192 5 | 1.8159749723196192 1.8159749723196192 0.0000000000000000 6 | Cu 7 | 1 8 | Direct 9 | 0.0000000000000000 0.0000000000000000 0.0000000000000000 10 | 11 | 0.00000000E+00 0.00000000E+00 0.00000000E+00 12 | -------------------------------------------------------------------------------- /examples/copper/fs.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs.png -------------------------------------------------------------------------------- /examples/copper/fs_matplotlib_screenshot.jpg: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_matplotlib_screenshot.jpg -------------------------------------------------------------------------------- /examples/copper/fs_matplotlib_screenshot.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_matplotlib_screenshot.png -------------------------------------------------------------------------------- /examples/copper/fs_mayavi_screenshot.jpg: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_mayavi_screenshot.jpg -------------------------------------------------------------------------------- /examples/copper/fs_mayavi_screenshot.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_mayavi_screenshot.png -------------------------------------------------------------------------------- /examples/copper/fs_uc_mayavi_screenshot.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_uc_mayavi_screenshot.png -------------------------------------------------------------------------------- /examples/copper/fs_xcrysden_screenshot.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/QijingZheng/VASP_FermiSurface/3fcc498f5a1179063159d7919bc9c9c4502690ba/examples/copper/fs_xcrysden_screenshot.png -------------------------------------------------------------------------------- /fs.py: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env python 2 | 3 | import numpy as np 4 | from ase.io import read 5 | import os 6 | import sys 7 | import argparse 8 | 9 | 10 | def find_fermi_level(band_energies, kpt_weight, 11 | nelect, occ=None, sigma=0.01, nedos=100, 12 | soc_band=False, 13 | nmax=1000): 14 | ''' 15 | Locate ther Fermi level from the band energies, k-points weights and number 16 | of electrons. 17 | 1.0 18 | Ne = \sum_{n,k} --------------------------- * w_k 19 | ((E_{nk}-E_f)/sigma) 20 | e + 1 21 | 22 | 23 | Inputs: 24 | band_energies: The band energies of shape (nspin, nkpts, nbnds) 25 | kpt_weight: The weight of each k-points. 26 | nelect: Number of electrons. 27 | occ: 1.0 for spin-polarized/SOC band energies, else 2.0. 28 | sigma: Broadening parameter for the Fermi-Dirac distribution. 29 | nedos: number of discrete points in approximately locating Fermi level. 30 | soc_band: band energies from SOC calculations? 31 | nmax: maximum iteration in finding the exact Fermi level. 32 | ''' 33 | 34 | if band_energies.ndim == 2: 35 | band_energies = band_energies[None, :] 36 | 37 | nspin, nkpts, nbnds = band_energies.shape 38 | 39 | if occ is None: 40 | if nspin == 1 and (not soc_band): 41 | occ = 2.0 42 | else: 43 | occ = 1.0 44 | 45 | if nbnds > nedos: 46 | nedos = nbnds * 5 47 | 48 | kpt_weight = np.asarray(kpt_weight, dtype=float) 49 | assert kpt_weight.shape == (nkpts,) 50 | kpt_weight /= np.sum(kpt_weight) 51 | 52 | emin = band_energies.min() 53 | emax = band_energies.max() 54 | e0 = np.linspace(emin, emax, nedos) 55 | de = e0[1] - e0[0] 56 | 57 | # find the approximated Fermi level 58 | nelect_lt_en = np.array([ 59 | np.sum(occ * (band_energies <= en) * kpt_weight[None, :, None]) 60 | for en in e0 61 | ]) 62 | ne_tmp = nelect_lt_en[nedos//2] 63 | if (np.abs(ne_tmp - nelect) < 0.05): 64 | i_fermi = nedos // 2 65 | i_lower = i_fermi - 1 66 | i_upper = i_fermi + 1 67 | elif (ne_tmp > nelect): 68 | for ii in range(nedos//2-1, -1, -1): 69 | ne_tmp = nelect_lt_en[ii] 70 | if ne_tmp < nelect: 71 | i_fermi = ii 72 | i_lower = i_fermi 73 | i_upper = i_fermi + 1 74 | break 75 | else: 76 | for ii in range(nedos//2+1, nedos): 77 | ne_tmp = nelect_lt_en[ii] 78 | if ne_tmp > nelect: 79 | i_fermi = ii 80 | i_lower = i_fermi - 1 81 | i_upper = i_fermi 82 | break 83 | 84 | ############################################################ 85 | # Below is the algorithm used by VASP, much slower 86 | ############################################################ 87 | # find the approximated Fermi level 88 | # x = (e0[None, None, None, :] - band_energies[:, :, :, None]) / sigma 89 | # x = x.clip(-100, 100) 90 | # dos = 1./sigma * np.exp(x) / (np.exp(x) + 1)**2 * \ 91 | # kpt_weight[None, :, None, None] * de 92 | # ddos = np.sum(dos, axis=(0,1,2)) 93 | # 94 | # nelect_from_dos_int = np.sum(ddos[:nedos/2]) 95 | # if (np.abs(nelect_from_dos_int - nelect) < 0.05): 96 | # i_fermi = nedos / 2 - 1 97 | # i_lower = i_fermi - 1 98 | # i_upper = i_fermi + 1 99 | # elif (nelect_from_dos_int > nelect): 100 | # for ii in range(nedos/2, -1, -1): 101 | # nelect_from_dos_int = np.sum(ddos[:ii]) 102 | # if nelect_from_dos_int < nelect: 103 | # i_fermi = ii 104 | # i_lower = i_fermi 105 | # i_upper = i_fermi + 1 106 | # break 107 | # else: 108 | # for ii in range(nedos/2, nedos): 109 | # nelect_from_dos_int = np.sum(ddos[:ii]) 110 | # if nelect_from_dos_int > nelect: 111 | # i_fermi = ii 112 | # i_lower = i_fermi - 1 113 | # i_upper = i_fermi 114 | # break 115 | 116 | # Locate the exact Fermi level using bisectioning 117 | e_lower = e0[i_lower] 118 | e_upper = e0[i_upper] 119 | lower_B = False 120 | upper_B = False 121 | for ii in range(nmax): 122 | e_fermi = (e_lower + e_upper) / 2. 123 | 124 | z = (band_energies - e_fermi) / sigma 125 | z = z.clip(-100, 100) 126 | F_nk = occ / (np.exp(z) + 1) 127 | N = np.sum(F_nk * kpt_weight[None, :, None]) 128 | # print ii, e_lower, e_upper, N 129 | 130 | if (np.abs(N - nelect) < 1E-10): 131 | break 132 | if (np.abs(e_upper - e_lower / (np.abs(e_fermi) + 1E-10)) < 1E-14): 133 | raise ValueError("Cannot reach the specified precision!") 134 | 135 | if (N > nelect): 136 | if not lower_B: 137 | e_lower -= de 138 | upper_B = True 139 | e_upper = e_fermi 140 | else: 141 | if not upper_B: 142 | e_upper += de 143 | lower_B = True 144 | e_lower = e_fermi 145 | 146 | if (ii == nmax - 1): 147 | raise ValueError("Cannot reach the specified precision!") 148 | 149 | return e_fermi, F_nk 150 | 151 | 152 | def get_brillouin_zone_3d(cell): 153 | """ 154 | Generate the Brillouin Zone of a given cell. The BZ is the Wigner-Seitz cell 155 | of the reciprocal lattice, which can be constructed by Voronoi decomposition 156 | to the reciprocal lattice. A Voronoi diagram is a subdivision of the space 157 | into the nearest neighborhoods of a given set of points. 158 | 159 | https://en.wikipedia.org/wiki/Wigner%E2%80%93Seitz_cell 160 | https://docs.scipy.org/doc/scipy/reference/tutorial/spatial.html#voronoi-diagrams 161 | """ 162 | 163 | cell = np.asarray(cell, dtype=float) 164 | assert cell.shape == (3, 3) 165 | 166 | px, py, pz = np.tensordot(cell, np.mgrid[-1:2, -1:2, -1:2], axes=[0, 0]) 167 | points = np.c_[px.ravel(), py.ravel(), pz.ravel()] 168 | 169 | from scipy.spatial import Voronoi 170 | vor = Voronoi(points) 171 | 172 | bz_facets = [] 173 | bz_ridges = [] 174 | bz_vertices = [] 175 | 176 | # for rid in vor.ridge_vertices: 177 | # if( np.all(np.array(rid) >= 0) ): 178 | # bz_ridges.append(vor.vertices[np.r_[rid, [rid[0]]]]) 179 | # bz_facets.append(vor.vertices[rid]) 180 | 181 | for pid, rid in zip(vor.ridge_points, vor.ridge_vertices): 182 | # WHY 13 ???? 183 | # The Voronoi ridges/facets are perpendicular to the lines drawn between the 184 | # input points. The 14th input point is [0, 0, 0]. 185 | if(pid[0] == 13 or pid[1] == 13): 186 | bz_ridges.append(vor.vertices[np.r_[rid, [rid[0]]]]) 187 | bz_facets.append(vor.vertices[rid]) 188 | bz_vertices += rid 189 | 190 | bz_vertices = list(set(bz_vertices)) 191 | 192 | return vor.vertices[bz_vertices], bz_ridges, bz_facets 193 | 194 | 195 | def get_primitive_cell_3d(cell): 196 | """ 197 | Get the vertices, lines and facets of the primitive cell. 198 | """ 199 | cell = np.asarray(cell, dtype=float) 200 | assert cell.shape == (3, 3) 201 | 202 | dx, dy, dz = np.mgrid[0:2, 0:2, 0:2] 203 | dxyz = np.c_[dx.ravel(), dy.ravel(), dz.ravel()] 204 | px, py, pz = np.tensordot(cell, [dx, dy, dz], axes=[0, 0]) 205 | points = np.c_[px.ravel(), py.ravel(), pz.ravel()] 206 | 207 | lines = [] 208 | faces = None 209 | 210 | for ii in range(len(points)): 211 | for jj in range(ii): 212 | if np.abs(dxyz[ii] - dxyz[jj]).sum() == 1: 213 | lines.append(np.vstack([points[ii], points[jj]])) 214 | 215 | return points, lines, faces 216 | 217 | 218 | class ebands3d(object): 219 | ''' 220 | ''' 221 | 222 | def __init__(self, eigvenval='EIGENVAL', efermi=None, kmesh=[], 223 | symprec=1E-5, 224 | poscar=None, 225 | kpoints=None): 226 | ''' 227 | Init 228 | ''' 229 | 230 | self._fname = eigvenval 231 | # the directory containing the input file 232 | self._dname = os.path.dirname(eigvenval) 233 | if self._dname == '': 234 | self._dname = '.' 235 | if poscar is None: 236 | self.poscar = self._dname + '/POSCAR' 237 | if kpoints is None: 238 | self.kpoints = self._dname + '/KPOINTS' 239 | 240 | # read bands, k-points of the irreducible Brillouin Zone 241 | self.read_eigenval() 242 | # set the Fermi energy 243 | self.set_efermi(efermi) 244 | # set the k-points mesh 245 | self.set_kmesh(kmesh) 246 | # read POSCAR 247 | self.atoms = read(self.poscar) 248 | # create the grid to ir map 249 | self.ir_kpts_map(symprec=symprec) 250 | # 251 | self.get_fermi_ebands3d() 252 | 253 | def interpolate_ebands3d(self, mesh): 254 | ''' 255 | Interpolate the band energies in the primitive cell using zero-padding fft. 256 | ''' 257 | 258 | # Somehow, only numpy irfft/rfft support fft interpolation. 259 | from numpy.fft import irfftn, rfftn, irfft, rfft 260 | 261 | mesh = np.asarray(mesh, dtype=int) 262 | assert mesh.shape == (3,), "Invalid dimension of new mesh!" 263 | nx, ny, nz = self.kmesh 264 | nnx, nny, nnz = mesh 265 | 266 | ebands3d_uc_interp = [] 267 | for ispin in range(self.nspin): 268 | uc_tmp = [] 269 | for b3d in self.fermi_ebands3d_uc[ispin]: 270 | b0 = irfft(rfft(b3d), nnz, axis=2) * nnz / nz # 1. interpolate along the z-axis 271 | b1 = np.swapaxes(b0, 1, 2) # 2. interpolate along the y-axis 272 | b2 = irfft(rfft(b1), nny, axis=2) * nny / ny 273 | b3 = np.swapaxes(b2, 0, 2) # 2. interpolate along the x-axis 274 | b4 = irfft(rfft(b3), nnx, axis=2) * nnx / nx 275 | b5 = np.swapaxes( 276 | np.swapaxes(b4, 0, 2), 1, 2 277 | ) 278 | 279 | # b3d_interp = irfftn( 280 | # rfftn(b3d), s=mesh, axes=(0,1,2) 281 | # ) * np.prod(mesh) / np.prod(self.kmesh) # only keep the real part, the imaginary part is supposed to be small 282 | uc_tmp.append(b5.copy()) 283 | ebands3d_uc_interp.append(uc_tmp) 284 | 285 | self.kmesh = mesh 286 | self.fermi_ebands3d_uc = ebands3d_uc_interp 287 | 288 | 289 | def get_fermi_ebands3d(self): 290 | ''' 291 | For those bands that corss the Fermi level, unfold the band energies on 292 | the irreducible BZ onto the whole reciprocal primitive cell. 293 | ''' 294 | 295 | # band energies of the k-points within the primitive cell 296 | self.fermi_ebands3d_uc = [] 297 | # band energies of the k-points within the Brillouin Zone 298 | # self.fermi_ebands3d_bz = [] 299 | 300 | # nx, ny, nz = self.kmesh 301 | for ispin in range(self.nspin): 302 | uc_tmp = [] 303 | # bz_tmp = [] 304 | for iband in self.fermi_xbands[ispin]: 305 | # the band energies of the k-points within primitive cell 306 | etmp = self.ir_ebands[ispin, self.grid_to_ir_map, iband] 307 | etmp.shape = list(reversed(self.kmesh)) # the first axis of etmp corresponds to the last cell direction 308 | # x-index in the "grid" of ir_kpts_map runs fastest 309 | etmp = np.swapaxes(etmp, 0, 2) # swap axes to turn back the order 310 | # # make the band energies periodic in the primitive cell 311 | # etmp = np.tile(etmp, (2,2,2))[:nx+1, :ny+1, :nz+1] 312 | uc_tmp.append(etmp) 313 | 314 | # # the band energies of the k-points within Brillouin Zone 315 | # btmp = np.tile(etmp, (2,2,2)) 316 | # s = btmp.shape 317 | # btmp.shape = (btmp.size) 318 | # # set the band energies of the k-points outside BZ to a large 319 | # # one so that the energy isosurface will not extent outside 320 | # # beyond the BZ. 321 | # btmp[np.logical_not(self.bz_in_kgrid_2uc)] = self.emax + 100. 322 | # btmp.shape = s 323 | # bz_tmp.append(btmp) 324 | 325 | self.fermi_ebands3d_uc.append(uc_tmp) 326 | # self.fermi_ebands3d_bz.append(bz_tmp) 327 | 328 | # periodic band energies, mesh size +1 329 | # self.kmesh = [nx+1, ny+1, nz+1] 330 | 331 | def to_bxsf(self, prefix='ebands3d', ncol=6): 332 | ''' 333 | Output the ebands3d as Xcrysden .bxsf format. 334 | ''' 335 | 336 | with open('{:s}.bxsf'.format(prefix), 'w') as out: 337 | out.write("BEGIN_INFO\n") 338 | out.write(" # Launch as: xcrysden --bxsf ebands3d.bxsf\n") 339 | out.write(" Fermi Energy: {:12.6f}\n".format(self.efermi)) 340 | out.write("END_INFO\n\n") 341 | 342 | out.write("BEGIN_BLOCK_BANDGRID_3D\n") 343 | out.write(" band_energies\n") 344 | out.write(" BANDGRID_3D_BANDS\n") 345 | 346 | # the number of bands that corss the Fermi level 347 | number_fermi_xbands = sum([len(xx) for xx in self.fermi_xbands]) 348 | out.write(" {:d}\n".format(number_fermi_xbands)) 349 | # number of data-points in each direction (i.e. nx ny nz for 3D gr) 350 | out.write(" {:5d}{:5d}{:5d}\n".format(*(x for x in self.kmesh))) 351 | # origin of the bandgrid. 352 | # Warning: origin should be (0,0,0) (i.e. Gamma point) 353 | out.write(" {:16.8f}{:16.8f}{:16.8f}\n".format(0.0, 0.0, 0.0)) 354 | # Reciprocal lattice vector 355 | out.write( 356 | '\n'.join([" " + ''.join(["%16.8f" % xx for xx in row]) 357 | for row in self.atoms.cell.reciprocal()]) 358 | ) 359 | 360 | # write the band grid for each band, the values inside a band grid 361 | # are specified in row-major (i.e. C) order. This means that values 362 | # are written as: 363 | # 364 | # C-syntax: 365 | # for (i=0; i