Solver

This article is stub, waiting for the new class Numerov.cpp

Brief Explanation of Schroedinger.cpp

trap_array

Numerical integration of a function using the Trapezoidal rule.

f_sol_Numerov

solve_Numerov

solve_Numerov solves the pointed potential using the Numerov algorithm and renormalizing the output wavefunction to 1. To do this it must try the solutions for different energies. The natural solution to the second degree differential equation is the exponential. But my boundary conditions impose 0 at both beginning and end of the wavefunction, so you have to try until you find such solution by finding where the exponential solution changes sign.

bisec_Numerov

Applies a bisection algorith to the numerov method to find the energy that gives the non-trivial (non-exponential) solution with the correct boundary conditions (@param wavefunction[0] == @param wavefunction[@param nbox] == 0)