function?enkf（NMLT_0T_f）
%[UaUf]=enkf（NMLT_0T_f）
format?long;


%?u=zeros（L+1dT+1）;

dT=（T_f-T_0）/N?;%?T分N份
Tn=T_0:dT:T_f;
%lT=length（Tn）;
x_comp_length=100;
dx?=?x_comp_length/L;
x?=?0:dx:x_comp_length;
x?=?x‘;
dt?=dT/M;
for?n=1:N
?t?=?Tn（n）:dt:Tn（n+1）;
end
c?=?1;
%?display（‘this?number?should?be?less?than?1?for?stability:‘）
mu?=?c*dt/dx;



%?define?the?mesh?in?time
%?display（‘this?number?should?be?less?than?1?for?stability:‘）

%?for?i=1:lT;
%??for?j=1:10:L;
%?	temp=?normrnd（01）;%?隨機誤差
%?????for?r=1:10;
%?????????u（1j+r）=?a（1j+r）+temp;????

%?????????u（1100）=u（11）;
%?????end
%??end
P=10;

u=zeros（L+1M+1）;
v=zeros（L+1N+1）;
Uf=zeros（L+1P）;
Ua=zeros（L+1P）;
%ini_u=zeros（L+1P）;
%?Df=zeros（4N+1）;
%?a?=?normrnd（01L+11）;
%?a（L+