Hi. I am trying to get eigenvalues of a Hamiltonian operator. The problem is, the operator is a 2x2 matrix containing derivatives, and I don't know how to implement that. The equation is given by HcPsi = Energija*Psi, where Psi is a vector of length 2, dependent of z, and Hc is a 2x2 matrix containing first and second derivatives with respect to z, and functions of z, and Energija. The idea is to use NDEigensystem to get eigenvalues of Hc, which are the values of Energija, and eigenvectors which are the values of components of Psi. The code included is very long, and most of it are definitions of various functions and constants. Hc is defined in the last line. Thank you!
Milan Vicentijevic
Attachments:
