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');