function E = gafit(x0)
load gaussfit.dat
x = gaussfit(:,1);
y = gaussfit(:,2);
E = sum( ( x0(1)*exp(-x0(2)*x.^2)-y).^2);