求教四个方程的常微分方程组求解问题

yxlnbu

大家好，我有个四个方程的常微分方程组，定义如下，该方程可以求解，可是解基本上NaN。我想请大家帮我看看，是什么问题。我想得到的是时间从e-6—e-4这一段的数值
谢谢您了
function xdot=dianci(t,x)
A=0.01;
c=3e10;
ro=8.924;
Y0=69e10;
K=36;
rw=0.05642;
Rs=0.004;
Rc=4.73e14;
Lc=1.542e-18;
U=13.34256384;
br=-( -4.6861.*x(1).^5 + 72.6026.*x(1).^4 - 447.6495.*x(1).^3 + 1376.6343.*x(1).^2-2121.6799.*x(1)+1323.3641);
mr=(396.51-413.32.*x(1)+192.46.*x(1).^2-46.45.*x(1).^3+5.69.*x(1).^4-0.28.*x(1).^5);
L=(4*pi/c^2)*(log(8*x(1)/rw)-1.75);
dl=(4*pi/c^2)*(log(8*x(1)/rw)-0.75);%dl/dr
M=(mr.*2.*pi./c.^2);
dm=((2.*pi./c.^2).*br.*x(1));%dm/dr
xdot=[ x(2);%r
      br.*x(3).*x(4)./(ro.*A.*c^2)+(x(4)^2./(4.*pi.*A.*x(1)))*dl-(Y0.*(1+K.*x(2)./x(1)).^0.45)./(ro.*x(1));%vr
      (1/(M^2-L*Lc))*(L*(Rc.*x(3)+x(4).*x(2).*dm+13.33)-M*(Rs.*x(4)+x(3).*x(2).*dm+x(4)*x(2).*dl));%i1
      (1/(M^2-L*Lc))*(Lc*(Rs.*x(4)+x(3).*x(2).*dm+x(4).*x(2).*dl)-M.*(Rc.*x(3)+x(4).*x(2).*dm+U));%i2
      ];%函数方程

[t y]=ode45(@dianci,[0,1],[1.6;0;0;0])%解法及初值