?????????????? ?????? ? ????????? ?????????? ??????? ?????

call sgefa(C,vetvey,vetvey,ipvtC,info)
call sgedi(C,vetvey,vetvey,ipvtC,det,workC,1)
deallocate(ADJ,NEXT,VECTOR,VISIT,LIST)
end subroutine derevo
!******* подпрограмма выделения слова в строке **************
subroutine addword( words, st, wbegin, wend, neww )
integer wbegin, wend, neww
character (len = *) words(*) ! Перенимающий размер массив
character (len = *) st ! Строка, перенимающая длину
neww = neww + 1
words(neww) = st(wbegin : wend)
wbegin = 0
end
!******* подпрограмма инвертирование матрицы ***************
subroutine inverse4(A,n)
real:: A(n,n)
do i=1,n
if (A(i,i)==0.) then
print*,a(i,i),n,i; stop 'деление на ноль!!!'
endif
A(i,i)=1.0/A(i,i)
do 5 j=1,n
if (j-i) 1,5,1
1 A(j,i)=A(j,i)*A(i,i)
do 4 k=1,n
if(k-i) 2,4,2
2 A(j,k)=A(j,k)-A(j,i)*A(i,k)
if (j-n) 4,3,4
3 A(i,k)=- A(i,i)*A(i,k)
4 continue
5 continue
enddo
do j=1,n-1
A(n,j)=- A(n,n)*A(n,j)
end do
return
end subroutine inverse4
92