Matlab - Error: DImensions must agree

%Creamos un vector w que contiene N frecuencias aleatorias distribuidas segun lorentziana.

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

a=0;

b=1;

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

c=0.5;

w=c*tan(G*pi-pi/2);

%K= constante, ponemos 3

K=3;

%Calculamos unas condiciones iniciales para las thetha(i) que la vamos a llamar

%xo(dentro del bucle de tiempo)

z=0;

v=2*pi;

xo=z+(v-z)*rand(1,N)

%Tenemos que calcular R que va entre 0 y 1 y el angulo medio phi.

dt=0.5; %intervalos de tiempo

R=0;

phi=0;

theta=xo;

for t=2:1:202      %200 (t)time steps

  rx=0;

  ry=0;

    TO=(t-1)*dt; %tiempo inicial

    T=t*dt;      %tiempo final dentro del intervalo

    %parámetros ode45: tiempo integracion y opciones

    tspan=[TO   T];

    options=odeset('RelTol',1e-6);

    %funcion que queremos integrar

  -->  f=@(t,theta)( w + (K*R)*sin(phi-theta) );

    %ode45 integra funcion

    [t,theta]=ode45(f,tspan,xo,options);

    yo=theta;

    rx=(1/N)*cos(yo);

    RX=sum(rx);

    ry=(1/N)*sin(yo);

    RY=sum(ry);

  %bucle N

R=sqrt(RX^2+RY^2);

phi=atan2(RY,RX);

plot(T,R);

hold on

end  %bucle t

[t,y] = ode45(@(t,y) cos(t.^2), [0.9 55], 0);

[~,w] = findpeaks(y);

figure

plot(t,y,'-b')

hold on

plot(t(w),y(w),'*r')

fprintf('tiene  %i máximos locales en este intervalo  \n',length(w))