function r=fitexp2(d,X);
% r is sum(residual²)
x=X(:,1);
y=X(:,2);
g=gaussian(x,d(1),d(2),d(3));
h=gaussian(x,d(4),d(5),d(6));
ymod=g+h;
r=sum((ymod-y).^2);
#plot(x,y,'o',x,ymod);
#drawnow;
