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

ganquanlife

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
( I+ n7 l0 U# T  s! @, K; vn=size(g,1);
if n>=3
6 r/ D" ?7 S/ V' h1 s" G$ |% n    a=sum(g);
    b=sum(a==2);
    if b==n
7 I/ m; l& Y5 [  [        fprintf('本图为圈，无割点。n')
        nc=0;
    end
else
- W- w+ i8 _4 f3 O    [w,k]=dfs3(g);
2 J! b; q) ]% L- q    %nc=[];
) m" f. r) f2 I3 q' a! q/ w$ q    nc=isncf(w,k);
    n=size(g,1);
4 Q! j; s1 o. x6 z9 E    for i=1:n
        for j=1:n
            if w(i,j)>1
; U7 Q4 S: c& H/ Z5 H( b1 ^$ G                if k(i)>k(j)
3 P4 I& ]2 k1 X: U                    g(i,j)=2;
" Y0 b; m% A( X- b, n. C2 q                else
* W& ~( e+ A; [, `6 A                    g(i,j)=3;
  Z( L6 G+ b- ^" Q5 d+ ^0 ^                end
            end
. F; _" J5 D6 r) v1 G- r1 y        end
- R) y$ z9 J5 }' l! E4 j1 B: q    end
9 q+ t9 \' B1 e6 X3 r* l& p" s   
' }. B( a) I# ?5 S    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)]);
    end
   
    for i=1:n
3 ~6 }# W9 J2 y* n        for j=1:n
& q, C9 W7 _$ {9 @0 \: _+ T* V0 a            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
5 _. f# a  ^0 S            end
' |3 k. c, g8 w        end
    end
end
0 ^4 d6 k# q/ n7 w6 Cend


function nc=isncf(w,k)
# B5 C% r5 O7 o/ F& `9 \1 [' c: K    nc=[];
    t=zeros(size(w));
. X3 D% q, C  p2 E. q/ K3 N# p, r    n=size(w,1);
: m. K' O. v2 I! k/ [4 J' l    a=find(w~=0);
    for i=1:length(a)
        d(i)=w(a(i));
  D& W0 D: C0 E" M: h/ A- I5 ?% m" O        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
        else
  P* o0 h& D7 R5 p+ A" o5 g& R            t(i)=floor(a(i)/n);
4 H' D# Y6 i/ d+ O        end
        t1(i)=mod(a(i),n);
3 S0 e5 h# J5 C        if t1(i)==0
6 G; w8 d" o9 V2 ]& Q' m) C            t1(i)=n;
        end
    end
    [b,c]=sort(d);
2 X- A% o) {5 N. I5 u    p=[1];pc=0;
. E& I. w8 G* c; N    for i=1:length(a)
/ F0 P7 F: k% _& v, t$ ]  d5 Y# b        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
1 \) H9 t% m- V! S$ t* y) _3 I: E            t(t1(c(i)),t(c(i)))=3;
( r& _0 a6 t, s3 {        end
/ {4 W; z  `& u  d        if pc==0
& `+ a# [8 _7 {6 z3 k, ?            tc=isempty(setdiff([1:n],p));
            if tc
# f& o) n9 C1 G7 N- Z& @  p2 v2 C                t0=sum(t(1,==3);
8 t' J! \2 \8 T! T" X( h                if  t0>=2
  _3 K: v9 @: ]$ A9 u$ u                    nc=union(nc,1);
                end
+ l' G, b8 S1 h7 {' ^                break;
) T- Y, S* H8 b7 R1 L            end
3 F/ N1 ~8 I' J: D* O# s- n        end
* V; U0 n6 s! s! s7 J' p* d6 n    end
   
        
end

ganquanlife

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