clear all
[x1,x2]=meshgrid(–1:0.05:4, –1:0.05:6);
f=–2*x1+x2;
g1=x1+2*x2–6;
g2=x1–2;
g3=–x1;
g4=–x2;
cla reset
axis auto
xlabel(‘x1’),ylabel(‘x2’)
title(‘Exercise 8.39’)
hold on
cv1=[0:0.04:0.6];
contour(x1,x2,g1,cv1,‘g’);
cv2=[0:0.01:0.02];
contour(x1,x2,g1,cv2,‘k’);
cv3=[0:0.02:0.25];
contour(x1,x2,g2,cv3,‘g’);
cv4=[0:0.01:0.01];
contour(x1,x2,g2,cv4,‘k’);
cv5=[0:0.025:0.3];
grid off;
hold off;