Matlab - ODE45 tspan

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

a=0;

b=1;

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

c=0.5;

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

%fValor de la constante K

K=3;

%Calculamos 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(N,1);

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

%evolucionan en el tiempo la fase theta.

dt=0.5; %intervalos de tiempo

R=1;

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

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

    q=sin(phi-theta)

    %funcion que queremos integrar

    f=@(t,theta)(w+(K*R*q));

    %ode45 integra funcion

    [t,theta]=ode45(f,[0 0.5],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