1.6.m

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

r=randn(1,200)
a=0
for j=5.992:0.00008:6.008
   a=a+1
   r(a)=j
end
c=[1,-30,360,-2160,6480,-7776]
p=polyval(c,r)
for j=1:200
p1(j)=r(j)^5-30*r(j)^4+360*r(j)^3-2160*r(j)^2+6480*r(j)-7776
p2(j)=(r(j)-6)^5
end
subplot(2,1,1);plot(p);
hold on
subplot(2,1,2);plot(p1,'r');
hold on
subplot(2,1,1);plot(p2,'k');
hold on
subplot(2,1,2);plot(p2,'k');