Tout(k,1)=k*T;

Nu=2;

vd2=0;

u_real=zeros(Nr,2);

u_piao=zeros(Nr,2);

Q_cell=cell(Tsim,Tsim);

for i=1:1:Tsim

   for j=1:1:Tsim

      if i=j

         Q_cell(i,j)=q;

      else

         Q_cell(i,j)=zeros(Nx,Nx);

         end

      end

end