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

ganquanlife

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
if n>=3
    a=sum(g);
8 m6 w7 H. {5 E    b=sum(a==2);
/ T) \  |8 m" _  d    if b==n
5 Y( s: W! T; e, y        fprintf('本图为圈，无割点。n')
        nc=0;
    end
else
, ]3 C6 E6 W+ W+ d4 _5 P1 V, V    [w,k]=dfs3(g);
  t7 H  r) M- c    %nc=[];
    nc=isncf(w,k);
# F1 J! r( g1 o% k, |3 |    n=size(g,1);
    for i=1:n
        for j=1:n
            if w(i,j)>1
% j% t0 r1 \$ d  q1 ?# I                if k(i)>k(j)
3 |8 L0 V7 b6 T# Z" d7 B1 Z                    g(i,j)=2;
                else
                    g(i,j)=3;
                end
3 G; j- x4 f& f; I- P3 Q* s            end
! t! d8 \- E0 w; E0 I; m* U+ B5 s2 {        end
    end
* g& V) X+ n# a2 M: E; u   
8 [' Q. r6 b% B9 e1 j* }: P' o$ d4 q    for i=1:n
* ?( Q5 l7 ]5 m5 R        f1=find(g(i,==2);
        f2=find(g(i,==3);
        f=union(f1,f2);
, ?8 T; ~3 D: L# t4 H) K) n9 V$ x        l(i)=min([k(f) k(i)]);
: a/ y& \( h/ K$ m8 K9 V3 L    end
   
    for i=1:n
        for j=1:n
$ K4 K) r% Y! D  E. G+ H7 q% L            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
            end
2 A7 B& ?6 W1 D        end
    end
end
end
2 ?% y+ I) v" `/ _7 L1 W

3 t/ I  s( S# ?, `$ P: b' o" d4 Ffunction nc=isncf(w,k)
* p: J) f% p& i) k7 p6 n- b    nc=[];
) {8 l! F! a6 R( z9 S    t=zeros(size(w));
* n' W& c0 L9 j% G6 w" F    n=size(w,1);
+ ]0 R  i; C9 z( A* c5 m! O    a=find(w~=0);
    for i=1:length(a)
) F+ T5 D  o7 |5 \% _        d(i)=w(a(i));
/ }9 P' [' P, U8 a        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
; r& C8 B0 [. p9 J% U% ?        else
8 V& G% U' o1 c            t(i)=floor(a(i)/n);
3 T. w7 i0 q# V$ _  `        end
0 }* n+ D$ V$ @1 A# y: q5 e        t1(i)=mod(a(i),n);
        if t1(i)==0
7 b, H6 E9 B9 R            t1(i)=n;
/ E- `" P  ?& ]9 e/ a        end
    end
    [b,c]=sort(d);
    p=[1];pc=0;
% O4 D7 A2 Z8 [6 @: x3 q& Y    for i=1:length(a)
/ P$ G! p/ K. Y: _( I5 g$ B$ ^        if k(t1(c(i)))<k(t(c(i)))
" I' Z- m9 ?+ }# w" d            p=union(p,t(c(i)));
- w/ o/ L' }+ R6 s+ h            t(t1(c(i)),t(c(i)))=3;
8 Q. j/ k' j( N        end
        if pc==0
            tc=isempty(setdiff([1:n],p));
            if tc
' V+ E" [& I1 U" M; W; J. ^                t0=sum(t(1,==3);
6 P% l3 g! F8 C& C7 R# K' y                if  t0>=2
+ A0 V- r- n" q* m$ W% V2 T                    nc=union(nc,1);
                end
/ ?- b) ]# v7 D4 L                break;
+ q, ~+ x& W) E# L% @            end
        end
1 b+ J, V4 p/ y% n9 k% V    end
! Y+ J$ Y8 @: Y6 D) E3 f   
        
1 [- t+ s- s- Hend

ganquanlife

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