Matlab - Error: matrix dimension must agree

>> clear all

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

c=0;

 b=1;

 G=c+(b-c)*rand(N,1);

 w=(1./(-2*G+sqrt(4*G.^2-1)).^(1/3))+(-2*G+sqrt(4*G.^2-1)).^(1/3);

 disp(w);

 Introduce el numero de osciladores :5

   1.5029 + 0.0000i

   1.0103 - 0.2490i

   1.0904 - 0.7531i

   1.0917 - 0.7588i

   1.6155 + 0.0000i

% 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)

c=0;

b=1;

G=c+(b-c)*rand(N,1);

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

%frecuencia parabolica en funcion de G(w)

%a=1;%parabolica entre -1 y 1

w=1./((1-2*G+2*sqrt(-G+G.^2))).^(1/3)+(1-2*G+2*sqrt(-G+G.^2)).^(1/3)

disp(w);