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

ganquanlife

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
+ |$ k: R& f. q; T( Z) uif n>=3
    a=sum(g);
2 F& t( c) N' B# |1 N2 a* X    b=sum(a==2);
    if b==n
        fprintf('本图为圈，无割点。n')
        nc=0;
# ?; @- C9 y, d, e2 {( K! j: g/ H2 s    end
else
    [w,k]=dfs3(g);
9 z5 T# b9 M9 K) R    %nc=[];
6 b& \$ \6 {) q5 p# c    nc=isncf(w,k);
) B6 F" X/ M# U% G( B% {4 @    n=size(g,1);
6 h( t* o6 p# C$ Q  `- g( `, i    for i=1:n
        for j=1:n
5 ]! D1 j* N/ j: B0 f/ t1 s6 d+ K- q- C            if w(i,j)>1
+ V/ P, I% [: n3 M0 G8 y8 z                if k(i)>k(j)
$ c7 W$ g6 ^( z6 E7 l  t  t: m" H                    g(i,j)=2;
                else
                    g(i,j)=3;
& W/ z. j: F3 o: m, J! D                end
            end
        end
+ S( ~! ^' \7 Y0 g7 L# _$ ?) G8 m    end
   
& j9 z6 q4 E) M    for i=1:n
% s; \. }/ o) }( }" R; e        f1=find(g(i,==2);
6 t! ]% E6 |# ^& ^9 L" V        f2=find(g(i,==3);
# z; J2 f2 ?, T* p# t9 B        f=union(f1,f2);
        l(i)=min([k(f) k(i)]);
    end
5 ~+ H  o7 {, ~   
    for i=1:n
        for j=1:n
! M* G+ S. Q1 R) n3 R* F            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
0 a- H# K- {/ X, q# |            end
        end
    end
- G+ y* R1 t; p. N" cend
end

7 Z3 h' x' d  Q, x1 g7 ^1 V
: n4 S1 i. O) h: H5 K, Ffunction nc=isncf(w,k)
    nc=[];
1 M4 }( g& }! T  W# A1 P! X, ?: z    t=zeros(size(w));
    n=size(w,1);
    a=find(w~=0);
    for i=1:length(a)
        d(i)=w(a(i));
8 c: c  M4 e5 [9 D        if a(i)/n>floor(a(i)/n)
+ `  S6 V- J9 ?% }9 F* u          t(i)=floor(a(i)/n)+1;
        else
            t(i)=floor(a(i)/n);
! }$ M7 h+ v' g" r% g% g        end
2 d5 d- ~/ M7 B8 J& g: x        t1(i)=mod(a(i),n);
  m2 [* }; G4 a* I) [( ]) j2 j7 u        if t1(i)==0
+ @2 a% B" z  s: s5 i$ @            t1(i)=n;
        end
    end
    [b,c]=sort(d);
    p=[1];pc=0;
( L3 j* S. V- r    for i=1:length(a)
. X( I" D8 ]) q0 c, M        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
0 _/ s9 P+ w4 `! b! l            t(t1(c(i)),t(c(i)))=3;
3 d% H/ x  K1 t3 `" f; v: i/ c        end
' m4 K9 @, X3 S# [: j; H        if pc==0
            tc=isempty(setdiff([1:n],p));
( L) e* d5 S8 k' D# x1 }6 m( `1 n            if tc
                t0=sum(t(1,==3);
4 q$ ^/ r' C9 H# v% ~                if  t0>=2
8 N6 ]; P6 z$ j/ Y& h' O                    nc=union(nc,1);
                end
; h! o2 H: D* Q4 B3 M9 g, o                break;
            end
  v" ~  ?/ z; c: P        end
9 d3 Z6 M( t: W9 s; H  L    end
' i& S& h7 U9 R, ~5 t! o   
        
end

ganquanlife

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