接上一篇

function GET(m)

i=1
for x=-1:0.01:1
H(1)=x
H(2)=(3x^2-1)/2
for i=1:m
H(i+2)=((2i+1)xH(i+1)-i*H(i))/(i+1);
end
i=1
subplot(3,3,i)
plot(x,H(i+2),’-b.’,‘linewidth’,3)
hold on
subplot(3,3,i+1)
plot(x,H(i+3),’-r.’,‘linewidth’,3)
hold on
subplot(3,3,i+2)
plot(x,H(i+4),’-g.’,‘linewidth’,3)
hold on
end
figure(gcf)