@@ -80,8 +80,8 @@ function [globalPar, cPars] = ...
%first estimate input functions
%take function with maximum response as IVF surrogate
- [m mi]=max(centers);
- [fm mj]=max(m);
+ [m, mi]=max(centers);
+ [~, mj]=max(m);
cay=centers(mi(mj),:)';
w=ones(size(cay));