function?[Lcc]=yima（Mne）
%y=1:2^n;
%Y=daoweixu（yn）;
%LL=-ones（12^n）;
%for?i=1:2^n
?%???LL（i）=M（Y（i））;
%end
w0=1-e;
w1=e;
W0=e;
W1=1-e;
LL=daoweixu（Mn）;
for?i=1:2^n
????if?（LL（i）==0）
????????LL（i）=w0/w1;
????else
????????LL（i）=W0/W1;
????end
end
[~P]=polarize（ne）;
L=-ones（2^nn+1）*0.1;
L（:n+1）=LL;
u=-ones（12^n）;
S=-ones（2^nn）;???????????%%%%%%%%%???L?S?V?初始化為-1的矩陣
V=-ones（12^n-2）;
for?j=1:n?????????????%指數項矩陣右上角初始賦值0
????for?i=1:2^（j-1）
????????S（ij）=0;
????end
end
l=1;
while（l<=2^n）
????if?（P（l）~=1）
????????u（l）=0;
????else
????????j=1;
????????i=l;
????????if?（l>1）?????????????????%%%%%%%%
????????????S（l1）=u（l-1）;
????????end???????????????????????%%%%%%%%
????????if?（l>2）
????????????if?（mod（l2）==0）
????????????????p=l-2;????????????????
????????????????for?t=1:p
????????????????????V（t）=u（t）;
????????????????end
????????????else
????????????????p=l-1;
????????????????for?t=1:p
????????????????????V（t）=u（t）;
????????????????end
????????????end
????????end
????????while（j????????????if?（rem（fix（（i-1）/（2^（j-1）））2）==1）
????????????????if?（L（i-2^（j-1）j+1）==-0.1）
???????????????????????%%%%%%%%%%%%%%%%%%%%%子程序2???這里?
????????????????????????????if?（l>=3）????????????????????????
????????????????????????????????p=fix（（p+1）/2）;
????????????????????????????????for?q=1:p
????????????????????????????????????V（q）=mod（V（2*q-1）+V（2*q）2）;
????????????????????????????????end
????????????????????????????????i=i-2^（j-1）;
????????????????????????????????j=j+1;
????????????????????????????????if（S（ij）==-1）
????????????????????????????????????S（ij）=V（p）;
????????????????????????????????%else
?????????????????????????????????%???if?（V（p）~=S（ij））
??????????????????????????????????%??????k=10000000
???????????????????????????????????%?end
????????????????????????????????end
????????????????????????????????i;
????????????????????????????????j;
????????????????????????????????V;
????????????????????????????else
????????????????????????????????i=i-2^（j-1）;
????????????????????????????????j=j+1;
????????????????????????????end
????????????????????elseif?（L（i-2^（j-1）j+1）~=-0.1&&L（ij+1）==-0.1）
?????????????????????????????%%%%%%%%%%%%%%%%%%%%%子程序4???????????????????????????
????????????????????????????if?（l>=3）
????????????????????????????????p=fix（（p+1）/2）;
????????????????????????????????for?q=1:p
????????????????????????????????????V（q）=V（2*q）;
????????????????????????????????end
????????????????????????????????j=j+1;
????????????????????????????????if（S（ij）==-1）
????????????????????????????????????S（ij）=V（p）;
???????????????????????????????%?else
????????????????????????????????%????if?（V（p）~=S（ij））
?????????????????????????????????%???????k=10000000
???????????????????????????????????%?end
????????????????????????????????end
????????????????????????????????i;
????????????????????????????????j;
????????????????????????????????V;
??????????????