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

ganquanlife

function [nc] = ncutf(g)
5 ?( J9 m  f. J9 I% y7 ]6 z# D%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
/ X) L0 u% Q$ b- T! X* T7 aif n>=3
( @$ J7 J. }2 K    a=sum(g);
9 w/ e6 A8 C+ e    b=sum(a==2);
    if b==n
        fprintf('本图为圈，无割点。n')
% Z2 R, w7 ?" X, }! O        nc=0;
    end
else
9 A; Z- G+ R  D0 s2 U    [w,k]=dfs3(g);
    %nc=[];
    nc=isncf(w,k);
    n=size(g,1);
* j2 g  s: V6 [5 a7 h    for i=1:n
( t: E7 n; V( a) w# m, z& {$ s; k" s" y! Y        for j=1:n
  R. O. [- v) u1 H' \3 z3 H: G1 a            if w(i,j)>1
                if k(i)>k(j)
) q- N- l% J9 l; j                    g(i,j)=2;
                else
. ^7 s& Q2 T: ]+ ^$ A                    g(i,j)=3;
8 d, {( u- q$ h% t2 ~. s                end
            end
        end
    end
" y3 E) Q+ F# A! G. X( w' G  K) s   
    for i=1:n
  u6 B+ F, h; U8 U9 Q0 |  w        f1=find(g(i,==2);
8 h9 g/ ?8 j; [  t4 J; F        f2=find(g(i,==3);
        f=union(f1,f2);
        l(i)=min([k(f) k(i)]);
( A$ f5 }5 c$ K2 \, q    end
+ Y; h% v0 h$ Q4 B/ q6 _8 C( O4 \   
5 S: [. ~3 l3 i+ F0 G9 V7 n9 [6 h, @" Q    for i=1:n
0 c% Y* T. C  j2 |6 L' `. f        for j=1:n
            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
  R# d; ?8 P9 C0 P' M# n            end
        end
    end
end
  Y9 G' q$ q' ~8 Zend
& w  k9 A. w% ?8 _8 [) d6 l
+ w6 H+ ]$ Y/ a% x  ~  Q6 a
function nc=isncf(w,k)
    nc=[];
  D! o* d; m0 T0 T$ @0 ^/ k. g    t=zeros(size(w));
$ F9 P4 z6 d+ k& A    n=size(w,1);
    a=find(w~=0);
    for i=1:length(a)
, F3 n3 y! I5 f# ^5 m) k! H        d(i)=w(a(i));
        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
: m) g( L5 L- ~, S/ ]+ ?8 e        else
/ x6 @! o6 D1 p( I; T! u1 }            t(i)=floor(a(i)/n);
$ b1 A6 _) l) o: A' o1 P        end
- I* r: k6 l! u/ Q        t1(i)=mod(a(i),n);
        if t1(i)==0
            t1(i)=n;
% V$ a8 H$ ?+ ^5 B9 e, C4 D  ], f        end
2 h. w6 m, F& r# O    end
) M& L0 U9 J" _7 x8 X7 L" j* r    [b,c]=sort(d);
    p=[1];pc=0;
; \/ R+ A/ G' n( {% Z1 ?    for i=1:length(a)
        if k(t1(c(i)))<k(t(c(i)))
            p=union(p,t(c(i)));
& x  U: q0 v  s9 S6 Y2 V2 p/ C4 l            t(t1(c(i)),t(c(i)))=3;
        end
        if pc==0
. }( V( Q: s: N( \            tc=isempty(setdiff([1:n],p));
/ ^8 A9 u5 B2 N6 p- B, j            if tc
                t0=sum(t(1,==3);
! W. L! R4 Y, n4 {' H9 h                if  t0>=2
% c" H, u. M# j! d                    nc=union(nc,1);
                end
                break;
            end
        end
) O8 T3 H6 o; e' P    end
   
" C  {  a) O! Q$ S3 D# e        
end

ganquanlife

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