AA(I,J) = -B |
DD(I,J) = B+B+A*REL |
PXX = X(I+1,J)-2.*X(I,J)+X(I-1,J) |
QXX = Y(I+1,J)-2.*Y(I,J)+Y(I-1,J) |
PYY = X(I,J+1)-2.*X(I,J)+X(I,J-1) |
QYY = Y(I,J+1)-2.*Y(I,J)+Y(I,J-1) |
PXY = X(I+1,J+1)-X(I+1,J-1)-X(I-1,J+1)+X(I-1,J-1) |
QXY = Y(I+1,J+1)-Y(I+1,J-1)-Y(I-1,J+1)+Y(I-1,J-1) |
RX(I,J) = A*PXX+B*PYY-C*PXY |
RY(I,J) = A*QXX+BCa*QYY-C*QXY |
END DO |
END DO |
Continued from Previous Page |