%clear all
fem1 =  ModelUtil.model('Model');
mesh1 = mphxmeshinfo(fem1,'solname', 'sol1');
xxx = mesh1.nodes.coords(1,:);
yyy = mesh1.nodes.coords(2,:);
soltranspose  = fem1.sol('sol1').getU(1);
soltranspose2 = mphgetu(fem1,'solnum',1)
solu = soltranspose';
[XX,YY] = meshgrid(linspace(-1,+1,100),linspace(-1,1,100));
ZZ = griddata(xxx,yyy,solu,XX,YY);
surf(XX,YY,ZZ)
plot3(xxx,yyy,solu,'.');

pause
solu = soltranspose2';
[XX,YY] = meshgrid(linspace(-1,+1,100),linspace(-1,1,100));
ZZ = griddata(xxx,yyy,solu,XX,YY);
surf(XX,YY,ZZ)
plot3(xxx,yyy,solu,'.');