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

ganquanlife

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
2 S4 B! D- M% S7 o$ g# {1 N& k$ Hn=size(g,1);
if n>=3
! @4 x* S, e) i# n3 |    a=sum(g);
    b=sum(a==2);
+ A9 S  ]" S6 A( X) U3 m    if b==n
/ y! p: O$ m1 V' T% z        fprintf('本图为圈，无割点。n')
5 R" u  d9 J7 d% g0 D2 `4 M        nc=0;
- K9 O( X: w3 N$ G& \% ~7 p    end
else
6 x6 @+ y% ?$ }2 C+ K: h- Q    [w,k]=dfs3(g);
) w- W9 B- z$ S% u    %nc=[];
    nc=isncf(w,k);
    n=size(g,1);
# {# |: J  I8 P' B9 @: Y8 V    for i=1:n
; J3 z* X$ c- X: o) ~) Q- v' q6 H        for j=1:n
) q3 s; l# z+ J) W( o/ k7 b/ y            if w(i,j)>1
                if k(i)>k(j)
                    g(i,j)=2;
  q1 z3 h2 Q# E0 }$ z# C7 Y                else
                    g(i,j)=3;
                end
            end
5 q- L* r0 s5 X        end
    end
   
# m4 h3 `( {8 J9 O* q  p    for i=1:n
        f1=find(g(i,==2);
        f2=find(g(i,==3);
        f=union(f1,f2);
4 c5 L! N* P1 d/ g' w        l(i)=min([k(f) k(i)]);
# r0 y  m2 ~1 E% ]1 e3 T* g/ x    end
   
    for i=1:n
: V! @! U! z5 A" v& V        for j=1:n
            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
            end
) Q& ]/ A9 S, v  q; n: d% h        end
    end
4 A) j7 C2 B/ Z+ t/ S1 p7 x5 Tend
end
  U# q9 f& t. r! ^2 ^1 ?8 q

function nc=isncf(w,k)
    nc=[];
* V# U' [. ]5 L5 n. h    t=zeros(size(w));
! q9 G- s0 g0 S3 P    n=size(w,1);
    a=find(w~=0);
; y" n4 G5 }0 k) s    for i=1:length(a)
3 Y4 K* z) g) g0 F) ~        d(i)=w(a(i));
7 X0 r( q  e- p        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
        else
            t(i)=floor(a(i)/n);
        end
& V; \4 F' Q! C) T0 d        t1(i)=mod(a(i),n);
        if t1(i)==0
            t1(i)=n;
        end
    end
    [b,c]=sort(d);
    p=[1];pc=0;
4 E# t& D$ ^+ P" S% m7 [- D    for i=1:length(a)
3 X- }" V4 l! o+ n6 @        if k(t1(c(i)))<k(t(c(i)))
+ r/ z& x" D$ l7 m7 ~: t. w0 x            p=union(p,t(c(i)));
            t(t1(c(i)),t(c(i)))=3;
        end
        if pc==0
% ^+ V& \. Q% w2 y; u. t3 G            tc=isempty(setdiff([1:n],p));
3 }5 X4 J4 [. ^$ W8 i7 w            if tc
0 L1 Q2 l% D# {. f                t0=sum(t(1,==3);
; N3 @' z7 |4 p1 {2 h! B                if  t0>=2
                    nc=union(nc,1);
                end
                break;
            end
        end
    end
   
        
end

ganquanlife

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