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

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

dy1=y2;
dy2=-y1;

dy=[dy1; dy2];
