T = readtable("https://jp.mathworks.com/matlabcentral/answers/uploaded_files/1716431/1%E5%88%97%E7%9B%AE%E3%81%8C%EF%BD%98%E5%BA%A7%E6%A8%99%E3%80%812%E5%88%97%E7%9B%AE%E3%81%8C%EF%BD%99%E5%BA%A7%E6%A8%99%E3%81%A7%E3%81%99.txt");
websave("fit_ellipse.zip", "https://jp.mathworks.com/matlabcentral/mlc-downloads/downloads/submissions/3215/versions/1/download/zip");
unzip("fit_ellipse.zip");
fit_ellipse(T.Var1, T.Var2, axes) % fit_ellipse - finds the best fit to an ellipse for the given set of points.
hold on;
scatter(T.Var1, T.Var2, 5);