1.7.m

来自「学习使用各种数值方法计算数值代数问题。其中包括计算矩阵的特征值问题」· M 代码 · 共 10 行

c=[1,-30,360,-2160,6480,-7776]
r=roots(c)
axis ('equal')
for j=1:5
    if (imag(r(j))==0)
        r(j)=r(j)+0.0000001*i
    end
plot(r(j),'*')
hold on
end