% https://en.wikipedia.org/wiki/Matrix_representation_of_conic_sections
clear
x=[-0.3616  1.1616  1.1616 -0.3616   0.4000]';
y=[-0.1707  0.7707 -0.9323  1.5323  -0.9042]';
disp('      A         B         C         D         E        F')
disp(N)
A=N(1); B=N(2); C=N(3); D=N(4); E=N(5); F=N(6);
disp ('      xc        yc')
disp([xc,yc])
AQ = [A B/2 D/2;B/2 C E/2;D/2 E/2 F];
A33 = AQ(1:2,1:2);  % det(A33)>0 & det(AQ)<>0
K = -det(AQ)/det(A33);
a = Lambda(1) ; b = Lambda(2);
disp('axes length ')
disp('semi-major,  semi-minor')
disp([a, b])
axis([-1.5 2 -1.5 2])
hold on
plot(xc+[1,-1]*V(1,1)*a,yc+[1,-1]*V(2,1)*a,'--')
plot(xc+[1,-1]*V(1,2)*b,yc+[1,-1]*V(2,2)*b,'--')
plot(x,y,'o','MarkerFaceColor','r');grid

方程组求解就是了