function [y] = Rabi_xN_Func(n, v, t) offset = v(1)/n; %Read parameters. Order: amplitude, center, width for i=1:n amplitude(i) = v(2+4*(i-1)); freq(i) = v(3+4*(i-1)); dec(i) = v(4+4*(i-1)); phs(i) = v(5+4*(i-1)); end y=0; for i=1:n y = y + offset + amplitude(i).*cos(2*pi*freq(i)*t + phs(i)).*exp((-t/dec(i))); end