使用最小二乘法的椭圆

x = [1.02 .95 .87 .77 .67 .56 .44 .30 .16 .01]';
y = [0.39 .32 .27 .22 .18 .15 .13 .12 .13 .15]';

% function Ax = b
N = numel(x);
A = [x.^2 x.*y y.^2 x y];
b = [-1*ones(N,1)];

%coeffs = [a b c d e]';
coeffs = A\b; 

a = coeffs(1)
b = coeffs(2);
c = coeffs(3);
d = coeffs(4);
e = coeffs(5);

plot(x,y, 'bo')
hold on

% prepare contour plot
xmin = min(x);
xmax = max(x);
ymin = min(y);
ymax = max(y);
%set deltas 
deltay = .001;
deltax = .001;
%create arrays X and Y
[X,Y] = meshgrid(xmin:deltax:xmax,ymin:deltay:ymax);
% evaluate quadratic form 
Z = a*X.^2 + b*X.*Y + c*Y.^2 + d*X + e*Y +f;
% plot where Z =0. line
contour(X,Y,Z,[0,0])
title('Planetary Orbit');

我试着做一个椭圆,但代码正在做一个抛物线。我想这是我设置等高线图的方式。如果有人能给我一些需要改变的方向,我将不胜感激。