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

ganquanlife

function [nc] = ncutf(g)
  I. X5 \1 y( K0 I$ y%求割点的算法 g为邻接矩阵 nc为割点的集合
8 r0 t( ]8 v+ L+ I( ln=size(g,1);
8 k- }3 x' p' o2 a, ?& {if n>=3
( a; {2 u: x: [' R9 ~    a=sum(g);
( I7 F5 z/ C* ?& a    b=sum(a==2);
! t2 V" z  B# r    if b==n
        fprintf('本图为圈，无割点。n')
        nc=0;
    end
+ Q& W# K) G+ b& t* b( p: c" S/ L/ belse
1 i2 c: L1 f9 _- W    [w,k]=dfs3(g);
" x: q( T/ M: r    %nc=[];
6 K8 y. E# _+ ]6 j( [" {/ Z- X5 K    nc=isncf(w,k);
    n=size(g,1);
. a' `$ T( ]  ~1 @, P& \" s    for i=1:n
$ k8 C# ]4 A0 ?. O        for j=1:n
            if w(i,j)>1
                if k(i)>k(j)
3 t. M7 d* J/ u8 B& r4 ^! x3 J' T4 ?                    g(i,j)=2;
( x, V, C5 D* L/ Q                else
                    g(i,j)=3;
  H) n5 t. A* w( @& Q: e! _% _; I. [5 `                end
            end
        end
    end
   
    for i=1:n
' B4 f, q( _) T3 ^$ u2 X        f1=find(g(i,==2);
+ ^0 U2 L- h* c; R        f2=find(g(i,==3);
6 v# H, i8 j( e) e        f=union(f1,f2);
        l(i)=min([k(f) k(i)]);
    end
   
& |! \3 t% }, t2 |    for i=1:n
8 |2 V# i2 X* W9 P! X5 R8 c: J- [        for j=1:n
% t( N  M- Q  s3 M4 V0 _            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
7 t! p+ {9 U, f                nc=union(i,nc);
0 A7 @* q7 R" F( Z3 S            end
- v/ P5 `, T" O8 g/ k: ^        end
    end
end
end


function nc=isncf(w,k)
    nc=[];
    t=zeros(size(w));
( s2 R3 Z3 E: Q4 l0 _5 h  z& A    n=size(w,1);
    a=find(w~=0);
! m' H7 a7 s  T9 e+ b    for i=1:length(a)
% d4 P1 y4 Z/ Q0 y        d(i)=w(a(i));
        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
7 t# R3 }: z9 z        else
            t(i)=floor(a(i)/n);
6 {4 q9 ]/ X0 I& G        end
1 r8 I7 R9 J/ M" }" F        t1(i)=mod(a(i),n);
# T. K9 P  H) N; ?. ]4 l9 u8 r" r( W  N        if t1(i)==0
7 B5 L: e  J6 n7 Q. ]            t1(i)=n;
        end
6 T7 l. j' m. L/ A5 O    end
6 O) ]3 ?* E7 l8 Q    [b,c]=sort(d);
% S7 L2 K- J. H/ ?& F" o4 K  I! `    p=[1];pc=0;
    for i=1:length(a)
        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
            t(t1(c(i)),t(c(i)))=3;
3 W' P. e8 e+ X+ K: X: x8 p- k9 X2 o        end
) ~6 q; O% n- V  C; U+ p        if pc==0
            tc=isempty(setdiff([1:n],p));
% p3 E/ d$ c8 n! O/ T            if tc
                t0=sum(t(1,==3);
                if  t0>=2
) V4 u2 q/ l! S                    nc=union(nc,1);
3 G1 `1 K9 s0 M; o1 Q                end
                break;
$ t4 o1 k' C8 n; Z& O5 X            end
        end
    end
; N/ O" c. ?9 c7 U( M1 m- }0 g   
        
5 a: p6 K7 e0 u, E; Bend

ganquanlife

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