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

ganquanlife

function [nc] = ncutf(g)
9 l7 |- @3 \! S- u/ |8 v& k%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
if n>=3
4 G1 T2 Q. L' i& b" G2 A' B    a=sum(g);
    b=sum(a==2);
# P" n) ]- u" k/ S& h9 {9 c4 _    if b==n
        fprintf('本图为圈，无割点。n')
8 n% g1 |7 Z1 X% w) k        nc=0;
    end
else
    [w,k]=dfs3(g);
    %nc=[];
    nc=isncf(w,k);
    n=size(g,1);
    for i=1:n
        for j=1:n
            if w(i,j)>1
! o/ g2 x- K2 V$ B                if k(i)>k(j)
                    g(i,j)=2;
                else
                    g(i,j)=3;
0 J2 b5 q; T1 ~! y- Y9 z9 n$ T                end
& ~- B$ }( I: r$ q" [            end
! F* e! j7 ~0 K1 q9 ?% f4 t        end
    end
   
. ?) W' D0 f/ I. u/ w; }2 o8 e    for i=1:n
        f1=find(g(i,==2);
        f2=find(g(i,==3);
        f=union(f1,f2);
        l(i)=min([k(f) k(i)]);
% v5 m& F" L. c8 y7 L: k- }    end
   
* j* X9 Z: H7 L8 C    for i=1:n
2 |* w( _% P$ g, c, P        for j=1:n
0 ~) p; h* _( L. J  h; i            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
+ M) Z2 z* \: y# N! F4 m8 V7 W            end
        end
. p  v, C9 l" H" @    end
end
end

7 V8 R/ H, E! T; c
! U+ T9 _) C7 \2 l2 B. L' mfunction nc=isncf(w,k)
2 b9 ~9 _6 |6 D: D3 M) ^  f3 c    nc=[];
: O: p3 n" a6 X, P% V    t=zeros(size(w));
    n=size(w,1);
    a=find(w~=0);
2 m% R& x+ X( s- _7 W* W* v- ?    for i=1:length(a)
        d(i)=w(a(i));
        if a(i)/n>floor(a(i)/n)
% G* Z; A( j- y          t(i)=floor(a(i)/n)+1;
        else
            t(i)=floor(a(i)/n);
        end
* h% S1 q, t3 @4 |7 |        t1(i)=mod(a(i),n);
        if t1(i)==0
0 W" J" k% j8 l- E0 u  d            t1(i)=n;
) ^; f8 c  N. p8 c: a/ {        end
$ S0 H! U) n+ `' h" d    end
0 @8 k, \' ~- k. N( ]2 T8 V. \    [b,c]=sort(d);
    p=[1];pc=0;
, E. b# Q0 }! w9 _& N$ n    for i=1:length(a)
        if k(t1(c(i)))<k(t(c(i)))
# ~+ a  V% w  h: V% r$ t            p=union(p,t(c(i)));
            t(t1(c(i)),t(c(i)))=3;
2 L. u6 f. L" e: f% X        end
2 i. r# F3 J- E' O7 T5 R        if pc==0
            tc=isempty(setdiff([1:n],p));
( }4 I! b* {8 ]            if tc
7 R1 O$ p: ~+ g- [/ f                t0=sum(t(1,==3);
8 o1 q3 c" q' _' q3 }( w3 t7 z                if  t0>=2
                    nc=union(nc,1);
                end
" u6 z* Y) I" d- d- }5 f1 M                break;
6 ~+ ~9 `% s/ R4 [. b            end
7 f  p3 t$ s& }: E* j1 b% U  [+ l        end
0 h5 c, V, U5 t& P! B    end
$ a" q+ t  d3 I4 W3 c   
        
0 H6 G- A9 [% m  Y, T& e; @end

ganquanlife

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