运行时间历程程序会提示如下错误

dj15175771160

function  fuhuanjeffcottzhuanzixitong_solution
clear all
clc   
n=5000;
p=0;
wj=2*pi*n/60;
uo=0.018;       %外油膜粘度,pa.s
ui=0.018;       %内油膜粘度,pa.s
co=0.000114;    %外油膜间隙,m
ci=0.000093;    %内油膜间隙,m
Rro=0.018830;   %浮环外径,m
Rri=0.014093;   %浮环内径,m
Rj=0.014;       %轴颈半径,m
L=0.018;        %轴承宽度,m
wr=wj/(1+(uo/ui)*(ci/co)*(Rro/Rri)^3);   %浮环的转速
x0=[-0.1 0 -0.1 0 -0.1 0 -0.1 0];     %初值,另一组初值[0 0.1 0 0.1];
options=odeset('RelTol',1.e-7,'AbsTol',[1.e-8 1.e-8 1.e-8 1.e-8 1.e-8 1.e-8 1.e-8 1.e-8] );    %控制精度

figure          % 1时间历程
T=2*pi;         %一个周期
T1=2*pi/wj;
[t,X]=ode23s('fuhuanforce',[0:T1/100:14000*T1],x0,options);      
subplot(121)   
plot(t,X(:,3))
xlabel('t')
ylabel('Yj')       %  grid on
title('时间历程x-t')

subplot(122)   
plot(t,X(:,7))
xlabel('t')
ylabel('Yr')       %  grid on
title('时间历程x-t')