Comparing discretization techniques
In Periodic problems and plane-wave discretizations we saw how simple 1D problems can be modelled using plane-wave basis sets. This example invites you to work out some details on these aspects yourself using a number of exercises. The solutions are given at the bottom of the page.
For this example we consider the discretization of
\[ H = - \frac12 Δ + V(x) \quad \text{with $V(x) = \cos(x)$}\]
on $[0, 2π]$ with periodic boundary conditions. The $\cos(x)$ takes the role of a lattice-periodic potential. We will be interested in the smallest eigenvalues of this discretized Hamiltonian. Of note the boundary condition matters: The spectrum we will get is different from e.g. the spectrum of $H$ on $\mathbb{R}$.
Finite differences
We approximate functions $ψ$ on $[0, 2\pi]$ by their values at grid points $x_k = 2\pi \frac{k}{N}$, $k=1, \dots, N$. The boundary conditions are imposed by $ψ(x_0) = ψ(x_N), ψ(x_{N+1}) = ψ(x_1)$. We then have
\[ \big(Hψ\big)(x_k) \approx \frac 1 2 \frac{-ψ_{k-1} + 2 ψ_k - ψ_{k+1}}{δx^2} + V(x_k) ψ(x_k)\]
with $δx = \frac{2π}{N}$.
This can be put in matrix form in the following way:
# Finite differences Hamiltonian -1/2 Delta + V on [0, 2pi] with periodic bc.
# Pass it a function V.
using LinearAlgebra
function build_finite_differences_matrix(Vfunction, N::Integer)
δx = 2π/N
# Finite-difference approximation to -½Δ
T = 1/(2δx^2) * Tridiagonal(-ones(N-1), 2ones(N), -ones(N-1))
# The type Tridiagonal is efficient, but to establish the periodic boundary conditions
# we need to add elements not on the three diagonals, so convert to dense matrix
T = Matrix(T)
T[1, N] = T[N, 1] = -1 / (2δx^2)
# Finite-difference approximation to potential: We collect all coordinates ...
x_coords = [k * δx for k=1:N]
V = Diagonal(Vfunction.(x_coords)) # ... and evaluate V on each of the x_coords
T + V
end;Show that the finite-difference approximation of -½Δ is indeed an approximation of the second derivative. Obtain an estimate of the first eigenvalue of $H$. Hint: Take a look at the eigen function from LinearAlgebra.
Plane waves method
In this method, we expand states on the basis
\[ e_G(x) = \frac{1}{\sqrt{2\pi}} e^{iGx} \qquad \text{for $G=-N,\dots,N$}.\]
Show that
\[ \langle e_G, e_{G'}\rangle = ∫_0^{2π} e_G^\ast(x) e_{G'}(x) d x = δ_{G, G'}\]
and (assuming $V(x) = \cos(x)$)
\[ \langle e_G, H e_{G'}\rangle = \frac 1 2 \left(|G|^2 \delta_{G,G'} + \delta_{G, G'+1} + \delta_{G, G'-1}\right).\]
What happens for a more general $V(x)$?
Code this and check the first eigenvalue agrees with the finite-difference case. Compare accuracies at various basis set sizes $N$.
Using DFTK
We now use DFTK to do the same plane-wave discretization in this 1D system. To deal with a 1D case we use a 3D lattice with two lattice vectors set to zero.
using DFTK
a = 2π
lattice = a .* [[1 0 0.]; [0 0 0]; [0 0 0]];Define Hamiltonian: Kinetic + Potential
terms = [Kinetic(),
ExternalFromReal(r -> cos(r[1]))] # r is a vector of size 3
model = Model(lattice; n_electrons=1, terms, spin_polarization=:spinless); # One spinless electronEcut defines the number of plane waves by selecting all those $G$, which satisfy the relationship $½ |G|^2 ≤ \text{Ecut}$.
Ecut = 500
basis = PlaneWaveBasis(model; Ecut, kgrid=(1, 1, 1))PlaneWaveBasis discretization:
architecture : DFTK.CPU()
num. mpi processes : 1
num. julia threads : 1
num. DFTK threads : 1
num. blas threads : 2
num. fft threads : 1
Ecut : 500.0 Ha
fft_size : (135, 1, 1), 135 total points
kgrid : MonkhorstPack([1, 1, 1])
num. red. kpoints : 1
num. irred. kpoints : 1
Estimated memory usage (per MPI process):
single ψ : 3.94KiB
single ρ : 1.05KiB
peak memory (SCF) : 44.2KiB
Discretized Model(custom, 1D):
lattice (in Bohr) : [6.28319 , 0 , 0 ]
[0 , 0 , 0 ]
[0 , 0 , 0 ]
unit cell volume : 6.2832 Bohr
num. electrons : 1
spin polarization : spinless
temperature : 0 Ha
terms : Kinetic()
ExternalFromReal(Main.var"##363".var"#3#4"())We now seek the ground state using the self-consistent field algorithm.
scfres = self_consistent_field(basis; tol=1e-4)
scfres.energiesEnergy breakdown (in Ha):
Kinetic 0.2090845
ExternalFromReal -0.7441493
total -0.535064852288On this simple linear (non-interacting) model, the SCF converges in one step. The ground state energy of is simply the lowest eigenvalue; it should match the smallest eigenvalue of $H$ computed above.
Plotting
We can also get the first eigenvector (in the plane wave basis) and plot it
using Plots
ψ_fourier = scfres.ψ[1][:, 1]; # first k-point, all G components, first eigenvector
# Transform the wave function to real space
ψ = ifft(basis, basis.kpoints[1], ψ_fourier)[:, 1, 1]
# Eigenvectors are only defined up to a phase. We fix it by imposing that psi(0) is real
ψ /= (ψ[1] / abs(ψ[1]))
plot(real(ψ); label="")Again this should match with the result above.
Look at the Fourier coefficients of ψ_fourier and compare with the result above.
The DFTK Hamiltonian
We can ask DFTK for the Hamiltonian
E, ham = energy_hamiltonian(basis, scfres.ψ, scfres.occupation; ρ=scfres.ρ)
H = ham.blocks[1]
typeof(H)DFTK.DftHamiltonianBlockThis is an opaque data structure, which encodes the Hamiltonian. What can we do with it?
using InteractiveUtils
methodswith(typeof(H), supertypes=true)9-element Vector{Method}:- Array(block::HamiltonianBlock) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:74
- *(H::HamiltonianBlock, ψ) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:63
- (Matrix)(block::HamiltonianBlock) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:75
- eltype(block::HamiltonianBlock) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:64
- size(block::HamiltonianBlock, i::Integer) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:65
- size(block::HamiltonianBlock) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:66
- mul!(Hψ::AbstractArray, H::DFTK.DftHamiltonianBlock, ψ::AbstractArray) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/terms/Hamiltonian.jl:137
- PreconditionerNone(::HamiltonianBlock) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/eigen/preconditioners.jl:21
- PreconditionerTPA(ham::HamiltonianBlock; kwargs...) in DFTK at /home/runner/work/DFTK.jl/DFTK.jl/src/eigen/preconditioners.jl:45
This defines a number of methods. For instance, it can be used as a linear operator:
H * DFTK.random_orbitals(basis, basis.kpoints[1], 1)63×1 Matrix{ComplexF64}:
-0.047560152954042444 + 0.05625609756786979im
-0.030476614320948267 + 0.04559753267556719im
0.04290726519236661 - 0.15222332456148094im
0.11680421404613703 - 0.5248399676452765im
-0.01472246737360134 + 0.15909839111953974im
0.5524474882176351 + 1.131405896003262im
-0.9518701129244916 - 1.6793157312953104im
-1.3513351698770542 + 1.9707951692049095im
-2.585522827194829 + 0.8790274741145313im
-2.1662103048266155 + 5.570026666941963im
-1.0406747007799073 + 8.322511078640348im
-6.623424081318346 - 1.117844323627062im
6.6171326793500596 + 3.8344474114334703im
-0.6718907317100934 - 7.424711576668564im
-3.3568719624966707 + 9.906070662177797im
8.267594790232605 + 6.637710364631848im
-2.820114491218241 - 4.651030786289675im
19.478675389489723 - 10.806711487905156im
-21.933451642350754 + 1.0198514015496314im
-17.646777983309605 + 39.7793311803421im
-25.183677962115876 - 26.99697781533597im
1.0737352706420913 + 5.880151106135673im
2.0876481665919746 - 10.432487611234544im
-4.9985631352735025 - 22.36321411774483im
-26.391933609783383 - 41.64229235433832im
-1.8854696509804671 + 14.125941351508056im
21.72503767681319 - 6.749588197963241im
-2.2946674094164683 + 2.241101507008985im
56.605355971341424 + 16.36304805582955im
-42.57740848583217 + 44.21609192468564im
-1.7513596778527782 + 7.651064256376149im
-45.469320295922124 + 48.542711183106306im
39.32131427577507 - 4.880821716290422im
20.64938627014884 + 5.065581273983556im
2.387408534029374 - 55.46824817400782im
67.58984033341093 - 76.98038163329413im
-14.167463743878244 + 43.066120079543786im
-10.861117311686346 + 28.144425641737662im
15.027698940110648 + 25.80205467058822im
32.03995479284555 - 5.128336643080557im
7.288095204696067 + 29.53434850283691im
-8.911634372577263 - 22.213365336499237im
-0.31254527889141104 + 8.52600975028564im
-6.56322743013501 + 17.67831532981195im
2.8050002800164053 - 0.34659257156884105im
11.681523862773302 - 13.420435350521545im
-22.14021331597608 - 7.300254975609509im
15.299996467842611 - 1.6276437654297673im
-2.192225385572011 - 4.429909020557902im
-2.6736426525516865 - 8.116918985402872im
-6.500306136021729 + 3.1398666589761417im
2.9824689690398904 - 7.270731722753138im
-1.9114366081548866 + 2.252747210713393im
10.377987621631037 - 10.071029365521664im
4.432587070797425 - 0.5005241260720572im
-0.8149431479059134 + 0.7384108588894804im
-1.4126207244033995 - 1.0989416910121121im
-1.3442144909943075 + 1.8144848927261574im
-0.5040629596967884 + 0.8784410874407214im
1.9411901935744504 - 1.763792492515834im
-0.1794547386710232 - 0.12922834677999276im
0.0036705600088862145 - 0.14060387980124608im
-0.12580653934766317 - 0.03564154607265484imWe can also get its full matrix representation:
Array(H)63×63 Matrix{ComplexF64}:
9.38036e-19+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im
0.5+2.37729e-18im 0.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im
6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 2.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im
-2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 4.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im
9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 8.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im
5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 12.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im
1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 18.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im
-6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 24.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im
6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 32.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im
5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 40.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im
-5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 50.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im
1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 60.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im
-3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 72.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im
-4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 84.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im
6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 98.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im
1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 112.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im
3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 128.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im
4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 144.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im
6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 162.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im
-6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 180.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im
-2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 200.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im
-4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 220.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im
-1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 242.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im
1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 264.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im
-4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 288.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im
-2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 312.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im
1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 338.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im
6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 364.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 1.25324e-17+7.3538e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im
-7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 392.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 2.40812e-18-1.91415e-17im 1.25324e-17+7.3538e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im
0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 420.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im 9.85983e-18+4.45332e-18im 2.40812e-18-1.91415e-17im 1.25324e-17+7.3538e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im
1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 450.0+0.0im 0.5-1.91936e-18im 4.3417e-18+1.90195e-17im 9.85983e-18+4.45332e-18im 2.40812e-18-1.91415e-17im 1.25324e-17+7.3538e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im -2.32476e-18-1.03486e-17im
-2.32476e-18-1.03486e-17im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 480.5+0.0im -6.21078e-18-8.87719e-18im 4.3417e-18+1.90195e-17im 9.85983e-18+4.45332e-18im 2.40812e-18-1.91415e-17im 1.25324e-17+7.3538e-18im 2.94354e-19-1.42818e-17im -6.80457e-18-1.65903e-18im -5.05513e-18+9.40828e-18im -1.24221e-17+9.18236e-18im 2.3641e-17+8.69902e-18im 7.46908e-18-4.89191e-18im -3.96529e-18-1.63017e-17im -6.5294e-18-5.68244e-18im -2.09005e-17+1.30029e-17im 1.30228e-18-1.05794e-17im 7.50769e-18+1.8538e-17im 7.79391e-18-1.28996e-18im -4.55526e-18-1.45669e-17im -2.07003e-17+3.24918e-17im 2.76834e-17+2.04324e-17im 1.26093e-17-3.72139e-17im -1.9652e-17-2.7462e-17im -5.50336e-18+2.58911e-17im 1.50263e-17+4.32004e-18im 2.08933e-18-9.22629e-18im -2.1533e-17+1.12123e-17im -1.02584e-17-3.88425e-20im 5.20786e-18+2.01367e-18im 1.19815e-17-3.24554e-18im -4.04378e-18+5.31047e-18im -5.18221e-18+5.03797e-20im
-8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 1.25324e-17-7.3538e-18im -4.44267e-19-4.10411e-18im 9.85983e-18-4.45332e-18im 8.49259e-18-1.66418e-17im -6.21078e-18+8.87719e-18im 480.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im
1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 1.25324e-17-7.3538e-18im -4.44267e-19-4.10411e-18im 9.85983e-18-4.45332e-18im 8.49259e-18-1.66418e-17im 0.5+2.37729e-18im 450.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im
-8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 1.25324e-17-7.3538e-18im -4.44267e-19-4.10411e-18im 9.85983e-18-4.45332e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 420.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im
-7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 1.25324e-17-7.3538e-18im -4.44267e-19-4.10411e-18im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 392.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im
6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 1.25324e-17-7.3538e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 364.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im
1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 2.94354e-19+1.42818e-17im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 338.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im
-2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.80457e-18+1.65903e-18im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 312.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im
-4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im -5.05513e-18-9.40828e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 288.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im
1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im -1.24221e-17-9.18236e-18im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 264.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im
-1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im 2.3641e-17-8.69902e-18im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 242.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im
-4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 7.46908e-18+4.89191e-18im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 220.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im
-2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.96529e-18+1.63017e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 200.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im
-6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -6.5294e-18+5.68244e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 180.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im
6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im -2.09005e-17-1.30029e-17im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 162.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im
4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.30228e-18+1.05794e-17im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 144.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im
1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 7.50769e-18-1.8538e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 128.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im
1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 3.42423e-18+1.3351e-17im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 112.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im
6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im -4.55526e-18+1.45669e-17im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 98.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im
-4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -7.79391e-18-1.28996e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 84.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im
-3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im 2.76834e-17-2.04324e-17im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 72.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im
1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im 1.26093e-17+3.72139e-17im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 60.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im
-5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.9652e-17+2.7462e-17im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 50.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im
5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im -5.50336e-18-2.58911e-17im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 40.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im
6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im 1.50263e-17-4.32004e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 32.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im
-6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im 2.08933e-18+9.22629e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 24.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im
1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im -2.1533e-17-1.12123e-17im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 18.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im
5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im -1.02584e-17+3.88425e-20im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 12.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im
9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im 5.20786e-18-2.01367e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 8.0+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im
-2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 1.19815e-17+3.24554e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 4.5+0.0im 0.5-1.91936e-18im 6.84621e-18-4.62374e-18im
6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im -4.04378e-18-5.31047e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 2.0+0.0im 0.5-1.91936e-18im
0.5-1.91936e-18im 6.84621e-18-4.62374e-18im -2.10986e-18-7.33721e-18im 9.53855e-19+2.05757e-17im 5.30139e-18+9.96275e-18im 1.15239e-17-1.06252e-18im -6.09606e-18-2.01355e-17im 6.73916e-18-9.6349e-18im 5.08456e-18-5.8805e-19im -5.68698e-18+2.3003e-19im 1.48741e-17+1.10137e-17im -3.44396e-18-3.19427e-17im -4.12338e-17-4.82849e-18im 6.96319e-18+3.06423e-17im 1.41552e-18-4.45966e-18im 1.06844e-17+1.07046e-17im 4.92777e-18+5.57465e-18im 6.87119e-18-1.51362e-19im -6.53743e-18-2.77989e-18im -2.28337e-18-2.65376e-18im -4.23963e-18+3.77063e-18im -1.07755e-17-3.20506e-18im 1.64216e-17-3.33489e-18im -4.8805e-18-2.16383e-18im -2.56782e-17-4.87146e-18im 1.35838e-17+5.69649e-18im 6.55612e-18+3.83412e-18im -7.57355e-18-5.37897e-18im -8.89745e-18+9.6317e-18im 1.13146e-17+6.44402e-18im -8.98722e-18-7.74616e-18im -5.18221e-18-5.03797e-20im 1.13146e-17-6.44402e-18im 0.0+6.78716e-18im -7.57355e-18+5.37897e-18im 6.55612e-18-3.83412e-18im 1.35838e-17-5.69649e-18im -2.56782e-17+4.87146e-18im -4.8805e-18+2.16383e-18im 1.64216e-17+3.33489e-18im -1.07755e-17+3.20506e-18im -4.23963e-18-3.77063e-18im -2.28337e-18+2.65376e-18im -6.53743e-18+2.77989e-18im 6.87119e-18+1.51362e-19im 4.92777e-18-5.57465e-18im 3.61585e-19-1.24616e-17im 1.41552e-18+4.45966e-18im 6.5791e-18-3.57518e-18im -4.12338e-17+4.82849e-18im -3.44396e-18+3.19427e-17im 1.48741e-17-1.10137e-17im -5.68698e-18-2.3003e-19im 5.08456e-18+5.8805e-19im 6.73916e-18+9.6349e-18im -6.09606e-18+2.01355e-17im 1.15239e-17+1.06252e-18im 5.30139e-18-9.96275e-18im 9.53855e-19-2.05757e-17im -2.10986e-18+7.33721e-18im 6.84621e-18+4.62374e-18im 0.5+2.37729e-18im 0.5+0.0imCompare this matrix Array(H) with the one you obtained in Exercise 3, get its eigenvectors and eigenvalues. Try to guess the ordering of $G$-vectors in DFTK.
Increase the size of the problem, and compare the time spent by DFTK's internal diagonalization algorithms to a full diagonalization of Array(H). Hint: The @belapsed and @benchmark macros (from the BenchmarkTools package) are handy for this task. Note that there are some subtleties with global variables (see the BenchmarkTools docs for details). E.g. to use it to benchmark a function like eigen(H) run it as (note the $):
using BenchmarkTools
@benchmark eigen($H)Solutions
Exercise 1
If we consider a function $f : [0, 2π] → \mathbb{R}$, to first order
\[f(x + δx) = f(x) + δx f'(x) + O(δx^2)\]
therefore after rearrangement
\[f'(x) = \frac{f(x + δx) - f(x)}{δx} + O(δx).\]
Similarly
\[f''(x) = \frac{f'(x + δx) - f'(x)}{δx} + O(δx),\]
such that overall
\[f''(x) \simeq \frac{f(x + 2δx) - f(x + δx) - f(x + δx) + f(x)}{δx^2} = \frac{f(x + 2δx) - 2f(x + δx) + f(x)}{δx^2}\]
In finite differences we consider a stick basis of vectors
\[\left\{ e_i = (0, …, 0, \underbrace{δx}_\text{$i$-th position}, 0, …, 0) \middle| i = 1, … N \right\}.\]
Keeping in mind the periodic boundary conditions (i.e. $e_0 = e_N$) projecting the Hamiltonian $H$ onto this basis thus yields the proposed structure.
We start off with $N = 100$ to obtain
Hfd = build_finite_differences_matrix(cos, 100)
L, V = eigen(Hfd)
L[1:5]5-element Vector{Float64}:
-0.5351640695081414
0.3429576842587494
0.8532206682511969
2.053679933442362
2.078512110605267This is already pretty accurate (to about 4 digits) as can be estimated looking at the following convergence plot:
function fconv(N)
L, V = eigen(build_finite_differences_matrix(cos, N))
first(L)
end
Nrange = 10:10:100
plot(Nrange, abs.(fconv.(Nrange) .- fconv(200));
yaxis=:log, legend=false, mark=:x, xlabel="N", ylabel="Absolute error")Exercise 2
We note that
\[\langle e_G, e_{G'}\rangle = ∫_0^{2π} e_G^\ast(x) e_{G'}(x) d x = 1/2π ∫_0^{2π} e^{i(G'-G)x} d x\]
Since $e^{iy}$ is a periodic function with period $2\pi$, $\int_0^{2\pi} e^{i m y} = \delta_{0,m}$. Therefore if $G≠G'$ we have that $\langle e_G, e_{G'}\rangle = 0$, while $G=G'$ implies $\langle e_G, e_{G'}\rangle = 1$. In summary:
\[\langle e_G, e_{G'}\rangle = δ_{G, G'}\]
Next fo $V(x) = \cos(x)$ we obtain
\[\langle e_G, H e_{G'}\rangle = \frac 1 2 ∫_0^{2π} e_G^\ast(x) H e_{G'}(x) d x\]
We start by applying the Hamiltonian to a plane-wave:
\[H e_{G'}(x) = - \frac 1 2 (-|G|^2) \frac 1 {\sqrt{2π}} e^{iG'x) + cos(x) \frac 1 {\sqrt{2π}} e{iG'x}\]
Then, using the result of the first part of the exercise and the fact that $cos(x) = \frac 1 2 \left(e{ix} + e{-ix}\right)$, we get:
\[\begin{align*} ⟨ e_G, H e_{G'}⟩ &= \frac 1 2 G^2 δ_{G, G'} + \frac 1 {4π} \left(∫_0^{2π} e^{ix ⋅ (G'-G+1)} d x + ∫_0^{2π} e^{ix ⋅ (G'-G-1)} d x \right) \\ &= \frac 1 2 \left(|G|^2 \delta_{G,G'} + \delta_{G, G'+1} + \delta_{G, G'-1}\right) \end{align*}\]
In case a more general $V(x)$ was employed, this potential still has to be periodic over $[0, 2\pi]$ to fit our setting. Assuming sufficient regularity in $V$ we can employ a Fourier series:
\[V(x) = \sum_{G=- \infty}^{\infty} \hat{V}_G e_G(x)\]
where
\[\hat{V}_G = \frac{1}{\sqrt{2π}} ∫_0^{2π} V(x) e^{-iGx} dx = ∫_0^{2π} V(x) e_G^\ast dx .\]
Note that one can change of this as a change of basis from the position basis to the plane-wave basis.
Based on this expansion
\[\begin{align*} ⟨ e_G, V e_{G'} ⟩ &= \left\langle e_G, ∑_{G''} \hat{V}_{G''} \, e_{G'+G''} \right\rangle \\ &= \sum_{G''=-\infty}^\infty \hat{V}_{G''} ⟨ e_G, e_{G'+G''} ⟩ \\ &= \sum_{G''=-\infty}^\infty \hat{V}_{G''} \, δ_{G-G', G''} ⟩ \\ &= \hat{V}_{G-G'} \end{align*}\]
and therefore
\[⟨ e_G, H e_{G'} ⟩ = \frac 1 2 |G|^2 \delta_{G,G'} + \hat{V}_{G-G'},\]
i.e. essentially the Fourier transform of $V$ determines the contribution to the matrix elements of the Hamiltonian.
Exercise 3
The Hamiltonian matrix for the plane waves method can be found this way:
# Plane waves Hamiltonian -½Δ + cos on [0, 2pi].
function build_plane_waves_matrix_cos(N::Integer)
# Plane wave approximation to -½Δ
Gsq = [float(i)^2 for i in -N:N]
# Hamiltonian as derived in Exercise 2:
1/2 * Tridiagonal(ones(2N), Gsq, ones(2N))
end;Then we check that the first eigenvalue agrees with the finite-difference case, using $N = 10$:
Hpw_cos = build_plane_waves_matrix_cos(10)
L, V = eigen(Hpw_cos)
L[1:5]5-element Vector{Float64}:
-0.5350648522878154
0.34336012839908336
0.8536343543207997
2.0565044112660984
2.081227363352144We look at the convergence plot to compare the accuracy for various numbers of plane-waves $N$:
function fconv(N)
L = eigvals(build_plane_waves_matrix_cos(N))
first(L)
end
Nrange = 2:10
plot(Nrange, abs.(fconv.(Nrange) .- fconv(200)); yaxis=:log, legend=false,
ylims=(1e-16,Inf), ylabel="Absolute error", xlabel="N", mark=:x)Notice how compared to exercise 1 the considered basis size $N$ is much smaller, indicating that plane-wave methods more quickly lead to accurate solutions than finite-difference methods.
Exercise 4
For efficiency reasons the data in Fourier space is not ordered increasingly with $G$. Therefore to plot the Fourier space representation sensibly, we need to sort by ascending values of the $G$ vectors first. For this we extract the Fourier vector of each plane-wave basis function in the index order:
coords_G_vectors = G_vectors_cart(basis, basis.kpoints[1]) # Get coordinates of first and only k-point
# Only keep first component of each vector (because the others are zero for 1D problems):
coords_Gx = [G[1] for G in coords_G_vectors]
p = plot(coords_Gx, real(ψ_fourier); label="real part", xlims=(-10, 10))
plot!(p, coords_Gx, imag(ψ_fourier); label="imaginary part")The plot is symmetric about the zero (confirming that the orbitals are real) and only takes peaked values, which corresponds to the expected result for a cosine potential.
Exercise 5
To figure out the ordering we consider a small basis and build the Hamiltonian:
basis_small = PlaneWaveBasis(model; Ecut=5, kgrid=(1, 1, 1))
ham_small = Hamiltonian(basis_small)
H_small = Array(ham_small.blocks[1])
H_small[abs.(H_small) .< 1e-12] .= 0 # Drop numerically zero entries
H_small7×7 Matrix{ComplexF64}:
0.0+0.0im 0.5+1.00424e-16im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.5-1.00424e-16im
0.5-1.00424e-16im 0.5+0.0im 0.5+1.00424e-16im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im
0.0+0.0im 0.5-1.00424e-16im 2.0+0.0im 0.5+1.00424e-16im 0.0+0.0im 0.0+0.0im 0.0+0.0im
0.0+0.0im 0.0+0.0im 0.5-1.00424e-16im 4.5+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im
0.0+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im 4.5+0.0im 0.5+1.00424e-16im 0.0+0.0im
0.0+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.5-1.00424e-16im 2.0+0.0im 0.5+1.00424e-16im
0.5+1.00424e-16im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.0+0.0im 0.5-1.00424e-16im 0.5+0.0imThe equivalent version using the build_plane_waves_matrix_cos function is N=3 (both give rice to a 7×7 matrix).
Hother = build_plane_waves_matrix_cos(3)7×7 LinearAlgebra.Tridiagonal{Float64, Vector{Float64}}:
4.5 0.5 ⋅ ⋅ ⋅ ⋅ ⋅
0.5 2.0 0.5 ⋅ ⋅ ⋅ ⋅
⋅ 0.5 0.5 0.5 ⋅ ⋅ ⋅
⋅ ⋅ 0.5 0.0 0.5 ⋅ ⋅
⋅ ⋅ ⋅ 0.5 0.5 0.5 ⋅
⋅ ⋅ ⋅ ⋅ 0.5 2.0 0.5
⋅ ⋅ ⋅ ⋅ ⋅ 0.5 4.5By comparing the entries we find the ordering is 0,1,2,...,-2,-1, which can also be found by inspecting
first.(G_vectors(basis_small, basis_small.kpoints[1]))7-element Vector{Int64}:
0
1
2
3
-3
-2
-1Both matrices have the same eigenvalues:
maximum(abs, eigvals(H_small) - eigvals(Hother))1.7763568394002505e-15and in the eigenvectors we find the same rearrangements in the entries of the eigenvectors of both matrices, matching the DFTK ordering of is 0,1,2,...,-2,-1.
eigvecs(Hother)[:, 1]7-element Vector{Float64}:
-0.008461095314567344
0.08520425570524764
-0.42353481192309905
0.7915641575577494
-0.42353481192310566
0.08520425570524905
-0.008461095314567471eigvecs(H_small)[:, 1]7-element Vector{ComplexF64}:
-0.7915641575577493 + 1.5898485476322146e-16im
0.42353481192310877 - 1.7013307112974292e-16im
-0.08520425570525011 + 5.1339563909648e-17im
0.008461095314567624 - 6.797609502050578e-18im
0.008461095314567069 + 3.3988047510250303e-18im
-0.08520425570524359 - 1.711318796988132e-17im
0.4235348119230965 + 0.0imNotice, that eigenvectors are only defined up to a phase, so the sign may globally be inverted between the two eigenvectors.
Exercise 6
We benchmark the time needed for a full diagonalization (instantiation of the Array plus call of eigen) versus the time needed for running the SCF (i.e. iterative diagonalization using plane waves).
using Printf
for Ecut in 200:200:1600
basis_time = PlaneWaveBasis(model; Ecut, kgrid=(1, 1, 1))
t_eigen = @elapsed eigen(Array(Hamiltonian(basis_time).blocks[1]))
t_scf = @elapsed self_consistent_field(basis_time; tol=1e-6, callback=identity);
@printf "%4i eigen=%8.4fms scf=%8.4fms\n" Ecut 1000t_eigen 1000t_scf
end 200 eigen= 2.2855ms scf= 4.7796ms
400 eigen= 4.1037ms scf= 2.6723ms
600 eigen= 6.1691ms scf= 2.8989ms
800 eigen= 9.6372ms scf= 2.9078ms
1000 eigen= 11.7182ms scf= 2.8268ms
1200 eigen= 13.9919ms scf= 3.0490ms
1400 eigen= 16.6610ms scf= 3.1162ms
1600 eigen= 19.7095ms scf= 2.9872ms