# Splitting 2nd order DE into two 1st order DE's

1. Apr 25, 2010

### markov4

1. The problem statement, all variables and given/known data

1.a) Taking the Schrodinger second order differential equation given, split it into two first order differential equations for numerical solving.
We're given the relevant constants in a table, such as $$\hbar$$, $$\alpha$$, etc. 'z' is taken to be an indepedent variable which we incorporate into our initial conditions (which we must state).
3. The attempt at a solution

Have i done this right? Is it right to say that $$\Psi(z)$$=vz ? Since the $$\Psi(z)$$ that is in the equation, the one multiplied by ..+[U(z)-En]$$\Psi(z)$$ has to be set to something too (?).

#### Attached Files:

• ###### 15r1sm0.jpg
File size:
53 KB
Views:
5
2. Apr 25, 2010

### gabbagabbahey

No, just leave $\Psi(z)$ as $\Psi(z)$

It's true that $d\Psi=vdz$, but until you solve for $v(z)$ you have no way of integrating $vdz$.

3. Apr 25, 2010

### markov4

^^But for the second part of the question, we're required to "Write a program that uses an inbuilt numerical solver (such as Matlab's ode45) to numerically solve for $$\Psi(z)$$ for z1<z<z2..." so we have to choose our own initial conditions to incorporate into our program, along with using ode45 in matlab.

So then does that mean $$\Psi(z)$$=vz still? I don't understand why you would integrate it though? Unless i've completely missed something in the numerical solving part of this course.

4. Apr 25, 2010

### Gigasoft

v depends on z, so you can't integrate $$\frac{d\Psi}{dz}=v$$ to get $$\Psi=vz$$.

5. Apr 25, 2010

### markov4

^^So when i put these coupled equations into matlab, what do i define $$\Psi(z)$$ to be?

6. Apr 25, 2010

### gabbagabbahey

It's one of the functions you are trying to solve your system of equations for...you simply leave it as $\Psi(z)$

7. Apr 25, 2010

### markov4

^^I see then. So i'm assuming that ODE45 in matlab will accomodate for leaving $$\Psi(z)$$ in the equation? I'm reading up on it now (ODE45) but i don't think i understand it.

8. Apr 26, 2010

### markov4

Ok after doing some reading on ode45, i think i get it. So yes i'll leave Psi(z) as it is, and implementing it into my matlab code, it will be solved according to the initial conditions i set. Just need to try it out in matlab now!