fem1 = mphload('1d_helm_6eigs.mph')
mesh1 = mphxmeshinfo(fem1,'solname', 'sol1');
xxx = mesh1.nodes.coords(1,:);
soltranspose = mphgetu(fem1,'solnum',1)
solu = soltranspose';
plot(solu)
pause;  display('Press a key to continue...')
%  the solution is ordered from left to right
plot(xxx,solu)
%  the mesh is NOT ordered from left to right meaning that the vectors for
%  the solution and the mesh are NOT SYNCED.