!........................................................................
!
! Subroutine :
!
! Description:
!
! Arguments  :          
!
! Mod/Commons:    
!
! Calls      :      
!
! Author     : 
!
! Modified   :     
!             
!........................................................................
!
! $Id$
!
! $Log$
! Revision 1.2  2003/04/30 17:14:51  cesrulib
! dlr's changes since last import
!
! Revision 1.1.1.1  2002/12/13 19:23:28  cesrulib
! import bmadz
!
!
!........................................................................
!
  SUBROUTINE GJDET_sc(A,B,X,N)

  use precision_def

  implicit none

  integer i,n,k,j,m

  real(rp) A(n,n),B(n),X(n)
  real(rp) const

  DO 1 I=1,N
  DO 2 K=1,N
  IF(K==I)GO TO 2
100     FORMAT(1X,2I4,E12.5,5X,E12.5)
  CONST=-A(I,K)/A(I,I)
  DO 3 J=1,N
  A(J,K)=A(J,K)+CONST*A(J,I)
  IF(J==I)A(J,K)=0.
3  CONTINUE
  B(K)=B(K)+CONST*B(I)
2  CONTINUE
  CONST=A(I,I)

  A(1:n,I)=A(1:n,I)/CONST
  A(I,I)=1.
  B(I)=B(I)/CONST
1  CONTINUE
  DO 5 M=1,N
5  X(M)=B(M)
  RETURN
  END