% This function is used to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m; %沿圆周方向均匀划分m格
delta_n=2/n; %沿轴向均匀划分n格
for j=1:n+1 %对各节点赋压力初值为0
for i=1:m+1
P(i,j)=0;
Q(i,j)=0;
end
end
S=0;
T=0;
for j=2:n %第一次计算各内部节点的压力值
for i=2:m
A=(1+e*cos((i+1/2-1)*delta_m))^3;
B=(1+e*cos((i-1/2-1)*delta_m))^3;
C=((1/ratio)*delta_m/delta_n)^2*(1+e*cos((i-1)*delta_m))^3;
D=C;
E=A+B+C+D;
F=delta_m*e*(cos((i+1/2-1)*delta_m)-cos((i-1/2-1)*delta_m));
P(i,j)=(A*P(i+1,j)+B*P(i-1,j)+