%x=-1.2:.01:-.8;
%z=-0.2:.01:.2;
 x=0:.1:.2;
 z=0:.01:2;
 [X,Z] = meshgrid(x,z);
%f=((1-t)-1).^2     +   ((1-2.*t).*(1-t)-mu).^2;
f=(Z-1).^2     +   (X.*Z-mu).^2;
figure(1);
contour(x,z,f,20);
figure(2);
mesh(X,Z,f);
