@fffm_84 The savemesh function saves the coordinates of the mesh, and Vh(i,j) is showing the connectivities of the mesh. So they are different things. If you want to save the data (fespace), use the ffmatlib.idp. I added below my edited version with the load data function added.