1.3.m

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

r=randn(1,200)
a=0
for j=5.992:0.00008:6.008
   a=a+1
   r(a)=j
end
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