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

ganquanlife

function [nc] = ncutf(g)
%求割点的算法 g为邻接矩阵 nc为割点的集合
n=size(g,1);
if n>=3
9 c3 M, M+ I6 W3 x7 a+ n5 Q' w  _2 S    a=sum(g);
* w1 p- J1 |. y) e! G/ w7 q6 S    b=sum(a==2);
6 S1 e1 I; w# j; `( ~$ m! i# a    if b==n
        fprintf('本图为圈，无割点。n')
        nc=0;
    end
/ E( Z- C, X. q+ U$ p+ }else
: a- g# ^, N0 T% E- b    [w,k]=dfs3(g);
# [5 [+ M! o( k7 N9 d/ M4 P3 ~, i/ k    %nc=[];
    nc=isncf(w,k);
- p( |+ X2 V3 e2 A$ q5 P    n=size(g,1);
    for i=1:n
1 u% K0 P% j4 K. w6 \6 o  }& |9 {; m        for j=1:n
" ~$ k& J/ p* }( Q! n( v            if w(i,j)>1
                if k(i)>k(j)
$ L/ I" `* k' s5 f  k5 c                    g(i,j)=2;
                else
                    g(i,j)=3;
                end
            end
        end
, ?0 @3 i& u, c7 G    end
   
$ c+ ~3 X0 N: ]& I( b    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)]);
0 ?/ Y! W- Y7 o$ E8 k! p    end
5 O4 B+ N. R, A0 W" r; V4 w/ q) g   
; b4 _: l/ F+ h- r) D% [& {1 V* f    for i=1:n
- Q/ L# j0 x6 X& \9 ^6 d7 P" e        for j=1:n
3 T+ p; m! a$ Z/ e' U            if g(i,j)==3 & k(i)>1&l(j)>=k(i)
                nc=union(i,nc);
* h( g* [( J7 ^* S& M# W            end
        end
6 U! T# G+ H# t  v) j0 c1 X% H    end
6 l" M3 a* |4 Z7 V: k+ [end
end
5 i9 Z* M- z8 h2 X9 Z( l
3 {) y  ?" g* P) G- H) h
4 d+ q; X' w- [4 d% r" ~function nc=isncf(w,k)
6 z2 }5 c4 q/ k1 C! v+ l% M6 Z  o) E4 `6 e    nc=[];
9 p9 L1 t$ Y1 R9 Q    t=zeros(size(w));
4 Z6 l* \8 X& z5 d+ W% U- \    n=size(w,1);
    a=find(w~=0);
    for i=1:length(a)
* l* \. L. f1 N+ S, d8 h) b        d(i)=w(a(i));
  P) P1 E8 J6 j" y4 m        if a(i)/n>floor(a(i)/n)
          t(i)=floor(a(i)/n)+1;
        else
6 S8 Z  e' a, E: i: D            t(i)=floor(a(i)/n);
        end
6 _: t9 A  @6 \# c6 D4 A        t1(i)=mod(a(i),n);
        if t1(i)==0
            t1(i)=n;
        end
% Y5 A% K' ?2 S3 C    end
% n& Q  R2 t* B1 x# S& e! R    [b,c]=sort(d);
, N; W1 h$ w% z" V    p=[1];pc=0;
/ q5 o# ]+ f# N- o; n  Y4 `  Z* X    for i=1:length(a)
+ ?0 {4 \5 [5 j        if k(t1(c(i)))<k(t(c(i)))
$ Y, ?' C* D3 F& P, H. v' A- U            p=union(p,t(c(i)));
% }- u  {+ U  y; i            t(t1(c(i)),t(c(i)))=3;
% i  Y- r2 c5 n        end
# m" ?& U$ e6 Y' L2 J4 V        if pc==0
            tc=isempty(setdiff([1:n],p));
, b* ^# d/ _; ~- j8 Z            if tc
                t0=sum(t(1,==3);
                if  t0>=2
                    nc=union(nc,1);
3 e" n9 g7 M. V: y! `                end
                break;
            end
        end
7 O$ |8 @9 K' X2 u3 W2 y. V  j    end
$ M, O  l1 [& h" a! K- g9 W1 j   
5 c3 ]8 t. {9 K& ^        
end

ganquanlife

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