
clear all;

[x,y]=ellipse_data;
V=[x.^2 y.^2];
t=ones(size(x));
c=V\t;
format long
c

clf;
hold on;
plot(x,y,'ro');
theta_test=0:0.01:2*pi;
a=sqrt(1/c(1)) 
b=sqrt(1/c(2))
plot(a*cos(theta_test),b*sin(theta_test),'k-');
hold off;