diff options
Diffstat (limited to 'hw8/10-3-new.jl')
| -rw-r--r-- | hw8/10-3-new.jl | 32 |
1 files changed, 0 insertions, 32 deletions
diff --git a/hw8/10-3-new.jl b/hw8/10-3-new.jl index fb819b8..42e7162 100644 --- a/hw8/10-3-new.jl +++ b/hw8/10-3-new.jl @@ -8,41 +8,9 @@ xmax = 6.5 N = 2000 x = collect(range(xmin, stop = xmax, length = N + 1)) # one extra, want 1000 spaces dx = x[2] - x[1] -println("dx = ", dx) -p = 0.0 -V0 = 2 * p^2 / (2.0 * m) -sig = 0.25 -x0 = 2.0 # a bit away to not interfere with the walls - -# setup the potential barrier -# V = 0 .* x -# # println("V = ", length(V)) -# for i in 1:length(V) -# V[i] = abs(x[i])^n -# end - -# make the initial wavefunction -psi0 = exp.(-((x[2:end-1] .+ x0) .^ 2) / (2 * sig^2)) .* exp.(1.0im * p * x[2:end-1]) C = sum(abs2.(psi0) .* dx) psi0 = psi0 / sqrt(C) -# plot(x[2:end-1], abs.(psi0), xlabel = "x", ylabel = "Re(psi)", title = "Initial wavefunction", legend = :none) - -# # make the Hamiltonian -# A = (hbar^2 / (2 * m * dx^2)) -# # main diagonal -# d = A * ones(N - 1) + V[2:end-1] # add potential -# # upper and lower diagonals -# dl = (A / -2.0) * ones(N - 2) -# du = (A / -2.0) * ones(N - 2) -# H = Tridiagonal(dl, d, du) - -# # diagonalize matrix -# E, psi = eigen(H) -# psi1 = psi[:, 1] - -# print the energy levels -# println("Energy levels: ", E[1:10]) function get_energy_levels_for_n(n, pot_plot) V = 0 .* x |