syms?x?f;
syms?l;

data=zeros（110）;
y=zeros（110）;
for?p=1:10
????data（p）=rand（1）*100
????y（p）=rand（1）*100
end
m=50;
l（1）=1;
clear?figure;
clear?f;
for?n=2:m
l（n）=exp（x）/factorial（n-1）*diff（x^（n-1）*exp（-x）xn-1）;
end
for?i=1:m
eval（[‘vec_‘num2str（i）‘=subs（l（i）xdata）‘]）
end
vec_1=ones（1length（data））;
%?vec_2=data;
%?vec_3=ones（1length（data））;
%length（x）=length（x1）
matrix=zeros（mm）;
yx=zeros（1m）;
out=zeros（mlength（data））