主题:请问各位大牛,怎么求复数矩阵的行列式
我想求一个一般复数矩阵的行列式,请问大家该怎么编程啊,就是算det()
我用的是先把一般复数矩阵的本征值求出来,然后让所有本征值相乘的方法
可是算出来有问题,大家帮我看看
CALL DEVCCG(m/2-1,Ax,m/2-1,evalx,vrx,m/2-1)
do i=1,m/2-1
write(222,*) i,evalx(i)
enddo
det=(1.0,0.0)
do i=1,m/2-1
DET=det*evalx(i)
enddo