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. MODEL CREATION You may build a model using any of: nlorz('defaults') builds a 2 Lorentzians model nlorz(n) builds an n Lorentzians model nlorz and nlorz('gui') shows a Dialogue to enter the number of Lorentzians. nlorz([ ... ]) use a 3n+1 values parameters for n Lorentzians Reference: http://en.wikipedia.org/wiki/Lorenztian_function MODEL EVALUATION 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: Nov. 26, 2018 See also iFunc, iFunc/fits, iFunc/plot, ngauss (c) E.Farhi, ILL. License: EUPL.