matlab 数学建模 图像处理

I = imread('ellipse.BMP');
Ig = rgb2gray(I);
imshow(Ig);
hold on
conicP = ginput(15);
plot(conicP(:,1)', conicP(:,2)', 'xr');     %drawing sample points
a0 = [1 1 1 1 1 1];
f = @(</