quantum states does not "feel" the dirac delta function

[giladsofer]

While using quantumstates along with the potential V=dirac(x), it seems that the energies eigenstates computed by quantumstates are exactly the same as that of V=0. This remains the same for several delta functions, located at different positions. It seems that the delta function potential is 'invisible' to quantumstates. Any ideas on how to fix this? (assuming I want to use the dirac(x) function MATLAB offers, and not write an ad hoc delta potential myself)

[nickhale]

Unfortunately there is currently no support for solving ODEs (or ODE eigenvalue problems) with delta functions. If I recall; we started looking at it, but it would be quite a lot of effort to implement robustly. We should at least throw a warning/error if someone attempts this, however.

[trefethen]

As @nickhale says, one can't just write it down and solve it as one would like. On the other hand, one could fake the delta function by making it a boundary with a jump condition on the derivative. This should work fine though without maximal Chebfun elegance.

[giladsofer]

As @nickhale says, one can't just write it down and solve it as one would like. On the other hand, one could fake the delta function by making it a boundary with a jump condition on the derivative. This should work fine though without maximal Chebfun elegance.

Thank you - this could definitely solve my problem. How can I implement the jump boundary condition on the derivative using quantumstates? (if I cannot, is there something else I can use?)