Added files

Run Mag_bias.m to check. Now, it does not converge correctly. I check
the discrete system dynamic computation and residual computation, no big
problem found. There is error of system dynamic computation, but not
big. I will update later.