1.4.m

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

r=rand(1,200)*0.016+5.992
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');