matlab ode45求解齿轮动力学,MATLAB ode45求解轴承动力学微分方程组，显示积分公差无法满足...

function fun=ballbearings6(t,y)

global Mri  delta_r theta_2 alpha_2 Cbl alpha_1 Cbr delta_l theta_1 Chol  Chor g  Fx Fy Fz B1 B2 Mhol  Mhor  Khxl  Khyl  Khzl  Khxr  Khyr  Khzr  Ix  Iy  Iz  l1  l2  l  D  f01  v01  Ns   Nb  wc  d    Ke  alpha01  A   deltaep Fcr Fcl  Flx  Fly  Flz  FILx  FILy1  FILy2 Frx Fry Frz FIRx FIRy1 FIRy2

Mri = 3.5;

Cbl = 2.9911*10^3;   %轴承的阻尼

Cbr = 2.9911*10^3;

Chol = 2.9911*10^3;

Chor = 2.9911*10^3;

g = 10;

Fx= 100*sin(wc*t) ;  %径向力

Fy= 0;

Fz= 0;    %轴向力

Mhol = 4.2;   %外圈和轴承座的质量总和

Mhor = 4.2;

Khxl = 2.9911*10^7;   %轴承座刚度

Khyl = 2.9911*10^7;

Khzl = 2.9911*10^7;

Khxr = 2.9911*10^7;

Khyr = 2.9911*10^7;

Khzr = 2.9911*10^7;

Ix = 0.5177;    % 转动惯量

Iy = 0.5177;

Iz = 0