ajarifi
commented
3 months ago \xi may need to be defined earlier in the potential.
@doc raw"""
HarmonicOscillator(k=1.0, m=1.0, ℏ=1.0)
k is the force constant, m is the mass of particle and \hbar is the reduced Planck constant (Dirac's constant).
""" HarmonicOscillator
@doc raw"""
V(model::HarmonicOscillator, x)
V(x)
= \frac{1}{2} k x^2
= \frac{1}{2} m \omega^2 x^2
= \frac{1}{2} \hbar \omega \xi^2,where \omega = \sqrt{k/m} is the angular frequency and \xi = \sqrt{\frac{m\omega}{\hbar}}x.
""" V(model::HarmonicOscillator, x)
ohno
HarmonicOscillator(α=1.0, m=1.0, ℏ=1.0). α should be replaced by k. Like this: HarmonicOscillator(k=1.0, m=1.0, ℏ=1.0).