a=randi(10,1,4)
I=1:length(a);
J=[I;I-1];J=J(:)'
str=sprintf('a(%d)*t.^(%d)+',J)
str=['@(t) ',str(1:end-1)]
P=eval(str)
clear a
P(2)
P(1:5)
