figure;hold on
xv = linspace(5e5,1e7,100);
yv = linspace(0.6, 200, 100);
[X, Y] = meshgrid(xv, yv);
Z = (0.037 .* (X.^0.8) .* Y )./ (1 + 2.443 .* (X .^ -0.1) .* (Y .^ (2/3) - 1));
h=surf(log10(X), log10(Y), real(log10(Z)),'facecolor','interp','EdgeColor','none');
xlabel('Re'),ylabel('Pr'),zlabel('Nu');
zl=[2 6];
Xc=[X(1,1) X(1,end) X(end,1) X(end,end)];
Yc=[Y(1,1) Y(1,end) Y(end,1) Y(end,end)];
Zc=[Z(1,1) Z(1,end) Z(end,1) Z(end,end)];
Z0=[zl(1) zl(1) zl(1) zl(1)];
for i=1:4;
hline(i)=plot3(log10([Xc(i) Xc(i)]),log10([Yc(i) Yc(i)]),[real(log10(Zc(i))) Z0(i)],'k')
end
view(3)
grid on
set(h,'facealpha',0.7)
set(gca,'zlim',zl)
set(hline,'linewidth',2)
set(hline(4),'linestyle','--')
Best Answer