iFit/nlorz

 y = nlorz(p, x, [y]) : multiple Lorentzians

   iFunc/nlorz multiple Lorentzian fitting function
     y = sum p(i)*exp(-0.5*((x-p(i+1))/p(i+2)).^2) + p(end);

   This expression assumes that the Amplitudes are independent from the Widths.

 to initiate n Lorentzian use: nlorz(n)
 will result in an iFunc model of n Lorentzian functions.

 Reference: http://en.wikipedia.org/wiki/Lorenztian_function

 input:  p: multiple Lorentzian model parameters (double)
            p = [ Amplitude1 Centre1 HalfWidth1 ... BackGround ]
          or 'guess'
         x: axis (double)
         y: when values are given and p='guess', a guess of the parameters is performed (double)
 output: y: model value
 ex:     y=nlorz([1 0 1 0.5 2 0.5 0], -10:10); or plot(nlorz(3))

 Version: Aug. 22, 2017
 See also iFunc, iFunc/fits, iFunc/plot, ngauss
 (c) E.Farhi, ILL. License: EUPL.