function [chi y]=e2fit(x1,r,img,xs,x2)

x=[x1 x2];

r0=x(1)+i*x(2);
r1=x(3)+i*x(4);
a0=x(5)+i*x(6);
b0=x(7)+i*x(8);
a1=x(9)+i*x(10);
b1=x(11)+i*x(12);


y=abs(real((r-r0).*a0)+i*real((r-r0).*b0)) + abs(real((r-r1).*a1)+i*real((r-r1).*b1));

chi=mean((img(:)-y(:)).^2);

%disp([log10(chi)-2.7944 xs-x]);