I have a time iteration function looked on a 2D surface like this.
Since the numbers wee very small i.e. hbar=6.6260700404e-34./(2*pi), my professor told me to use our own "scaled unites" during the calculation. ie. instead of using $1e-9$ (meter), scale the unit and use $1e-6$ ($\mu m$) during the calculation. e.t.c.
It worked very well when I was only changing the space parameter, (m) However, when I adjust the space scale and time scale together. I encounter a problem, where the function stopped being evolving smoothly.
I checked my substitution, and they were correct. But what happened here? Can we really save time of calculation by using scaled units?