Вопрос для умных! Реализация метода координат на фортране

Fortran 90
Сижу и пытаюсь написать прогу на фортране. Ну вот никак не могу въехать в методы гиперплоскостей. Как там распараллелить?
Задание вот такое:
DO 30 I=1,N
DO 30 J=I,N
DO 30 K=J,N
A(I,J)=B(J,K)+C(I,J)

Решить нужно методом координат.
Кто-нибудь может чем-нибудь помочь. Может просто подсказать про этот метод. Застрял я окончательно, а сдавать уже... уже сегодня надо.
Если кто шаритт в этом, помогите хоть чем-то, пожалуйста.

Кое-что написал. Вот, посмотрите пожалуйста:
DIMENSION B(100,100),C(100,100),A(100,100),AA(100,100)
DIMENSION TMP (100)
S=50
N=10
DO 10 I=1,N
DO 10 J=I,N
DO 10 K=J,N
A(I,J)=I
B(J,K)=J*2-I
C(I,J)=K-7+I*J
10 CONTINUE
CSBEG
DO 30 I=1,N
DO 30 J=I,N
DO 30 K=J,N
A(I,J)=B(J,K)+C(I,J)
30 CONTINUE
CSEND
CPBEG
DO 20 I=1,N
#CALL EADDVV(TMP,B(1,I),C(I,1),N)
AA(1,I)=TMP (I)
20 CONTINUE
CPEND
IFLAG=0
DO 40 I=1,N
DO 40 J=1,N
IF (A(I,J).NE.AA(I,J)) IFLAG=1
40 CONTINUE
50 FORMAT (' IFLAG= ',I1)
WRITE(6,50) IFLAG
STOP
END