u=@(x,y) 3*exp(-((x-1).^2+y.^2)/10).*cos(x-1).*sin(y); [q2,me2,toGlobal2]=fc_meshtools.simplicial.getMeshArrays('2D',2); [q1,me1,toGlobal1]=fc_meshtools.simplicial.getMeshArrays('2D',1); U2=fc_meshtools.eval(u,q2); figure(1) fc_graphics4mesh.plot(q2,me2,U2,'plane',true) shading interp axis image;axis off colorbar figure(2) fc_graphics4mesh.plotmesh(q2,me2,'color','LightGray') hold on fc_graphics4mesh.plot(q1,me1,U2(toGlobal1),'LineWidth',2,'Plane',true) caxis([min(U2) max(U2)]) axis image;axis off colorbar

