# Homework Help: Solving Differential Equation numerically

1. Mar 7, 2016

### aaaa202

1. The problem statement, all variables and given/known data
I am supposed to write a script that can solve the Schrödinger equation on a nonuniform grid.

2. Relevant equations
Finite element approximation to the second derivative as in:

3. The attempt at a solution
I have defined a grid x = [x1,x2,x3,....,x4] with nonuniform spacing and a potential V(x). Specifically the nonuniform spacing is such that the first 10 points have a spacing of 10-10 and the others have 10-9.
I have then solved the Schrödinger equation by finding eigenvectors of:
H = -ħ2/2m D + V

The problem is that these have some kind of weird oscillatory behaviour, which stems from the discontinuity in the distance between the grid points. Do anyone have an idea what I could be doing wrong?

File size:
7.6 KB
Views:
66
2. Mar 7, 2016

### Simon Bridge

How did you construct D?

3. Mar 7, 2016

4. Mar 9, 2016

### aaaa202

I think if you use an approximation for the second derivative based on a nonuniform mesh, will this create some problems for the eigenvalue equation? Is it still valid that the eigenfunctions of the second derivative is just the eigenvectors of the matrix that represents it?