Tutorial - 1D finite square well

This m-script demonstrates a split-operator method by calculating wave functions and an energy spect

Description

Simulate a time dependent solution of Schrodinger's equation and calculate energy spectrum of a 1D finite square well using split-operator method.

Reference: M. D. Feit, J. A. Fleck, Jr., and A. Steiger, "Solution of the Schrodinger Equation by a Spectral Method", Journal of Computational Physics 47, 412-433 (1982).

Acknowledgements

The author wishes to acknowledge the following in the creation of this submission:
Eigen Energy Solver for Schroedinger Equation

• MATLAB release: MATLAB 7 (R14)

[quang khai le]

??? Error: File: D:\Programming\TimeDependentSchrodinger.m Line: 136 Column: 18
"identifier" expected, "(" found.

[Narupon Chattrapiban]

I've run the file again on my computer and it's still OK. I'm not sure if it's because of the version of Matlab you used.

[Narupon Chattrapiban]

I had a chance to use this script in matlab6.5 and I can see the problem now.

The function fzero in matlab6.5 and matlab7 are not the same.

For matlab6.5
You can change the script on ~ line 134
from

g = inline('tan(z) - sqrt((z0/z).^2 - 1)','z','z0');

to

g = inline('tan(x) - sqrt((P1/x).^2 - 1)',1);

and change

zn(nrn) = fzero(@(z) g(z,z0),crss_n(nrn));

to

zn(nrn) = fzero(g,crss_n(nrn),[],z0);

[Hyn Ahn]

interesting code!!

[Shams]

Useful ! Thanks for sharing ! Cheers & Good Luck

[William]