Define LibPub kronreduction(smat,row)= Func : :If dim(smat)[1]<row Then : Goto end :EndIf : :Local tmat,i,j :newMat(dim(smat)[1],dim(smat)[2])→tmat : :For i,1,dim(smat)[1] : For j,1,dim(smat)[2] : smat[i,j]-((smat[i,row]*smat[row,j])/(smat[row,row]))→tmat[i,j] : EndFor :EndFor : :Return tmat :Lbl end :EndFunc Kron Reduction 이라는 것을 알게되어, 개략적인 기능만을 구현하여 보았습니다. 간단한 함수입니다만, 행렬을 어떻게 처리하는게 좋을지...