Matlab - Lorentziana

% N numero osciladores con los que trabajamos

N=input(' Introduce el numero de osciladores :');

%Creamos un vector ramdon de frecuencias iniciales w para N osciladores

%a y b intervalos de la distribucion acumulada G(w)

a=0;

b=1;

v=a+(b-a)*rand(N);

disp(v);

G=v(1,:);

%Ahora entonces como tenemos numero aleatorios entre 0 y 1 , el valor de la

%frecuencia lorentziana en funcion de G(w)=tan^(-1)(w/a)/pi que es la integral de la g(w) lorentziana, si despejo la w:

c=0.5;

w=c*tan(G*pi);

disp(w)

n=input(' Introduce el numero de osciladores :');

x=rand(1,n);

s=10;

m=100;

corte=10*s;

f=s*tan(pi*x)+m;

hist((f(abs(f-m)<corte)),n^(1/2))