Unable to plot a wavefunction

2 次查看(过去 30 天)
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?


David Goodmanson
David Goodmanson 2018-1-4
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 个评论
Sergio Manzetti
Sergio Manzetti 2018-1-4
Thanks David, this is excellent! Indeed E=1 was set to test. I will try this out as you suggest! Best wishes


更多回答(0 个)


Help CenterFile Exchange 中查找有关 Annotations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by