function dy = vdpp(t, y, flag,u)
    dy = zeros(2,1);
    dy(1) = y(2);
    dy(2) = u*(1-y(1)^2)*y(2)-y(1);
end