雷诺方程matlab算法

function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)Fx=0;Fy=0;epsilon=0.6;%偏心率=0.6L=0.03;d=0.06;%轴颈直径dlamb