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

ganquanlife

function [nc] = ncutf(g)
7 H2 M: [( w; D1 E%求割点的算法 g为邻接矩阵 nc为割点的集合
+ D! D2 Q" D' k  }. s. Jn=size(g,1);
" f: |7 @9 E$ N0 bif n>=3
( ^4 K6 n" K: W; _7 G+ a/ u& F    a=sum(g);
    b=sum(a==2);
    if b==n
        fprintf('本图为圈，无割点。n')
        nc=0;
" O/ v& \0 u7 e2 ?6 M' w' ^    end
else
3 s/ ?  y3 W/ s: Z$ w    [w,k]=dfs3(g);
    %nc=[];
- t. y% @6 R( F6 Q- z4 Q! d" o; o    nc=isncf(w,k);
    n=size(g,1);
    for i=1:n
0 j4 z& F7 }$ t6 `' n+ |  c        for j=1:n
            if w(i,j)>1
7 ~0 m9 N! J+ q/ {                if k(i)>k(j)
6 @7 |5 i4 _# M! Q; C7 m5 |  k! {                    g(i,j)=2;
; q, T7 B2 L5 v0 U4 s                else
7 s  `! G" w# `0 a% `" f                    g(i,j)=3;
! e6 R# q" l% H* x8 z* t3 T                end
            end
        end
2 X- m: T/ h1 w/ ^1 ^    end
/ |, N3 {9 @$ P2 ~4 H  u   
- S9 s! K+ x% Z    for i=1:n
9 c2 F- ?6 }4 b) p  a; \        f1=find(g(i,==2);
: W* D6 m& x: \, \0 p. _        f2=find(g(i,==3);
/ C3 k' M7 D* G3 P4 n; e        f=union(f1,f2);
$ E; ^5 d) T7 @1 k- A1 w$ A' W, u3 y        l(i)=min([k(f) k(i)]);
2 D1 Y& G6 o+ R  u  H    end
0 }' L" U$ ]) \, E1 t; d   
    for i=1:n
: G4 x) ]7 K# Q        for j=1:n
6 V2 C* p, |9 w; P) i. V* Q            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
  o( L2 J1 u: A8 G                nc=union(i,nc);
            end
        end
, G0 o& I3 W7 O! `, j    end
0 e  O, [- L( D: Mend
end


function nc=isncf(w,k)
: e! z, y5 n1 ~6 u5 ~    nc=[];
8 _4 M; x- A( z3 J+ a    t=zeros(size(w));
    n=size(w,1);
) s1 V0 I$ v: O. ~" @6 ~+ N( y    a=find(w~=0);
    for i=1:length(a)
        d(i)=w(a(i));
4 A; l5 k4 C5 B  O. S) I        if a(i)/n>floor(a(i)/n)
! O( s, o4 b: N& N4 ^! f* V          t(i)=floor(a(i)/n)+1;
+ C  M3 v; t% k2 J        else
            t(i)=floor(a(i)/n);
) R& \6 o; F1 `9 ]% c        end
: g3 \, X7 n8 \3 S8 k; _        t1(i)=mod(a(i),n);
        if t1(i)==0
" t* S1 O8 z$ h. U8 o4 s            t1(i)=n;
        end
    end
' ?0 y3 Q* ~' A; g% G6 m  U" P, b  Z    [b,c]=sort(d);
# W, h4 C: U; T& a6 [  h# j4 F    p=[1];pc=0;
- y9 E0 J# f) O, \7 H8 I    for i=1:length(a)
        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
+ V/ _& ]5 O; l2 Z+ [$ S& j( o            t(t1(c(i)),t(c(i)))=3;
1 f6 k9 \! M9 o! A0 V0 J, {        end
6 K1 k. n3 l) Z, l4 S# J        if pc==0
$ J% z& q8 R9 M- d, [- X1 j9 q            tc=isempty(setdiff([1:n],p));
7 X) h4 g( q, H: r" v5 B3 i            if tc
6 U9 Y' d9 t! v                t0=sum(t(1,==3);
' B1 T* H) P6 ]0 o! z( R                if  t0>=2
5 u9 Y* I# a5 `$ p9 `                    nc=union(nc,1);
2 ^2 Q3 |, m  J, G2 T                end
: _: R# H( d0 V7 _                break;
            end
( ?3 e$ X8 ^, j3 d/ N! S        end
    end
/ m( _0 F: I) S% J   
        
end

ganquanlife

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