Real [int,int] data(nt,m):// nt: mesh of Time , m: degree of freedom (numbers of nodes in the mesh)
// m=Vh.ndof;
for(t=0;I<t;t++) //t : times
Data(t,:)=u^n+1[];
Real [int,int] data(nt,m):// nt: mesh of Time , m: degree of freedom (numbers of nodes in the mesh)
// m=Vh.ndof;
for(t=0;I<t;t++) //t : times
Data(t,:)=u^n+1[];