function sse=myfit(params,x,y)
A = params(1);
B = params(2);
C = params(3);
Y = A.*(1-exp(-B.*x))+C;
err = Y - y;
sse=sum(err.^2);