Correcting fit functions to use delays, in some rare conditions non-offseted time values were still used

matlab/fitSingle.m

@@ -26,7 +26,7 @@ function fv=fitSingle(globalPar, par, cax)
     if alpha==1/tau
         fv=cax1.*e0;
         if k2==1/tau
-            r0=k1.*cax1.*cax.*ek;
+            r0=k1.*cax1.*cax1.*ek;
         else
             r0=kt*(cax1.*e0-(e0-ek)/(k2-1/tau));
         end
@@ -54,4 +54,4 @@ function fv=fitSingle(globalPar, par, cax)
 
 end
 
-%par=ones(2,
+%par=ones(2,

matlab/fitSingleZero.m

@@ -13,7 +13,7 @@ function fv=fitSingleZero(globalPar, cax)
     
     cax1=cax-dt;
     if tau==1/alpha
-        fv=cax1.*exp(-cax/tau);
+        fv=cax1.*exp(-cax1/tau);
     else
         fv=tau*(exp(-alpha*cax1)-exp(-cax1/tau))/(1-alpha*tau);
     end
@@ -25,4 +25,4 @@ function fv=fitSingleZero(globalPar, cax)
 
 end
 
-%par=ones(2,
+%par=ones(2,