* * $Id: sfeqn.inc,v 1.1.1.1 1996/02/15 17:49:04 mclareni Exp $ * * $Log: sfeqn.inc,v $ * Revision 1.1.1.1 1996/02/15 17:49:04 mclareni * Kernlib * * * * sfeqn.inc * IF(K .LE. 0) GOTO 299 DO 220 L = 1, K B(1,L) = A(1,1)*B(1,L) 220 CONTINUE IF(N .EQ. 1) GOTO 299 DO 243 L = 1, K DO 232 I = 2, N IM1 = I-1 S21 = - B(I,L) DO 231 J = 1, IM1 S21 = DOTF(A(I,J),B(J,L),S21) 231 CONTINUE B(I,L) = - A(I,I)*S21 232 CONTINUE NM1 = N-1 DO 242 I = 1, NM1 NMI = N-I S22 = - B(NMI,L) DO 241 J = 1, I NMJP1 = N - J+1 S22 = DOTF(A(NMI,NMJP1),B(NMJP1,L),S22) 241 CONTINUE B(NMI,L) = - S22 242 CONTINUE 243 CONTINUE 299 CONTINUE