用matlab程序求多项式a(s)=(s+3)(s+4)(s+6)的展开式

可以使用符号函数的expand()展开多项式，例如：

>> syms s
>> a = expand((s+3)*(s+4)*(s+6))
 
a = s^3 + 13*s^2 + 54*s + 72

你是要输出每个i对应s的值吧。程序如下
s=ones(1,20);
y1=1;
y2=1;
y3=1;
fori=0:1:20
a1=y1;a2=y2;a3=y3;
y1=(100*a2+150*a3)*0.5;
y2=(0.09*a1)*0.5;
y3=(0.2*a2)*0.5;
s(i+1)=y1+y2+y3;
plot(i,s(i+1),'.')
holdon
vpa(s(i+1),6)
end
gridon
axis([020120650000])

结果s的值：

ans =

125.145

ans =

15.3795

ans =

282.589

ans =

76.8403

ans =

641.015

ans =

268.264

ans =

1468.22

ans =

819.938

ans =

3394.03

ans =

2340.38

ans =

7913.3

ans =

6411.35

ans =

18594.8

ans =

17096.3

ans =

44002.1

ans =

44742.3

ans =

104775.0

ans =

115521.0

ans =

250844.0

ans =

295284.0

ans =

603387.0