function [rrsrevrd,yfitrev,residualrevrd] = fitrevrd(CrevnoxErd,xhigh,yhigh,handle) nox=1; yfitrevrd= (CrevnoxErd(1)./[1+10.^(nox*(-CrevnoxErd(2)+xhigh)/0.06)]); residualrevrd = yhigh-yfitrevrd; rrsrevrd = sum(residualrevrd.^2); set(gcf, 'DoubleBuffer','on'); set(handle,'ydata',yfitrevrd) drawnow pause(.08)