%%
%% Listing 7.13 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 [du, flag, data] = cvrdrhs(t, u, data)
global Lap2D;
global beta;

du = -beta*(Lap2D*u) + u - u.^3;
flag = 0;
end

