%Perturbation_VanderPol_fun1.m
function dy=Perturbation_VanderPol_fun1(t,y,epsilon)

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

dy1=y2;
dy2=-y1-(1-cos(t)^2)*sin(t);

dy=[dy1; dy2];
