function [rrsfor,yfitfor,residualfor] = fitfor(CfornoxEox,xhigh,yhigh,handle)
nox =1;
yfitfor= CfornoxEox(1)./[1+10.^(nox*(-CfornoxEox(2)+xhigh)/0.06)];
residualfor = yhigh-yfitfor;
rrsfor = sum(residualfor.^2);
% ssefor = sum(errorfor.^2);

% ssrfor = sum((yfitfor-mean(yfitfor)).^2);
% sstfor = ssefor + ssrfor;
% r2 = ssrfor/sstfor;

set(gcf, 'DoubleBuffer','on');
set(handle,'ydata',yfitfor)
drawnow
pause(.08)