求助关于求图割点的代码，哪里错了？怎么改正啊？

ganquanlife

function [nc] = ncutf(g)
8 h# q3 h  w- _# C* j! M%求割点的算法 g为邻接矩阵 nc为割点的集合
  M" y* |5 N4 Z" h+ z0 M# Ln=size(g,1);
3 I( ]& J4 w4 w  nif n>=3
    a=sum(g);
5 l$ p  T$ i9 z. O) n' X    b=sum(a==2);
4 m; `: K8 ^2 `$ V  j6 ?    if b==n
0 o4 l/ j8 P/ z" x6 _3 {        fprintf('本图为圈，无割点。n')
% Z4 _& }. x9 \& m* h- ?        nc=0;
) _! f' s6 r' k2 n8 T+ Y    end
else
6 j2 A% v, {" U0 f/ T    [w,k]=dfs3(g);
    %nc=[];
6 F3 Y' g: k! @8 z    nc=isncf(w,k);
    n=size(g,1);
0 V$ S9 F% w* }" e- ^$ h2 m    for i=1:n
        for j=1:n
+ Q  f& F7 Z2 W& ~" N. P            if w(i,j)>1
                if k(i)>k(j)
" E: C# N: d0 `! {9 |& p                    g(i,j)=2;
+ @( W. |' m6 K; P* f                else
                    g(i,j)=3;
3 |7 h0 z- b! a% N/ P6 N                end
% T6 t, V9 ^1 I& ?            end
        end
    end
0 R& R; T" t. }- {$ ~  G$ k   
; o$ o, k* y* a* ~    for i=1:n
8 m% L( M! K/ g0 s        f1=find(g(i,==2);
  Y& C6 ^" y( d; u, U$ y& }, o        f2=find(g(i,==3);
        f=union(f1,f2);
        l(i)=min([k(f) k(i)]);
    end
   
0 H2 p5 C( P/ R4 p8 N$ |* x2 h    for i=1:n
% q/ E- V+ f% }1 J. \& p        for j=1:n
            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
2 ]: I0 K$ _5 Y6 x9 P                nc=union(i,nc);
            end
& g' U5 `' x' T& c9 f        end
/ h7 z' Q& c4 Y4 u    end
end
end

2 I$ X" D* I& H) Z; J  F
function nc=isncf(w,k)
    nc=[];
: N$ p" B2 z0 m    t=zeros(size(w));
    n=size(w,1);
    a=find(w~=0);
3 v# ~7 E; x5 B# h    for i=1:length(a)
- Z+ ^- T5 Y" ?' n        d(i)=w(a(i));
: r( O6 F  C5 I8 v        if a(i)/n>floor(a(i)/n)
- L' Q6 |. `  Z: E/ \3 b          t(i)=floor(a(i)/n)+1;
' K; s4 A1 ~2 T. ~        else
            t(i)=floor(a(i)/n);
        end
        t1(i)=mod(a(i),n);
        if t1(i)==0
, {; g6 E& s7 d* X6 v            t1(i)=n;
        end
% X! L- q1 x5 d4 \6 k$ i    end
    [b,c]=sort(d);
) O" a$ \* U4 Y3 a    p=[1];pc=0;
) y, z0 I& ?0 S/ d    for i=1:length(a)
4 Q* F& i- M2 Q/ y6 t( v* G        if k(t1(c(i)))<k(t(c(i)))
2 ~' y8 T) X, t' p; {! c0 h/ z            p=union(p,t(c(i)));
            t(t1(c(i)),t(c(i)))=3;
7 {$ n  V) ?/ L% u8 E7 \        end
( G$ C& x' o" n# M        if pc==0
            tc=isempty(setdiff([1:n],p));
* X% x( w0 W  O            if tc
, J, |1 c1 Q' L7 l- e1 J7 f! y                t0=sum(t(1,==3);
                if  t0>=2
                    nc=union(nc,1);
6 y' D7 H% t; g7 H5 c6 \! R                end
$ N+ Y5 ?( ]8 X. Q: m0 b' I                break;
            end
        end
    end
% B3 N9 J% B4 \3 _! Y" h   
        
6 D0 u8 ], h9 Eend

ganquanlife

如果没法改正，另外编一个求割点的代码也行阿