%SingularPerturbation_Example_Phase_fun.m

function dy=SingularPerturbation_Example_Phase_fun(t,y,epsilon)

y1=y(1);
y2=y(2);

dy1=-y1+y2;
dy2=(-atan(y1+y2))/epsilon;

dy=[dy1; dy2];