Adding ploting of IVF fits in fitCentersIVF

matlab/fitCentersIVF.m

@@ -102,7 +102,7 @@ function [globalPar, cPars] = ...
     %ofile=fullfile(path,sprintf('CIFinal_%d_%d.txt',nclass,realizationId));
     %writematrix(fci,ofile,'Delimiter','tab');
     
-    globalPar=[fp(1:3) 0 0 0];
+    globalPar=[fp(1:4) 0 0];
     
     cPars=reshape(qfit,[],nclass)';
     %reset k1 for those where BVF is nearly one
@@ -111,8 +111,8 @@ function [globalPar, cPars] = ...
     %cPars(cPars(:,3)>1,1)=0;
     
     
-    %opath=fullfile(path,patientID);
-    %plotFunctions(opath,code,globalPar,cPars,cax,centers);
+    opath=fullfile(path,patientID);
+    plotFunctions(opath,code,globalPar,cPars,cax,centers);
     
     
 end