! Calculeaza x din: S*x =y
! Prin substitutie inapoi
! Matricea L - data de Cholesky; Vectorul 'y'->'b' dat de Forward
! 'x' stocat in 'b'
!
subroutine back(n, s,b)
!common n
real(4) s(n*(n+1)/2), b(n)
real(4) sik, sii

! i=n:
b(n) =b(n)/s(loca(n,n))

A_DO: do i =n-1,1,-1
xi =0.
	do k =i+1, n
	xk =b(k)
	sik =s(loca(i,k))
	xi =xi +sik*xk
	enddo
yi =b(i)

sii =s(loca(i,i))
xi =(yi -xi)/sii
b(i) =xi
ENDDO A_DO

END subroutine back

