clear all; close all;
load gaussfit.dat;
x = gaussfit(:,1);
y = gaussfit(:,2);
plot(x,y,'o')
x0 = [1 1];
coeff = fminsearch('gafit',x0);
xga = linspace(-3,3,100);
a = coeff(1); 
b = coeff(2);
yga = a*exp(-b*xga.^2);
plot(x,y,'o',xga,yga,'k')