nclear, clf, hold
off
nx1=0:0.1:2; y1=-2:0.1:2; [x,y]=meshgrid(x1,y1); [f1,f2]=funnonlin(x,y);
nfigure(1)
nsubplot(1,2,1)
nmesh(f1,'linewidth',2), hold on,
mesh(f2,'linewidth',2),
naxis([min(x1)
max(x1) min(y1) max(y1) -10 10]);
nxlabel('x'); ylabel('y');
zlabel('z'); grid on; hold off,
nsubplot(1,2,2)
n[c,h]=contour(x,y,f1,'-r','linewidth',2);
clabel(c,h);
nhold on
n[c,h]=contour(x,y,f2,'linewidth',2);
clabel(c,h);
nhold off
naxis([min(x1)
max(x1) min(y1) max(y1)]); xlabel('x');
ylabel('y'); grid on;
legend('f1','f2‘,2)
n
nx2=0:0.1:20;
y2=-2:0.1:20;
[x,y]=meshgrid(x2,y2);
[f1,f2]=funnonlin(x,y);
nfigure(2)
nsubplot(1,2,1)
nmesh(f1), hold on,
mesh(f2), axis([min(x2) max(x2)
min(y2) max(y2) -10 10]);
nxlabel('x');
ylabel('y'); zlabel('z'); grid
on; hold off,
nsubplot(1,2,2)
n[c,h]=contour(x,y,f1,'-r','linewidth',2);
clabel(c,h);
nhold on
n[c,h]=contour(x,y,f2,'linewidth',2);
clabel(c,h);
nhold off
naxis([min(x1)
max(x1) min(y1) max(y1)]); xlabel('x'); ylabel('y'); grid on; legend('f1','f2')
n