
function ZE = test(SE,scale,xoffset,yoffset,sigma)

pts = xoffset + (-50:50)';

ZE = scale * exp(-(pts/sigma).^2) + yoffset;

figure(2), clf, hold on
plot(SE)
plot(ZE,'r')

title(sprintf('Sum-Squared Difference = %f',sum((SE-ZE).^2)))
