1.2.m

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

r=randn(1,200)*0.016+5.992
s=zeros(40,6)
a=0
for j=1:5:196
    a=a+1
    k=[r(j),r(j+1),r(j+2),r(j+3),r(j+4)]
    s(a,:)=poly(k)
end
b=0;
for b=1:6
plot(s(:,b))
hold on
end