A=zeros(i+1,1);for j=0:i A(j+1,1)=q1*r1^(i-j)*q2*r2^j;end C{i+1，1}=A;end
a = cell(60,1);for i = 1:60 a{i,1} =eval(['data' num2str(i)]);end