N-R法求解雅可比矩阵出错 Matrix is singular to working precision. 矩阵奇异

diudiu1

Jacobicopy=Jacobi;   
4 U% o7 T6 p* I- W5 W# z8 @9 G   detaPQUcopy=detaPQU;
4 }' J1 q1 n) v' Q% M7 n   detaFEcopy=detaFE;
( `# y1 r- W4 s$ Q    for j=num:-1:2
( v: C9 l: ^- H: O2 |4 |) W        if flagofPV(j)==1
            Jacobicopy(2*(j-1),:)=[];
            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
' o+ q! S/ M) x1 C% U" H        end
    end   
8 M, u3 O' _& u# k$ }* ~, o2 J6 H   
6 ]" d9 E# }' A3 z+ d7 e detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
invofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy
( U6 w) I( J/ t4 |. y) m; i( |; b
, {' ]" ^( P: ]) e$ _+ G; V
------------------------------------------------------------------------------------------------

. Z- d" q. Y$ ^1 f% O# wJacobimodified =

1 _9 Y. ?6 ?  `" x3 b   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
    1.4851   77.6767   -0.8299    3.1120         0         0         0
    3.1120    0.8299   25.8604    4.7915         0  -31.7460         0
   -0.8299    3.1120    1.6223   45.6152         0         0         0
  -63.4921         0         0         0   63.4921         0         0
, `2 W! R+ G! }6 t         0         0  -31.7460         0         0   31.7460         0
         0         0         0  -31.7460         0         0         0

Warning: Matrix is singular to working precision.
> In yyg at 170

% a' U, f' j* m* v: D" q1 P# ginvofJ =

& u5 f0 R) A* y( ?: B! ~   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf