function y=gafit(x0)

load gaussfit.dat
x=gaussfit(:,1);
y=gaussfit(:,2);

sum=0;
for j=1:length(x)
   sum=sum+ ( x0(1)*exp(-x0(2)*x(j)^2) ...
      -y(j)  )^2 ;
end
y=sum;