matlab油膜压力程序,编的关于雷诺方程求解油膜压力的程序无法运行

程序是基于活塞环-缸套润滑模型的雷诺方程求解，采用有限差分法迭代求解，但程序无法在command window界面运行，求大家帮忙看一下

ita1=0.8*10^-6;

ita2=0.4*10-6;

ita=sqrt((0.8*10^-6)^2+(0.4*10^-6)^2);

Vr1=(ita1/ita)^2;

Vr2=(ita2/ita)^2;

e1=0.00001;

b=0.0015;

n=81;

k=0;

r=0.0575;

l=0.215;

q=890;

u=0.0025;

H=zeros(360,163);

Qs=zeros(360,163);

qs=zeros(360,163);

Qc=zeros(360,163);

Qx=zeros(360,163);

x=zeros(360,163);

for i=1:1:163

k=k+1;

x(i)=0.0015*(i-n-1)/n;

H0=2;

H(i)=H0+e1*x(i)^2/ita/b^2;

if 0<=H(i)<=5

qs(i)=1.899*H(i).^0.98.*exp(-0.92*H(i)+0.05*H(i).^2);

else

qs(i)=1.126.*exp(-0.25*H(i));

end