%phase_limit_linnolin_fun2.m

function dy=phase_limit_linnolin_fun2(t,y)

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

dy=[dy1; dy2];