%phase_limit_linnolin_fun1.m

function dy=phase_limit_linnolin_fun1(t,y)

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

dy=[dy1; dy2];