matlab ode45三体问题,小白急求~~关于ode45不能解的问题

ode45本身是变步长的呀，怎么设置最小步长？我改用其它ode函数计算都会出现这个警告Warning: Failure at t=5.373174e-001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.908935e-015) at time t.，用ode4函数算了也出现这个问题，缩小tspan也不行，现贴上我的程序请高手指点。

clear

clc

miu=0.0178;

L=201.6e-3;

R=210e-3;

LD=L/(2*R);

C=0.699e-3;

pofai=C/R;

mg=104860;

m=mg/9.8;

ome=450;

pa=2*miu*ome/(pofai.^2);

Kxx =    6.453356740612708;

Kxy =   -4.133413328965316;

Kyx =1.146014052781378; Kyy = 7.846056548809618;

Dxx =1.145752328144106;Dxy =1.200350375824454;

Dyx =1.200350375824454;Dyy = 2.277131377326176;

KXxx = 4.310137875942596;KXxy =1.621686159880267;KXyy=-32.99946362265302;

DXxx =