1 | AA(I,J) = -B |
2 | DD(I,J) = B+B+A*REL |
3 | PXX = X(I+1,J)-2.*X(I,J)+X(I-1,J) |
4 | QXX = Y(I+1,J)-2.*Y(I,J)+Y(I-1,J) |
5 | PYY = X(I,J+1)-2.*X(I,J)+X(I,J-1) |
6 | QYY = Y(I,J+1)-2.*Y(I,J)+Y(I,J-1) |
7 | PXY = X(I+1,J+1)-X(I+1,J-1)-X(I-1,J+1)+X(I-1,J-1) |
8 | QXY = Y(I+1,J+1)-Y(I+1,J-1)-Y(I-1,J+1)+Y(I-1,J-1) |
9 | RX(I,J) = A*PXX+B*PYY-C*PXY |
10 | RY(I,J) = A*QXX+BCa*QYY-C*QXY |
11 | END DO |
12 | END DO |
13 | Continued from Previous Page |