/* Steven Andrews, 1/00 */
/* See documentation called Quantum doc */
/* Copyright 2003 by Steven Andrews.  Permission is granted
   for non-commercial use of and modifications to the code. */

float psiHO(float x,int n,float k,float m);
complex Bracket(float *psib,float *psik,float *xr,int n);
float NormalKet(float *psi,float *xr,int n);
complex Hamiltonian(float *psib,float *v,float *psik,float *xr,int n,float mass);
complex Dipole(float *psib,float *psik,float *xr,int n,float q);
void EQMFschrod(phptr u,void *k,phptr dudt);
int TKFplotpsi(float t,phptr u,void *tkptr);
float *TimeEvolve(float *psi,float *ur,float m,float *xr,int n,float Dt);
int TKFdipcorr(float t,phptr u,void *tkptr);
void DipCorr(float *psib,float *psik,float *ur,float m,float q,float *xr,int n,float Dt,float *dip);
float EigenketAHO(float *xr,float *ur,float *psi,int n,float k,float m);
sptr DipTransform(float *dip,float dt,int nt,float wmin);