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

function [C,A,g]=formc(ta,he)
maxv1=max (ta);maxv2=max (he);maxv= max (maxv1,maxv2);
g = make _graph ('cep1',1,maxv,ta,he);
g1=make _graph ('cep2',0,maxv,ta,he);
a=graph _2_mat (g,'node-arc');
A0=full (a); arc = arc_number (g); node = node_number (g);
nc=connex (g);
if nc==1 then
t= min _weight _tree (node,g1);
yzl=maxv; hordsum=setdiff (1:arc,t)
at(1:node-1, 1:arc)=A0(1:node-1, 1:arc);
size_hord=size(hordsum);
for i=1:size_hord(2) at(node-1+i,hordsum(i))=1; end;
else
inow=0;kd=0;
for ks=1:nc
ns=con_nodes(ks,g); size_ns=size(ns); kyz=size_ns(2);
for j=1:kyz-1 yz0(kd+j)=ns(j); end; kd=kd+kyz-1;
k=0;vetv= zeros(arc);
for i=1:kyz
i1=ns(i);
for j=1:arc
jot=j;
if A0(i1,j)==- 1. then
k=k+1; vetv(k)=j;
end;
end;
end;
gpod1=subgraph (ns,'nodes',g1); node1 = node_number (gpod1);
t= min _weight _tree (node1,gpod1); size_t= size(t);
for ii=1:size_t(2)
tnew(ii)=vetv(t(ii));
end;
hord=setdiff (vetv',tnew');
for j=1:k-ii hordsum(inow+j)=hord(j); end; inow=inow+k-ii;
end;
hordsum=hordsum'; yzl=yz0'; size_yz0=size(yz0');
for j=1:size_yz0(2)
for i=1:node
if i==yz0(j) then
at(j,1:arc)=A0(i,1:arc);
end;
end;
end;
size_hord=size(hordsum); nhord=size_hord(2);
for i=1:nhord at(size_yz0(2)+i,hordsum(i))=1; end;
end;
A=at'; C=at^(- 1);
endfunction ;
82