Unable to plot a wavefunction

4 views (last 30 days)
Hi, I wonder if any of you can try to plot this and see if they also get a weird bar?
syms h g x C t m e p c
h = 1.0545718E-34
g = 5.344285879E-28
E = 1
Z = - (exp(-(x*(g*1i + (- 2*g^2 + E)^(1/2)))/h)*(g*1i - (- 2*g^2 + E)^(1/2)))/(2*(E - 2*g^2)^(1/2)) + (exp(-(x*(g*1i - (- 2*g^2 + E)^(1/2)))/h)*(g*1i + (- 2*g^2 + E)^(1/2)))/(2*(E - 2*g^2)^(1/2))
fplot(real(Z), [ -4, 4 ])
If so, why?

Accepted Answer

David Goodmanson
David Goodmanson on 4 Jan 2018
Hi Sergio,
This function appears to have scaling issues. If you plot
fplot(real(Z), [-2*h, 2*h])
fplot(imag(Z), [-2*h, 2*h])
you can see the function. I believe this addresses the question, but there are other things going on too. Possibly you were using E=1 just for test purposes, but compared to that value, 2*g^2 is infinitesimal. Also the function itself is an increasing exponential which is not good for a wave function unless there is some kind of hard boundary.
  1 Comment
Sergio Manzetti
Sergio Manzetti on 4 Jan 2018
Thanks David, this is excellent! Indeed E=1 was set to test. I will try this out as you suggest! Best wishes

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!