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

ganquanlife

function [nc] = ncutf(g)
+ j. `. E& ?% G% v( }%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
if n>=3
1 W) m/ D$ `# U& d& C# _    a=sum(g);
! r: R' K- b  V+ Y. b' U    b=sum(a==2);
    if b==n
        fprintf('本图为圈，无割点。n')
        nc=0;
9 R% C) Y* T( f7 E" }6 v) S    end
else
: k& n! S) m* [) `, u# T9 |    [w,k]=dfs3(g);
    %nc=[];
    nc=isncf(w,k);
- T) P, @, g0 C: |. ]( c! R) o    n=size(g,1);
+ W! _9 m) ?( D8 y* {    for i=1:n
        for j=1:n
* r2 g# I( ~) @2 S- `            if w(i,j)>1
                if k(i)>k(j)
                    g(i,j)=2;
                else
                    g(i,j)=3;
3 W7 w5 L: |5 @, ~" j2 l5 v                end
            end
        end
    end
   
    for i=1:n
. O! J) `! I- x: w, s- O( ?1 |: i" p        f1=find(g(i,==2);
  H4 h* Y! y" r5 {8 d1 {* K        f2=find(g(i,==3);
        f=union(f1,f2);
# c* N, ]% X! e) }! G        l(i)=min([k(f) k(i)]);
$ s7 ]2 p+ {2 M: c7 M. f    end
+ k) }3 A2 B, F% b# o9 N   
    for i=1:n
        for j=1:n
            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
2 F+ p- l9 ]% y/ T* Y) a                nc=union(i,nc);
( K+ U! E( o+ e" q9 I8 v' d            end
( H# O9 l/ F3 l' y% i! F! g        end
: H$ D+ z) ^+ G6 j2 X8 d$ c    end
! u/ b  _8 X9 X: f; L, _end
0 F7 U; ^; b- v, E& a3 h' Cend
, [& B" b- U( H: z' H- V9 o

' O; f1 |+ V, n4 V6 wfunction nc=isncf(w,k)
" a- b8 i: W. ^9 Y    nc=[];
- m  L) Z* `5 Q3 Z8 }9 m/ r    t=zeros(size(w));
! `! `! W) S0 f- }& m4 E" H    n=size(w,1);
+ D% ?0 @9 ~3 @/ R! e    a=find(w~=0);
    for i=1:length(a)
, }. x7 ~$ K7 e( g, q! I8 ?2 G  {        d(i)=w(a(i));
6 o/ d$ s, H- P+ ]3 f! w5 D        if a(i)/n>floor(a(i)/n)
; m2 g# v2 X. t          t(i)=floor(a(i)/n)+1;
        else
. w( X# ^* b. F5 S# W            t(i)=floor(a(i)/n);
/ q% i- h6 Y& F        end
+ ~  r/ ]- w' B$ x; i        t1(i)=mod(a(i),n);
2 z: h6 g0 S2 `  x# P0 K. b4 [5 |        if t1(i)==0
            t1(i)=n;
        end
9 |! q) K6 {& i3 N4 X3 J2 r5 ^    end
    [b,c]=sort(d);
0 I% m  ^6 b$ D$ Z; c    p=[1];pc=0;
    for i=1:length(a)
5 e/ }! g6 c* u* b; s$ v        if k(t1(c(i)))<k(t(c(i)))
& M1 i7 c  a2 d            p=union(p,t(c(i)));
# Y. [% Y' h0 n# T8 Y9 W; M* k            t(t1(c(i)),t(c(i)))=3;
        end
5 v1 V+ }1 Z6 R5 m! {. o) s        if pc==0
8 L% O7 l- F3 B' Y( S  L- L            tc=isempty(setdiff([1:n],p));
1 F9 j) [2 ^3 Y5 \& [            if tc
                t0=sum(t(1,==3);
                if  t0>=2
  X. R7 o- S0 o: Y3 ~$ E                    nc=union(nc,1);
                end
                break;
            end
/ k' V; {% D6 d) [8 [        end
- M# P( a% a8 V$ a9 {* P    end
   
        
end

ganquanlife

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