clcx0=-1.682:.01:1.682;[x,y]=meshgrid(x0);z=7.163514-3.094629*x+2.15714*y+1.566561*(x.^2)+1.904204*(y.^2);mesh(x,y,z);hold on;contour(x,y,z,100);