用ode45解微分方程遇到的实际问题

        最近在用ode45解微分方程数值解，试图复现论文中的图。一般来说说微分方程（组）只要按照响应的条件去撰写好对应的回调函数即可，基本没什么难度，但对于本文遇到的的这个问题，可能还需要一些技巧去实现解法，这篇文章就来说说我们其中遇到的几个问题。

       一、问题提出和简单分析：

       方程的条件和初值如下：

        常系数和初始变量：        

  模型参数：Ms=1.6x10^6；a=1000；alpha=0.001；,mu0=4πx10^7 ；c=0.1；gamma1=7x10^-18；gamma1_p=-1x10^-25 ；gamma2 =-3.3x10^-30 ；gamma2_p=2.1x10^-38 ；E=2.02x10^11 ；ksi=2000 ；H=40 ；

        自变量σ区间为[0,300]MPa 注意：1MPa=1×106 Pa