I was going to use this code to fit Gaussians to a known curve by minimizing the known - fit residual. I tried for 3 and 4 Gaussians, meaning 9 or 12 variables. This newer code fails giving the error:
Too many input arguments.
Error in ==> fminsearchbnd>@(x,varargin)fun(xtransform(x),varargin{:}) at 233
intrafun = @(x, varargin) fun(xtransform(x), varargin{:});
Error in ==> fminsearch at 205
fv(:,1) = funfcn(x,varargin{:});
Error in ==> fminsearchbnd at 264
[xu,fval,exitflag,output] = fminsearch(intrafun,x0u,options,varargin);
The original code by John D'Errico works perfectly.