ydot=eq_2disc02b(t,y,mu,K,F0);

function ydot=eq_2disc02b(t,y,mu,K,F0);
%Equation of the double disc dynamo chaotic attractor for a "TOY" model.
%Volobuev, D. "Toy" dynamo to describe the long-term solar activity cycles. Solar Physics 238, 421-430 (2006)
%
%Example of MatLab 5.x-6.x command line call:
%
%options=odeset;[T,Y] = ode45(@eq_2disc02b,1:1:120,[-3 0 -8 ],options,1,2,-11);figure;plot(T,Y);legend('I_1','I_2','\omega')
%
%mu=1;K=2; F0=-11;
A=mu*(K^2-K^-2);%constants, see paper, pp.423-424
dt=.06;%dt=1;%this is a time scale multiplier; dt=.06 is approximately  time in years, in paper (Figure 2), dt=1
%Equation 1, p.423
ydot(1)=dt*(-mu*y(1)+y(3)*y(2));
ydot(2)=dt*(-mu*y(2)+(y(3)-A)*y(1));
ydot(3)=dt*(F0-y(1)*y(2));%
ydot=ydot(:);