options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); 是解微分方程时的选项设置，'RelTol',1e-4,是相对误差设置，'AbsTol',[1e-4 1e-4 1e-5]是绝对误差设置
