%%
%% Listing 7.9 w książce "Obliczenia inżynierskie i naukowe", PWN 2011
%% (C) Piotr Krzyżanowski, 2011
%%
%% Możesz używać tego kodu zgodnie z licencją Creative Commons BY-SA, 
%% zob. http://creativecommons.org/licenses/by-sa/3.0/pl/
%%

function dY = vanderpoleps(t, Y, epsilon)
dY = [Y(2); ...
        epsilon*(1-Y(1)^2)*Y(2) - Y(1)];
end
