線形の連立方程式の解法 こちらのつづき 掃き出し法に近いもので、Gauss-Jordan法 拡大係数行列を変形させてみる N<-3 A<-matrix(runif(N*(N+1)),N,N+1) C<-A for(i in 1:N){ C[i,]<-C[i,]/C[i,i] if(i != N){ for(j in (i+1):N){ C[j,]<-C[j,]-C[j,i]*C[i,]…