function ell_test()

I = imread('ellipse.png');
x0 = 95; y0 = 125;
[J,p] = ell(I, 100, 40, x0, y0, 0, 256, 512, 1.5);
subplot(231); imshow(I); hold on; plot(x0, y0, 'r+'); % center of transform
subplot(232); imshow(J);
I = iell(J, p);
subplot(233); imshow(I);

I = imread('ellipse45.png');
x0 = 138; y0 = 148;
[J,p] = ell(I, 100, 40, x0, y0, pi/4, 256, 512, 1.5);
subplot(234); imshow(I); hold on; plot(x0, y0, 'r+'); % center of transform
subplot(235); imshow(J);
I = iell(J, p);
subplot(236); imshow(I);