流体力学中怎么用fortran语言写差分方程啊
发布网友
发布时间:2022-05-02 03:33
我来回答
共2个回答
热心网友
时间:2022-06-27 18:23
下面的语句用的中心差分,把一阶导数、二阶导数都写成了差分格式。
DO I=2,N-1
HX(I)=(H1(I+1)-H1(I-1))/(2*DX)
ENDDO
HX(1)=0
HX(N)=0
DO I=2,N-1
H2X(I)=(H1(I+1)-2*H1(I)+H1(I-1))/(DX*DX)
ENDDO
H2X(1)=0
H2X(N)=0
热心网友
时间:2022-06-27 18:23
下面的语句用的中心差分,把一阶导数、二阶导数都写成了差分格式。
DO I=2,N-1
HX(I)=(H1(I+1)-H1(I-1))/(2*DX)
ENDDO
HX(1)=0
HX(N)=0
DO I=2,N-1
H2X(I)=(H1(I+1)-2*H1(I)+H1(I-1))/(DX*DX)
ENDDO
H2X(1)=0
H2X(N)=0
希望对你有所帮助