function yd= ode(t,y);
global a b c r;

yd= [y(2); (a*r^2*exp(2*i*t)+b*r*exp(i*t)+c)*y(1)];

