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

ganquanlife

function [nc] = ncutf(g)
6 w; H* J" f5 F  j%求割点的算法 g为邻接矩阵 nc为割点的集合
) R4 _" c# |2 j* x* w7 qn=size(g,1);
% [: i( N% h: ?, M7 J) k* N4 Sif n>=3
; O( z' f2 D0 i4 i) {% e# N) J& Y    a=sum(g);
    b=sum(a==2);
    if b==n
0 q# a1 z& P5 X        fprintf('本图为圈，无割点。n')
5 Q9 R1 ~; }) I8 F: i8 h: u        nc=0;
' Z6 `: ~5 b, A! y    end
else
0 Z' [/ W+ c" `- d! }    [w,k]=dfs3(g);
' c& [* R) F2 `    %nc=[];
    nc=isncf(w,k);
! Y  c8 C& X9 f0 [9 X* y* W5 }) |) I    n=size(g,1);
+ M/ Z/ G/ L# g    for i=1:n
' b; U: \& T; a6 K        for j=1:n
6 T- f; i* {# ~            if w(i,j)>1
                if k(i)>k(j)
. k/ q1 g9 E, ^" E                    g(i,j)=2;
                else
                    g(i,j)=3;
                end
            end
( z! D' `  i, T2 a        end
    end
   
    for i=1:n
        f1=find(g(i,==2);
8 M# A1 j$ B$ m+ l* C        f2=find(g(i,==3);
. ~: Y& c$ S+ M$ ~0 B3 _; i        f=union(f1,f2);
! r( d$ N/ }, n" O. q6 n+ [7 z; t        l(i)=min([k(f) k(i)]);
    end
   
    for i=1:n
        for j=1:n
            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
5 @! N( u- `& g' x            end
7 }( z8 \' A$ i. w( c        end
# _! h& n0 Y  y3 ]% K    end
end
* W8 P  u$ l0 N5 }8 Bend

3 ~* P' Z8 W$ x. V0 i+ b: f
% k/ t8 B) h8 R0 x# t! [function nc=isncf(w,k)
  i7 l0 B, C+ |: a  T0 `    nc=[];
8 \2 o7 L' w7 N    t=zeros(size(w));
+ ?4 N/ K0 O. m4 N    n=size(w,1);
1 q& e5 e6 x7 {/ M    a=find(w~=0);
    for i=1:length(a)
        d(i)=w(a(i));
        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
        else
            t(i)=floor(a(i)/n);
        end
) i4 k2 j: R9 f( e- _( }6 d        t1(i)=mod(a(i),n);
% _3 d$ p5 d$ n        if t1(i)==0
            t1(i)=n;
; z9 F* `# w: t0 T& z        end
7 k( U4 s, r4 x# y  o0 k: k  m/ h    end
4 E0 y& t7 W- _  ]4 L  q    [b,c]=sort(d);
    p=[1];pc=0;
- Y/ e' {( \& C6 O    for i=1:length(a)
; E1 \& \- V- J- v* h) B        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
( `; `2 ~: L5 S! l            t(t1(c(i)),t(c(i)))=3;
  v0 G* T+ M6 z+ M8 a6 r! X2 w        end
4 ^& a& @* i2 k+ v; j% s        if pc==0
            tc=isempty(setdiff([1:n],p));
9 F8 ~* D* }. H: ?! n            if tc
  n' E7 q! [8 U0 a' V2 A. @                t0=sum(t(1,==3);
                if  t0>=2
                    nc=union(nc,1);
                end
                break;
            end
. H5 `& B9 ~- V        end
+ q/ N, z) _0 Z- _2 i6 c2 O) Z    end
4 q6 @; ~. p, D8 M   
        
end

ganquanlife

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