matlab雷诺曲线,关于用有限差分法解雷诺方程

这是我编写的程序，希望对你有用处！

function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)

Fx=0;

Fy=0;

epsilon=0.6;

L=0.03;

d=0.06;

lambda=L/d;

m=40;

n=30;

phi1=0;

phi2=2*pi;

delta_phi=(phi2-phi1)/m;

delta_lambda=2/n;

k=1;

P0=zeros(m+1,n+1);

A=zeros(m+1,n+1);

B=zeros(m+1,n+1);

C=zeros(m+1,n+1);

D=zeros(m+1,n+1);

E=zeros(m+1,n+1);

F=zeros(m+1,n+1);

H=zeros(m+1,n+1);

for i=1:1:m+1

%      theta(i)=(i-1)*delta_phi;

for j=1:1:n+1

H(i,j)=1+epsilon*cos((i-1)*delta_phi);

end

end

S=0;

T=0;

ERR=1e-3;

G