こちらのつづき ピボット操作というのは入っていないが、掃き出し法での式変形をもう少し再現できるようにしてみる N<-4 M<-5 A<-matrix(runif(N*M),N,M) C<-A L<-min(N,M-1) for(i in 1:L){ C[i,]<-C[i,]/C[i,i] for(j in 1:N){ if(j != i){ C[j,]<-C[j,]-C…

パッケージに "rootSolve"というものがあるらしい こちらのコメントで教えていただいたもの 詳しくは次を実行 vignette("rootSolve") 中身は以下の通り rootSolve uniroot.all : to solve for all roots of one (nonlinear) equation multiroot : to solve n…