%phase_limitcycle_fun.m

function dy=phase_limitcycle_fun(t,y)

dy1=y(1)+y(2)-y(1)*(y(1)^2+y(2)^2);
dy2=-2*y(1)+y(2)-y(2)*(y(1)^2+y(2)^2);

dy=[dy1; dy2];