Gauss-Jordan法
- こちらのつづき
- ピボット操作というのは入っていないが、掃き出し法での式変形をもう少し再現できるようにしてみる
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[j,i]*C[i,] } } } C if(M-N==1){ B<-A[,-M] x<-C[,M] B%*%x-A[,M] }