next up previous contents
Next: Debugging Up: Schrödinger energy solver Previous: rtbisp480()

func_Schrodinger()

double func_Schrodinger(double E,
                        int match, double V[], double r[], double dr[], int N)

Return value:

Input:

Design a routine func_Schrodinger() with the above prototype which calls upon schint() to integrate Schrödinger's equation from k1=0 to k2=N using Psi=0 and Psip=1 as initial conditions, and then returns the final value of Psi. This routine is most appropriate for the file ``physics.c''. (Don't forget to put the prototype in ``physics.h''!)

Hint: You will have to define an internal variable array Psiout[] to send into schint(). You should accomplish this dynamically in func_Schrodinger() just as you did for T[] in schint().



Tomas Arias
Sun Feb 11 18:46:55 EST 2001