Finally, we finish all fundamental code for simulating few body real-time dynamics code, with a internal plot and data processing function only works for closed system (non-superoperator form).
For simulating open system with Lindbladian operators, please take mesolve result out and do analysis using self-made functions.
Important update!!
we have good data saving and loading class, named by `qsave`, will make your life easier when doing simulating data saving, loading, and undoing.