h1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
h2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
h3 = surf(np,nd,Ua3,'FaceColor',[0,0,1],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8);
zdiff