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

diudiu1

Jacobicopy=Jacobi;   
6 o8 }: Z/ I: [6 @   detaPQUcopy=detaPQU;
: e3 R3 s* k1 u* |3 q. a   detaFEcopy=detaFE;
    for j=num:-1:2
        if flagofPV(j)==1
            Jacobicopy(2*(j-1),:)=[];
* ^" A" X1 }3 e' w7 C: l  f' _            Jacobicopy(:,2*(j-1))=[];  
# J5 R0 F# p& T. @            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
0 R# s1 _! a& ?2 W  K        end
    end   
  ?; |, D( q/ H. H   
( ~' f- O! S4 C! B# h detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
/ _# A- D6 U! rinvofJ=inv(Jacobicopy)
! ?! L& U* d$ f& D. U detaFEcopy=invofJ*detaPQUcopy
% ?$ N+ c  |1 w4 M3 u

) c/ ]% Z. g; w, ?3 u5 o. y9 c$ N------------------------------------------------------------------------------------------------
, x. f7 ^; P) S3 A
Jacobimodified =

; k' C* ?. H5 Y. y" V   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
$ M6 c4 M$ l2 K  -63.4921         0         0         0   63.4921         0         0
         0         0  -31.7460         0         0   31.7460         0
         0         0         0  -31.7460         0         0         0
2 Y$ e4 d4 d2 o7 M
Warning: Matrix is singular to working precision.
. x. f! @/ i- Y2 T4 P& e. t+ H> In yyg at 170
& v( K/ A4 E, s' L: y% {: `; @. D
invofJ =
/ N+ D* H0 v$ Q5 ?
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
% c9 g/ ]$ _! ~$ N+ g+ p   Inf   Inf   Inf   Inf   Inf   Inf   Inf
! x) |& ~5 H/ B4 L- b   Inf   Inf   Inf   Inf   Inf   Inf   Inf
1 Y: Q  u: w4 i* T( n7 F   Inf   Inf   Inf   Inf   Inf   Inf   Inf
- {# d) a5 o! [& G1 e& N   Inf   Inf   Inf   Inf   Inf   Inf   Inf